| 查看: 3044 | 回复: 2 | |||
njutsjh木虫 (小有名气)
|
[求助]
关于用Broyden法解非线性方程组的问题 已有1人参与
|
|
我用Broyden法解非线性方程组的时候，Matlab需要计算很长时间才能出结果，是不是代码有不合理的地方，该怎么修改，谢谢大家！ 代码如下： function =Broyden(func,u,TOL,m) % 应用Broyden方法求非线性方程组法f(x)=0。 % func为待求解的非线性方程组。 % u为初始近似向量。 % TOL为误差容限。 % m为最大迭代次数。 N=size(func); n=N(1); var=findsym(func); A=subs(jacobian(func),var,u); v=subs(func,var,u); H=inv(A); s=-H*v; u=u+s; for k=1:m w=v; v=subs(func,var,u); y=v-w; z=-H*y; p=-s'*z; if(p==0) disp('Method failed!'); return; end C=p*eye(n)+(s+z)*s'; H=C*H/p; s=-H*v; u=u+s; disp(vpa(u',6)); if(max(abs(s))<TOL) break; end end u=vpa(u,6); |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
27博士申请（市政/环境）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助如何用matlab求非线性方程组的解
已经有8人回复
- 求解决matlab解非线性方程组出现的问题
已经有6人回复
- MATLAB解非线性方程组
已经有9人回复
- 请教高手，如何解非线性方程组！！！
已经有13人回复
- matlab 解非线性偏微分方程组
已经有6人回复
- 求用牛顿迭代法解非线性方程组的matlab程序
已经有4人回复
- 一个关于有约束非线性方程组的求解问题
已经有7人回复
- 求教 如何用matlab解带约束的非线性方程组
已经有4人回复
- 求助解复杂非线性方程组的好的方法
已经有24人回复
- Matlab用牛顿法求解非线性方程组问题
已经有8人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 用fortran程序遗传算法解非线性方程组
已经有7人回复
- 牛顿迭代法解含对数的非线性方程组
已经有11人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】非线性方程组 的表示问题
已经有6人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3916
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24222
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2015-01-08 09:43:57
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
【答案】应助回帖
感谢参与，应助指数 +1
| 试试这个，你源程序里面太多符号运算，还有就是你的Jacobi矩阵是调用http://www.pcs.cnu.edu/~bbradie/matlab/systems/broyden.m |
3楼2015-01-08 10:52:33
