数值分析上机实验最小二乘法.docVIP

数值分析实验报告五 最小二乘法 题目 设有如下数据 -3 -2 -1 0 1 2 3 -1.76 0.42 1.2 1.34 1.43 2.25 4.38 用三次多项式拟合这组数据，并绘出图形。 方法 最小二乘法 程序 M文件： syms x f; xx=input(请输入插值节点 as [x1,x2...]\n); ff=input(请输入插值节点处对应的函数值 as [f1,f2...]\n); m=input(请输入要求的插值次数m=); n=length(xx); for i=1:(m+1) syms fai x; fai=x^(i-1); for j=1:n x=xx(j); H(i,j)=eval(fai); end end A=ff*(H)*inv(H*(H)); syms x; f=0; for i=1:(m+1) f=f+A(i)*x^(i-1); end f plot(xx,ff,*) hold on ezplot(f,[xx(1),xx(n)]) 结果 save and run 之后： 请输入插值节点 as [x1,x2...] [-3 -2 -1 0 1 2 3] 请输入插值节点处对应的函数值 as [f1,f2...] [-1.76 0.42 1.2 1.34 1.43 2.25 4.38] 请输入要求的插值次数m=3 f = 133/