欢迎来到天天文库
浏览记录
ID:38268832
大小:87.25 KB
页数:5页
时间:2019-05-25
《列选主元矩阵求逆C源程序》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、列选主元矩阵求逆C源程序列选主元矩阵求逆C源程序/**Copyright(C)2011,北京航天时代激光导航技术有限责任公司**Allrightsreserved.**程序功能:求矩阵的逆矩阵***当前版本:V1.0**作者:谢页华**手机:18611242620**邮箱:hitboy_0@163.com**完成日期:2011-08-06*/#include"math.h"#include"stdio.h"#include"stdlib.h"intMatrix_Athwart(intn,double*A){in
2、ti,j,k;doubletemp,*B;B=(double*)malloc(sizeof(double)*n*n);for(i=0;i3、;}}if(temp==0)return(0);//MessageBox("矩阵不可逆。");if(k!=i){for(j=i;j4、k)-(*(A+j*n+i))*(*(A+i*n+k))/(*(A+i*n+i));}for(k=0;k=0;i--){for(k=i;k>=0;k--){for(j=0;j5、+j)=*(B+k*n+j)-(*(A+k*n+i+1))/(*(A+(i+1)*n+i+1))*(*(B+(i+1)*n+j));}}}for(i=0;i6、tn,double*Matrix_A){inti,j;if(Matrix_Athwart(n,Matrix_A)!=1){printf("矩阵不可逆!");列选主元矩阵求逆C源程序}else{printf("矩阵求逆的结果是:");for(i=0;i7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
3、;}}if(temp==0)return(0);//MessageBox("矩阵不可逆。");if(k!=i){for(j=i;j4、k)-(*(A+j*n+i))*(*(A+i*n+k))/(*(A+i*n+i));}for(k=0;k=0;i--){for(k=i;k>=0;k--){for(j=0;j5、+j)=*(B+k*n+j)-(*(A+k*n+i+1))/(*(A+(i+1)*n+i+1))*(*(B+(i+1)*n+j));}}}for(i=0;i6、tn,double*Matrix_A){inti,j;if(Matrix_Athwart(n,Matrix_A)!=1){printf("矩阵不可逆!");列选主元矩阵求逆C源程序}else{printf("矩阵求逆的结果是:");for(i=0;i7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
4、k)-(*(A+j*n+i))*(*(A+i*n+k))/(*(A+i*n+i));}for(k=0;k=0;i--){for(k=i;k>=0;k--){for(j=0;j5、+j)=*(B+k*n+j)-(*(A+k*n+i+1))/(*(A+(i+1)*n+i+1))*(*(B+(i+1)*n+j));}}}for(i=0;i6、tn,double*Matrix_A){inti,j;if(Matrix_Athwart(n,Matrix_A)!=1){printf("矩阵不可逆!");列选主元矩阵求逆C源程序}else{printf("矩阵求逆的结果是:");for(i=0;i7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
5、+j)=*(B+k*n+j)-(*(A+k*n+i+1))/(*(A+(i+1)*n+i+1))*(*(B+(i+1)*n+j));}}}for(i=0;i6、tn,double*Matrix_A){inti,j;if(Matrix_Athwart(n,Matrix_A)!=1){printf("矩阵不可逆!");列选主元矩阵求逆C源程序}else{printf("矩阵求逆的结果是:");for(i=0;i7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
6、tn,double*Matrix_A){inti,j;if(Matrix_Athwart(n,Matrix_A)!=1){printf("矩阵不可逆!");列选主元矩阵求逆C源程序}else{printf("矩阵求逆的结果是:");for(i=0;i7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
7、rix_A;Matrix_A=(double*)malloc(sizeof(double)*n*n);printf("请输入待求逆矩阵:");for(i=0;i8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
8、atrix_Result_Output(n,Parameter_Input(n));return(0);}
此文档下载收益归作者所有