C語言解3元1次方程組 用初中學(xué)的最基本的聯(lián)合消元法
沒學(xué)過線性代數(shù)，但是很多算法都和矩陣相關(guān)，所以就硬著頭皮學(xué)。
最近就想自己能不能先寫個算線性方程組的程序呢？后來就想了這么個方法，暫時只能算3元的，任意元的接下來繼續(xù)想。有太多硬編碼，希望有興趣的讀者可以給點修改建議！
#include "stdafx.h"http://VS2010需要
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
double x[3];//存放解x,y,z，使用數(shù)組便于輸出
//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個方程聯(lián)合消去z，得到x，y的兩個方程，存入B
double B[2][4];
//B的兩個方程聯(lián)合消去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);//調(diào)試用
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個參數(shù)
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后則依次算出其他數(shù)字
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];
//輸出結(jié)果
for(int i=0;i<3;i++)
printf("x[%d]=%lf\t",i,x[i]);
printf("\n");
system("pause");
return 0;
}
2.截圖 
相關(guān)文章
Qt 使用QDialog實現(xiàn)界面遮罩的示例(蒙版)
界面遮罩在很多時候都可以用到，例如彈窗，本文主要介紹了Qt 使用QDialog實現(xiàn)界面遮罩的示例，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2022-04-04
這篇文章主要介紹了C++設(shè)計模式之建造者模式Builder的相關(guān)資料，具有一定的參考價值，感興趣的小伙伴們可以參考一下2018-03-03
這篇文章主要介紹了Qt如何實現(xiàn)在窗口上繪制任意大小的矩形，并且通過邊角的拖曳按鈕可改變矩形大小，感興趣的小伙伴可以跟隨小編一起學(xué)習(xí)一下2022-06-06
C語言實現(xiàn)數(shù)據(jù)結(jié)構(gòu)和雙向鏈表操作
這篇文章主要介紹了C語言實現(xiàn)數(shù)據(jù)結(jié)構(gòu)雙向鏈表操作,需要的朋友可以參考下2017-03-03
這篇文章主要介紹了C++11強類型枚舉的相關(guān)資料，幫助大家更好的理解和學(xué)習(xí)c++11，感興趣的朋友可以了解下2020-08-08
這篇文章主要為大家詳細(xì)介紹了DOS簡易版C語言貪吃蛇，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2019-11-11
c/c++中，函數(shù)可以傳遞的參數(shù)有三種形式,值、引用和指針。以下分別對這三種形式進行了介紹，需要的朋友可以過來參考下2013-07-07
這篇文章主要為大家詳細(xì)介紹了opencv實現(xiàn)圖形輪廓檢測，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2021-04-04
這篇文章主要介紹了C++內(nèi)聯(lián)函數(shù)，文章主要圍繞C++內(nèi)聯(lián)函數(shù)的相關(guān)資料展開詳細(xì)內(nèi)容,需要的朋友可以參考一下，希望對大家有所幫助2021-11-11
C++探索構(gòu)造函數(shù)私有化會產(chǎn)生什么結(jié)果
C++的構(gòu)造函數(shù)的作?：初始化類對象的數(shù)據(jù)成員。即類的對象被創(chuàng)建的時候，編譯系統(tǒng)對該對象分配內(nèi)存空間，并?動調(diào)?構(gòu)造函數(shù)，完成類成員的初始化。構(gòu)造函數(shù)的特點：以類名作為函數(shù)名，?返回類型2022-05-05