可视化GUI的非线性方程组的牛顿求解法

可视化（GUI） 非线性方程组 牛顿求解法 matlab ecjtu

function x=newton_solve（Fvx0e） 
max1=20;

if （~ischar（F） || ~ischar（v）） %定义类型
    error（‘Parameter F and v should be char type!‘）; 
end 
if （~isnumeric（x0） || ~isnumeric（e）） 
    error（‘Parameter x0 and e should be numeric type!‘）; 
end 
F=sym（F）;v=sym（v）;
dF=jacobian（Fv）;
for index=1:numel（v） 
    %把x分量化，把x1变成x（1），把x2变成x（2）
    F=subs（Fv（index）[‘X（‘ num2str（index） ‘）‘]）;
    dF=subs（dFv（index）[‘X（‘ num2str（index） ‘）‘]）;
end
%dFF

X=x0‘;
P=X;
Y=eval（F‘）
for k=1:max1
    X=P
    J=eval（dF）
    J\Y‘
    Q=P-（J\Y‘）‘;
    X=Q;
    Z=eval（F‘）;
    QP
    err=norm（Q-P）;
    P=Q;
    Y=Z;
    k
    if （err        break
    end
end
X=P;
eval（F‘）
x=P

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         726  2012-06-20 18:11  linfeng1.m
     文件        1287  2012-06-20 17:28  wangqiu.m
     文件         336  2012-06-20 18:17  新建文本文档.txt
     文件      385024  2012-06-19 15:03  课程设计.doc