>> x=-5:0.01:5;
>> y=x.^5+x.^4+4*random('Normal',0,1,1,length(x));
>> p=polyfit(x,y,5);
>> f=polyval(p,x);
>> plot(x,y,':r',x,f,'b');
polyfit 和polyval就是做这个的 你如果是要用在计算方法导论这类课的作业的话恐怕不合适
附图一张局部的 好像不清楚
