duoxiaongshinihelizi.m

来自「多项式拟合的一个小例子」· M 代码 · 共 22 行

 x=[0 1 2 3 4 5]; 

 y=[0 20 60 68 77 110]; 

 newx=0:0.05:5; % 以新阵列形成更小增量以利回归计算及绘图 

for n=2:9 

f(:,n)=polyval(polyfit(x,y,n),newx)'; 

 plot(newx,f(:,n),x,y,'o') 

title(['Poly. regression, deg=',int2str(n)]) 

 xlabel('Time'), ylabel('Temp'), grid 

 pause % 每次要暂停看清楚图再执行下一步骤 

 end