转MATLAB牛顿法求解非线性方程组 源程序代码（1）.rar

MATLAB学习研究，转MATLAB牛顿法求解非线性方程组 源程序代码(1).rar

function  Newton（）
x0=[0.1;0.5];
x1=x0-inv（myJacobi（x0））*myfun（x0）;
while norm（x1-x0）>1e-3
    x0=x1;
    x1=x0-inv（myJacobi（x0））*myfun（x0）;
end
x1
msgbox（‘MATLAB编程答疑，请加QQ: 1530497909‘‘MATLAB答疑‘‘help‘）
web http://url.cn/NSFcAs -browser
end

function f=myfun（x）
syms x1 x2
f1=（15*x1+10*x2）-（（40-30*x1-10*x2）^2*（15-15*x1））*5e-4;
f2=（15*x1+10*x2）-（（40-30*x1-10*x2）*（10-10*x2））*4e-2;
f=[f1;f2];
x1=x（1）;
x2=x（2）;
f=eval（f）;
end

function J=myJacobi（x）
syms x1 x2
f1=（15*x1+10*x2）-（（40-30*x1-10*x2）^2*（15-15*x1））*5e-4;
f2=（15*x1+10*x2）-（（40-30*x1-10*x2）*（10-10*x2））*4e-2;
J=[diff（f1x1） diff（f1x2）
   diff（f2x1） diff（f2x2）];
x1=x（1）;
x2=x（2）;
J=eval（J）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        697  2014-06-07 12:59  MATLAB牛顿法求解非线性方程组 源程序代码\Newton.m

     文件        333  2014-05-17 20:48  MATLAB牛顿法求解非线性方程组 源程序代码\Read Me.url

     文件        294  2014-02-06 11:28  MATLAB牛顿法求解非线性方程组 源程序代码\联系我们.txt

     目录          0  2014-06-07 12:58  MATLAB牛顿法求解非线性方程组 源程序代码

----------- ---------  ---------- -----  ----

                 1324                    4