欢迎来到天天文库
浏览记录
ID：57687280
大小：35.50 KB
页数：18页
时间：2020-09-01
《C语言实现矩阵的LU分解、施密特正交化、Givens分解、Householder分解.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、C语言实现矩阵的LU分解、施密特正交化、Givens分解、Householder分解ByKim.Wang，UCAS#include#include#include#defineHS10#defineLS10intn,m;floata[HS][LS],bc[HS][LS];voidgivens(){floatfm,sc,cos,sin,r[HS][LS],q[HS][LS],swap[HS][LS],p[HS][LS];intih,jh,i,
2、j,kh,iw;for(i=0;i3、*r[i][j]+r[j][j]*r[j][j]);cos=r[j][j]/fm;sin=r[i][j]/fm;p[i][i]=p[j][j]=cos;p[i][j]=-sin;p[j][i]=sin;for(ih=0;ih4、[jh]=swap[ih][jh];for(ih=0;ih5、)printf("%0.5f",q[i][j]);printf("");}printf("Givens约减矩阵R的结果如下：");for(i=0;i6、][LS],p[HS][LS];for(i=0;i7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
3、*r[i][j]+r[j][j]*r[j][j]);cos=r[j][j]/fm;sin=r[i][j]/fm;p[i][i]=p[j][j]=cos;p[i][j]=-sin;p[j][i]=sin;for(ih=0;ih4、[jh]=swap[ih][jh];for(ih=0;ih5、)printf("%0.5f",q[i][j]);printf("");}printf("Givens约减矩阵R的结果如下：");for(i=0;i6、][LS],p[HS][LS];for(i=0;i7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
4、[jh]=swap[ih][jh];for(ih=0;ih5、)printf("%0.5f",q[i][j]);printf("");}printf("Givens约减矩阵R的结果如下：");for(i=0;i6、][LS],p[HS][LS];for(i=0;i7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
5、)printf("%0.5f",q[i][j]);printf("");}printf("Givens约减矩阵R的结果如下：");for(i=0;i6、][LS],p[HS][LS];for(i=0;i7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
6、][LS],p[HS][LS];for(i=0;i7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
7、iw][k];em=sqrt(em);u[i]=r[i][k]-em;}elseu[i]=r[i][k];}for(udm=0,ie=k;ie8、;is
8、;is
此文档下载收益归作者所有
