版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (5058)
>基金申请 (547)
>导师招生 (386)
>文献求助 (357)
>虫友互识 (322)
>休闲灌水 (199)
>论文投稿 (161)
>招聘信息布告栏 (160)
>博后之家 (144)
>硕博家园 (114)
>考博 (82)
>考研 (70)
>教师之家 (64)
>公派出国 (48)
>环境 (47)
>能源 (45)

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 foolishstar 的 5 个金币

foolishstar

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2359.1
散金: 375
红花: 1
帖子: 266
在线: 92.9小时
虫号: 434695
注册: 2007-08-19
性别: GG
专业: 工程热物理相关交叉领域

[交流] 【求助，5金币】matlab赛德尔迭代求解非线性方程组

请见附件，运行结果出错（计算出来有复数结果）
有高手看看怎么回事吗？

[ Last edited by csfn on 2008-12-29 at 20:10 ]

回复此楼

» 猜你喜欢

国社科又开始会评了，不知道这次命运如何已经有11人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
有多少人是今天查系统知道结果的？已经有19人回复
基金未中，这种答复是模板吗？已经有7人回复
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识已经有7人回复
面上意见出来了已经有23人回复
能否申诉？已经有8人回复

1楼 2008-03-06 18:27:02

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

RE：matlab赛德尔迭代求解非线性方程组

★ ★ ★
lei0736(金币+3,VIP+0):多谢

可以参考一下C代码，（VC++6.0编译通过）

---------------------------------------------------------------------------------------
               b[1]=m_b21;
   b[2]=m_b22;
   b[3]=m_b23;
   a[1][1]=1;
   a[1][2]=0.5;
   a[1][3]=0.5;
   a[2][1]=0.5;
   a[2][2]=1;
   a[2][3]=0.5;
   a[3][1]=0.5;
   a[3][2]=0.5;
   a[3][3]=1;
   X2[0]=m_X20;
   Y2[0]=m_Y20;
   Z2[0]=m_Z20;
   wcx2=m_wcx2;
   X2[1]=(b[1]-a[1][2]*Y2[0]-a[1][3]*Z2[0])/a[1][1];
   Y2[1]=(b[2]-a[2][1]*X2[1]-a[2][3]*Z2[0])/a[2][2];
   Z2[1]=(b[3]-a[3][1]*X2[1]-a[3][2]*Y2[1])/a[3][3];
   i=1;
   while(fabs(X2-X2[i-1])+fabs(Y2-Y2[i-1])+fabs(Z2-Z2[i-1])>wcx2)
   {  i=i+1;
      X2=(b[1]-a[1][2]*Y2[i-1]-a[1][3]*Z2[i-1])/a[1][1];
      Y2=(b[2]-a[2][1]*X2-a[2][3]*Z2[i-1])/a[2][2];
      Z2=(b[3]-a[3][1]*X2-a[3][2]*Y2)/a[3][3];
   }
   m_x_2=X2;
   m_y_2=Y2;
   m_z_2=Z2;
   UpdateData(FALSE);

赞一下(10人)