| 查看: 622 | 回复: 4 | |||
| 当前主题已经存档。 | |||
小木虫:)荣誉版主 (著名写手)
|
[交流]
【求助】请高手帮小虫编一个关于非线性方程组牛顿法的程序（已解决）
|
||
|
各位数值计算的高手，小虫一个朋友让我编一个程序，用牛顿迭代法数值求解一个非线性方程组 一式，(（a1-x2）/a1)^2+(x3/b1)^2-1=0 二式，(（a1-x2-a）/a2)^2+(x3/b2)^2-1=0 三式，x1+x2=h0 其中x1,x2,x3为未知数，其它字母皆为常数。 由于小虫电脑上还没有安装MATLAB和Fortran，所以调试成了大问题，害怕编错了，还是想向各位求助，能否帮小虫编一个计算这个非线性方程组的程序啊？ [ Last edited by formleaf on 2009-10-28 at 06:54 ] |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有26人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
haixing2008
荣誉版主 (文坛精英)
- 应助: 22 (小学生)
- 贵宾: 2.252
- 金币: 32603.8
- 散金: 2020
- 红花: 44
- 沙发: 1
- 帖子: 25068
- 在线: 429.1小时
- 虫号: 540264
- 注册: 2008-04-06
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 信息科学
2楼2009-10-27 08:26:50
小木虫:)
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
|
俺的电脑真的没有装MATLAB，朋友明天就要要了，哪位高手帮小虫调试一下这个程序吧 clear clc X=[0.5;0.5;0.5]; %x1=0.5;x2=0.5;x3=0.5;赋予初值 a1=1;a2=2;a=1;b1=1;b2=3;h0=3; %赋予常数值 err=1; %赋予err初值 while err>0.2 %循环条件要求总误差小于0.2 X0=X; %把X赋给X0 FD=[0,2*(X0(2)-a1)/(a1*a1),2*X0(3)/b1;0,2*(X0(2)+a-a1)/(a2*a2),2*X0(3)/b2;1,1,0]; F=[((a1-X0(2))/a1)^2+(X0(3)/b1)^2-1;((a1-X0(2)-a)/a2)^2+(X0(3)/b2)^2-1;X0(1)+X0(2)-h0]; X=X0-inv(FD)*F; err=(X-X0)'*(X-X0); end X 不胜感激 [ Last edited by 小木虫  on 2009-10-27 at 22:58 ] |
3楼2009-10-27 21:41:52
4楼2009-10-28 00:00:32
小木虫:)
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
5楼2009-10-28 00:06:37
50