最小二乘法C语言编程.docxVIP

下载本文档

24
0
约2.85千字
约 5页
2017-08-09 发布于浙江
举报
版权申诉

最小二乘法C语言编程.docx

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

最小二乘法C语言编程

#include hidef.h /* common defines and macros */#include derivative.h /* derivative-specific definitions */float a[3][3]={0};float a8[]={1,2,3,4,5,6,7,8,9,10,11,12,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27} ;float a9[]={147,155,146,146,144,146,145,144,143,146,146,149,147,148,150,152,155,157,160,163,168,173,175,177,180,199,203} ; float b[3][1]={0}; float b3[3][1]; float b4[3][3]; float b1[3][3]={0}; float a1[3][3]={0}; float x,x1,x2,x3,x4,x5,x6,x7,x8,x9,c,*p,*pp,*ppp; int i,j,ii,jj,d1=0;void main(void) { /////////////////////最小二乘初始化////////////////////////// x=0;p=a9;for(i=0;i10;i++)x+=*(p++);pp=*(b+0);*(pp++)=x; x=0; p=a8; for(ii=0;ii3;ii++) { if(ii==0){ pp=*(a+0); *(pp++)=10; for(i=0;i10;i++) x1=x+=*(p++); *(pp++)=x; p=a8; for(i=0;i10;i++) { x=*(p++); x2+=x*x; } *(pp++)=x2; } if(ii==1){ pp=*(a+1); *(pp++)=x1; *(pp++)=x2; p=a8; x1=0; for(i=0;i10;i++) { x=*(p++); x1+=x*x*x; } *(pp++)=x1; } if(ii==2) { pp=*(a+2); *(pp++)=x2; *(pp++)=x1; x3=0; p=a8; for(i=0;i10;i++) { x=*(p++); x2=x*x*x*x; x3+=x2; } *(pp++)=x3; } }pp=*(b+1); p=a9; ppp=a8; x3=0; for(i=0;i10;i++) { x=*(p++); x2=*(ppp++); x3+=x*x2; } *pp=x3; pp=*(b+2); p=a9; ppp=a8; x3=0; for(i=0;i10;i++) { x=*(p++); x2=*(ppp++); x3+=x*x2*x2; } *pp=x3; ///////////////////////求矩阵的模//////////////////////////////// p=*(a+0); c=*p; p=*(a+1); p++; for(i=1;i3;i++) {for( j=1;j3;j++) {switch (d1) { case 0:x=x1= *(p++); d1++;break; case 1:x4=x5= *(p++); d1++;break; case 2:x2= *(++p); d1++;break; case 3:x3= *(++p); break; } } }x7=c*(x1*x3-x5*x2); p=*(a+1); c=*p; p=*(a+0); p++; x1=*(p++); x5=*p; x8=c*(x1*x3-x5*x2); p=*(a+2); c=*p; x9=c*(x1*x4-x5*x); c=1/(x7-x8+x9);/////////////////求伴随阵////////////////// p=*(a+0); x1=*(p++); x2=*(p++); x3=*(p++); p=*(a+1); x4=*(p++); x5=*(p++); x6=*(p++); p=*(a+2); x7=*(p++); x8=*(p++); x9=*(p++); p=*(b1+0); *(p++)=(x5*x9-x6*x8); *(p++)=-1*(x2*x9-x3*x8); *(p++)=(x2*x6-x3*x5); p=*(b1+1); *(p++)=-1*(x4*x9-x6