青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

隨筆-48  評論-259  文章-1  trackbacks-0

#include "iostream.h"
#include"fstream.h"
#include "math.h"

#define N 100

void Jacobi()
{
  int i,j,k,g;
  double a[N][N],b[N],x[N],X[N];
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++)
    istrm>>a[i][j];
    istrm>>b[i];
  }
 
  for(i=0;i<n;i++)
  {
   b[i]/=a[i][i];
   for(j=0;j<n;j++)
   {  if(j==i) continue;
    a[i][j]/=a[i][i];
   }
   a[i][i]=0;
  }
  cout<<"設置初始向量:";
  for(i=0;i<n;i++)
   cin>>x[i];
  int m;
  cout<<"輸入跌代的次數:";
   cin>>m;
  for(g=0;g<m;g++)
  { 
   for(i=0;i<n;i++)
   {
      double sum=0;
      for(j=0;j<n;j++)
   {
    if(j==i) continue;
    sum+=a[i][j]*x[j];
   }
   X[i]=b[i]-sum;
 
   }
   cout<<g<<" ";
   for(k=0;k<n;k++)
   {
    x[k]=X[k];
         cout<<X[k]<<" ";
   }
   cout<<endl;
 
  }
}
void GaussSeidel()
{
  int i,j,k;
  double a[N][N],b[N],x[N];
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++)
    istrm>>a[i][j];
    istrm>>b[i];
  }
  for(i=0;i<n;i++)
  {
   b[i]/=a[i][i];
   for(j=0;j<n;j++)
   {  if(j==i) continue;
    a[i][j]/=a[i][i];
   }
   a[i][i]=0;
  }
  cout<<"設置初始向量";
  for(i=0;i<n;i++)
   cin>>x[i];
  int m;
  cout<<"跌代的次數:";
   cin>>m;
  for(int g=0;g<m;g++)
  {
     for(i=0;i<n;i++)
   {
      double sum=0;
      for(j=0;j<n;j++)
   {
    if(j==i) continue;
    else  sum+=a[i][j]*x[j];
   
   }
   x[i]=b[i]-sum;
 
   }
  cout<<g<<" ";
   for(k=0;k<n;k++)
   {
   
         cout<<x[k]<<" ";
   }
   cout<<endl;
 
  }
}
void SOR()
{
  int i,j,k,cnt=0;
  double a[N][N],b[N],x[N],X[N];
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++)
    istrm>>a[i][j];
    istrm>>b[i];
  }
  for(i=0;i<n;i++)
  {
   b[i]/=a[i][i];
   for(j=0;j<n;j++)
   {  if(j==i) continue;
    a[i][j]/=a[i][i];
   }
   a[i][i]=0;
  }
  cout<<"設置初始向量";
  for(i=0;i<n;i++)
   cin>>x[i];
  double e,m=1;
  cout<<"輸入允許的誤差:";
  cin>>e;
  int w;
  cout<<"輸入松弛因子:";
  cin>>w;
 
  while(m>e)
  {
  cnt++;
     for(i=0;i<n;i++)
   {
         double sum=0;
         for(j=0;j<n;j++)
   {
     if(j==i) continue;
     if(j>i)  sum+=a[i][j]*x[j];
     if(j<i)  sum+=a[i][j]*X[j];
   
   }
      X[i]=x[i]+w*(b[i]-sum-x[i]);
 
  }cout<<cnt<<" ";
     for(k=0,m=0;k<n;k++)
   {
    m+=(X[k]-x[k])*(X[k]-x[k]);
    x[k]=X[k];
         cout<<x[k]<<" ";

   }
     m=sqrt(m);
     cout<<endl;
  }
  cout<<"跌代的次數為"<<cnt<<endl;
}
void Gauss()
{
  int i,j,k;
  float a[N][N],b[N],x[N];
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++)
  istrm>>a[i][j];
   istrm>>b[i];
  }
  for(i=1;i<n;i++)
   for(j=i;j<n;j++)
   { float q=a[j][i-1];
     for(k=i-1;k<n;k++) a[j][k]+=-q*a[i-1][k]/a[i-1][i-1];
  b[j]+=-q*b[i-1]/a[i-1][i-1];
   }
   x[n-1]=b[n-1]/a[n-1][n-1];
   for(i=n-2;i>=0;i--)
   {
     float sum=0;
  for(j=n-1;j>i;j--) sum+=x[j]*a[i][j];
        x[i]=(b[i]-sum)/a[i][i];
   }
  for(i=0;i<n;i++)
    cout<<x[i]<<" ";
}
void SanJiao()
{
  int i,j,r;
  float a[N][N],b[N],l[N][N],u[N][N];
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++)
  istrm>>a[i][j];
   istrm>>b[i];
  }
 for(i=0;i<n;i++)
  for(j=0;j<n;j++)
  {
    if(i<j) l[i][j]=0;
    if(i==j) l[i][j]=1;
    else u[i][j]=0;
  }
  u[0][0]=a[0][0];
   for(i=1;i<n;i++)
  {
    u[0][i]=a[0][i];
 l[i][0]=a[i][0]/u[0][0];
  }
  for(r=1;r<n;r++)
  {
    for(i=r;i<n;i++)
 {
   u[r][i]=a[r][i];
   for(j=0;j<=r-1;j++)
   {  
    u[r][i]-=l[r][j]*u[j][i];
   
   }

 }
 for(i=r+1;i<n;i++)
 {
   l[i][r]=a[i][r];
   for(j=0;j<=r-1;j++)
   {l[i][r]-=l[i][j]*u[j][r];}
   l[i][r]/=u[r][r];
 }
  }
  for(i=0;i<n;i++)
 {  cout<<"|";
  for(j=0;j<n;j++)
   {
    cout<<l[i][j]<<" ";
   }
  cout<<"|";
  cout<<" "<<"|";
  for(j=0;j<n;j++)
   {
     cout<<u[i][j]<<" ";
   }
    cout<<"|"<<endl;
 }
}
void Sqrt()
{
    int i,j,k;
  double temp;
  double a[10][10],l[10][10]={0};
  int n;
  cout<<"輸入行列式的階數:";
  cin>>n;
  ifstream istrm("data.txt");
  for(i=1;i<=n;i++)
    for(j=1;j<=n;j++)
  istrm>>a[i][j];
  for(j=1;j<=n;j++)
  {temp=0;
   for(k=1;k<j;k++) temp+=(l[j][k]*l[j][k]);
    l[j][j]=sqrt(a[j][j]-temp);
 for(i=j+1;i<=n;i++)
 {
  temp=0;
   for(k=1;k<j;k++) temp+=(l[j][k]*l[i][k]);
   l[i][j]=(a[i][j]-temp)/l[j][j];
  
 }
  }
  for(i=1;i<=n;i++)
 {  cout<<"|";
  for(j=1;j<=n;j++)
   {
    cout<<l[i][j]<<" ";
   }
  cout<<"|"<<endl;
  
 
  }
}
void main()
{
  int n,m;
  cout<<"***********提示菜單***************"<<endl;
  cout<<"*******1-線性方程的直接解法*******"<<endl;
  cout<<"*******2-線性方程的間接解法*******"<<endl;
  cout<<"*******0-退出程序*****************"<<endl;
  cout<<"**********************************"<<endl;
  while(1)
  {
    cout<<"輸入你要進行的操作:";
 cin>>n;
 if(n==0) break;
 switch(n)
 {
 case 1:cout<<"***********************"<<endl;
     cout<<"***1-高斯消去法********"<<endl;
     cout<<"***2-矩陣三角分解法****"<<endl;
     cout<<"***3-平方根分解法******"<<endl;
     cout<<"***0-返回上一層********"<<endl;
     cout<<"***********************"<<endl;
     cin>>m;
     switch(m)
     {
     case 1:Gauss();break;
     case 2:SanJiao();break;
     case 3:Sqrt();break;
     case 0:break;
     default:cout<<"輸入錯誤請重新輸入"<<endl;      
     }
     break;
 case 2:cout<<"***********************"<<endl;
     cout<<"***1-雅可比迭代法******"<<endl;
     cout<<"***2-高斯-塞德爾迭代法*"<<endl;
     cout<<"***3-松弛迭代法********"<<endl;
     cout<<"***0-返回上一層********"<<endl;
     cout<<"***********************"<<endl;
     cin>>m;
     switch(m)
     {
     case 1:Jacobi();break;
     case 2:GaussSeidel();break;
     case 3:SOR();break;
           case 0:break;
     default:cout<<"輸入錯誤請重新輸入"<<endl;
     }break;
 default:cout<<"輸入錯誤請重新輸入"<<endl;
 }
  }

}

