x = -1021;y=11+25x^2;f1=polyfitxy3;f=poly2symf1;y1=polyvalf1x;x2=linspace-1110;y2=interp1xyx2;plotxyr-xy1b-hold onplotx2y2klegend数据点3次拟合曲线3次多项式插值xlabelXylabelYf = poly2strf1x代码逐句解析

x = -1:0.2:1;：定义一个从-1到1，步长为0.2的向量x。
y=1./(1+25.*x.^2);：计算y的值，其中y=1/(1+25x^2)。
f1=polyfit(x,y,3);：用3次多项式拟合x和y，得到拟合后的多项式系数。
f=poly2sym(f1);：将多项式系数转换为符号表达式。
y1=polyval(f1,x);：用拟合后的多项式系数计算x对应的y值。
x2=linspace(-1,1,10);：生成一个从-1到1的等差数列，共10个数。
y2=interp1(x,y,x2);：用3次多项式插值算法，将x、y插值成x2对应的y值。
plot(x,y,'r*-',x,y1,'b-')：绘制x,y的散点图和3次多项式拟合曲线。
hold on：保持当前图形，以便在同一图形中添加其他图形元素。
plot(x2,y2,'k')：绘制3次多项式插值曲线。
legend('数据点','3次拟合曲线','3次多项式插值')：添加图例，标注数据点、3次拟合曲线和3次多项式插值曲线。
xlabel('X'),ylabel('Y')：添加坐标轴标签。
f = poly2str(f1,'x')：将多项式系数转换为字符串