注意!!!syms属性不能直接用,我们先要去mathtool商店下载对应插件,下载方法首次执行报错后有指示。
Matlab
clc;clear;close all;
xi=[37;38;39;40;41;42;43];
yi=[3.40;3.00;2.10;1.53;1.80;1.90;2.90];
A=[xi.^0,xi.^1,xi.^2]; %生成方程组系数矩阵,以向量生成矩阵
AA=A'*A; %生成法方程ATA
a=inv(AA)*(A'*yi) %求解系数,inv为求逆矩阵
syms x
px=a(1)+a(2)*x+a(3)*x.^2;
xx=linspace(37,43,100); %在37和43之间平均取100个点
yy=subs(px,x,xx); %求100点的函数值
figure
hold on %曲线图合并在一起
plot(xi,yi,'o') %画出数据点
plot(xx,yy,'-') %拟合曲线