zyj8119

木蟲 (著名寫手)

[交流] 【求助】程序運行已經(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ù)此樓

» 猜你喜歡

數(shù)一英一271專碩（085401）求調(diào)劑，可跨已經(jīng)有8人回復(fù)
295求調(diào)劑已經(jīng)有4人回復(fù)
0856求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有10人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑已經(jīng)有6人回復(fù)
289求調(diào)劑已經(jīng)有12人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有10人回復(fù)
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩已經(jīng)有5人回復(fù)

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

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

【求助】程序運行但沒達到目的？已經(jīng)有10人回復(fù)
【求助】此程序也是沒有問題了，但是依然運行不出結(jié)果已經(jīng)有7人回復(fù)
【求助】gaussian03運行程序出現(xiàn)拒絕訪問的錯誤已經(jīng)有5人回復(fù)
【求助】此FORTRAN90程序沒有錯誤，但是運行不出來。已經(jīng)有4人回復(fù)
【求助】程序運行后出現(xiàn)的錯誤 access violation 已經(jīng)有15人回復(fù)
【求助】為什么我這個程序編譯都沒有問題，就是運行不出來？【已完結(jié)】已經(jīng)有9人回復(fù)
【求助】安裝ansys時，運行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

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

我用VC6運行了你的程序，采用的也是你給出的例子文件，沒有問題，可以輸出結(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的回帖

查看全部 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的回帖

zyj8119

木蟲 (著名寫手)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

smart0193

木蟲 (小有名氣)

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

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

zyj8119(金幣+5): 謝謝提醒。。。 2011-04-07 16:55:09

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

贊一下

回復(fù)此樓

5樓2011-04-07 16:01:45

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

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

[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +3	大饞小子 2026-03-28	3/150	2026-03-28 20:10 by wxiongid
[考研] 一志愿華北電力大學(xué)能動專碩，293，求調(diào)劑 +3	15537177284 2026-03-23	5/250	2026-03-28 16:11 by xxxsssccc
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +6	nnnnnnn5 2026-03-25	6/300	2026-03-28 13:33 by 果果媽咪
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +4	Zuhui0306 2026-03-25	4/200	2026-03-28 13:07 by 唐沐兒
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂呵呵的追夢人
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4	傾____蕭 2026-03-27	5/250	2026-03-27 11:52 by zhshch
[碩博家園] 招收生物學(xué)/細胞生物學(xué)調(diào)劑 +3	IceGuo 2026-03-26	4/200	2026-03-27 05:35 by user003
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 282求調(diào)劑 +3	wcq131415 2026-03-24	3/150	2026-03-25 12:16 by userper
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