| 查看: 2620 | 回复: 8 | |||
z65668396金虫 (正式写手)
|
[求助]
Matlab用牛顿法求解非线性方程组问题
|
|
我的问题是 比如这个方程组 x1.^2-10x1+x2.^2+A=0 x1*x2.^2+x1-10x2+8=0 其中A是在一定范围变化的 例如从-5~5； 用牛顿法解这个方程组，绘制方程组的解随A变化的曲线 如果A是常量 牛顿法很简单 比如A=8 程序如下： 牛顿迭代M文件 function s=NewtonIterate(x,eps) %Newton迭代法求非线性方程组 %x为迭代初值，eps为允许误差值 if nargin==1 eps=1.0e-6; elseif nargin<1 error return end x1=fx1(x);%非线性方程组 x2=-dfx1(x);%非线性方程组导数 x3=inv(x2); x0=x3*x1'; while norm(x0)>=eps %循环迭代 x=x0'+x; x1=fx1(x); x2=-dfx1(x); x3=inv(x2); x0=x3*x1'; end s=x0'+x; return 非线性方程组M文件 function y=fx1(x) y(1)=x(1)*x(1)-10*x(1)+x(2)*x(2)+8; y(2)=x(1)*x(2)*x(2)+x(1)-10*x(2)+8; y=[y(1) y(2)]; 偏导M文件 function y=dfx1(x) y(1)=2*x(1)-10; y(2)=2*x(2); y(3)=x(2)*x(2)+1; y(4)=2*x(1)*x(2)-10; y=[y(1) y(2);y(3) y(4)]; 如果A是变量，就要改动非线性方程组的M函数， function y=fx1(x) A=-5：5; y(1)=x(1)*x(1)-10*x(1)+x(2)*x(2)+A; y(2)=x(1)*x(2)*x(2)+x(1)-10*x(2)+8; y=[y(1) y(2)]; 运行报错 不知道哪位大神能够解答一下 非常感谢 |
回复此楼» 猜你喜欢
291求调剂
已经有7人回复
-
300求调剂
已经有8人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这个方程组怎么解啊
已经有14人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- Matlab：Newton-Raphson 解线非线性性方程组的问题
已经有12人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- matlab调试问题
已经有3人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 关于matlab线性规划 背包问题的 解法 求救呀~~方程知道 但不知道如何解答····
已经有4人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- 纯新手，全方面求助关于MATLAB........
已经有9人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 请问牛顿法和高斯牛顿法有什么不同？非常感谢
已经有6人回复
- 怎么才能做出三个9的线性方程啊？？！！愁死我了
已经有11人回复
- matlab求解非线性方程组
已经有16人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求一个mathmetica的程序
已经有10人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】求解非线性方程
已经有8人回复
- 【讨论】关于线性试验中线性方程的标准
已经有8人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
2楼2013-09-05 14:04:15
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
3楼2013-09-05 14:18:05
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
4楼2013-09-05 14:28:32
z65668396
金虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1120.8
- 散金: 1682
- 红花: 3
- 帖子: 469
- 在线: 369.6小时
- 虫号: 1774509
- 注册: 2012-04-23
- 性别: GG
- 专业: 机械动力学
5楼2013-09-05 15:44:21
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
6楼2013-09-05 16:06:05
z65668396
金虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1120.8
- 散金: 1682
- 红花: 3
- 帖子: 469
- 在线: 369.6小时
- 虫号: 1774509
- 注册: 2012-04-23
- 性别: GG
- 专业: 机械动力学
7楼2013-09-05 16:21:48
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2013-10-26 23:26:29
z65668396
金虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1120.8
- 散金: 1682
- 红花: 3
- 帖子: 469
- 在线: 369.6小时
- 虫号: 1774509
- 注册: 2012-04-23
- 性别: GG
- 专业: 机械动力学
9楼2013-10-27 09:52:25
5