fortran90源程序最小二乘法解多项式拟合问题.docVIP

下载本文档

169
0
约2.37千字
约 3页
2019-07-27 发布于江西
举报
版权申诉

fortran90源程序最小二乘法解多项式拟合问题.doc

1、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

最小二乘法解多项式拟合问题 doubleprecision a(100,100),c(100,100),d(100),alpha(100) doubleprecision gs(100,100),temp,q,delt(100),x(100),y(100) integer?gsn open(100,file=in.txt,status=old) read(100,*)?n,m read(100,*)?(x(i),i=1,n),(y(j),j=1,n) close(100) do?1?i=1,n ?do?2?j=1,m+1 910 ??a(i,j)=x(i)**(j-1) 2 ?continue 1 continue open(200,file=process.txt,status=unknown) write(200,*)?a(m,n):????[,m:,n,???????????n:,(m+1),] do?3?i=1,n ?write(200,35)?(a(i,j),j=1,(m+1)) 3 continue do?4?i=1,m+1 ?do?7?j=1,m+1 920 ??c(i,j)=0 ??do?8?k=1,n ???c(i,j)=c(i,j)+x(k)**(i+j-2) 8 ??continue 7 ?continue 4 continue do?9?i=1,m+1 ?d(i)=0 ?do?11?k=1,n 929 ??d(i)=d(i)+y(k)*(x(k)**(i-1)) 11 ?continue 9 continue write(200,*)?c(n,n):????[,n:,(m+1),] do?12?i=1,m+1 ?write(200,35)?(c(i,j),j=1,m+1) 12 continue write(200,*)?d(n):????[,n:,] write(200,35)?(d(i),i=1,m+1) close(200) loop=0 940 gsn=m+1 open(300,file=equation.txt,status=unknown) write(300,*)?The?added?matrix?is: do?13?i=1,gsn ?write(300,35)?(c(i,j),j=1,gsn),d(i) 13 continue do?14?i=1,gsn ?do?15?j=1,gsn ??gs(i,j)=c(i,j) 15 ?continue 950 ?gs(i,gsn+1)=d(i) 14 continue do?16?k=1,gsn+1 ?do?17?i=k+1,gsn+1 ??if(dabs(gs(k,k)).lt.dabs(gs(i,j)))?then ???do?18?j=1,gsn+1 ????temp=gs(k,j) ????gs(k,j)=gs(i,j) ????gs(i,j)=temp 18???????continue ??endif 17?????continue ???????temp=gs(k,k) ?do?19?j=K,gsn+1 ??gs(k,j)=gs(k,j)/temp 19 ?continue ???????do?21?i=k,gsn+1 ??temp=gs(i,j) ??do?22?j=k,gsn+1 ???gs(i,j)=gs(i,j)-temp*gs(k,j) 22? ??continue 21 ?continue ?write(300,*)?------------------------------ ?loop=loop+1 ?write(300,*)?loop:,loop ?do?23?i=1,m ??write(300,35)?(gs(i,j),j=1,gsn+1) 23 ?continue 16 continue alpha(gsn)=gs(gsn,gsn+1) do?24?k=gsn-1,1,-1 ?temp=0.0d0 ?do?25?j=k+1,gsn ??temp=temp+gs(k,j)*alpha(j) 25 ?continue 24 continue write(300,*)?------------------------------ write(300,35)?(alpha(i),i=1,gsn) q=0.0d0 do?26?i=1,n ?delt(i)=-y(i) ?do?27?j=1,m+1 ??delt(i)=delt(i)+(x(i