例2-64.m

来自「这是一个MATLAB7.0基础与提高例题的所有源码」· M 代码 · 共 19 行

>> x=0:0.1:1;
>> y=[-0.447 1.973 3.16 6.43 7.07 7.58 7.68 9.78 9.43 9.32 11.2];
>> n=2;
>> p=polyfit(x,y,n)
>> xi=linspace(0,1,100);
>> z=polyval(p,xi);
>> plot(x,y,'o',x,y,xi,z,':')
>> xlabel('x'),ylabel('y'),title('二次曲线拟合')
>> pd=polyder(p)
>> z=polyval(pd,xi);
>> plot(xi,z)     
>> xlabel('x'),ylabel('dy/dx')
>> title('曲线微分')
>> dy=diff(y)./diff(x);
>> xd=x(1:max(size(x))-1);
>> plot(xd,dy)
>> title('用diff进行近似差分')
>> ylabel('dy/dx轴'),xlabel('x轴')