| 查看: 459 | 回复: 2 | |||
bakyym木虫 (小有名气)
|
[求助]
关于牛顿法的问题。急急急 已有1人参与
|
| 利用Newton-Raphson求解方程式时，得到的解跳出了定义的解的范围该怎么办。怎样获得一种安全性比较高的求根方法？ |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有289人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求用牛顿迭代法解非线性方程组的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人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68181.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-06-25 18:07:32
bakyym
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2503.2
- 帖子: 66
- 在线: 27.6小时
- 虫号: 2240593
- 注册: 2013-01-14
- 性别: GG
- 专业: 高分子物理与高分子物理化
|
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.所以不知道该怎么办 |
3楼2014-06-25 18:18:40
50