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

zyj8119

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲(chóng)號(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不能打開(kāi)這個(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不能打開(kāi)這個(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

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

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

回復(fù)此樓

» 猜你喜歡

環(huán)境工程 085701，267求調(diào)劑已經(jīng)有5人回復(fù)
275求調(diào)劑已經(jīng)有17人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有15人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有9人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有10人回復(fù)
南京大學(xué)化學(xué)調(diào)劑已經(jīng)有9人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有13人回復(fù)
考研調(diào)劑已經(jīng)有6人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有13人回復(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

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲(chóng)號(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的回帖

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

zyj8119

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲(chóng)號(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

木蟲(chóng) (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 3224.1
帖子: 200
在線: 76.8小時(shí)
蟲(chóng)號(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)題是不是由于你的輸入文件還在打開(kāi)著？在執(zhí)行程序的時(shí)候，要關(guān)閉它

贊一下(2人)

回復(fù)此樓

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

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

smart0193

木蟲(chóng) (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 3224.1
帖子: 200
在線: 76.8小時(shí)
蟲(chóng)號(hào): 1120220

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

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +7	WW.' 2026-03-23	9/450	2026-03-29 12:42 by fmesaito
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落睿可思
[考研] 320分，材料與化工專(zhuān)業(yè)，求調(diào)劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 一志愿廈門(mén)大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 322求調(diào)劑 +5	舊吢 2026-03-24	5/250	2026-03-28 13:26 by Iveryant
[考研] 277跪求調(diào)劑 +5	1915668 2026-03-27	9/450	2026-03-28 09:58 by zhshch
[考研] 359求調(diào)劑 +4	王了個(gè)楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚(yú)的貓
[考研] 336材料求調(diào)劑 +7	陳瀅瑩 2026-03-26	9/450	2026-03-27 00:20 by wxiongid
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 材料考研求調(diào)劑 +3	Dendel 2026-03-23	6/300	2026-03-26 17:51 by fmesaito
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 26考研-291分-廈門(mén)大學(xué)（085601）-柔性電子學(xué)院材料工程專(zhuān)業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 296求調(diào)劑 +4	汪��？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。�！
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 材料專(zhuān)碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060