untitledquxiannihe2.m

数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序

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];
plot(t,y,'r*')
legend('实验数据(ti,yi)')
xlabel('t'),ylabel('y'),
title('数据散点图')

syms a1 a2 a3
t=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];
fi=a1.*t.^2+a2.*t+a3

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];
fi=a1.*t(i).^2+a2.*t(i)+a3;
fy=fi-y;fy2=fy.^2;J=sum(fy.^2)

Ja1=diff(J,a1);Ja2=diff(J,a2);Ja3=diff(J,a3);
Ja11=simple(Ja1),Ja21=simple(Ja2),Ja31=simple(Ja3)

A=[2992 36992 487696;272 29992 36992;32 272 29992];
B=[771842/25 68447/25 2973/10];
C=B/A,f=poly2sym(C)

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];
n=length(xi);
f=13.9206*ti^2-15.1955*ti-207.6085;
t=1:0.5:16;
F=13.9206*t^2-15.1955*t-207.6085;
fy=abs(f-y);fy2=fy.^2;EW=max(fy),