215215

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
散金: 134
紅花: 1
帖子: 125
在線: 53.5小時
蟲號: 860062
注冊: 2009-09-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】程序修改！已有4人參與

PROGRAM PENDULUM
!
! Main Program for a driven pendulum under damping solved with
! the fourth-order Runge-Kutta algorithm.
! Parameters: Q, B,and W (omega_0).
!
   PARAMETER (N=1000,L=100,M=1)
   DIMENSION  Y(2,N)
      COMMON/CONST/Q,B,W
!
  PI = 4.0*ATAN(1.0)
  H  = 3.0*PI/L
  Q  = 0.5
  B  = 0.9
  W  = 2.0/3.0
  Y(1,1) = 0.0
  Y(2,1) = 2.0
!
! Using the Runge-Kutta algorithm to integrate the
! equation
!
  DO 100 I = 1, N-1
T  = H*I
Y1 = Y(1,I)
Y2 = Y(2,I)
DK11 = H*G1(Y1,Y2,T)
DK21 = H*G2(Y1,Y2,T)
DK12 = H*G1((Y1+DK11/2.0),(Y2+DK21/2.0),(T+H/2.0))
DK22 = H*G2((Y1+DK11/2.0),(Y2+DK21/2.0),(T+H/2.0))
DK13 = H*G1((Y1+DK12/2.0),(Y2+DK22/2.0),(T+H/2.0))
DK23 = H*G2((Y1+DK12/2.0),(Y2+DK22/2.0),(T+H/2.0))
DK14 = H*G1((Y1+DK13),(Y2+DK23),(T+H))
DK24 = H*G2((Y1+DK13),(Y2+DK23),(T+H))
Y(1,I+1) = Y(1,I)+(DK11+2.0*(DK12+DK13)+DK14)/6.0
Y(2,I+1) = Y(2,I)+(DK21+2.0*(DK22+DK23)+DK24)/6.0
!
! Bring theta back to the region [-pi,pi]
!
IF(ABS(Y(1,I+1)).GT.PI) THEN
   Y(1,I+1) = Y(1,I+1)-2.0*PI*ABS(Y(1,I+1))/Y(1,I+1)
      END IF
100 CONTINUE
!
WRITE(6,999)(Y(1,I),Y(2,I),I=1,N,M)
STOP
  999 FORMAT (2F16.8)
   END
!
  FUNCTION G1 (Y1,Y2,T)
  COMMON /CONST/Q,B,W
   G1=Y2
  RETURN
  END
!
FUNCTION G2 (Y1,Y2,T)
COMMON /CONST/Q,B,W
  G2 = -Q*Y2-SIN(Y1)+B*COS(W*T)
RETURN
END
————————————————————————————————————

以上程序如果把B=0.9修改為B在區(qū)間【1，1.5】，要想繼續(xù)運行的話該怎么修改此程序��？？？？？請賜教！

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

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

寫了一個fortran90的小程序，編譯通不過，請大俠幫忙已經(jīng)有59人回復(fù)
【求助】求高人幫我修改程序已經(jīng)有5人回復(fù)
【求助】幫忙改寫展寬程序已經(jīng)有11人回復(fù)
【原創(chuàng)】vasp的POSCAR轉(zhuǎn)換成cif格式程序修改版已經(jīng)有7人回復(fù)
【求助】修改MATLAB支持向量機核函數(shù)程序遇到的疑問？求助已經(jīng)有5人回復(fù)

nevergiveup!

1樓 2010-04-21 18:25:53

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
jjdg(金幣+2):辛苦了 2010-04-26 17:29

不是很明白你什么意思，這個程序中的參數(shù)包括了B，你修改了，重新編譯就可以了，這個會有什么問題嗎？

贊一下(2人)

回復(fù)此樓

2樓2010-04-26 14:36:37

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

holmescn

金蟲 (正式寫手)

程序強帖: 37
應(yīng)助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
jjdg(金幣+2):辛苦了 2010-04-26 17:29
wangen994(金幣+2):活動期間額外獎勵 2010-04-27 09:56

