plot3([1,4,7],[2,5,8],[3,6,9],[1,4,7],[2,5,8],[3,6,9],'b.'册拦消,'州知衡兆markersize',25)grid on
下面是我编的一个简单程序:交流学习一下...
clear, close all
x=[0,3,5,7,9,11,12,13,14,15] y=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6]
plot(x,y)
xi=0:0.1:15
yi_cubic=interp1(x,y,xi,'cubic')
plot(x,y,'ro',xi,yi_cubic)
pp=csape(x,y,'second') v=ppval(pp,xi)
v%x坐标每改变0.1时的y坐标
T=(ppval(pp,0.1)-ppval(pp,0))/0.1%求零点的导数值
angle=atan(T)*180/pi%得到零点的切线角度约为23.7843°
s=v(130:151)%将13≤x≤15的y值赋给数组s
ss=min(s)%求出最小值为0.9721
下面是插值完成后历拦的曲线,x坐标每改变0.1时的y坐标可笑烂伏碰携以由数组v读出,斜率使用近似求法求出,如果还有更好的方法请分享~~
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)