C语言解3元1次方程组 用初中学的最基本的联合消元法
没学过线性代数，但是很多算法都和矩阵相关，所以就硬着头皮学。
最近就想自己能不能先写个算线性方程组的程序呢？后来就想了这么个方法，暂时只能算3元的，任意元的接下来继续想。有太多硬编码，希望有兴趣的读者可以给点修改建议！
#include "stdafx.h"//VS2010需要
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
double x[3];//存放解x,y,z，使用数组便于输出
//x=1,y=2,z=3
//double A[3][4]=
//{
// 1,1,1,-6,
// 1,1,-1,0,
// 1,-1,1,-2
//};
//x=2,y=5,z=10
double A[3][4]=
{
1,2,3,-42,
2,-1,5,-49,
-1,3,-3,17
};
//A的第一个方程分别与2、3个方程联合消去z，得到x，y的两个方程，存入B
double B[2][4];
//B的两个方程联合消去y，得到x
double C[1][4];
void CombineB(int i)//A[0]与A[1]或A[2]消去z
{
double m0=abs(A[i][2]);
int b=A[0][2]*A[i][2]>0?-1:1;
double mi=b*abs(A[0][2]);
//printf("m0=%lf,mi=%lf\n",m0,mi);//调试用
for(int j=0;j<4;j++)
{
B[i-1][j]=A[0][j]*m0+A[i][j]*mi;
}
}
void CombineC()//B[0]与B[1]消去y
{
double m0=abs(B[1][1]);
int b=B[0][1]*B[1][1]>0?-1:1;
double m1=b*abs(B[0][1]);
for(int j=0;j<4;j++)
{
C[0][j]=B[0][j]*m0+B[1][j]*m1;
}
}
//显示4个参数
void Show(double N[][4],int d1)
{
char W[3]={'x','y','z'};
for(int i=0;i<d1;i++)
{
for(int j=0;j<3;j++)
{
printf("%lf * %c + ",N[i][j],W[j]);
}
printf("%lf = 0\n",N[i][3]);
}
printf("\n\n");
}
int _tmain(int argc, _TCHAR* argv[])
{
Show(A,3);
CombineB(1);
CombineB(2);
Show(B,2);
CombineC();
Show(C,1);
//得到x后则依次算出其他数字
x[0]=(0-C[0][3])/C[0][0];
x[1]=(0-B[0][3]-B[0][0]*x[0])/B[0][1];
x[2]=(0-A[0][3]-A[0][0]*x[0]-A[0][1]*x[1])/A[0][2];
//输出结果
for(int i=0;i<3;i++)
printf("x[%d]=%lf\t",i,x[i]);
printf("\n");
system("pause");
return 0;
}
2.截图 
相关文章
函数是一组一起执行一个任务的语句。每个 C 程序都至少有一个函数，即主函数 main() ，所有简单的程序都可以定义其他额外的函数。本文就为大家详细讲讲C语言中函数的相关知识点，希望有所帮助2022-07-07
这篇文章主要介绍了论C++的lambda是函数还是对象，对于有捕获的lambda，其等价于对象。对于没有任何捕获的lambda，其等价于函数，下面来看看具体的相关内容,需要的朋友可以参考一下2022-02-02
getchar()是在输入缓冲区顺序读入一个字符(包括空格、回车和Tab) ，getchar()是stdio.h中的库函数，它的作用是从stdin流中读入一个字符，本文给大家介绍getchar的使用，感兴趣的朋友跟随小编一起看看吧2022-12-12
今天小编就为大家分享一篇opencv3/C++ FLANN特征匹配方式，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-12-12
C++如何将一个vector内容赋值给另一个vector,及swap与assign区别
在本文中,我们将主要介绍5种将一个vector内容赋值给另一个vector的方式,顺便讨论下swap与assign的区别,如有错误或未考虑完全的地方,望不吝赐教2023-08-08
我们常用的x86结构都是小端模式，而大部分DSP，ARM也是小端模式，不过有些ARM是可以选择大小端模式。所以对于上面的maxHeight是应该以小端模式来存放，具体情况请看下面两表2013-10-10
这篇文章主要介绍了C++设计模式之状态模式,本文讲解了什么是状态模式、状态模式的使用场合、状态模式的实现代码等内容,需要的朋友可以参考下2014-10-10
C++在多线程中使用condition_variable实现wait
这篇文章主要介绍了C++中的condition_variable中在多线程中的使用，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习吧2022-09-09
c++之std::get_time和std::put_time
std::get_time和std::put_time是C++中用于日期和时间的格式化和解析的函数,它们都包含在<iomanip>头文件中,std::get_time用于从输入流中解析日期时间字符串,而std::put_time则用于将std::tm结构格式化为字符串2024-10-10
C,语言常用的排序方法有很多种。比如说冒泡排序,直接交换排序,直接选择排序,直接插入排序,二分插入排序,快速排序,归并排序等等,下面这篇文章主要给大家介绍了关于c语言实现几种常用的排序算法,需要的朋友可以参考下2021-06-06
最新评论