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

返回列表

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號: 664177

[交流] 【求助】程序運(yùn)行已經(jīng)無錯(cuò)誤，但是輸入文件有問題，求教。

CODE:

#include

#include

#include

#include

using namespace std;

class splin1

{

private:

        int n,m;

        double *x,*y,*dy,*ddy,*t,*z,*dz,*ddz,integ;

        double *s;

public:

        splin1 (int nn,int mm)

        {

                n=nn;m=mm;

                x=new double[n];

                y=new double[n];

                dy=new double[n];

                ddy=new double[n];

                s=new double[n];

                t=new double[n];

                z=new double[m];

                dz=new double[m];

                ddz=new double[m];

        }

        void input();

        void interp();

        void output();

        ~splin1()

        {delete[] x,y,dy,ddy,t,z,dz,ddz,s;}

};

void splin1::input()

{

        int k;

        char str1[20];

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

        cin>>str1;

        ifstream fin(str1);

        if(!fin)

        {cout<<"\n不能打開這個(gè)文件"<
        for(k=0;k
        {fin>>x[k];fin>>y[k];}

        fin>>dy[0]>>dy[n-1];

        for(k=0;k>t[k];

        fin.close();

}

void splin1::interp()

{

        int i,j;

        double h0,h1,alpha,beta;

        s[0]=dy[0];dy[0]=0.0;

        h0=x[1]-x[0];

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

        {

                h1=x[j+1]-x[j];

                alpha=h0/(h0+h1);

                beta=(1.0-alpha)*(y[j]-y[j-1])/h0;

                beta=3.0*(beta+alpha*(y[j+1]-y[j])/h1);

                dy[j]=-alpha/(2.0+(1.0-alpha)*dy[j-1]);

                s[j]=(beta-(1.0-alpha)*s[j-1]);

                s[j]=s[j]/(2.0+(1.0-alpha)*dy[j-1]);

                h0=h1;

        }

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

                dy[j]=dy[j]*dy[j+1]+s[j];

        for(j=0;j<=n-2;j++)s[j]=x[j+1]-x[j];

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

        {

                h1=s[j]*s[j];

                ddy[j]=6.0*(y[j+1]-y[j])/h1-2.0*(2.0*dy[j]+dy[j+1])/s[j];

        }

        h1=s[n-2]*s[n-2];

        ddy[n-1]=6.*(y[n-2]-y[n-1])/h1+2.*(2.*dy[n-1]+dy[n-1])/s[n-2];

        integ=0.0;

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

        {

                h1=0.5*s[i]*(y[i]+y[i+1]);

                h1=h1-s[i]*s[i]*s[i]*(ddy[i]+ddy[i+1])/24.0;

                integ=integ+h1;

        }

        for(j=0;j<=m-1;j++)

        {

                if(t[j]>=x[n-1])i=n-2;

                else

                {

                        i=0;

                        while(t[j]>x[j+1])i=i+1;

                }

                h1=(x[j+1]-t[j])/s[i];

                h0=h1*h1;

                z[j]=(3.0*h0-2.0*h0*h1)*y[i];

                z[j]=z[j]+s[i]*(h0-h0*h1)*dy[i];

                dz[j]=6.0*(h0-h1)*y[i]/s[i];

                dz[j]=dz[j]+(3.0*h0-2.0*h1)*dy[i];

                ddz[j]=(6.0-12.0*h1)*y[i]/(s[i]*s[i]);

                ddz[j]=ddz[j]+(2.0-6.0*h1)*dy[i]/s[i];

                h1=(t[j]-x[i])/s[i];

                h0=h1*h1;

                z[j]=z[j]+(3.0*h0-2.0*h0*h1)*y[i+1];

                z[j]=z[j]-s[i]*(h0-h0*h1)*dy[i+1];

                dz[j]=dz[j]-6.0*(h0-h1)*y[i+1]/s[i];

                dz[j]=dz[j]+(3.0*h0-2.0*h1)*dy[i+1];

                ddz[j]=ddz[j]+(6.0-12.0*h1)*y[i+1]/(s[i]*s[i]);

                ddz[j]=ddz[j]-(2.0-6.0*h1)*dy[i+1]/s[i];

        }

}

void splin1::output()

{

        int i;

        char str2[20];

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

        cin>>str2;

        ofstream fout(str2);

        if(!fout)

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

                fout<
                    <
                cout<
                        <
        }

        fout<
        cout<
        <
        for(i=0;i
        {

                fout<
                <
                cout<
                <
        }

        fout<
        fout<
        fout.close();

}

int main()

{

        splin1 solution(12,8);

        solution.input();

        solution.interp();

        solution.output();

}

輸入文件為：
0.52  5.28794
8.0 13.8400
17.95 20.2000
28.65 24.9000
50.65 31.1000
104.6 36.5000
156.6 36.6000
260.7 31.0000
364.4 20.9000
468.0 7.80000
507.0 1.50000
520.0 0.200000

1.86548 -0.046115

4.0  14.0  30.0  60.0  130.0  230.0 450.0  515.0
結(jié)果出現(xiàn)這個(gè)問題：[Session started at 2011-04-06 17:46:04 +0800.]
GNU gdb 6.3.50-20050815 (Apple version gdb-1469) (Wed May  5 04:36:56 UTC 2010)
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB.  Type "show warranty" for details.
This GDB was configured as "x86_64-apple-darwin".tty /dev/ttys000
Loading program into debugger…
Program loaded.
run
[Switching to process 1601]
Running…

輸入文件名:splin1.txt

不能打開這個(gè)文件splin1.txt

Debugger stopped.
Program exited with status value:1.(gdb)

[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 354求調(diào)劑 +6	Tyoumou 2026-03-18	9/450	2026-03-21 20:47 by lbsjt
[考研] 333求調(diào)劑 +5	87639 2026-03-21	7/350	2026-03-21 19:31 by ColorlessPI
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 350求調(diào)劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao

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

[交流] 【求助】程序運(yùn)行已經(jīng)無錯(cuò)誤，但是輸入文件有問題，求教。

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

我這里很正常啊。你單獨(dú)跟蹤一下，看看那里出問題了。附圖：

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

我這里很正常啊。你單獨(dú)跟蹤一下，看看那里出問題了。附圖：