| 3 | 1/1 | 返回列表 |
| 查看: 393 | 回復(fù): 2 | |||
bakyym木蟲 (小有名氣)
|
[求助]
關(guān)于牛頓法的問題。急急急 已有1人參與
|
| 利用Newton-Raphson求解方程式時(shí),得到的解跳出了定義的解的范圍該怎么辦。怎樣獲得一種安全性比較高的求根方法? |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

木蟲 (小有名氣)
|
SUBROUTINE NEWTON_RAPHSON(I,J,LAMBDA,F) IMPLICIT NONE COMMON/CONF1/RX,RY COMMON/CONF3/SEMIA,SEMIB COMMON/CONF4/THETA COMMON/RAD/NBLA,NBLB INTEGER,PARAMETER:: MAX=20,N=200 INTEGER I,J,M DOUBLE PRECISION,PARAMETER::TOL=1.0D-6 DOUBLE PRECISION RX(N),RY(N) DOUBLE PRECISION NBLA(N),NBLB(N) DOUBLE PRECISION THETA(N) DOUBLE PRECISION SEMIA(N),SEMIB(N) DOUBLE PRECISION LAMBDA,RXIJ,RYIJ DOUBLE PRECISION DX,F,DF2,DF LAMBDA=SEMIA(I)/(SEMIA(I)+SEMIA(J)) DO M=1,100 CALL FUNCD(I,J,LAMBDA,1,F,DF,DF2) DX=DF/DF2 LAMBDA=LAMBDA-DX IF(ABS(DX)<TOL)THEN IF((0.0-LAMBDA)*(LAMBDA-1.0)<0)THEN PAUSE 'LAMBDA JUMP OUT THE BRACKETS ' END IF GOTO 800 END IF END DO 具體的程序是這樣的,函數(shù)是已知的,具體就是求df=0在λ{(lán)0,1}范圍內(nèi)的解。有時(shí)求出的λ值大于1或者小于0.所以不知道該怎么辦 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料,環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +3 | 17501029541 2026-03-23 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|