C++实现复数矩阵求逆 matlab inv

一、引言

之前偶尔一次有用到将matlab转为C++语言的需求，其中matlab有一个inv函数可以非常方便的求矩阵的逆，甚至是复数矩阵。而C++中没有类似的函数。在csdn上有一个matlab2c的库的博客（github地址第80行开始），但是只有实数矩阵求逆的代码，而我又在百度上搜到一篇文献中写到将复数矩阵转为两个实数矩阵然后在进行求解。然后我就将两个代码一结合，就实现了复数矩阵求逆。经验证与matlab的inv函数结果一致。

二、原理

2.1 实数矩阵求逆

这里直接给出matlab2c中求实数矩阵的代码,其中函数inv的两个参数：a代表实数矩阵，用double二维数组存储、num代表这个矩阵的行数/列数。

void swap(double* a, double* b);  //声明子程序
double** inv(double** a, int num)
{
	int* is, * js, i, j, k;
	int n = num;
	double temp, fmax;
	double** tp = new double* [num];
	for (int i = 0; i < num; i++) tp[i] = new double[num];
	for (i = 0; i < n; i++) {
		for (j = 0; j < n; j++) {
			tp[i][j] = a[i][j];
		}
	}
	is = new int[n];
	js = new int[n];
	for (k = 0; k < n; k++)
	{
		fmax = 0.0;
		for (i = k; i < n; i++) {
			for (j = k; j < n; j++)
			{
				temp = fabs(tp[i][j]);//找最大值
				if (temp > fmax)
				{
					fmax = temp;
					is[k] = i; js[k] = j;
				}
			}
		}

		if ((fmax + 1.0) == 1.0)
		{
			delete[] is;
			delete[] js;
			return NULL;
		}
		if ((i = is[k]) != k)
			for (j = 0; j < n; j++)
				swap(&tp[k][j], &tp[i][j]);//交换指针
		if ((j = js[k]) != k)
			for (i = 0; i < n; i++)
				swap(&tp[i][k], &tp[i][j]);  //交换指针
		tp[k][k] = 1.0 / tp[k][k];
		for (j = 0; j < n; j++)
			if (j != k)
				tp[k][j] *= tp[k][k];
		for (i = 0; i < n; i++)
			if (i != k)
				for (j = 0; j < n; j++)
					if (j != k)
						tp[i][j] = tp[i][j] - tp[i][k] * tp[k][j];
		for (i = 0; i < n; i++)
			if (i != k)
				tp[i][k] *= -tp[k][k];
	}
	for (k = n - 1; k >= 0; k--)
	{
		if ((j = js[k]) != k)
			for (i = 0; i < n; i++)
				swap(&tp[j][i], &tp[k][i]);
		if ((i = is[k]) != k)
			for (j = 0; j < n; j++)
				swap(&tp[j][i], &tp[j][k]);
	}
	delete[] is;
	delete[] js;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
		}
	}
	return tp;
}
void swap(double* a, double* b)
{
	double c;
	c = *a;
	*a = *b;
	*b = c;
}

2.2 复数矩阵求逆

在那篇文献中，我们可以看到将复数矩阵C拆分为实部矩阵A+虚部矩阵Bi，则