把Program到END中間的代碼放到一個subroutine里，這個subrutine要有一個參數(shù)B。然后在主程序里用do循環(huán)調(diào)用。

具體實現(xiàn)請查閱Fortran相關(guān)書籍。

贊一下(3人)

回復(fù)此樓

3樓2010-04-26 17:06:38

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

215215

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
散金: 134
紅花: 1
帖子: 125
在線: 53.5小時
蟲號: 860062
注冊: 2009-09-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by snoopyzhao at 2010-04-26 14:36:37:
不是很明白你什么意思，這個程序中的參數(shù)包括了B，你修改了，重新編譯就可以了，這個會有什么問題嗎？

我就是想運行一下B屬于【1，1.5】時的程序，上邊是B=0.9，就是怎么修改此程序改成是運行B屬于【1，1.5】時的程序，因為我不懂FORTRAN，希望高手能給我修改編譯好此程序讓我運行一下而已，謝謝�。�！

贊一下

回復(fù)此樓

nevergiveup!

4樓2010-04-27 13:39:01

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

215215

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
散金: 134
紅花: 1
帖子: 125
在線: 53.5小時
蟲號: 860062
注冊: 2009-09-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

jjdg:你要多少步長呢？ 2010-04-27 16:35

引用回帖:

Originally posted by holmescn at 2010-04-26 17:06:38:
把Program到END中間的代碼放到一個subroutine里，這個subrutine要有一個參數(shù)B。然后在主程序里用do循環(huán)調(diào)用。

具體實現(xiàn)請查閱Fortran相關(guān)書籍。

我不懂FORTRAN，希望高手能幫我調(diào)整好此程序讓我運行一下知道結(jié)果而已（就是把b=0.9修改為b屬于區(qū)間【1，1.5】運行出結(jié)果），謝謝！

贊一下(1人)

回復(fù)此樓

nevergiveup!

5樓2010-04-27 13:41:28

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
jjdg(金幣+1):感謝交流！ 2010-04-27 16:35

引用回帖:

（就是把b=0.9修改為b屬于區(qū)間【1，1.5】運行出結(jié)果）

修改為多少？1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 還是多少？

贊一下(2人)

回復(fù)此樓

6樓2010-04-27 15:29:57

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

215215

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
散金: 134
紅花: 1
帖子: 125
在線: 53.5小時
蟲號: 860062
注冊: 2009-09-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by snoopyzhao at 2010-04-27 15:29:57:

修改為多少？1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 還是多少？

ok,我把原話給你發(fā)上來你看吧，
Modify the program to obtain a bifurcation diagram(w-b plot at a fixed “西塔角” )with q=0.5and w=2/3 and in the region of b屬于【1，1.5】。

贊一下

回復(fù)此樓

nevergiveup!

7樓2010-04-28 12:22:18

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
resonant(金幣+1):感謝交流：-） 2010-04-28 19:57

你只要修改 b 的值，比如 b = 1.2，然后編譯程序，運行之就好了，這個似乎沒有什么難度吧。

如果你不喜歡每改一個值都編譯，那么就像 holmescn 說的那樣，把它改成一個帶參數(shù)的，呵呵……

贊一下(2人)

回復(fù)此樓

8樓2010-04-28 15:28:26

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

215215

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 388.2
散金: 134
紅花: 1
帖子: 125
在線: 53.5小時
蟲號: 860062
注冊: 2009-09-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by snoopyzhao at 2010-04-28 15:28:26:
你只要修改 b 的值，比如 b = 1.2，然后編譯程序，運行之就好了，這個似乎沒有什么難度吧。

如果你不喜歡每改一個值都編譯，那么就像 holmescn 說的那樣，把它改成一個帶參數(shù)的，呵呵……

b這里不是一個確定的值，呵呵，我也不清楚怎么回事！呵呵呵

贊一下

回復(fù)此樓

nevergiveup!

9樓2010-04-28 22:07:10

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

這讓人怎么改呢？呵呵……

贊一下(1人)

回復(fù)此樓

10樓2010-04-29 11:24:01

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

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

[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 生物學(xué)調(diào)劑招人！�。� +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客