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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab 求解非线性方程组 牛顿解法
21/1返回列表
查看: 910  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

早春绿茶

铜虫 (初入文坛)

[求助] matlab 求解非线性方程组 牛顿解法

function Y=f(x,y,z)
Y=[x^2+y^2+z^2-1;
    2*x^2+y^2-4*z;
    3*x^2-4*y+z^2];


function y=JF(x,y,z)
f1='x^2+y^2+z^2-1';
f2='2*x^2+y^2-4*z';
f3='3*x^2-4*y+z^2';
df1x=diff(sym(f1),'x');
df1y=diff(sym(f1),'y');
df1z=diff(sym(f1),'z');
df2x=diff(sym(f2),'x');
df2y=diff(sym(f2),'y');
df2z=diff(sym(f2),'z');
df3x=diff(sym(f3),'x');
df3y=diff(sym(f3),'y');
df3z=diff(sym(f3),'z');
这个是网上求三个非线性方程组的部分程序。。我的模型中想法是采用matlab用牛顿迭代法将非线性方程组线性化，但是如果划分网格为50的话，会有50个非线性方程，中间46都是一样的形式，请教一下高手怎么输入方程组以及怎么构建贾克比迭代的系数矩阵呢？？
回复此楼

» 收录本帖的淘帖专辑推荐

matlab

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-17 16:33:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

readytogo

荣誉版主 (文学泰斗)

瓜帅-催化剂

优秀区长优秀版主优秀版主优秀版主
持续跟进向高人学习
赞一下 回复此楼
才须学也，非学无以广才，非志无以成学。【科研生活区】、【碳一化工】（煤、天然气化工）欢迎您常来交流
2楼2012-09-18 08:52:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 早春绿茶 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 4/200 2026-08-01 21:11 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 4/200 2026-08-01 21:10 by Kddn1pB6pAzX
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 4/200 2026-08-01 21:06 by Kddn1pB6pAzX
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 21:02 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 20:57 by Kddn1pB6pAzX
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 12/600 2026-08-01 20:31 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 20:13 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 11/550 2026-08-01 20:03 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 8/400 2026-08-01 20:03 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 20:02 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 10/500 2026-08-01 19:57 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 19:57 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 19:53 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 19/950 2026-08-01 15:18 by 40862026
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
信息提示
关闭
请填处理意见
关闭