| 查看: 507 | 回复: 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 ] |
回复此楼» 猜你喜欢
11408 321分求调剂
已经有4人回复
-
301求调剂
已经有12人回复
-
085801 总分275 本科新能源 求调剂
已经有9人回复
-
085600 295分求调剂
已经有17人回复
-
289求调剂
已经有23人回复
-
291求调剂
已经有23人回复
-
304求调剂
已经有15人回复
-
285求调剂
已经有11人回复
-
材料调剂
已经有11人回复
-
292分，材料与化工，申请调剂
已经有19人回复
小木虫:)
无虫
- 应助: 0 (幼儿园)
- 贵宾: 2.542
- 金币: 9239.2
- 散金: 130
- 红花: 22
- 帖子: 1211
- 在线: 343.9小时
- 虫号: 800712
- 注册: 2009-06-29
- 性别: GG
- 专业: 基础物理学
- 管辖: 物理
5楼2009-10-28 00:06:37
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
