原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
#include stdio.h
#include malloc.h
#include math.h
void ScanMatrix(float *arr,int m,int n); //输入矩阵
void matrixmul(float*arr1,float*arr2,float*arr3,int m,int n,int k); //矩阵数乘函数
void printm(char c,float*arr,int m,int n); //矩阵打印函数
void PALU(float *arr,int m,int n); //PA=LU分解函数
void Gram_Schmint(float *arr,int m,int n); //QR分解函数
void Household(float *arr,int m,int n); //household分解函数
void Givens(float *arr,int m,int n); //givens分解函数
void divide(); //分解矩阵总函数
int main()
{
divide();
return 0;
}
void divide()
{
int m,n,method,flag=1;
char c=0;
while(c!=#)
{
printf(请输入mXn矩阵维度m:);
scanf(%d,m);
printf(请输入mXn矩阵维度n:);
scanf(%d,n);
float *arr= (float *)malloc(sizeof (float) * m*n);
ScanMatrix(arr,m,n);
while(flag)
{
printf(请输入矩阵的分解方式:\n1-PA=LU,\n2-Gram-Schmit正交化A=QR,\n3-Household分解A=QU,\n4-Givens分解A=QU。\n5-结束分解.\n);
scanf(%d,method);
switch(method)
{
case 1:
PALU(arr,m,n);
break;
case 2:
Gram_Schmint(arr,m,n);
break;
case 3:
Household(arr,m,n);
break;
case 4:
Givens(arr,m,n);
break;
case 5:
flag=0;
break;
}
}
flag=1;
printf(输入#结束，否则继续重新输入矩阵开始分解:);
getchar();
c=getchar();
free(arr);
}
}
void PALU(float *arr,int m,int n)
{
int i,j,k,p,q,exsit=1,change1=0,change2=0;
float t,s;
float *arrt= (float *)malloc(sizeof (float)*m*n);
float *P= (float *)malloc(sizeof (float) * m*m);
float *L= (float *)malloc(sizeof (float) * m*m);
for(i=0;im;i++)
for(j=0;jn;j++)
{
*(arrt+i*n+j)=*(arr+i*n+j);
}
for(i=0;im;i++)
for(j=0;jm;j++)
{
*(P+i*m+j)=*(L+i*m+j)=0;
}
for(i=0;im;i++)
{
*(P+i*m+i)=*(L+i*m+i)=1;
}
for(i=0;iminexsit==1;i++)
{
if(*(arrt+i*n+i)==0) //判断主元是否为0；
{
for(j=i+1;jm+1;j++)//对主元下面的元素进行考察；
{
if(*(ar
文档评论（0）