Matlab的离散点曲线导数曲率数值模拟方法.doc

《Matlab的离散点曲线导数曲率数值模拟方法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、Matlab的离散点曲线导数曲率数值模拟方法曲率公式：曲率计算.JPG(2.31KB)2010-3-717:01模拟方法：复制内容到剪贴板代码:clc;clearall;closeall;x0=linspace(0,1);y0=sin(x0).*cos(x0);h=abs(diff([x0(2),x0(1)]));%模拟一阶导figure;boxon;holdon;ythe1=cos(x0).^2-sin(x0).^2;%理论一阶导yapp1=gradient(y0,h);%matlab数值近似plo

2、t(x0,ythe1,'.');plot(x0,yapp1,'r');legend('理论值','模拟值');title('模拟一阶导');%模拟二阶导figure;boxon;holdon;ythe2=(-4)*cos(x0).*sin(x0);%理论二阶导yapp2=2*2*del2(y0,h);%matlab数值近似plot(x0,ythe2,'.');plot(x0,yapp2,'r');legend('理论值','模拟值');title('模拟二阶导');%模拟曲率symsxyy=sin(x

3、)*cos(x);yd2=diff(y,2);yd1=diff(y,1);k=abs(yd2)/(1+yd1^2)^(3/2);k1=subs(k,x,x0);k2=abs(yapp2)./(1+yapp1.^2).^(3/2);figure;boxon;holdon;plot(x0,k1,'.');plot(x0,k2,'r');legend('理论值','模拟值','Location','NorthWest');title('模拟曲率');效果：未命名.bmp(595.01KB)2010-3-71

4、7:01未命名1.bmp(574.92KB)2010-3-717:01