 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 643 | 回復: 1 | ||
[求助]
求幫助！改Function中的可變系數(shù)進行循環(huán)計算的問題
|
|
沒有金幣了啊�。�！希望發(fā)揚雷鋒精神啊！求大神幫忙解決問題��！ 我想用Fortran 求解非線性方程組，程序是Monte Carlo法 程序由主程序，F(xiàn)unction, Subroutine三部分組成。 Function給出了非線性方程組，通過subroutine進行計算非線性方程組，主程序給出了初始計算值X1,X2,X3。 問題是， Function 中有個系數(shù)M是變化的， 在subroutine（或者主程序）中，先設(shè)定了初始值M0，每次計算結(jié)束了X后，由X計算新的M=f(x). 如果新Mc 和舊M差別在1.0E-05內(nèi)就結(jié)束計算， 否則重新賦予Function 系數(shù)M, 重新調(diào)用子程序循環(huán)計算X. DIMENSION X(N), Y( N) REAL NRND1 External F Double precision R, M, M0, Mc Double precision R=1.0E-4; d=30.0E-3; k=383.0; * h=40000.0; T=303.0 A= B K= 1 R= 1.0E0 M0=450.0 40 M=M0 Z= F(X, N, M) 10 IF (A. GE. EPS) THEN L= L+ 1 DO 20 I= 1, N 20 Y(I) = -A+ 2.0* A* NRND1(R) + X(I) Z1= F(Y, N, M) K= K+ 1 IF ( Z1. GE. Z) THEN IF (K .GT. M) THEN K= 1 A= A/2.0 ENDIF GOTO 10 ELSE K= 1 DO 30 I= 1, N 30 X(I) = Y( I) Z= Z1 IF ( Z. GE. EPS) GOTO 10 ENDIF ENDIF Mc=T+(1000.0*X(1)-M)*(d/k+1/h)/R IF(ABS(M-Mc) .GE. EPS) THEN M0=(M+Mc)/2 GOTO 40 ENDIF END 上面程序是對subroutine的修改。倒數(shù)第6行直接調(diào)用X1對新的Mc計算，并賦值給M0，然后返回行40。 但是貌似根本就沒有把新的M傳遞給方程F(X,N,M)進行循環(huán)計算。最后只有第一次的結(jié)果。 難道計算得到的X1不能直接調(diào)用，用于方程計算么？(嘗試過把X1直接調(diào)出主程序計算倒數(shù)第6行的方程，方程結(jié)果根本就不對) 此子程序有否錯誤呢？ 謝謝！ |
回復此樓» 猜你喜歡
303求調(diào)劑
已經(jīng)有6人回復
-
306求調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑
已經(jīng)有6人回復
-
275求調(diào)劑
已經(jīng)有5人回復
-
求調(diào)劑
已經(jīng)有7人回復
-
考研調(diào)劑
已經(jīng)有4人回復
-
285求調(diào)劑
已經(jīng)有6人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有7人回復
-
318求調(diào)劑
已經(jīng)有6人回復
-
305分求調(diào)劑（食品工程）
已經(jīng)有5人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【求助】怎么用spss計算兩個變量間的相關(guān)系數(shù)
已經(jīng)有5人回復
- 利用matlab求一階變系數(shù)微分方程組的數(shù)值解。
已經(jīng)有7人回復
- 求助�。徃迦颂岢龅年P(guān)于檢測限的問題，不懂是什么意思，求幫助。。。
已經(jīng)有14人回復
- 變系數(shù)偏微分方程matlab求解
已經(jīng)有6人回復
- 求回歸方程參數(shù)的最優(yōu)化算法matlab代碼
已經(jīng)有17人回復
- 用1stOpt進行復數(shù)曲線擬合時，結(jié)果出錯。急求教�。�！
已經(jīng)有5人回復
- 1stopt 擬合隱函數(shù)的四參數(shù)
已經(jīng)有7人回復
- 麻煩大家?guī)兔纯催@個審稿人什么意思（關(guān)于數(shù)據(jù)分析和”化學計量系數(shù)？“）
已經(jīng)有5人回復
- 哪位蟲友幫幫忙，我的comsol計算時出現(xiàn)如下錯誤
已經(jīng)有9人回復
- 兩條不同的氮氣吸附脫附曲線，求幫助~~
已經(jīng)有18人回復
- 關(guān)于用化學吸附儀做化學脈沖實驗計算金屬分散度的一些問題！求幫助�。�
已經(jīng)有11人回復
- 【求助】FLUENT UDF 氣液兩相流 汽蝕 汽化潛熱(問題詳見內(nèi)容)
已經(jīng)有21人回復
- 請問，參比制劑溶出曲線的變異系數(shù)偏大，是怎么回事��？
已經(jīng)有9人回復
- matlab中系數(shù)是矩陣的積分
已經(jīng)有5人回復
- 如何計算變位齒輪的齒廓重疊干涉系數(shù)
已經(jīng)有4人回復
- 【求助】迎風角變化下對長方體的風阻系數(shù)的關(guān)系，請求幫助！
已經(jīng)有7人回復
- 【求助/交流】生物統(tǒng)計學問題，變異系數(shù)
已經(jīng)有11人回復
- 【求助】MUSIC的MD模塊計算甲烷在IRMOF-1上的擴散系數(shù)
已經(jīng)有11人回復
- 【求助】MATLAB計算NRTL活度系數(shù)
已經(jīng)有7人回復
- 【求助】常微分方程組求解中系數(shù)與某變量值關(guān)聯(lián)的問題
已經(jīng)有12人回復
- 【求助】如何從粒度分布曲線上求算變異系數(shù)CV值？
已經(jīng)有4人回復
學蟲dj
金蟲 (正式寫手)
- 應(yīng)助: 18 (小學生)
- 金幣: 892.4
- 散金: 63
- 帖子: 688
- 在線: 89.9小時
- 蟲號: 1829966
- 注冊: 2012-05-22
- 性別: GG
- 專業(yè): 應(yīng)用地球物理學
2樓2014-05-26 16:40:35
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lernen8 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
|
[考研] 考研化學學碩調(diào)劑，一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3 | 123阿標 2026-03-17 | 3/150 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|