你好,你这里还有指数啊,所以类似的,我也定义一个函数
x = -pi:0.1:pi; y = exp(sin(x)); p = polyfit(x,log(y),3);%取n=3 yfit = exp(p(end))*exp(polyval([p(1:end-1),0],x)); plot(x,y,'r-', x,yfit,'b-') legend('y','yfit')