gielcheen

Constant n = 10000;
Function y=(6/pi^2)*(Sum(i=1:n)((1/pi^2)*exp(-n^2*pi^2*a*x/0.000001));
Data ;
0 0.98501
0.21 29.86622
0.42 32.67416
1 36.97943
4 43.80644
7 48.87115
10 51.80471
14 53.65742
18 57.08069
22 59.32729
26 60.3525
30 61.04827
35 62.75492
42 65.76456
49 68.92578
60 74.15131
這是我寫(xiě)的程序，運(yùn)行不出來(lái)，請(qǐng)問(wèn)有什么錯(cuò)誤嗎？

贊一下

回復(fù)此樓

13樓2012-02-27 20:31:04

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

gielcheen(金幣+10): ★有幫助 2012-02-28 14:04:02
小雨萌萌(金幣+2): 3Q 2012-02-29 09:44:41

1：左右括號(hào)不匹配，改正后可以運(yùn)行
2：僅有一個(gè)參數(shù)，擬合效果不好，能確認(rèn)公式正確嗎？n=10000時(shí)僅計(jì)算一次化多少時(shí)間？簡(jiǎn)化下吧，或改成積分形式。

贊一下

回復(fù)此樓

14樓2012-02-28 11:57:44

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

gielcheen

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 16 (小學(xué)生)
金幣: 8255
散金: 758
紅花: 15
帖子: 676
在線(xiàn): 560.1小時(shí)
蟲(chóng)號(hào): 883648
注冊(cè): 2009-10-25
專(zhuān)業(yè): 生物醫(yī)用高分子材料

引用回帖:

10樓: Originally posted by gielcheen at 2012-02-27 20:24:38:
0 0.98501
0.21 29.86622
0.42 32.67416
1 36.97943
4 43.80644
7 48.87115
10 51.80471
14 53.65742
18 57.08069
22 59.32729
26 60.3525
30 61.04827
35 62.75492
42 65.76456
49 68.92578
60 ...

Parameter a;
Variable x,y;
Function y=1-(6/pi^2)*(Sum(n=1:999)((1/pi^2)*exp(-n^2*pi^2*a*x/0.000001)));
data;
0 0.98501
0.21 29.86622
0.42 32.67416
1 36.97943
4 43.80644
7 48.87115
10 51.80471
14 53.65742
18 57.08069
22 59.32729
26 60.3525
30 61.04827
35 62.75492
42 65.76456
49 68.92578
60 74.15131
我也沒(méi)辦法簡(jiǎn)化了，n=999已經(jīng)是極限了，n=1000時(shí)就報(bào)告說(shuō)函數(shù)出錯(cuò)，您看您能幫我寫(xiě)一下程序嗎？

贊一下

回復(fù)此樓

15樓2012-02-28 14:03:40

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

gielcheen(金幣+50): ★★★很有幫助非常感謝，但是編程沒(méi)有看懂。況且n是到正無(wú)窮，如果只到1200，還是不太好！ 2012-02-28 16:44:11

用下列編程模式代碼(n= 1200)：

CODE:

Parameter a;

StartRange = [1E-100,1e-10];

Variable  x,y;

StartProgram [Pascal];

Procedure MainModel;

var i, n: integer;

    td: double;

Begin

    for i := 0 to DataLength - 1 do begin

        td := 0;

        for n := 1 to 1200 do

            td := td + exp(-n^2*pi^2*a*x[i]/0.000001);

        y[i] := 1/pi^2*td;

    end;

End;

EndProgram;

data;

0        0.98501

0.21        29.86622

0.42        32.67416

1        36.97943

4        43.80644

7        48.87115

10        51.80471

14        53.65742

18        57.08069

22        59.32729

26        60.3525

30        61.04827

35        62.75492

42        65.76456

49        68.92578

60        74.15131

參數(shù) 最佳估算
-------------------- -------------
a 1.75957508897304E-14

====== 輸出結(jié)果 =====

文件：數(shù)據(jù)文件 - 1
No 目標(biāo)y 計(jì)算y
1 0.98501 121.585420370805
2 29.86622 119.487562839931
3 32.67416 117.454445286095
4 36.97943 112.156261269926
5 43.80644 90.7639936726885
6 48.87115 76.4042344209492
7 51.80471 66.3661814206772
8 53.65742 57.0694131575808
9 57.08069 50.6002923996484
10 59.32729 45.8466101914636
11 60.3525 42.1938718181288
12 61.04827 39.2848251790096
13 62.75492 36.3697677962097
14 65.76456 33.1972842183255
15 68.92578 30.7310720717936
16 74.15131 27.7666907960382

贊一下

回復(fù)此樓

16樓2012-02-28 14:55:18

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

小雨萌萌(金幣+3): 3Q 2012-02-29 09:45:00

你可以把代碼中的1200換成12000或更多，只是太化時(shí)間。代碼已給出了自己慢慢算吧。

贊一下

回復(fù)此樓

17樓2012-02-28 20:15:44

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

gielcheen

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 16 (小學(xué)生)
金幣: 8255
散金: 758
紅花: 15
帖子: 676
在線(xiàn): 560.1小時(shí)
蟲(chóng)號(hào): 883648
注冊(cè): 2009-10-25
專(zhuān)業(yè): 生物醫(yī)用高分子材料

引用回帖:

樓: Originally posted by dingd at 2012-02-28 20:15:44:
你可以把代碼中的1200換成12000或更多，只是太化時(shí)間。代碼已給出了自己慢慢算吧。

能不用數(shù)代替，用正無(wú)窮算嗎？

贊一下

回復(fù)此樓

18樓2012-02-29 08:11:03

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

gielcheen(金幣+20): ★★★很有幫助 2012-02-29 09:15:57

“正無(wú)窮”只是個(gè)概念上的東西，數(shù)值計(jì)算時(shí)必須有具體值的。要是你的公式隨著n的增加而收斂（趨于穩(wěn)定），此時(shí)“正無(wú)窮”就可以用這個(gè)n值代替，否則如果發(fā)散說(shuō)明你的公式有問(wèn)題。

贊一下

回復(fù)此樓

19樓2012-02-29 08:33:27

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

gielcheen

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 16 (小學(xué)生)
金幣: 8255
散金: 758
紅花: 15
帖子: 676
在線(xiàn): 560.1小時(shí)
蟲(chóng)號(hào): 883648
注冊(cè): 2009-10-25
專(zhuān)業(yè): 生物醫(yī)用高分子材料

引用回帖:

樓: Originally posted by dingd at 2012-02-29 08:33:27:
“正無(wú)窮”只是個(gè)概念上的東西，數(shù)值計(jì)算時(shí)必須有具體值的。要是你的公式隨著n的增加而收斂（趨于穩(wěn)定），此時(shí)“正無(wú)窮”就可以用這個(gè)n值代替，否則如果發(fā)散說(shuō)明你的公式有問(wèn)題。

我把你給我的程序復(fù)制到1stopt里面了，但是無(wú)法計(jì)算，點(diǎn)了執(zhí)行計(jì)算毫無(wú)反應(yīng)。還有，能解釋一下各行的意思嗎？我在這方面是小白！

贊一下

回復(fù)此樓

20樓2012-02-29 09:15:49

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

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

[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 材料與化工專(zhuān)碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 0856專(zhuān)碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料專(zhuān)碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky