工程数值方法实验题

 matlab代码综合运用

clear
format long
data=load('a-t.txt');
t=data(:,1);
a=data(:,2);
scatter(t,a,'g.');
title('散点图及拟合曲线图');
xlabel('t-时间');
ylabel('a-加速度');
hold on
Q=[t.^0 t.^1 t.^2 t.^3 t.^4 t.^5];
k=(inv(Q'*Q))*(Q'*a);
fprintf('拟合方程系数（最低次到最高次）\n');
disp(k);
x=0:0.0001:3;

y=k(1)+k(2).*x+k(3).*x.^2+k(4).*x.^3+k(5).*x.^4+k(6).*x.^5;
plot(x,y,'-r');
fprintf('预测t=2.9时对应的a\n');
disp(y(29000)) 
gtext('拟合（2.9）=0.047405457695625');