| 查看: 583 | 回复: 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 ] |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
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
