dianzu2.m

来自「1、掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法」· M 代码 · 共 11 行

t=[20.5 32.5 51 73 95.7];
r=[765 826 873 942 1032];                  
a=polyfit(t,r,1); 
%R=[t' ones(5,1)]; % r1=t,r2=1
%a=R\r';% Ra=r
a1=a(1)
a2=a(2)
y=polyval(a,t);
plot(t,r,'k+',t,y,'r')