untitledquxiannihe.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 9 行

fprintf('i     t(i)      y(i)     y=a(1)*x^2+a(2)*x+a(3)     偏差\n');
for i=1:16
    t=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];
    y=[4.00,6.40,8.00,8.80,9.22,9.50,9.70,9.86,10.00,10.20,10.32,10.42,10.50,10.55,10.58,10.60];
    a=polyfit(t,y,2)
    error=y(i)-(a(1)*t(i)^2+a(2)*t(i)+a(3));
    fprintf('%2d       %2.1f       %5.2f       %10.4f       %10.5f\n',i,t(i),y(i),a(1)*t(i)^2+a(2)*t(i)+a(3),error);
end