版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

lernen8

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 57
帖子: 14
在線: 4.3小時(shí)
蟲號(hào): 2137826
注冊(cè): 2012-11-20
專業(yè): 鋼鐵冶金

[求助] 求幫助！改Function中的可變系數(shù)進(jìn)行循環(huán)計(jì)算的問(wèn)題

沒(méi)有金幣了�。。。∠Ｍl(fā)揚(yáng)雷鋒精神��！求大神幫忙解決問(wèn)題��！
我想用Fortran 求解非線性方程組，程序是Monte Carlo法
程序由主程序，F(xiàn)unction, Subroutine三部分組成。
Function給出了非線性方程組，通過(guò)subroutine進(jìn)行計(jì)算非線性方程組，主程序給出了初始計(jì)算值X1,X2,X3。

問(wèn)題是， Function 中有個(gè)系數(shù)M是變化的，
在subroutine（或者主程序）中，先設(shè)定了初始值M0，每次計(jì)算結(jié)束了X后，由X計(jì)算新的M=f(x).
如果新Mc 和舊M差別在1.0E-05內(nèi)就結(jié)束計(jì)算，
否則重新賦予Function 系數(shù)M, 重新調(diào)用子程序循環(huán)計(jì)算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

上面程序是對(duì)subroutine的修改。倒數(shù)第6行直接調(diào)用X1對(duì)新的Mc計(jì)算，并賦值給M0，然后返回行40。
但是貌似根本就沒(méi)有把新的M傳遞給方程F(X,N,M)進(jìn)行循環(huán)計(jì)算。最后只有第一次的結(jié)果。

難道計(jì)算得到的X1不能直接調(diào)用，用于方程計(jì)算么？(嘗試過(guò)把X1直接調(diào)出主程序計(jì)算倒數(shù)第6行的方程，方程結(jié)果根本就不對(duì))
此子程序有否錯(cuò)誤呢？
謝謝！

回復(fù)此樓

» 猜你喜歡

一志愿華中科技大學(xué)071000，求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)

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

【求助】怎么用spss計(jì)算兩個(gè)變量間的相關(guān)系數(shù) 已經(jīng)有5人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
求助�。徃迦颂岢龅年P(guān)于檢測(cè)限的問(wèn)題，不懂是什么意思，求幫助。。。已經(jīng)有14人回復(fù)
變系數(shù)偏微分方程matlab求解已經(jīng)有6人回復(fù)
求回歸方程參數(shù)的最優(yōu)化算法matlab代碼已經(jīng)有17人回復(fù)
用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教！��！已經(jīng)有5人回復(fù)
1stopt 擬合隱函數(shù)的四參數(shù) 已經(jīng)有7人回復(fù)
麻煩大家?guī)兔纯催@個(gè)審稿人什么意思（關(guān)于數(shù)據(jù)分析和”化學(xué)計(jì)量系數(shù)？“）已經(jīng)有5人回復(fù)
哪位蟲友幫幫忙，我的comsol計(jì)算時(shí)出現(xiàn)如下錯(cuò)誤已經(jīng)有9人回復(fù)
兩條不同的氮?dú)馕矫摳角€，求幫助~~ 已經(jīng)有18人回復(fù)
關(guān)于用化學(xué)吸附儀做化學(xué)脈沖實(shí)驗(yàn)計(jì)算金屬分散度的一些問(wèn)題！求幫助��！已經(jīng)有11人回復(fù)
【求助】FLUENT UDF 氣液兩相流汽蝕汽化潛熱(問(wèn)題詳見(jiàn)內(nèi)容) 已經(jīng)有21人回復(fù)
請(qǐng)問(wèn)，參比制劑溶出曲線的變異系數(shù)偏大，是怎么回事��？已經(jīng)有9人回復(fù)
matlab中系數(shù)是矩陣的積分已經(jīng)有5人回復(fù)
如何計(jì)算變位齒輪的齒廓重疊干涉系數(shù) 已經(jīng)有4人回復(fù)
【求助】迎風(fēng)角變化下對(duì)長(zhǎng)方體的風(fēng)阻系數(shù)的關(guān)系，請(qǐng)求幫助！已經(jīng)有7人回復(fù)
【求助/交流】生物統(tǒng)計(jì)學(xué)問(wèn)題，變異系數(shù) 已經(jīng)有11人回復(fù)
【求助】MUSIC的MD模塊計(jì)算甲烷在IRMOF-1上的擴(kuò)散系數(shù) 已經(jīng)有11人回復(fù)
【求助】MATLAB計(jì)算NRTL活度系數(shù) 已經(jīng)有7人回復(fù)
【求助】常微分方程組求解中系數(shù)與某變量值關(guān)聯(lián)的問(wèn)題已經(jīng)有12人回復(fù)
【求助】如何從粒度分布曲線上求算變異系數(shù)CV值？已經(jīng)有4人回復(fù)

1樓 2014-05-23 17:59:18

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

學(xué)蟲dj

金蟲 (正式寫手)

應(yīng)助: 18 (小學(xué)生)
金幣: 892.4
散金: 63
帖子: 688
在線: 89.9小時(shí)
蟲號(hào): 1829966
注冊(cè): 2012-05-22
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

看了半天沒(méi)看懂~樓主最好把代碼打包，問(wèn)題說(shuō)清楚~

贊一下

回復(fù)此樓

窮則獨(dú)善其事，達(dá)則兼濟(jì)天下

2樓2014-05-26 16:40:35

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lernen8 的主題更新