C語(yǔ)言解3元1次方程組 用初中學(xué)的最基本的聯(lián)合消元法
沒學(xué)過(guò)線性代數(shù)，但是很多算法都和矩陣相關(guān)，所以就硬著頭皮學(xué)。
最近就想自己能不能先寫個(gè)算線性方程組的程序呢？后來(lái)就想了這么個(gè)方法，暫時(shí)只能算3元的，任意元的接下來(lái)繼續(xù)想。有太多硬編碼，希望有興趣的讀者可以給點(diǎn)修改建議！
#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的第一個(gè)方程分別與2、3個(gè)方程聯(lián)合消去z，得到x，y的兩個(gè)方程，存入B
double B[2][4];
//B的兩個(gè)方程聯(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個(gè)參數(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)文章
C語(yǔ)言學(xué)習(xí)之函數(shù)知識(shí)總結(jié)
函數(shù)是一組一起執(zhí)行一個(gè)任務(wù)的語(yǔ)句。每個(gè)?C?程序都至少有一個(gè)函數(shù)，即主函數(shù)?main()?，所有簡(jiǎn)單的程序都可以定義其他額外的函數(shù)。本文就為大家詳細(xì)講講C語(yǔ)言中函數(shù)的相關(guān)知識(shí)點(diǎn)，希望有所幫助2022-07-07
這篇文章主要介紹了論C++的lambda是函數(shù)還是對(duì)象，對(duì)于有捕獲的lambda，其等價(jià)于對(duì)象。對(duì)于沒有任何捕獲的lambda，其等價(jià)于函數(shù)，下面來(lái)看看具體的相關(guān)內(nèi)容,需要的朋友可以參考一下2022-02-02
getchar()是在輸入緩沖區(qū)順序讀入一個(gè)字符(包括空格、回車和Tab) ，getchar()是stdio.h中的庫(kù)函數(shù)，它的作用是從stdin流中讀入一個(gè)字符，本文給大家介紹getchar的使用，感興趣的朋友跟隨小編一起看看吧2022-12-12
今天小編就為大家分享一篇opencv3/C++ FLANN特征匹配方式，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-12-12
C++如何將一個(gè)vector內(nèi)容賦值給另一個(gè)vector,及swap與assign區(qū)別
在本文中,我們將主要介紹5種將一個(gè)vector內(nèi)容賦值給另一個(gè)vector的方式,順便討論下swap與assign的區(qū)別,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2023-08-08
我們常用的x86結(jié)構(gòu)都是小端模式，而大部分DSP，ARM也是小端模式，不過(guò)有些ARM是可以選擇大小端模式。所以對(duì)于上面的maxHeight是應(yīng)該以小端模式來(lái)存放，具體情況請(qǐng)看下面兩表2013-10-10
這篇文章主要介紹了C++設(shè)計(jì)模式之狀態(tài)模式,本文講解了什么是狀態(tài)模式、狀態(tài)模式的使用場(chǎng)合、狀態(tài)模式的實(shí)現(xiàn)代碼等內(nèi)容,需要的朋友可以參考下2014-10-10
C++在多線程中使用condition_variable實(shí)現(xiàn)wait
這篇文章主要介紹了C++中的condition_variable中在多線程中的使用，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)吧2022-09-09
c++之std::get_time和std::put_time
std::get_time和std::put_time是C++中用于日期和時(shí)間的格式化和解析的函數(shù),它們都包含在<iomanip>頭文件中,std::get_time用于從輸入流中解析日期時(shí)間字符串,而std::put_time則用于將std::tm結(jié)構(gòu)格式化為字符串2024-10-10
C,語(yǔ)言常用的排序方法有很多種。比如說(shuō)冒泡排序,直接交換排序,直接選擇排序,直接插入排序,二分插入排序,快速排序,歸并排序等等,下面這篇文章主要給大家介紹了關(guān)于c語(yǔ)言實(shí)現(xiàn)幾種常用的排序算法,需要的朋友可以參考下2021-06-06