24小时热门版块排行榜

返回列表

当前主题已经存档。

小木虫:)

荣誉版主 (著名写手)

应助: 0 (幼儿园)
贵宾: 2.542
金币: 9239.2
散金: 130
红花: 22
帖子: 1211
在线: 343.9小时
虫号: 800712
注册: 2009-06-29
性别: GG
专业: 基础物理学
管辖: 物理

[交流] 【求助】请高手帮小虫编一个关于非线性方程组牛顿法的程序（已解决）

各位数值计算的高手，小虫一个朋友让我编一个程序，用牛顿迭代法数值求解一个非线性方程组

一式，(（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 ]

回复此楼

» 猜你喜欢

霸道做事，厚道做人

1楼 2009-10-26 22:26:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

haixing2008

荣誉版主 (文坛精英)

应助: 22 (小学生)
贵宾: 2.252
金币: 32603.8
散金: 2020
红花: 44
沙发: 1
帖子: 25068
在线: 429.1小时
虫号: 540264
注册: 2008-04-06
性别: GG
专业: 计算数学与科学工程计算
管辖: 信息科学

友情建议：“电脑上没有安装MATLAB和Fortran，所以调试成了大问题”
会吓人的，呵呵！还是安装一下吧。

自己可以先查查书，应该不难的。
编程还是得自己摸索一下，实在不行再问哦，不然学的东西肯定不多的，这是我的自身体会，呵呵

赞一下

回复此楼

平平淡淡才是真！

2楼2009-10-27 08:26:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

小木虫:)

荣誉版主 (著名写手)

应助: 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 ]