關於矩陣求聯立解

jhsdenny

【問題敘述】：矩陣求聯立解寫到一半，不過輸出格式一直沒辦法正確。
【程式語言】：C
【Compiler 軟體】： DEV-C++
【作業系統】： Win XP

【目前程式碼】：


【Compiler 訊息】：無誤
【目前執行結果】：     
                             
【期望可得到之執行結果】：
...

瀏覽完整內容，請先 註冊 或 登入會員

p12332145600

1. 
   
2. #include<stdio.h>
   
3. #define p printf
   
4. void main(void)
   
5. {
   
6.      int a[3][4]={{1,-2,4,12},{2,-1,5,18},{-1,3,-3,-8}};
   
7.      int i,j,x,tmp,divide1,divide2;
   
8.      for(i=0;i<3;i++){
   
9.           for(j=0;j<4;j++)
   
10.                p("%d ",a[i][j]);
    
11.      p("\n");     
    
12.      }
    
13. 
    
14.      for(i=0;i<3;i++)
    
15.      {
    
16.           for(j=i;j<4;j++)
    
17.           {
    
18.                if(i==j) //(0.0)(1.1)(2.2) 位置
    
19.                {     
    
20.                     tmp=a[i][j];     
    
21.                     a[i][j]=1;     //(0.0)(1.1)(2.2) 位置設1     
    
22.                     for(x=i+1;x<4;x++)
    
23.                     {
    
24.                          a[i][x]=a[i][x]/tmp;     
    
25.                     }
    
26.                     if(i==0)
    
27.                     {
    
28.                          divide1=a[i+1][j]-a[i+1][j]*2;     
    
29.                          divide2=a[i+2][j]-a[i+2][j]*2;
    
30.                          //求出差取反項
    
31.                     }
    
32.                     else if(i==1)
    
33.                     {
    
34.                          divide1=a[i-1][j]-a[i-1][j]*2;
    
35.                          divide2=a[i+1][j]-a[i+1][j]*2;
    
36.                     }
    
37.                     else if(i==2)
    
38.                     {
    
39.                          divide1=a[i-2][j]-a[i-2][j]*2;
    
40.                          divide2=a[i-1][j]-a[i-1][j]*2;
    
41.                     }
    
42.                }
    
43.                if(i==0)
    
44.                {
    
45.                     a[i+1][j]=a[i][j]*divide1+a[i+1][j];
    
46.                     a[i+2][j]=a[i][j]*divide2+a[i+2][j];
    
47.                     //求出反項乘那一行在和別行做相加
    
48.                }
    
49.                else if(i==1)
    
50.                {
    
51.                     a[i-1][j]=a[i][j]*divide1+a[i-1][j];
    
52.                     a[i+1][j]=a[i][j]*divide2+a[i+1][j];
    
53.                }
    
54.                else if(i==2)
    
55.                {
    
56.                     a[i-2][j]=a[i][j]*divide1+a[i-2][j];
    
57.                     a[i-1][j]=a[i][j]*divide2+a[i-1][j];
    
58.                }
    
59.           }
    
60.      }
    
61.      for(i=0;i<3;i++){
    
62.           for(j=0;j<4;j++)
    
63.                p("%d ",a[i][j]);     
    
64.      p("\n");     
    
65.      }
    
66. }

複製代碼

...

瀏覽完整內容，請先 註冊 或 登入會員

ponchi96

回復

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

p12332145600

可以幫忙翻譯一下他的題目嗎@@
我看不懂他input 跟 output的關聯是什麼

p12332145600

回復

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

ponchi96
1*X1 - 2*X2 + 4*X3 = 12
2*X1 - 1*X2 + 5*X3 = 18
-1*X1 + 3*X2 - 3*X3 = -8

把它的係數取出來  求 X1 X2 X3

...

瀏覽完整內容，請先 註冊 或 登入會員

jhsdenny

回復

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

p12332145600

非常感謝 大大精采的解說
將大大的程式 改變一下後 (改成輸入 解答用分數表示 使用gcd) 也越接近想要的輸出模式

目前正在研究無窮解 (找出一個一組解)


此外小弟還想請教 如果 變數的數目範圍(0<n<50)  然後 每行代表一個方程式(1<= m <=n)
...

瀏覽完整內容，請先 註冊 或 登入會員