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

登錄注冊

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

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

返回列表

lernen8

新蟲 (初入文壇)

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

[求助] fortran 求得的非線性解如何帶入其他式子里面計(jì)算

fortran用牛頓法計(jì)算了非線性方程組,但是需要把這個(gè)解用于其他方程計(jì)算驗(yàn)證后再循環(huán)求解。
但是貌似方程的解只能write出，但是不能直接調(diào)用到其他等式啊,計(jì)算出來的結(jié)果不對（永遠(yuǎn)只有一組循環(huán)值）。
（下面只給了主程序，其他的子程序都沒有給）

DIMENSION X(3) , Y(3) , A(3, 3) , B(3) , JS(3)
DOUBLE PRECISION X, Y, A, B, T, H, Tm, Tm0, q
Double precision R=1.0E-4; dm=30.0E-3; km=383.0;
* hw=40000.0; Tw=303.0
DATA X/ 0.1, 0.25, 0.25/
EPS= 1.0E-07
T= 0.1
H= 0.1
Tm0=550.0
5 Tm=Tm0
   CALL FS(X,Y,N,Tm)          （用于給出非線性方程組的子程序）
   CALL DNETN( 3, X, Y, EPS, FS, T, H, A, B, L, JS, Tm) （用于計(jì)算非線性方程組的子程序）
   WRITE( * , * )
   WRITE( * , 30) L
   WRITE( * , * )
   DO 10 I= 1, 3
10 WRITE( * , 20) I, X( I)
20 FORMAT( 1X, 'X( ', I2, ') = ', D15. 6)
30 FORMAT( 1X, 'L= ', I4)    WRITE( * , * )
   q=(1000.0*X(1)-Tm)/R                            （！�。。。【褪沁@個(gè)計(jì)算公式需要調(diào)用計(jì)算的方程組的一個(gè)解X(1)用于計(jì)算，但是貌似式子計(jì)算結(jié)果不對，實(shí)際并沒有能正確的循環(huán)計(jì)算，循環(huán)結(jié)果永遠(yuǎn)只有一組相同的數(shù)據(jù)）
Tmc=Tw+q*(dm/km+1/hw)
   If(abs(Tm-Tmc) .gt. 1.0E-2) then
      Tm0=(Tm+Tmc)/2                            （賦值給新的Tm0, 以用于5的循環(huán)計(jì)算）
      GOTO 5
   endif
   END

希望大蝦能幫忙看看！謝謝！

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工專碩調(diào)劑已經(jīng)有6人回復(fù)
086500 325 求調(diào)劑已經(jīng)有3人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有7人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有6人回復(fù)
環(huán)境工程調(diào)劑已經(jīng)有9人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有9人回復(fù)
287求調(diào)劑已經(jīng)有5人回復(fù)
22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) 已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有3人回復(fù)

1樓 2015-03-17 22:28:15

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

求助誤差函數(shù)的計(jì)算已經(jīng)有8人回復(fù)
求助：超越方程的解法已經(jīng)有13人回復(fù)
求fortran計(jì)算方法指導(dǎo) 14元非線性方程組之前8元的使用bfs法能解14元解不了呀求指導(dǎo) 已經(jīng)有4人回復(fù)
用mathematics解方程組已經(jīng)有19人回復(fù)
大家說說看流體力學(xué)和計(jì)算流體力學(xué)中數(shù)值分析的應(yīng)用已經(jīng)有9人回復(fù)
用fortran程序遺傳算法解非線性方程組已經(jīng)有7人回復(fù)
這個(gè)方程組怎么解啊已經(jīng)有14人回復(fù)
fortran中計(jì)算線性積分已經(jīng)有7人回復(fù)
非線性方程組的解法已經(jīng)有5人回復(fù)
求fortran解非線性方程組已經(jīng)有4人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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