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

zyj8119

木蟲 (著名寫手)

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

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

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è)問(wè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

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

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

回復(fù)此樓

» 猜你喜歡

0854電子信息求調(diào)劑已經(jīng)有7人回復(fù)
0854人工智能方向招收調(diào)劑已經(jīng)有3人回復(fù)
調(diào)劑已經(jīng)有4人回復(fù)
一志愿北化315 求調(diào)劑已經(jīng)有3人回復(fù)
有機(jī)合成求助已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有7人回復(fù)
收08調(diào)劑生已經(jīng)有6人回復(fù)
材料專碩找調(diào)劑已經(jīng)有5人回復(fù)
085602 289分求調(diào)劑已經(jīng)有5人回復(fù)

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

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

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

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

查看全部散金貼

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

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

zyj8119

木蟲 (著名寫手)

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

引用回帖:

Originally posted by smart0193 at 2011-04-07 20:18:40:
沒(méi)有，很快的，運(yùn)行，輸入文件名，直接就出結(jié)果了。
突然想起來(lái)，是不是防火墻的問(wèn)題呢？限制了程序?qū)ξ募脑L問(wèn)？或者通過(guò)程序來(lái)建立文件？

我檢查了服務(wù)里面的security center,是禁用狀態(tài)。。。

贊一下

回復(fù)此樓

9樓2011-04-07 22:05:19

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

查看全部 24 個(gè)回答

zyj8119

木蟲 (著名寫手)

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

引用回帖:

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ù)文件放在一個(gè)文件夾下面。。。

贊一下

回復(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.7小時(shí)
蟲號(hào): 1120220

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

我用VC6運(yùn)行了你的程序，采用的也是你給出的例子文件，沒(méi)有問(wè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

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

贊一下(2人)

回復(fù)此樓

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

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

zyj8119

木蟲 (著名寫手)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 求調(diào)劑 +5	林之夕 2026-03-24	5/250	2026-03-24 17:16 by dick_runner
[考研] 081700 調(diào)劑 267分 +9	迷人的哈哈 2026-03-23	9/450	2026-03-24 11:58 by 544594351
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 一志愿國(guó)科過(guò)程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5	jiajunX 2026-03-22	5/250	2026-03-23 18:07 by YMU施老師
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 311求調(diào)劑 +6	冬十三 2026-03-18	6/300	2026-03-22 20:18 by edmund7
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無(wú)際的草原
[考研] 求助 +5	夢(mèng)里的無(wú)言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438