最小二乘法程序说明及流程图.pdf

  1. 1、本文档共5页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 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)

xxj1658888 + 关注
实名认证
内容提供者

教师资格证持证人

该用户很懒,什么也没介绍

领域认证该用户于2024年04月12日上传了教师资格证

1亿VIP精品文档

相关文档