24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bakyym

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2503.2
帖子: 66
在线: 27.6小时
虫号: 2240593
注册: 2013-01-14
性别: GG
专业: 高分子物理与高分子物理化

[求助] 关于牛顿法的问题。急急急已有1人参与

利用Newton-Raphson求解方程式时，得到的解跳出了定义的解的范围该怎么办。怎样获得一种安全性比较高的求根方法？

回复此楼

» 猜你喜欢

推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有20人回复
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有257人回复
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）已经有0人回复
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求用牛顿迭代法解非线性方程组的matlab程序已经有4人回复
求大神帮忙matlab拟合函数求参数已经有9人回复
高手帮忙写个简单程序已经有6人回复
牛顿迭代法解方程。已经有11人回复
修正牛顿法和简化牛顿法是一样的吗已经有3人回复
高斯方程积分编程matlab或C语言已经有15人回复
关于共轭梯度法，类牛顿法中梯度符号判断的疑问已经有4人回复
求用牛顿迭代法求非线性方程组已经有4人回复
Matlab用牛顿法求解非线性方程组问题已经有8人回复
拟牛顿法求解复杂（带根号的）非线性方程组已经有17人回复
方程解的个数，求高手帮忙。急。已经有20人回复
帮忙解一个方程已经有9人回复
小小虫，求助《模式识别与人工智能》，非常感谢已经有11人回复
请问牛顿法和高斯牛顿法有什么不同？非常感谢已经有6人回复
非线性方程组matlab拟合已经有9人回复
重金求Matlab计算一个简单的问题已经有5人回复
求助svm在matlab中的应用已经有7人回复
【求助】帮忙求解四元方程组已经有11人回复
【求助】Levenberg-Marquardt法和LBFGS拟牛顿法，求C++版源码及使用方法已经有3人回复
急急急，求助大家，请帮帮忙求解一个方程组呀！已经有5人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼 2014-06-25 16:52:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

bakyym

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2503.2
帖子: 66
在线: 27.6小时
虫号: 2240593
注册: 2013-01-14
性别: GG
专业: 高分子物理与高分子物理化

引用回帖:

2楼: Originally posted by 月只蓝 at 2014-06-25 18:07:32
不知道你的具体程序是什么。
求解方程可以转换成最优问题，给根加上取值范围，等于是求解有约束的最优化问题。

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
具体的程序是这样的，函数是已知的，具体就是求df=0在λ{0,1}范围内的解。有时求出的λ值大于1或者小于0.所以不知道该怎么办