24小时热门版块排行榜

返回列表

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

bakyym

木虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有136人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有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.所以不知道该怎么办