C常用矩阵子函数
Mathcad2001-数学运算-向量和矩阵解读
其 中 的 “ Select a component to insert” 列 表 框 中 列 有 八 个 选 项 ， 常 用 的有Excel(输入Excel文件)、File Read or Write(读入数据文件或输出数据文件) 和Input Table(输入表)。前二项的使用 方法基本相同，用户只需逐次单击 “Next”按钮，便可完成数据文件的输 入。如选择“File Read or Write”项 后单击“Next”按钮，在下一对话框中 选择 “read from a file (读入数据
1.向量和矩阵
(1)创建向量和矩阵
在Mathcad2001中，根据线性代数的习 惯把单个变量称为标量，把包含多个变 量的一列变量称为向量，而把包含多列 的向量称为矩阵，向量和矩阵又合称为 数组。
创建向量和矩阵有下列几种方法：
(1)使用“Insert”菜单中的“Matrix” 命令，或单击“Matrix”工具面板中的 “ ”按钮，或按Ctrl＋M键，将打开 如图24所示的“Insert Matrix”对话框。
求矩阵的逆：若
， C 3 4 2 1 5 6
则
0.141 0.03 9 0.222
C 1 0.162 0.394
0.111
0.111 0.333 0.111
求矩阵各个向量对应的实部和虚部：
1 2
Re(A) 2 6
0 5
Im(B)4Fra bibliotek36 1
8
2
在Mathcad2001中，共提供与向量和矩阵有 关的内置函数39个，其中常用的有：
diag(v)：返回一个对角矩阵，对角元素 为向量v。 geninv(A)：返回矩阵A的逆矩阵。 rref(A)：返回矩阵A的阶形矩阵。 tr(M)：求矩阵斜对角线元素之和(迹)。 rank(A)：求矩阵A的秩。 eigenvec(M,z)：求矩阵M特征值z的特征 向量。 eigenvals(M)：求矩阵M的特征值。
c语言实现矩阵连乘问题
矩阵连乘问题描述：给定n 个矩阵12n {A ,A ,,A } ，其中i A 与1i A +是可乘的，1,2,,1i n =-。
考察这n 个矩阵的连乘积12n A A A 。
由于矩阵乘法满足结合律，所以计算矩阵的连乘可以有许多不同的计算次序。
这种计算次序可以用加括号的方式来确定。
若一个矩阵连乘积的计算次序完全确定，也就是说该连乘积已完全加括号，则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。
如何确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少?代码：#include<stdio.h>const int MAX=100;//静态变量，控制矩阵大小int n;//定义相乘矩阵 的个数int p[MAX+1],m[MAX][MAX],s[MAX][MAX];void matrixChain(){for(int i=1;i<=n;i++) //填主对角线d1{m[i][i]=0;}for(int r=2;r<=n;r++)//填次对角线dr(r=2~n){for(int i=1;i<=n-r+1;i++)//填次对角线的各个元素{int j=i+r-1;//计算次对角线dr 上第i 行的元素的列标 m[i][j]=m[i+1][j]+ p[i-1]*p[i]*p[j];//用计算Ai(Ai+1…Aj)的次数作为m[i][j]的初始值s[i][j]=i;//保存分界点for(int k=i+1;k<j;k++){//用m[i][k]和m[k+1][j]计算m[i][j]的新值 int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];if(t<m[i][j]){m[i][j]=t;s[i][j]=k;}}}}}void tb(int i,int j,int s[MAX][MAX]){//根据s[][]记录的各个子段的最优解，将其输出if(i==j) return;tb(i,s[i][j],s);tb(s[i][j]+1,j,s);printf("Multiply A%d,%dand A%d,%d\n",i,s[i][j],s[i][j]+1,j); }int main(){printf("Please input the number of n...\n");scanf("%d",&n);for(int i=0;i<=n;i++)scanf("%d",&p[i]);matrixChain();//调用函数tb(1,n,s);printf("%d\n",m[1][n]);//输出最优解printf("\n\n\n");for(int i=1;i<=n;i++)//输出矩阵的m[][]的值{for(int j=1;j<=n;j++){printf("%d\t",m[i][j]);}printf("\n");}printf("\n\n\n");for(int i=1;i<=n;i++)//输出矩阵的s[][]的值{for(int j=1;j<=n;j++){printf("%d\t",s[i][j]);}printf("\n");}return 0;}输出结果：。
c语言二次函数拟合,二次函数拟合算法
c语⾔⼆次函数拟合,⼆次函数拟合算法
⼆次函数拟合算法
原理：
在给定⼀组数据序列(x i，y i)，i=0,1,2…m，⽤⼆次多项式拟合这组数据时，设
p(x)=a0+a1x+a2x2，则根据拟合函数与数据序列的均⽅误差最⼩原则，可以得到⼆次多项式函数拟合的矩阵计算⽅程如下所⽰：(
m x i
m
i=1
x i2
m
i=1
x i
m
i=1
x i2
m
i=1
x i3
m
i=1
x i2
m
i=1
x i3
m
i=1
x i4
m
i=1
)(
a0
a1
a2
)= (
y i
m
i=1
x i y i
m
i=1
x i2y i
m
i=1
)
在我们的计算库伦效应实例中，Y即为每个Cycle对应的DischargeC/ChargeC的⽐值，X即
为每个Cycle对应的数字。
代码中定义的三个矩阵XX,AA,YY则分别对应原理公式中等式左边X系数矩阵，A系数矩阵
以及等式右边包含Y系数的矩阵。
具体步骤：
1：先将矩阵中需要的所有量计算出来，并且存放在XX,AA,YY三个矩阵中。
2：为了求得系数矩阵AA，我们需要先把XX矩阵求逆，然后与YY矩阵相乘。
函数MRinv
即为矩阵求逆的函数，返回时存放其逆矩阵。
3：得到系数矩阵AA之后，即得到了拟合好的⼆次函数，将此⼆次函数输出在Excel表中。
具体代码实现：步骤1对应代码：。
C语言常用算法归纳
C语言常用算法归纳应当掌握的一般算法一、基本算法：交换、累加、累乘二、非数值计算常用经典算法：穷举、排序（冒泡，选择）、查找（顺序即线性）三、数值计算常用经典算法：级数计算（直接、简接即递推）、一元非线性方程求根（牛顿迭代法、二分法）、定积分计算（矩形法、梯形法）四、其他：迭代、进制转换、矩阵转置、字符处理（统计、数字串、字母大小写转换、加密等）、整数各数位上数字的获取、辗转相除法求最大公约数（最小公倍数）、求最值、判断素数（各种变形）、数组元素的插入（删除）、二维数组的其他典型问题（方阵的特点、杨辉三角形）详细讲解一、基本算法1．交换（两量交换借助第三者）例1、任意读入两个整数，将二者的值交换后输出。
main(){ int a,b,t;scanf("%d%d",&a,&b);printf("%d,%d\n",a,b);t=a; a=b; b=t;printf("%d,%d\n",a,b);}【解析】程序中加粗部分为算法的核心，如同交换两个杯子里的饮料，必须借助第三个空杯子。
假设输入的值分别为3、7，则第一行输出为3，7；第二行输出为7，3。
其中t为中间变量，起到“空杯子”的作用。
注意：三句赋值语句赋值号左右的各量之间的关系！【应用】例2、任意读入三个整数，然后按从小到大的顺序输出。
main(){ int a,b,c,t;scanf("%d%d%d",&a,&b,&c);/*以下两个if语句使得a中存放的数最小*/if(a>b){ t=a; a=b; b=t; }if(a>c){ t=a; a=c; c=t; }/*以下if语句使得b中存放的数次小*/if(b>c) { t=b; b=c; c=t; }printf("%d,%d,%d\n",a,b,c);}2．累加累加算法的要领是形如“s=s+A”的累加式，此式必须出现在循环中才能被反复执行，从而实现累加功能。
c语言编写指针函数求n阶矩阵matrix最大元素与最小元素的值
当涉及到指针和矩阵操作时，我们可以编写一个函数来求解一个 n 阶矩阵中的最大元素和最小元素的值。
下面是用 C 语言编写的一个示例函数：```c#include <stdio.h>// 定义矩阵的维度#define N 3// 求解矩阵中的最大元素值int findMaxElement(int *matrix, int size) {int max = *matrix;for (int i = 1; i < size; i++) {if (*(matrix + i) > max) {max = *(matrix + i);}}return max;}// 求解矩阵中的最小元素值int findMinElement(int *matrix, int size) {int min = *matrix;for (int i = 1; i < size; i++) {if (*(matrix + i) < min) {min = *(matrix + i);}}return min;}int main() {int matrix[N][N] = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};int *ptr = &matrix[0][0]; // 获取矩阵的指针int max = findMaxElement(ptr, N * N);int min = findMinElement(ptr, N * N);printf("最大元素值为: %d\n", max);printf("最小元素值为: %d\n", min);return 0;}```在这个示例中，我们首先定义了一个 3x3 的矩阵 `matrix`，然后获取了矩阵的指针 `ptr`，并将其传递给 `findMaxElement` 和 `findMinElement` 函数来找到最大元素和最小元素的值。
矩阵分析
k − d +1
. .
1 ak Ck λi ∑k =0 m k −1
∑k =0 ak λi
m
k
p(λi ) p' (λi ) p′′(λi ) / 2! p(λi ) ... = ...
... ... ... ...
p(d −1) (λi ) . . p' (λi ) p(λi ) d×d
p( J 2 )
∑m ak λi k k =0 p(J p(Ji)= ‖ Σk=0makJik
1 ak Ck λi ∑k =0 m
k −1
... ... ... ... ... ... ...
1 (d −1)!
∑k =0 ak λi
m
k
ak Ckd −1λi ∑k =0 .
=PAP (PAP-1)k=PAP-1PAP-1. A=diag(J1,J2,J3)= A2=diag(J
1,J2,J3
. .PAP-1PAP-1=PAkP-1 PAP
J3
J1
J2
)2
=
J1 J2
J1 J2 J 3
0 1 1 1 0 0 1 0 −1
.用初等 .用初等 P -1=
0 1 0 1 −1 1 0 1 − 1
A-2E=
0 0 0 2 0 0 2 2 = 1 −1 1 1 1 1 − 1 −1 3 2 1 −1 1 0 1 1 0 0 1 0 − 1 1 0 = 1 − 1 1 0 1
矩阵运算——C语言实现
#include "stdafx.h"#include <stdlib.h>//结构体typedef struct{double ** mat;int row;int col;}Matrix;void InitialMatrix(Matrix *T, int row,int col); //只分配空间不初始化；void InitialMatrixZero(Matrix *T,int row, int col); //初始化为0void InitialMatrixRand(Matrix *T,int row, int col); //初始化为50以内随机正整数void InputMatrix(Matrix *T); //键盘输入矩阵void DestroyMatrix(Matrix *T); // 释放矩阵空间void PrintfMatrix(Matrix *T); //矩阵输出int AddMatrix(Matrix *A,Matrix *B,Matrix *C); // 矩阵加int MinusMatrix(Matrix *A,Matrix *B,Matrix *C); // 矩阵减int MultiMatrix(Matrix *A,Matrix *B,Matrix *C); //矩阵乘法double MeanMatrix(Matrix *T); //矩阵元素均值int SubMatrix(Matrix *T1,Matrix *T2,int BeginRow,int BeginCol,int EndRow,int EndCol); //求T1的子矩阵T2;测试程序参考界面如下所示：矩阵函数测试，请选择功能，输入对应的数字：***************************************************1：输入一个矩阵，求矩阵均值；2：产生一个随机数矩阵，求矩阵均值；3：输入两个个矩阵，求矩阵和；4：输入两个个矩阵，求矩阵差；5：输入两个个矩阵，求矩阵积；6：产生两个个随机数矩阵，求矩阵和；7：产生两个个随机数矩阵，求矩阵差；8：产生两个个随机数矩阵，求矩阵积；9：求矩阵的子阵，如矩阵的2-4行，1-3列的子阵；0：结束!#include <stdio.h>#include <stdlib.h>//结构体typedef struct {double **mat;int row;int col;} Matrix;//函数声明void InitialMatrix(Matrix *T, int row, int col);void InitialMatrixZero(Matrix *T, int row, int col);void InitialMatrixRand(Matrix *T, int row, int col);void InputMatrix(Matrix *T);void DestroyMatrix(Matrix *T); // 释放矩阵空间void PrintfMatrix(Matrix *T); //矩阵输出int AddMatrix(Matrix *A, Matrix *B, Matrix *C); // 矩阵加int MinusMatrix(Matrix *A, Matrix *B, Matrix *C); // 矩阵减int MultiMatrix(Matrix *A, Matrix *B, Matrix *C); //矩阵乘法double MeanMatrix(Matrix *T); //矩阵元素均值//int SubMatrix(Matrix *T1, Matrix *T2, int BeginRow, int BeginCol, int EndRow, int EndCol); //求T1的子矩阵T2;void NMatrix(void);//求逆矩阵//主函数int main(){int inputevent;//int i, j;int row, col;Matrix T;Matrix A;Matrix B;Matrix C;printf(" 矩阵函数测试，请选择功能，输入对应的数字：\n");printf(" ***************************************************\n\n"); printf("1：输入一个矩阵，求矩阵均值；\n");printf("2：产生一个随机数矩阵，求矩阵均值；\n");printf("3：输入两个个矩阵，求矩阵和；\n");printf("4：输入两个个矩阵，求矩阵差；\n");printf("5：输入两个矩阵，求矩阵积;");printf("\n6：产生两个随机数矩阵，求矩阵和；\n");printf("7：产生两个随机数矩阵，求矩阵差；\n");printf("8：产生两个随机数矩阵，求矩阵积;\n");printf("9:求矩阵的子阵,如矩阵的2-4行D,1-3列的子阵;\n");printf("10：输入一个方阵，求其逆矩阵\n");printf("0:结束!\n");printf("\n\n选择：");scanf("%d", &inputevent);while (inputevent != 0){if (inputevent == 1) {printf("您要输入的矩阵的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);InitialMatrixZero(&T, T.row, T.col);InputMatrix(&T);printf("\n矩阵的平均值为：\n%lf\n\n", MeanMatrix(&T));DestroyMatrix(&T);}if (inputevent == 2) {printf("输入您要产生的矩阵的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);InitialMatrixRand(&T, row, col);MeanMatrix(&T);PrintfMatrix(&T);printf("\n矩阵的平均值为：\n%lf\n\n", MeanMatrix(&T));DestroyMatrix(&T);}if (inputevent == 3) {printf("您要输入的矩阵A的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.col);InputMatrix(&B);InitialMatrix(&C, B.row, B.col);InitialMatrixZero(&C, C.row, C.col);if (AddMatrix(&A, &B, &C) == 1) {printf("\n\n矩阵的和为：A*B\n\n");PrintfMatrix(&C);}else printf("\n\n无法计算！\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 4) {printf("您要输入的矩阵A的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.col);InputMatrix(&B);InitialMatrix(&C, B.row, B.col);InitialMatrixZero(&C, C.row, C.col);if (MinusMatrix(&A, &B, &C) == 1) {printf("\n矩阵的差为：A-B=\n\n");PrintfMatrix(&C);}else printf("\n\n无法计算！\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 5) {int able;printf("您要输入的矩阵A的行数和列数e.g: 5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.col);InputMatrix(&B);InitialMatrix(&C, A.row, B.col);InitialMatrixZero(&C, C.row, C.col);able = MultiMatrix(&A, &B, &C);if (able == 1) {printf("\n积为：A*B\n\n");PrintfMatrix(&C);}else printf("\n\n无法计算！\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 6) {printf("您要产生的矩阵A的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf("A为：\n\n");PrintfMatrix(&A);printf("您要产生的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);printf("B为：\n\n");PrintfMatrix(&B);InitialMatrix(&C, A.row, A.col);InitialMatrixZero(&C, C.row, C.col);if (AddMatrix(&A, &B, &C) == 1) {printf("\n矩阵的和为：A+B=\n\n");PrintfMatrix(&C);}else printf("\n\n无法计算！\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 7) {printf("您要产生的矩阵A的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf("您要产生的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);InitialMatrix(&C, A.row, A.col);InitialMatrixZero(&C, C.row, C.col);if (MinusMatrix(&A, &B, &C) == 1) {printf("\n矩阵的差为：A-B=\n\n");PrintfMatrix(&C);}else printf("\n\n无法计算！\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 8) {printf("您要产生的矩阵A的行数和列数e.g:5,6: \n");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf("A为：\n\n");PrintfMatrix(&A);printf("您要产生的矩阵B的行数和列数e.g:5,6: \n");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);printf("B为：\n\n");PrintfMatrix(&B);InitialMatrix(&C, A.row, A.col);InitialMatrixZero(&C, C.row, C.col);if (MultiMatrix(&A, &B, &C) == 1) {printf("\n积为：A*B=\n\n");PrintfMatrix(&C);;}else printf("\n\n无法计算；\n\n");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);}if (inputevent == 9) printf("对不起，该函数尚在完善中\n\n");if (inputevent == 10) NMatrix();if (inputevent == 0)break;printf(" 矩阵函数测试，请选择功能，输入对应的数字：\n");printf(" ***************************************************\n\n"); printf("1：输入一个矩阵，求矩阵均值；\n");printf("2：产生一个随机数矩阵，求矩阵均值；\n");printf("3：输入两个个矩阵，求矩阵和；\n");printf("4：输入两个个矩阵，求矩阵差；\n");printf("5：输入两个矩阵，求矩阵积;");printf("\n6：产生两个随机数矩阵，求矩阵和；\n");printf("7：产生两个随机数矩阵，求矩阵差；\n");printf("8：产生两个随机数矩阵，求矩阵积;\n");printf("9:求矩阵的子阵,如矩阵的2-4行D,1-3列的子阵;\n");printf("10：输入一个方阵，求其逆矩阵\n");printf("0:结束!\n");printf("\n\n选择：");scanf("%d", &inputevent);}return 0;}//其他函数void InitialMatrix(Matrix *T, int row, int col){//printf("分配内存中......\n");int i;int succ = 1;//T=(Matrix *)malloc(sizeof(Matrix));T->row = row;T->col = col;T->mat = (double **)malloc(T->row * sizeof(double *));if (T->mat == NULL) {succ = 0;}else {for (i = 0; i < T->row; i++) {T->mat[i] = (double *)malloc(T->col * sizeof(double));if (T->mat[i] == NULL) {succ = 0;break;}}//if(succ==1)// printf("内存分配成功|；?\n");//else printf("内存分配失败；\n");}}void InitialMatrixZero(Matrix *T, int row, int col) { //printf("矩阵初始化为零中......\n");int i, j;for (i = 0; i < row; i++)for (j = 0; j < col; j++)T->mat[i][j] = 0;//printf("矩阵初始化为零矩阵成功；\n"); }void InitialMatrixRand(Matrix *T, int row, int col) { int i, j;for (i = 0; i < row; i++)for (j = 0; j < col; j++)(*T).mat[i][j] = rand() % 50;}void InputMatrix(Matrix *T) {printf("输入矩阵：\n");int i, j;for (i = 0; i < (*T).row; i++)for (j = 0; j < (*T).col; j++)scanf("%lf", &(*T).mat[i][j]);}void DestroyMatrix(Matrix *T){int i;for (i = 0; i < (*T).row; i++)free((*T).mat[i]);}void PrintfMatrix(Matrix *T){int i, j;for (i = 0; i < (*T).row; i++){for (j = 0; j < (*T).col; j++)printf("%lf ", (*T).mat[i][j]);printf("\n");}}int AddMatrix(Matrix *A, Matrix *B, Matrix *C){int i, j;if ((*A).row == (*B).row && (*A).col == (*B).col){for (i = 0; i < (*A).row; i++)for (j = 0; j < (*A).col; j++)(*C).mat[i][j] = (*A).mat[i][j] + (*B).mat[i][j];for (i = 0; i < (*A).row; i++)for (j = 0; j < (*A).col; j++)return 1;}else {printf("这两个矩阵不能相加！\n");return 0;}}int MinusMatrix(Matrix *A, Matrix *B, Matrix *C){int i, j;if ((*A).row == (*B).row && (*A).col == (*B).col){for (i = 0; i < (*A).row; i++)for (j = 0; j < (*A).col; j++)(*C).mat[i][j] = (*A).mat[i][j] - (*B).mat[i][j];return 1;}elseprintf("这两个矩阵不能相减！\n");return 0;}int MultiMatrix(Matrix *A, Matrix *B, Matrix *C){int i=0, j=0;int k = 0;if ((*A).col == (*B).row){for (i = 0; i < (*A).row; i++) {for (j = 0; j < (*B).col; j++)for(k=0;k <(A->col);k++)(*C).mat[i][j] += (*A).mat[i][k] * (*B).mat[k][j];}return 1;}elseprintf("这两个矩阵不能相乘！\n");return 0;}double MeanMatrix(Matrix *T) //矩阵元素均值{double mean;double sum = 0;int i, j;for (i = 0; i < (*T).row; i++)for (j = 0; j < (*T).col; j++)sum += (*T).mat[i][j];mean = sum / ((*T).row*(*T).col);return mean;}void NMatrix(void){#define M 20#define N 20int i,j,m,n;float y=1.0;float a[M][2 * M];float b[N][2 * N];float t, x;int k, T;printf("输入方阵的维数：\n"); //请输入方阵，即行和列相等的矩阵。
C语言计算逆矩阵
C语⾔计算逆矩阵花了4天写的，不过三天在重学线代。
1 #include<stdio.h>2 #include<stdlib.h> // 操作内存3 #include<math.h> // pow()函数,计算-1的n次⽅，可以不⽤这个函数，偷懒使⽤现成的45/*6显⽰矩阵7 matrix: 矩阵8 order: 阶数9*/10void showMatrix(float** matrix, int order)11 {12for (int i = 0; i < order; i++) {13for (int j = 0; j < order; j++) {14 printf(" %f ", matrix[i][j]);15 }16 printf("|\n");17 }18 }1920/*21交换两⾏22⼀开始想使⽤初等⾏变换计算，此函数可以删除，⽤不到23 x1:调换第⼀⾏24 x2:调换第⼆⾏25 order:矩阵阶数26 matrix:矩阵27*/28void replaceRow(int x1, int x2, int order, float **matrix)29 {30float temp;31 x1 -= 1;32 x2 -= 1;33for (int i = 0; i < order; i++) {34 temp = matrix[x1][i];35 matrix[x1][i] = matrix[x2][i];36 matrix[x2][i] = temp;37 }38 }3940/*41转置矩阵42 matrix: 矩阵43 order: 阶数44*/45void transposeMatrix(float** matrix, int order)46 {47float temp;48for (int i = 0; i < order; i++) {49for (int j = 0; j < i; j++) {50 temp = matrix[i][j];51 matrix[i][j] = matrix[j][i];52 matrix[j][i] = temp;53 }54 }55 }5657/*58获取除了某⾏某列的矩阵59 oldmatrix: 原本矩阵60 newmatrix: 新矩阵61 row: 要删除⾏62 col: 要删除列63 order: 阶数64*/65void get(float** oldmatrix, float** newmatrix, int row, int col, int order)66 {67// 删除了⼀⾏⼀列，所以新矩阵⾏列均⽐原矩阵少168int a = 0, b = 0;69int x, y, z = 0, w = 0;70// i,j循环原矩阵71for (int i = 0; i < order - 1; i++) {72for (int j = 0; j < order - 1; j++) {73// z,w代表⾏列的是否加⼀状态，防⽌多次加⼀，+1只需要1次74if (i >= row && z == 0) { a += 1; z = 1; }75if (j >= col && w == 0) { b += 1; w = 1; }7677 newmatrix[i][j] = oldmatrix[i+a][j+b];78 }79 a = 0;b = 0;80 z = 0;w = 0;81 }82 }8384/*85计算⾏列式86 matrix: 矩阵87 order: 阶数88*/89float calc(float** matrix, int order)90 {91// 递归求⾏列式值92float num=0;93int i, j;94if (order == 2) {95// 如果是⼆阶直接获取值96 num = matrix[0][0] * matrix[1][1] - matrix[0][1] * matrix[1][0];97 }98else {99// 创建更⼩⼆维数组100 order -= 1;101float** matrixlow = (float**)malloc(sizeof(float*) * order);102for (i = 0; i < order; i++) {103 matrixlow[i] = (float*)malloc(sizeof(float) * order);104 }105 order += 1;106107// 循环展开第⼀⾏108for (j = 0; j < order; j++) {109110get(matrix, matrixlow, 0, j, order);111// 此处开始递归，调⽤⾃⾝函数112 num += matrix[0][j] * pow(-1, j) * calc(matrixlow, order - 1); 113 }114115// 释放内存116for (i = 0; i < order-1; ++i)free(*(matrixlow + i));117 }118119return num;120 }121122123/*124主函数125*/126int main()127 {128int order; // 矩阵阶数129int i, j;130float det;131132// 获取矩阵阶133 printf("输⼊矩阵阶：");134 scanf("%d", &order);135136 printf("输⼊的阶是：%d\n\n", order);137138// 申请⼆维数组内存139float** matrix = (float**)malloc(sizeof(float*) * order);140for (i = 0; i < order; i++) {141 matrix[i] = (float*)malloc(sizeof(float) * order);142 }143144// 获取输⼊145for (i = 0; i < order; i++) {146for (j = 0; j < order; j++) {147 printf("位置：( %d , %d)请输⼊数据：",i + 1,j + 1);148 scanf("%f", &matrix[i][j]);149 }150 }151152// 计算并显⽰det153 det = calc(matrix, order);154 printf("\ndet值为：%f",det);155// 0不能做除数156if (det == 0)157 printf("\n矩阵接近奇异值，结果可能不准确！");158 printf("\n\n");159160// 显⽰输⼊矩阵161 printf("\n输⼊的矩阵是：\n\n");162 showMatrix(matrix, order);163164// 申请⼆维数组存储结果165float** Rmatrix = (float**)malloc(sizeof(float*) * order);166for (i = 0; i < order; i++) {167 Rmatrix[i] = (float*)malloc(sizeof(float) * order);168 }169170171// 开始计算172if (order == 2) {173// 阶数为⼆直接运⾏公式174float n = 1 / (matrix[0][0] * matrix[1][1] - matrix[0][1] * matrix[1][0]);175 Rmatrix[0][0] = n * matrix[1][1];176 Rmatrix[0][1] = -n * matrix[0][1];177 Rmatrix[1][0] = -n * matrix[1][0];178 Rmatrix[1][1] = n * matrix[0][0];179 }180else {181// 转置矩阵并显⽰182 transposeMatrix(matrix, order);183 printf("\n\n转置后为:\n\n");184 showMatrix(matrix, order);185186// 循环求i，j位的代数余⼦式187for (i = 0; i < order; i++) {188for (j = 0; j < order; j++) {189// 申请⼆维数组190 order -= 1;191float** matrixlow = (float**)malloc(sizeof(float*) * order);192for (int t = 0; t < order; t++) {193 matrixlow[t] = (float*)malloc(sizeof(float) * order);194 }195 order += 1;196197// 获取除了i，j⾏的值组成⾏列式198get(matrix, matrixlow, i, j, order);199// 计算⾏列式值除以det200 Rmatrix[i][j] = pow(-1, i + j) * calc(matrixlow, order - 1) / det; 201202// 释放内存203for (int t = 0; t < order-1; ++t)free(*(matrixlow + t));204 }205 }206207 }208209// 显⽰逆矩阵210 printf("\n\n逆矩阵为：\n\n");211 showMatrix(Rmatrix, order);212213214//// 释放⼆维数组215for (i = 0; i < order; ++i)free(*(matrix + i));216for (i = 0; i < order; ++i)free(*(Rmatrix + i));217218return0;219 }。