1、C語(yǔ)言解線性方程組一、實(shí)驗(yàn)名稱(chēng)：用C語(yǔ)言解線性方程組。二、實(shí)驗(yàn)原理：用初等行變化的方法將線性方程組的增廣矩陣化簡(jiǎn)成最簡(jiǎn)型；從最簡(jiǎn)型中可以看出線性方程組解的情況。三、實(shí)驗(yàn)工具：計(jì)算機(jī)和MicrosoftVisualC+6.0軟件。四、實(shí)驗(yàn)原程序代碼如下：#include<stdio.h>#include<math.h>#include<stdlib.h>#defineMAXR20#defineMAXC40floatMMAXC;floatM_BMAXRMAXC;voidf1(intm,intn)inti,j,r=0,c=0;floatx0;printf(&quo

2、t;n以下過(guò)程是將矩陣化為行階梯型n");printf("n輸出初始增廣矩陣!n");for(i=0;i<m;i+)for(j=0;j<n;j+)printf("%8.2f",M_Bij);/*原始增廣矩陣*/printf("n");printf("按任一健繼續(xù).n");getchar();while(c<n)for(i=r;i<m;i+)if(fabs(M_Bic)>=0.0001)break;if(i<m)if(i!=r)for(j=0;j<n;j+)x0=M

3、_Bij;M_Bij=M_Brj;M_Brj=x0;x0=M_Brc;for(j=0;j<n;j+)M_Brj/=x0;for(i=r+1;i<m;i+)x0=M_Bic;for(j=0;j<n;j+)M_Bij-=x0*M_Brj;r+;for(i=0;i<m;i+)for(j=0;j<n;j+)printf("%8.2f",M_Bij);printf("n");printf("按任一健繼續(xù).n");getchar();c+;if(r=m)break;printf("n矩陣的行階梯型如上!n&

4、quot;);printf("n以下過(guò)程是將行階梯型化為行最簡(jiǎn)型n");while(r)r-;for(j=0;j<n-1;j+)if(fabs(M_Brj)>=0.0001)break;c=j;for(i=0;i<r;i+)x0=M_Bic;for(j=0;j<n;j+)M_Bij-=x0*M_Brj;for(i=0;i<m;i+)for(j=0;j<n;j+)printf("%8.2f",M_Bij);printf("n");printf("按任一健繼續(xù).n");getchar

5、();floatMMAXC=0;intl1,l2;for(l1=0;l1<m;l1+)for(l2=0;l2<n;l2+)if(M_Bl1l2!=0)Ml2=1;break;if(Mn-1=1)printf("該方程無(wú)解!");return;for(l1=0;l1<m;l1+)for(l2=0;l2<n;l2+)if(M_Bl1l2!=0)Ml2=M_Bl1n-1;printf("輸出線性方程組的解:");for(l2=0;l2<n-1;l2+)printf("x%d(l2+1)=n",Ml2);prin

6、tf("行最簡(jiǎn)型即為所求!按任一健退出.n");getchar();voidmain(void)intp,q,m,n;printf("n請(qǐng)輸入線性方程組的行數(shù)m和列數(shù)n:");scanf("%d%d",&m,&n);printf("n現(xiàn)在輸入系數(shù)矩陣A和向量b:");for(p=1;p<=m;p+)printf("n請(qǐng)輸入a%d1-a%d%d系數(shù)和向量b%d:n",p,p,n,p);/*實(shí)現(xiàn)將每一行中的系數(shù)和向量一次性輸入，數(shù)之間用空格格開(kāi)，輸完后回車(chē)確定*/for(q=1

7、;q<=n+1;q+)/*將剛才輸入的數(shù)存入數(shù)組*/scanf("%f",&M_Bp-1q-1);f1(m,n+1);五、實(shí)驗(yàn)心得對(duì)于第一個(gè)實(shí)驗(yàn)，雖然實(shí)驗(yàn)原理很明了，但是用計(jì)算機(jī)用C語(yǔ)言來(lái)實(shí)現(xiàn)還是一件十分困難的事情，實(shí)驗(yàn)?zāi)軌蜻M(jìn)行的前提是必須理解實(shí)驗(yàn)的原理即初等行變換的精髓，理解實(shí)驗(yàn)原理以后就要明確實(shí)驗(yàn)的各個(gè)步驟如何去實(shí)現(xiàn)，想明白這些以后再?lài)L試編寫(xiě)程序，當(dāng)然其中會(huì)遇到很多困難，但是要有耐心有細(xì)心，就可以編寫(xiě)出自己的程序來(lái)完成這個(gè)實(shí)驗(yàn)。六、程序運(yùn)行情況：請(qǐng)輸入線性方程組的行數(shù)e和列數(shù)於44bl1巨4矢1數(shù)系1-4%al-呵3在輸T觀清L0 0 0 0 0 0 0

8、0 10 104.00-3.001.00-3.004.00-3.00-3.00-3.001.000.000.000.00-2.001.000.000.003.00-1.002.00-4.00-4.001.000.008.004.00-3.0012.00-24.00清輸入a21a24系數(shù)和向量b2：31-11-3請(qǐng)輸入a31a34系數(shù)和向量b3：L3011請(qǐng)輸入a41a44系數(shù)和向量b4：3-731-3以下近程是將矩陣化為行階梯型輸出初始增廣矩陣?-2.003.00-4.001.00-1.001.003.000.001.00-7.003.001.00接任一健繼續(xù)1.00-2.003.00-4.0

9、08.001.00-1.001.060.005.00-3.005.008.00-7.003.001.00接任一健繼續(xù)接任一健繼續(xù)1.00-2.003.00-4.004.006.001.00-1.001.00-3.000.000.001.000.006.000.000.000.008.000.00接任一健繼續(xù)1.002.003.00-4.004.000.001.001.001.00-3.000.000.001.000.00&.0B0.000.000.001.000.00按任一健繼續(xù)矩陣的行階梯型如上?以下過(guò)程是將行階梯型化為行最簡(jiǎn)型1.00-2.003.000.004.000.001.00-1.000.00-3.000.000.001.000.006.000.000.000.001.000.00接任一律繼續(xù).1.00-2.000.000.00-14.000.001.000.000.003.000.000.001.000.006.000.000.060.001.000.00接任一健繼續(xù)1.000.000.000.00-8.000.001.000.000.003.000.000.001.