posted on 2007-06-08 00:25 星夢情緣 閱讀(616) 評論(0)  編輯 收藏 引用 所屬分類: 數據結構的所有實現程序
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            欧美视频在线免费| 欧美国产精品劲爆| 国产欧美一区二区三区另类精品 | 蜜臀久久99精品久久久久久9 | 欧美激情1区2区| 宅男噜噜噜66国产日韩在线观看| 亚洲激情黄色| 欧美日韩国内自拍| 午夜国产不卡在线观看视频| 午夜欧美大尺度福利影院在线看 | 欧美夫妇交换俱乐部在线观看| 亚洲激情成人| 亚洲一区视频| 久久久999国产| 亚洲精品综合在线| 国产精品99久久久久久有的能看| 国产精品专区一| 欧美成人第一页| 欧美日韩一区三区| 久久久久国产精品一区| 美国成人毛片| 午夜一级久久| 欧美a级大片| 性欧美videos另类喷潮| 免费观看成人| 午夜精品视频一区| 嫩草影视亚洲| 久久黄色影院| 欧美日韩视频专区在线播放 | 一区二区三欧美| 欧美一级久久| 中文网丁香综合网| 久久精品麻豆| 午夜一区二区三区在线观看| 你懂的视频一区二区| 性欧美激情精品| 欧美欧美天天天天操| 久久久最新网址| 国产精品免费小视频| 亚洲七七久久综合桃花剧情介绍| 国产一区二区三区四区| 一区二区三区欧美在线| 最新国产成人av网站网址麻豆| 午夜精品久久久久久久99热浪潮| 999亚洲国产精| 麻豆精品91| 另类图片综合电影| 国产午夜精品理论片a级大结局 | 久久精品国产亚洲aⅴ| 欧美成人精品一区二区| 亚洲欧美日韩国产一区二区三区 | 亚洲综合电影| 欧美日韩国产精品一卡| 欧美国产精品va在线观看| 国产一区二区三区网站| 亚洲欧美日产图| 午夜久久久久| 国产精品视频免费一区| 亚洲美女一区| 亚洲亚洲精品三区日韩精品在线视频| 狂野欧美一区| 欧美激情亚洲| 亚洲精品小视频在线观看| 免费的成人av| 亚洲福利电影| 亚洲麻豆国产自偷在线| 国产精品美女午夜av| 亚洲美女精品成人在线视频| 亚洲美洲欧洲综合国产一区| 欧美国产亚洲精品久久久8v| 亚洲第一在线综合网站| 欧美丰满高潮xxxx喷水动漫| 亚洲成在线观看| 亚洲激情不卡| 欧美日韩精品一区二区在线播放| 亚洲精品国产精品国自产在线| 99av国产精品欲麻豆| 欧美性大战xxxxx久久久| 亚洲欧美精品suv| 久久亚洲一区| 最新亚洲一区| 国产精品99一区| 亚洲欧美日本精品| 欧美www视频在线观看| 亚洲国产片色| 国产精品国内视频| 久久都是精品| 亚洲国产精选| 亚欧成人在线| 精品成人国产| 欧美日韩视频第一区| 亚洲欧美久久久| 亚洲国产黄色| 欧美一级专区| 亚洲高清久久网| 国产精品久久7| 久久亚洲国产精品一区二区| 亚洲精品孕妇| 久久久久久91香蕉国产| 亚洲精品日韩一| 国产日韩欧美不卡在线| 暖暖成人免费视频| 亚洲欧美福利一区二区| 亚洲国产欧美日韩| 久久成人精品| 99精品国产一区二区青青牛奶| 国产伦精品一区二区三区在线观看 | 久久岛国电影| 日韩视频免费在线观看| 国产精品成人在线| 久久综合成人精品亚洲另类欧美| 亚洲美女免费视频| 免费在线观看成人av| 亚洲欧美精品| 日韩亚洲欧美精品| 狠狠色狠狠色综合日日91app| 欧美日韩一区二区三区在线视频| 久久久999国产| 亚洲在线观看| 99这里只有精品| 亚洲国产日韩综合一区| 久久这里只精品最新地址| 亚洲欧美日韩精品久久奇米色影视| 亚洲国产成人porn| 红桃视频国产精品| 国产女主播一区二区| 欧美日韩色一区| 欧美精品激情| 欧美激情视频一区二区三区免费 | 久久大逼视频| 亚洲综合电影一区二区三区| 亚洲日本中文字幕免费在线不卡| 免费精品视频| 麻豆91精品| 老司机免费视频一区二区| 欧美在线首页| 欧美在线亚洲一区| 欧美一区观看| 亚洲综合久久久久| 国产日韩精品在线| 国产精品爽爽爽| 国产精品一区二区三区成人| 欧美性色aⅴ视频一区日韩精品| 欧美国产免费| 欧美啪啪一区| 欧美日韩一二区| 国产精品久久久对白| 国产精品老牛| 国产日韩精品一区二区三区| 国产精品视频免费| 国产亚洲成av人在线观看导航| 国产欧美日韩激情| 国产亚洲成av人片在线观看桃| 国产一区二区三区黄| 国产综合久久| 尹人成人综合网| 亚洲欧洲精品天堂一级| a4yy欧美一区二区三区| 亚洲午夜视频| 欧美一区二区在线免费观看| 久久精品国产99| 米奇777在线欧美播放| 欧美激情第五页| 99精品黄色片免费大全| 亚洲欧美日韩另类| 久久嫩草精品久久久精品一| 免费在线观看成人av| 欧美日韩国产综合新一区| 国产精品久久毛片a| 狠狠久久婷婷| 日韩一级在线观看| 欧美影院一区| 亚洲春色另类小说| 亚洲视频狠狠| 玖玖玖国产精品| 欧美性猛交xxxx免费看久久久| 国产有码一区二区| 一本色道久久88精品综合| 久久av资源网站| 亚洲高清免费| 亚洲欧美视频在线观看视频| 老司机久久99久久精品播放免费| 欧美区日韩区| 黄色av日韩| 亚洲欧美激情视频| 亚洲夫妻自拍| 欧美一区二区观看视频| 欧美精选午夜久久久乱码6080| 国产精品午夜电影| 亚洲精品在线免费| 久久中文精品| 亚洲午夜精品在线| 欧美国产精品一区| 在线免费观看欧美| 香蕉久久夜色精品国产使用方法| 欧美freesex8一10精品| 午夜在线视频一区二区区别| 欧美精品精品一区| 亚洲国产精品一区二区尤物区| 欧美一区二区三区在线观看|