24小时热门版块排行榜    

查看: 813  |  回复: 1
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

早春绿茶

铜虫 (初入文坛)

[求助] 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都是一样的形式,请教一下高手怎么输入方程组以及怎么构建贾克比迭代的系数矩阵呢??
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

readytogo

荣誉版主 (文学泰斗)

瓜帅-催化剂

优秀区长优秀版主优秀版主优秀版主

持续跟进向高人学习
才须学也,非学无以广才,非志无以成学。【科研生活区】、【碳一化工】(煤、天然气化工)欢迎您常来交流
2楼2012-09-18 08:52:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 早春绿茶 的主题更新
信息提示
请填处理意见