| 查看: 449 | 回复: 2 | |||
bakyym木虫 (小有名气)
|
[求助]
关于牛顿法的问题。急急急 已有1人参与
|
| 利用Newton-Raphson求解方程式时，得到的解跳出了定义的解的范围该怎么办。怎样获得一种安全性比较高的求根方法？ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求用牛顿迭代法解非线性方程组的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
- 金币: 68171.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
5