當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

zyj8119

木蟲 (著名寫手)

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

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不能打開這個文件"<
        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不能打開這個文件"<
        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)這個問題：[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

不能打開這個文件splin1.txt

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

回復(fù)此樓

» 猜你喜歡

347求調(diào)劑已經(jīng)有4人回復(fù)
一志愿西交288化工專碩求調(diào)劑已經(jīng)有3人回復(fù)
一志愿中國科學(xué)院大學(xué)265求調(diào)劑已經(jīng)有4人回復(fù)
291求調(diào)劑已經(jīng)有13人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
322求調(diào)劑已經(jīng)有7人回復(fù)
0703化學(xué) 已經(jīng)有15人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有15人回復(fù)
313求調(diào)劑已經(jīng)有3人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有11人回復(fù)

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

求助vasp運(yùn)行過程中出現(xiàn)錯誤提示： WARNING: CNORMN: search vector ill defined 已經(jīng)有8人回復(fù)
求助matlab一程序運(yùn)行的問題，幫忙看看哪里不對已經(jīng)有4人回復(fù)
Perl 生物信息學(xué)編程求助，運(yùn)行出現(xiàn)問題，煩請幫忙查找一下錯誤，謝謝！已經(jīng)有3人回復(fù)
matlab程序運(yùn)行總是出現(xiàn)下面的錯誤，請各位高手指點(diǎn)一下，謝謝了！急急急已經(jīng)有5人回復(fù)
fortran 程序運(yùn)行中的錯誤已經(jīng)有5人回復(fù)
用opengl做了個工程，編譯無錯誤，但是運(yùn)行不起來。。。已經(jīng)有11人回復(fù)
【求助】運(yùn)行fortran程序出現(xiàn)錯誤 : error LNK2019 怎么回事啊已經(jīng)有8人回復(fù)
【求助】此程序沒有錯誤，怎么運(yùn)行不出來？已經(jīng)有6人回復(fù)
【求助】vc6.0運(yùn)行一個exe文件對圖像進(jìn)行canny邊緣檢測時就會出現(xiàn)以下錯誤，請教如何?/a> 已經(jīng)有5人回復(fù)

【求助】程序運(yùn)行但沒達(dá)到目的？已經(jīng)有10人回復(fù)
【求助】此程序也是沒有問題了，但是依然運(yùn)行不出結(jié)果已經(jīng)有7人回復(fù)
【求助】gaussian03運(yùn)行程序出現(xiàn)拒絕訪問的錯誤已經(jīng)有5人回復(fù)
【求助】此FORTRAN90程序沒有錯誤，但是運(yùn)行不出來。已經(jīng)有4人回復(fù)
【求助】程序運(yùn)行后出現(xiàn)的錯誤 access violation 已經(jīng)有15人回復(fù)
【求助】為什么我這個程序編譯都沒有問題，就是運(yùn)行不出來？【已完結(jié)】已經(jīng)有9人回復(fù)
【求助】安裝ansys時，運(yùn)行Keygen程序,但不會自動生成一個文件License.dat,為什么? 已經(jīng)有12人回復(fù)

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

查看全部散金貼

1樓 2011-04-06 20:39:34

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

smart0193

木蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 3224.1
帖子: 200
在線: 76.8小時
蟲號: 1120220

沒有，很快的，運(yùn)行，輸入文件名，直接就出結(jié)果了。
突然想起來，是不是防火墻的問題呢？限制了程序?qū)ξ募脑L問？或者通過程序來建立文件？

贊一下

回復(fù)此樓

7樓2011-04-07 20:18:40

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

查看全部 24 個回答

zyj8119

木蟲 (著名寫手)

引用回帖:

Originally posted by zyj8119 at 2011-04-06 20:39:34:
[code]#include<iostream>
#include<iomanip>
#include<fstream>
#include<cmath>
using namespace std;
class splin1
{
private:
int n,m;
double *x,*y,*dy,*ddy,*t,*z,*dz ...

我已經(jīng)把程序與數(shù)據(jù)文件放在一個文件夾下面。。。

贊一下

回復(fù)此樓

2樓2011-04-06 20:44:24

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

smart0193

木蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 3224.1
帖子: 200
在線: 76.8小時
蟲號: 1120220

★ ★ ★
zyj8119(金幣+1):謝謝參與
微塵、夢想(金幣+2): 謝謝參與應(yīng)助…… 2011-04-07 17:26:38

我用VC6運(yùn)行了你的程序，采用的也是你給出的例子文件，沒有問題，可以輸出結(jié)果。如下：
         0.52       5.28794 -6.27744e+066 2.84681e+066
            8          13.84 1.90782e+066 -6.58235e+065
         17.95          20.2 -5.40918e+065 1.66027e+065
         28.65          24.9 1.93591e+065 -2.87356e+064
         50.65          31.1 -7.10913e+064    4.6735e+063
         104.6          36.5 1.61149e+064 -1.44065e+063
         156.6          36.6 5.22698e+063 1.02188e+063
         260.7             31 -6.36429e+064 -2.34503e+063
         364.4          20.9 2.48875e+065 8.37238e+063
         468          7.8 -9.3144e+065 -3.11584e+064
         507          1.5 2.47047e+066 2.05615e+065
         520          0.2 -6.27744e+066 -2.89728e+066

-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132
-6.27744e+066 1.93191e+265 -9.23265e+198 2.94154e+132

6.20039e+068

你的問題是不是由于你的輸入文件還在打開著？在執(zhí)行程序的時候，要關(guān)閉它

贊一下(2人)

回復(fù)此樓

3樓2011-04-07 09:37:51

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

zyj8119

木蟲 (著名寫手)

引用回帖:

Originally posted by smart0193 at 2011-04-07 09:37:51:
我用VC6運(yùn)行了你的程序，采用的也是你給出的例子文件，沒有問題，可以輸出結(jié)果。如下：
0.52 5.28794 -6.27744e+066 2.84681e+066
8 13.84 1.90782e+066 ...

我原來用的是 xcode,現(xiàn)在用VC6.0編譯了，但是還是有問題，運(yùn)行到這一步就停止了：

贊一下

回復(fù)此樓

4樓2011-04-07 13:37:23

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

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

[考研] 材料學(xué)碩333求調(diào)劑 +15	北道巷 2026-03-24	15/750	2026-03-31 08:26 by yanflower7133
[考研] 085404 22408 315分 +3	zhuangyan123 2026-03-31	3/150	2026-03-31 07:57 by 機(jī)械調(diào)劑-北石化
[考研] 0856調(diào)劑 +3	曲聽筠 2026-03-30	3/150	2026-03-31 07:29 by 無際的草原
[考研] 一志愿大連理工大學(xué)材料求調(diào)劑 +6	Gymno 2026-03-30	6/300	2026-03-31 07:26 by 無際的草原
[考研] 材料科學(xué)與工程調(diào)劑 +10	深V宿舍吧 2026-03-30	11/550	2026-03-31 04:12 by fmesaito
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 21:07 by peike
[考研] 285求調(diào)劑 +6	AZMK 2026-03-29	9/450	2026-03-30 21:02 by dophin1985
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9	hanamiko 2026-03-27	9/450	2026-03-30 20:10 by 無際的草原
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +5	小小機(jī)靈蟲 2026-03-29	11/550	2026-03-30 15:02 by Wang200018
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +6	蘆lty 2026-03-25	7/350	2026-03-28 13:13 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671