C语言实现矩阵高斯消元法求解线性方程组

float** calculate2(int a, float** arr0,float**arr1)
{
	int i, j;
	// 将arr0的所有元素初始化为0
	for (i = 0; i < a; i++)
	{
		for (j = 0; j < 2*a; j++)
		{
			arr0[i][j] = 0;
		}
	}
	// 将arr1的值复制到arr0的相应位置，并在对角线上的位置添加1
	for (i = 0; i < a; i++)
	{
		for (j = 0; j < a; j++)
		{
			arr0[i][j] = arr1[i][j];
			if (i == j )
			{
				arr0[i][j + a] = 1;
			}
		}
	}
	// 测试用
	i = 0;
	cout << 'error';
	// 高斯消元法主循环
	for (j = 0; j < a; j++, i++)
	{
		int a = i;
		// 如果当前列的首个元素为0，则向下寻找非零元素并交换行
		if (arr0[i][j] == 0)
		{
			while (arr0[i][j] == 0)
				i++;
			// 交换两行的值
			int j2 = 0;
			for (; j2 < 2*a; j2++)
			{
				int tmp = arr0[a][j2];
				arr0[a][j2] = arr0[i][j2];
				arr0[i][j2] = tmp;
			}
			i = a;//更新i的值；
		}
		i++;
		// 将当前列的下方元素消为0
		for (; i < a; i++)
		{
			if (arr0[i][j] != 0)//不为0的情况
			{
				cout << 'error';
				float bi = arr0[i][j] / arr0[a][j];
				//变化该行
				int j1 = 0;
				for (; j1 < 2*a; j1++)
				{
					arr0[i][j1] -= arr0[a][j1] * bi;
				}
				//测试结果:成功			
			}
		}
		i = a;
	}
	outputjz(a,2*a,arr0);
	//测试用
	return NULL;

}

这段代码实现了一个名为'calculate2'的函数，该函数接受三个参数：

• a：矩阵的行数和列数
• arr0：一个二维浮点数组，用于存储计算结果
• arr1：一个二维浮点数组，存储输入的矩阵

该函数使用高斯消元法将输入的矩阵转换为上三角矩阵，并将结果存储在arr0中。函数中调用了一个名为outputjz的函数，用于输出结果，但该函数的具体实现未在代码中给出。

注意：

• 这段代码中包含一些测试代码，例如cout << 'error';，这些代码在实际使用中可以删除。
• 代码中缺少对outputjz函数的定义，因此无法确定其功能和输出方式。
• 这段代码仅供参考，实际使用时可能需要根据具体情况进行修改。