mm2003.m

《精通matlab7》“mastering matlab 7”的代码。

% mm2003.m

x = [0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1];
y = [-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2];
n = 2;
p = polyfit(x,y,n);
xi = linspace(0,1,100);
yi = polyval(p,xi);
pp = polyfit(x,y,10);
y10 = polyval(pp,xi); % evaluate 10th order polynomial
plot(x,y,'o',xi,yi,'--',xi,y10) % plot data
xlabel('x'), ylabel('y=f(x)')
title('Figure 20.3: 2nd and 10th Order Curve Fitting')