版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

[交流] 【求助】為什么我編譯了一個(gè)書本上的無錯(cuò)誤的C++程序，卻有這么多的錯(cuò)誤？

CODE:

#include

#include

#include

using namespace std;

class bint

{

private:

        int n;

        double **a,*b,*x,eps;

        double **p,*r,*e;

public:

        bint(int nn)

        {

                int i;

                n=nn;

                a=new double*[n];

                for(i=0;i
                p=new double*[n];

                for(i=0;i
                r=new double[n];

                e=new double[n];

                b=new double[n];

                x=new double[n];

        }

        void input();

        void a_bint();

        void gauss();

        void output();

        bint()

        {

                int i;

                for(i=0;i
                delete[] a;

                for(i=0;i
                delete[] p;

                delete[] r,e,b,x;

        }

};

void bint::input()

{

        int i,j;

        char str1[20];

        cout<<"\n輸入文件名:";

        cin>>str1;

        ifstream fin(str1);

        if(!fin)

        {cout<<"\n不能打開這個(gè)文件"<
        for(i=0;i
                for(j=0;j>a[i][j];

                for(i=0;i>b[i];

                fin.close();

}

void bint::a_bint()

{

        int i,j,k;

        double q,qq;

        i=100;

        for(k=0;k<=n-1;k++)

                for(j=0;j<=n-1;j++)p[k][j]=a[k][j];

                for(k=0;k<=n-1;k++)r[k]=b[k];

                gauss();

                for(k=0;k
                q=1.0+eps;

                while(q>=eps)

                {

                        if(i==0)return;

                        i=i-1;

                        for(k=0;k
                        {

                                e[k]=0.0;

                                for(j=0;j
                                        e[k]=e[k]+a[k][j]*x[j];

                        }

                        for(k=0;k<=n-1;k++)r[k]=b[k]-e[k];

                                for(k=0;k<=n-1;k++)

                                        for(j=0;j<=n-1;j++)p[k][j]=a[k][j];

                                        gauss();

                                        q=0.0;

                                        for(k=0;k<=n-1;k++)

                                        {

                                                qq=fabs(r[k])/(1.0+fabs(x[k]+r[k]));

                                                if(qq>q)q=qq;

                                        }

                                        for(k=0;k<=n-1;k++)x[k]=x[k]+r[k];

                }

}

void bint::gauss()

{

        int *js,l,k,i,j,is;

        double d,t;

        js=new int[n];

        l=1;

        for(k=0;k<=n-2;k++)

        {

                d=0.0;

                for(i=k;i<=n-1;i++)

                        for(j=k;j<=n-1;j++)

                        {

                                t=fabs(p[i][j]);

                                if(t>d){d=t;js[k]=j;is=i;}

                        }

                        if(d+1.0==1.0)l=0;

                        else

                        {if(js[k]!=k)

                        for(i=0;i<=n-1;i++)

                        {

                                t=p[i][k];

                                p[i][k]=p[i][js[k]];

                                p[i][js[k]]=t;

                        }

                        if(is!=k)

                        {

                                for(j=k;j<=n-1;j++)

                                {

                                        t=p[k][j];

                                        p[k][j]=p[is][j];

                                        p[is][j]=t;

                                }

                                    t=r[k];r[k]=r[is];r[is]=t;

                        }

                }

                        if(l==0)

                        {

                                delete[] js;

                                cout<<"\n系數(shù)矩陣奇異!無解."<
                                return;

                        }

                        d=p[k][k];

                        for(j=k+1;j<=n-1;j++)

                                p[k][j]=p[k][j]/d;

                        r[k]=r[k]/d;

                        for(i=k+1;i<=n-1;i++)

                                p[i][j]=p[i][j]-p[i][k]*p[k][j];

                        r[i]=r[i]-p[i][k]*r[k];

        }

}

            d=p[n-1][n-1];

            if(fabs(d)+1.0==1.0)

                        {

                delete[] js;

                                cout<<"\n系數(shù)矩陣奇異!無解."<
                            return;

                        }

            r[n-1]=r[n-1]/d;

            for(i=n-2;i>=0;i--)

                        {

                t=0.0;

                for(j=i+1;j<=n-1;j++)

                       t=t+p[i][j]*r[j];

                  r[i]=r[i]-t;

                        }

            js[n-1]=n-1;

          for(k=n-1;k>=0;k--)

          if(js[k]!=k)

                  {

               t=r[k];r[k]=r[js[k]];r[js[k]]=t;

                  }

       delete[] js;

    }

void bint::output()

{

        int i;

        char str2[20];

        cout<<"\n輸出文件名:";

        cin>>str2;

        ofstream fout(str2);

        if(!fout)

        {cout<<"\n不能打開這個(gè)文件"<
        fout<
        for(i=0;i
        {

                fout<
                        cout<
        }

        fout<
        fout.close();

}

void main()

{

        bint c(4);

        c.input();

        c.a_bint();

        c.output();

}

[ Last edited by zyj8119 on 2010-10-12 at 01:57 ]

回復(fù)此樓

» 猜你喜歡

304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)
330求調(diào)劑已經(jīng)有3人回復(fù)

好好學(xué)習(xí)，天天向上。

1樓 2010-10-12 01:50:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

ajian04:歡迎交流~ 2010-10-12 08:16:06

錯(cuò)誤是：
--------------------Configuration: 1 - Win32 Debug--------------------
Compiling...
bingtai.cpp
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(128) : error C2143: syntax error : missing ';' before 'return'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(139) : error C2501: 'd' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(139) : error C2065: 'p' : undeclared identifier
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(139) : error C2065: 'n' : undeclared identifier
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(139) : error C2109: subscript requires array or pointer type
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(139) : error C2109: subscript requires array or pointer type
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(140) : error C2143: syntax error : missing ';' before 'if'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(141) : error C2143: syntax error : missing ';' before '{'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(141) : error C2447: missing function header (old-style formal list?)
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(146) : error C2057: expected constant expression
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(146) : error C2466: cannot allocate an array of constant size 0
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(146) : error C2501: 'r' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(146) : error C2440: 'initializing' : cannot convert from 'int' to 'int []'
There are no conversions to array types, although there are conversions to references or pointers to arrays
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ';' before 'for'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ')' before ';'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ';' before '>='
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2501: 'i' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ';' before '>='
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ';' before '--'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2501: 'i' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2086: 'i' : redefinition
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2143: syntax error : missing ';' before '--'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(147) : error C2059: syntax error : ')'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(148) : error C2143: syntax error : missing ';' before '{'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(148) : error C2447: missing function header (old-style formal list?)
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(154) : error C2057: expected constant expression
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(154) : error C2466: cannot allocate an array of constant size 0
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(154) : error C2501: 'js' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(154) : error C2440: 'initializing' : cannot convert from 'int' to 'int []'
There are no conversions to array types, although there are conversions to references or pointers to arrays
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ';' before 'for'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ')' before ';'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ';' before '>='
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2501: 'k' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ';' before '>='
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ';' before '--'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2501: 'k' : missing storage-class or type specifiers
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2086: 'k' : redefinition
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2143: syntax error : missing ';' before '--'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(155) : error C2059: syntax error : ')'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(157) : error C2143: syntax error : missing ';' before '{'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(157) : error C2447: missing function header (old-style formal list?)
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(160) : error C2143: syntax error : missing ';' before 'delete'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(161) : error C2143: syntax error : missing ';' before '}'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(161) : error C2143: syntax error : missing ';' before '}'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(161) : error C2143: syntax error : missing ';' before '}'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(164) : error C2143: syntax error : missing ';' before '{'
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(164) : error C2447: missing function header (old-style formal list?)
Error executing cl.exe.

bingtai.obj - 47 error(s), 0 warning(s)

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí)，天天向上。

2樓2010-10-12 01:50:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

ajian04:加油~ 2010-10-12 08:16:14

引用回帖:

Originally posted by zyj8119 at 2010-10-12 01:50:50:
錯(cuò)誤是：
--------------------Configuration: 1 - Win32 Debug--------------------
Compiling...
bingtai.cpp
C:\Program Files\Microsoft Visual Studio\MyProjects\1\bingtai.cpp(128) : error C2143: sy ...

用VC++6.0編譯。

回復(fù)此樓

好好學(xué)習(xí)，天天向上。

3樓2010-10-12 01:55:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ajian04

禁蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
貴賓: 0.034
金幣: 48.8
散金: 795
紅花: 2
帖子: 637
在線: 39.3小時(shí)
蟲號(hào): 1007856
注冊(cè): 2010-04-28
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

zyj8119(金幣+3):大括號(hào)的位置有要求吧？ 2010-10-12 10:03:41

是你直接抄寫的嗎？遺漏了很多符號(hào)，還有可能是你的編譯環(huán)境不對(duì)，好好看看書上介紹的運(yùn)行環(huán)境吧

贊一下

回復(fù)此樓

4樓2010-10-12 08:15:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

木葉清風(fēng)

木蟲 (正式寫手)

應(yīng)助: 18 (小學(xué)生)
金幣: 4724.6
紅花: 9
帖子: 364
在線: 38.2小時(shí)
蟲號(hào): 276214
注冊(cè): 2006-09-02
專業(yè): 計(jì)算機(jī)軟件

★ ★
zyj8119(金幣+5):何謂連帶錯(cuò)誤？我的大括號(hào)的對(duì)數(shù)，感覺也是對(duì)的。 2010-10-12 09:55:25
ajian04(金幣+2):謝謝參與交流~ 2010-10-12 10:08:27

里面有些錯(cuò)誤是連帶的，可以參考下面的修改

CODE:

#include

#include

#include

using namespace std;

class bint

{

private:

        int n;

        double **a,*b,*x,eps;

        double **p,*r,*e;

public:

        bint(int nn)

        {

                int i;

                n=nn;

                a=new double*[n];

                for(i=0;i
                p=new double*[n];

                for(i=0;i
                r=new double[n];

                e=new double[n];

                b=new double[n];

                x=new double[n];

        }

        void input();

        void a_bint();

        void gauss();

        void output();

        ~bint()

        {

                int i;

                for(i=0;i
                delete[] a;

                for(i=0;i
                delete[] p;

                delete[] r,e,b,x;

        }

};

void bint::input()

{

        int i,j;

        char str1[20];

        cout<<"\n輸入文件名:";

        cin>>str1;

        ifstream fin(str1);

        if(!fin)

        {cout<<"\n不能打開這個(gè)文件"<
        for(i=0;i
                for(j=0;j>a[i][j];

        for(i=0;i>b[i];

        fin.close();

}

void bint::a_bint()

{

        int i,j,k;

        double q,qq;

        i=100;

        for(k=0;k<=n-1;k++)

                for(j=0;j<=n-1;j++)p[k][j]=a[k][j];

        for(k=0;k<=n-1;k++)r[k]=b[k];

        gauss();

        for(k=0;k
        q=1.0+eps;

        while(q>=eps)

        {

                if(i==0)return;

                i=i-1;

                for(k=0;k
                {

                        e[k]=0.0;

                        for(j=0;j
                                e[k]=e[k]+a[k][j]*x[j];

                }

                for(k=0;k<=n-1;k++)r[k]=b[k]-e[k];

                for(k=0;k<=n-1;k++)

                        for(j=0;j<=n-1;j++)p[k][j]=a[k][j];

                gauss();

                q=0.0;

                for(k=0;k<=n-1;k++)

                {

                        qq=fabs(r[k])/(1.0+fabs(x[k]+r[k]));

                        if(qq>q)q=qq;

                }

                for(k=0;k<=n-1;k++)x[k]=x[k]+r[k];

        }

}

void bint::gauss()

{

        int *js,l,k,i,j,is;

        double d,t;

        js=new int[n];

        l=1;

        for(k=0;k<=n-2;k++)

        {

                d=0.0;

                for(i=k;i<=n-1;i++)

                        for(j=k;j<=n-1;j++)

                        {

                                t=fabs(p[i][j]);

                                if(t>d){d=t;js[k]=j;is=i;}

                        }

                if(d+1.0==1.0)l=0;

                else

                {

                        if(js[k]!=k)

                        for(i=0;i<=n-1;i++)

                        {

                                t=p[i][k];

                                p[i][k]=p[i][js[k]];

                                p[i][js[k]]=t;

                        }

                        if(is!=k)

                        {

                                for(j=k;j<=n-1;j++)

                                {

                                        t=p[k][j];

                                        p[k][j]=p[is][j];

                                        p[is][j]=t;

                                }

                                t=r[k];r[k]=r[is];r[is]=t;

                        }

                }

                        if(l==0)

                        {

                                delete[] js;

                                cout<<"\n系數(shù)矩陣奇異!無解."<
                                        return;

                        }

                        d=p[k][k];

                        for(j=k+1;j<=n-1;j++)

                                p[k][j]=p[k][j]/d;

                        r[k]=r[k]/d;

                        for(i=k+1;i<=n-1;i++)

                                p[i][j]=p[i][j]-p[i][k]*p[k][j];

                        r[i]=r[i]-p[i][k]*r[k];

        }

d=p[n-1][n-1];

if(fabs(d)+1.0==1.0)

{

        delete[] js;

        cout<<"\n系數(shù)矩陣奇異!無解."<
                return;

}

r[n-1]=r[n-1]/d;

for(i=n-2;i>=0;i--)

{

        t=0.0;

        for(j=i+1;j<=n-1;j++)

                t=t+p[i][j]*r[j];

        r[i]=r[i]-t;

}

js[n-1]=n-1;

for(k=n-1;k>=0;k--)

if(js[k]!=k)

{

        t=r[k];r[k]=r[js[k]];r[js[k]]=t;

}

delete[] js;

        }

        void bint::output()

        {

                int i;

                char str2[20];

                cout<<"\n輸出文件名:";

                cin>>str2;

                ofstream fout(str2);

                if(!fout)

                {cout<<"\n不能打開這個(gè)文件"<
                fout<
                for(i=0;i
                {

                        fout<
                        cout<
                }

                fout<
                fout.close();

        }

        void main()

        {

                bint c(4);

                c.input();

                c.a_bint();

                c.output();

        }

贊一下(1人)

回復(fù)此樓

www.cvdelver.com

5樓2010-10-12 09:45:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

joy4210

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 3095.3
散金: 2
帖子: 219
在線: 114.7小時(shí)
蟲號(hào): 1096044
注冊(cè): 2010-09-11
專業(yè): 原子和分子物理

★
ajian04(金幣+1):歡迎交流~ 2010-10-12 10:08:37

是徐士良書上的吧

CODE:

  //3BINT.CPP

  //求解病態(tài)方程組

  #include  

  #include  

  #include  

  using namespace std;

  class  bint

  {

  private: 

           int n;

                   double  **a, *b, *x, eps;

                   double  **p, *r, *e;

  public:

               bint (int nn)

                   {

                           int i;

                           n = nn;

                   a = new double*[n];   //動(dòng)態(tài)分配內(nèi)存空間

                   for (i=0; i
               p = new double*[n];

                   for (i=0; i
               r = new double[n];

               e = new double[n];

                   b = new double[n];

                   x = new double[n];

                   }

               void input ();  //從文件讀入系數(shù)矩陣A以及常數(shù)向量B與eps

                   void a_bint ();    //執(zhí)行求解病態(tài)方程組的步驟

                   void gauss ();     //全選主元Gauss消去法

           void output ();   //輸出結(jié)果到文件并顯示

                   ~bint ()

                   {

                           int i;

                           for (i=0; i
                           delete [] a;

                           for (i=0; i
                           delete [] p;

                           delete [] r, e, b, x;

                   }

  };

  void bint::input ()    //從文件讀入系數(shù)矩陣A以及常數(shù)向量B與eps

  {

          int  i, j;

          char str1[20];

          cout <<"\n輸入文件名:  ";

          cin >>str1;

          ifstream  fin (str1);

          if (!fin)

          { cout <<"\n不能打開這個(gè)文件 " <
          for (i=0; i
                  for (j=0; j>a[i][j];

          for (i=0; i>b[i];           //讀入常數(shù)向量B

          fin.close ();

  }

  void bint::a_bint ()       //執(zhí)行求解病態(tài)方程組的步驟

  { 

          int i,j,k;

      double q,qq;

      i=100;

      for (k=0; k<=n-1; k++)

      for (j=0; j<=n-1; j++)  p[k][j]=a[k][j];

      for (k=0; k<=n-1; k++)  r[k]=b[k]; 

      gauss ();

      for (k=0; k
      q=1.0+eps;

      while (q>=eps)

      { 

                  if (i==0)  return;

          i=i-1;

                  for (k=0; k
                  {

                          e[k] = 0.0;

                          for (j=0; j
                                  e[k] = e[k] + a[k][j]*x[j];

                  }

          for ( k=0; k<=n-1; k++)  r[k]=b[k]-e[k];

          for ( k=0; k<=n-1; k++)

          for ( j=0; j<=n-1; j++)  p[k][j]=a[k][j];

          gauss ();

          q=0.0;

          for ( k=0; k<=n-1; k++)

          { 

                          qq=fabs(r[k])/(1.0+fabs(x[k]+r[k]));

              if (qq>q) q=qq;

          }

          for ( k=0; k<=n-1; k++)  x[k]=x[k]+r[k];

      }

  }

  void bint::gauss ()         //執(zhí)行全選主元Gauss消去法

  { 

          int *js,l,k,i,j,is;

      double d,t;

      js = new int[n];

      l=1;

      for (k=0; k<=n-2; k++)

      { 

                  d=0.0;

          for (i=k;i<=n-1;i++)

          for (j=k;j<=n-1;j++)

          { 

                          t=fabs(p[i][j]);

              if (t>d) { d=t; js[k]=j; is=i;}

          }

          if (d+1.0==1.0) l=0;

          else

          { if (js[k]!=k)

              for (i=0;i<=n-1;i++)

              { 

                  t=p[i][k]; 

                                  p[i][k]=p[i][js[k]]; 

                                  p[i][js[k]]=t;

              }

              if (is!=k)

              { 

                                  for (j=k;j<=n-1;j++)

                  { 

                      t=p[k][j]; 

                                          p[k][j]=p[is][j]; 

                                          p[is][j]=t;

                  }

                  t=r[k]; r[k]=r[is]; r[is]=t;

              }

          }

          if (l==0)

          { 

                          delete [] js;

                          cout <<"\n系數(shù)矩陣奇異！無解." <
              return;

          }

          d=p[k][k];

          for (j=k+1;j<=n-1;j++)

              p[k][j]=p[k][j]/d;

          r[k]=r[k]/d;

          for (i=k+1;i<=n-1;i++)

          { 

                          for (j=k+1;j<=n-1;j++)

                  p[i][j]=p[i][j]-p[i][k]*p[k][j];

              r[i]=r[i]-p[i][k]*r[k];

          }

      }

      d=p[n-1][n-1];

      if (fabs(d)+1.0==1.0)

      { 

                  delete [] js;

                  cout <<"\n系數(shù)矩陣奇異！無解." <
          return;

      }

      r[n-1]=r[n-1]/d;

      for (i=n-2;i>=0;i--)

      { 

                  t=0.0;

          for (j=i+1;j<=n-1;j++)

              t=t+p[i][j]*r[j];

          r[i]=r[i]-t;

      }

      js[n-1]=n-1;

      for (k=n-1;k>=0;k--)

        if (js[k]!=k)

        { 

                        t=r[k]; r[k]=r[js[k]]; r[js[k]]=t;

                }

    delete [] js;

  }

  void bint::output ()       //輸出結(jié)果到文件并顯示

  {

          int  i;

          char str2[20];

          cout <<"\n輸出文件名:  ";

          cin >>str2;

          ofstream fout (str2);

          if (!fout)

          { cout <<"\n不能打開這個(gè)文件 " <
          fout <
          for (i=0; i
          {

                  fout <
                  cout <
          }

          fout <
          fout.close ();

  }

  void main ()      //主函數(shù)

  {

          bint  c(4); 

          c.input ();         //從文件讀入系數(shù)矩陣A以及常數(shù)向量B與eps

          c.a_bint ();         //求解病態(tài)方程組

          c.output ();        //輸出結(jié)果到文件并顯示

  }

贊一下(1人)

回復(fù)此樓

6樓2010-10-12 09:50:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

ajian04:謝謝參與交流~ 2010-10-12 10:08:47

引用回帖:

Originally posted by joy4210 at 2010-10-12 09:50:07:
是徐士良書上的吧
[code]
  //3BINT.CPP
  //求解病態(tài)方程組
  #include
  #include
  #include
  using namespace std;
  class  bint
  {
  priv ...

是的，是徐士良的常用算法程序集（C++語言描述）

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí)，天天向上。

7樓2010-10-12 09:51:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

Originally posted by ajian04 at 2010-10-12 08:15:35:
是你直接抄寫的嗎？遺漏了很多符號(hào)，還有可能是你的編譯環(huán)境不對(duì)，好好看看書上介紹的運(yùn)行環(huán)境吧

是直接抄寫的，是不是大括號(hào)的位置也是有要求的？？？？？

贊一下

回復(fù)此樓

好好學(xué)習(xí)，天天向上。

8樓2010-10-12 10:01:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ajian04

禁蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
貴賓: 0.034
金幣: 48.8
散金: 795
紅花: 2
帖子: 637
在線: 39.3小時(shí)
蟲號(hào): 1007856
注冊(cè): 2010-04-28
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

★
zyj8119(金幣+2):可能我有些是中文情況下輸入的。 2010-10-12 10:57:59
余澤成(金幣+1):謝謝參與應(yīng)助！ 2010-10-12 17:44:04

引用回帖:

Originally posted by zyj8119 at 2010-10-12 10:01:17:

是直接抄寫的，是不是大括號(hào)的位置也是有要求的？？？？？

位置應(yīng)該沒有要求，只是不要在中文輸入法狀態(tài)下輸入就行了

贊一下(1人)

回復(fù)此樓

9樓2010-10-12 10:09:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zyj8119 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +7	Liwangman 2026-03-15	7/350	2026-03-18 20:08 by walc
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9	秋有木北 2026-03-14	9/450	2026-03-18 19:49 by macy2011
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[考研] 307求調(diào)劑 +3	冷笙123 2026-03-17	3/150	2026-03-18 09:55 by macy2011
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky

24小時(shí)熱門版塊排行榜

[交流] 【求助】為什么我編譯了一個(gè)書本上的無錯(cuò)誤的C++程序，卻有這么多的錯(cuò)誤？

» 猜你喜歡

[交流] 【求助】為什么我編譯了一個(gè)書本上的無錯(cuò)誤的C++程序，卻有這么多的錯(cuò)誤？