| 查看: 2726 | 回复: 5 | ||
杲1988木虫 (小有名气)
|
[求助]
matlab中利用LM算法解非线性方程组的问题 已有1人参与
|
|
劳烦各位看下：求解如下线性方程组，x(1)-0.7*sin(x(1))-0.2*cos(x(2))=0； x(2)-0.7*cos(x(1))+0.2*sin(x(2))=0 . 是利用改进的高斯牛顿算法，即LM算法求解，编制的三个m文件分别保存，LM算法主程序在lmm这个文件中，三个m文件如下： function [x,val,k]=lmm(Fk,JFk,x0) maxk=100; rho=0.55;sigma=0.4;muk=norm(feval(Fk,x0)); k=0;epsilon=1e-6;n=length(x0); while(k jfk=feval(JFk,x0); gk=jfk'*fk; dk=-(jfk'*jfk+muk*eye(n))\gk; if(norm(gk) while(m<20) newf=0.5*norm(feval(Fk,x0+rho^m*dk))^2; oldf=0.5*norm(feval(Fk,x0))^2; if(newf end m=m+1; end x0=x0+rho^mk*dk; muk=norm(feval(Fk,x0)); k=k+1; end x=x0; val=0.5*muk^2; function JF=JFk(x) JF=[1-0.7*cos(x(1)),0.2*sin(x(2)); 0.7*sin(x(1)),1+0.2*cos(x(2))]; function F = Fk(x) F(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2)); F(2)=x(2)-0.7*cos(x(1))+0.2*sin(x(2)); F=F(  ;命令窗调用格式如下： >> clear all; >> [x,val,k]=lmm(Fk,Jfk,[1,1]) 问题是，一运行，显示： ??? Input argument "x" is undefined. Error in ==> Fk at 3 F(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2)); 想知道是哪里未定义？还是调用格式不对？ 三个文件都没有报错，但调试时，显示 ??? Input argument "Fk" is undefined. Error in ==> lmm at 6 rho=0.55;sigma=0.4;muk=norm(feval(Fk,x0)); 求各位指点！ 对，版本是2011a,32位，XP系统。感激！！！ |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有249人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解非线性方程组解不开
已经有3人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 非线性方程组的解法
已经有5人回复
- Matlab：Newton-Raphson 解线非线性性方程组的问题
已经有12人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- 初学者，请教大家关于matlab解多元二次方程组
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
杲1988
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2174.5
- 散金: 5
- 帖子: 102
- 在线: 88.8小时
- 虫号: 2083787
- 注册: 2012-10-24
- 性别: GG
- 专业: 电机与电器
2楼2013-03-30 10:21:54
杲1988
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2174.5
- 散金: 5
- 帖子: 102
- 在线: 88.8小时
- 虫号: 2083787
- 注册: 2012-10-24
- 性别: GG
- 专业: 电机与电器
3楼2013-03-30 22:13:45
4楼2013-04-01 11:07:14
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2013-04-16 09:05:48
6楼2014-10-09 18:43:22
