- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
最小二乘法程序说明及流程图
陆韶琦 3110000441
程序说明:本程序用多项式拟合数据,程序会要求输入需要拟合的次数和数据点的个
数,数据文件应该保存在本程序运行时的 current folder 下,文件取名为“mytext.txt”
流程图: 输入需要拟合的次数
N ,和数据点的个数M
false
0i=N
+1 i+=1
false
0j=M
j+=1 C (i,j )=S(1,j)^(i-1)
(S(1,j)=0 时,C (i.j )=0)
T
A=C C,b=CY;
(Y 为数据点中yi 组成的
列向量)
通过列主元高斯消元法求解Ax=b
输出表达式
N
i
expr= a i + 1 x
i=0
对数据点S (1,1 )到S (1 ,M )区间按0.01 步长
作图
M 2
计算δ = δi ;
rms 1
δ =max |δ|.
max 1=i= i
end
程序代码:
%多项式最小二乘法拟合数据
N=input( please put in how many times the power will you overfit:);
M=input( how many couples of statistics are there in the table:);
%读入数据文件
f=fopen(mytxt.txt,r);
S=fscanf(f,%g,[M 2]);
fclose(f);
S=S;
%显示数据文件,确保正确输入
disp(S(x,y)=);
disp(S);
%建立多项式系数法方程组中间矩阵
C=zeros(N+1,M);
for i=1:N+1
for j=1:M
if S(1,j)==0
C(i,j)=0;
else
C(i,j)=S(1,j).^(i-1);
end
end
end
%建立法方程组
A=C*C;
Y=zeros(M,1);
for i=1:M
Y(i,1)=S(2,i);
end
b=C*Y;
%用列主元高斯消元法接法方程组
A=[A,b];
for i=1:N+1
max=abs(A(i,i));
for j=i+1:N+1
if abs(A(j,i))max
flag=j;
max=A(j,i);
end
end
for k=i:N+2
文档评论(0)