大家帮我看看用牛顿法解非线性方程组的matlab错在那里？
大家帮我看看用牛顿法解非线性方程组的matlab错在那里？总共三个函数，调用snow_whzh
function snow_whzh
syms x1 y1 fi1 x2 y2 fi2 x3 y3 fi3 t ti;
q=;%变量列表
Q=[fi1-pi*t-pi/4;
x1-cos(fi1); %fi1是已知变量，初值为pi/4
y1-sin(fi1);
x2-2*cos(fi2)-2*x1;
y2-2*sin(fi2);
y1-y2;
x3-x2-2*cos(fi2);
y3;
fi3;]; %Q是约束方程——非线性方程组
tn=input('过程分析的最大时间量tn=');
J=snow_jacobi(Q,q);
disp(q');
for ti=0:tn/10:tn
Q=subs(Q,t,ti)
q=snow_NR(Q,J,q);
disp(q');
end
function x=snow_NR(Q,J,x)
%用N-R法求解非线性方程组
%其中Q是约束方程向量，J是约束方程雅可比，%X是各个变量组成的向量
m=length(x);
for i=1:m
x0(i)=1-1/i;
end
x0=x0';
N=200; %最大迭代次数
k=1;
while k<N
xx=x0;
QQ=subs(Q,x(1),x0(1));
for i=2:m
QQ=subs(QQ,x(i),x0(i));
end
%disp(QQ);
%QQ=numeric(QQ);
JJ=subs(J,x(1),x0(1));
for i=2:m
JJ=subs(JJ,x(i),x0(i));
end
%JJ=numeric(JJ);
if abs(det(JJ))<eps %JACOBI矩阵奇异时停止计算
break ;
end
x0=x0-JJ\QQ;
x0=numeric(x0);
if norm(x0-xx)<eps%范数
break;
end
k=k+1;
end
x=x0;
function A=snow_jacobi(B,q)
%求约束方程的雅可比矩阵
n=length(B);
m=length(q);
for i=1:n
for j=1:m
A(i,j)=diff(B(i),q(j));
end
end 原帖由 tianya7071 于 2007-11-22 15:32 发表 http://www.chinavib.com/forum/images/common/back.gif
大家帮我看看用牛顿法解非线性方程组的matlab错在那里？
总共三个函数，调用snow_whzh
function snow_whzh
syms x1 y1 fi1 x2 y2 fi2 x3 y3 fi3 t ti;
q=;%变量列表
Q ...
何谓出错？请参照置顶帖的要求叙述清楚你的问题
大家帮我看看用牛顿法解非线性方程组的matlab错在那里
就是运行出现一下的错误，不知道为什么？??? Error using ==> <
Function '<' is not defined for values of class 'sym'.
Error in ==> D:\matlab6.5\work\牛顿迭代\snow_NR.m
On line 25==> if abs(det(JJ))<eps %JACOBI矩阵奇异时停止计算
Error in ==> D:\matlab6.5\work\牛顿迭代\snow_whzh.m
On line 18==> q=snow_NR(Q,J,q); 原帖由 tianya7071 于 2007-11-22 20:51 发表 http://www.chinavib.com/forum/images/common/back.gif
就是运行出现一下的错误，不知道为什么？
??? Error using ==> <
Function ' D:\matlab6.5\work\牛顿迭代\snow_NR.m
On line 25==> if abs(det(JJ)) D:\matlab6.5\work\牛顿迭代\snow_whzh.m
On line ...
符号变量不能进行大小比较，这是所有编程语言都遵循的规则。请按照出错提示自己调试一下、修改一下 怎么调试啊，能不能说的具体点啊 原帖由 tianya7071 于 2007-11-23 10:50 发表 http://www.chinavib.com/forum/images/common/back.gif
怎么调试啊，能不能说的具体点啊
请在置顶帖中找答案 符号变量不能比较
k－N<0 请将原问题公式帖一下，以便确定是否需要象你的代码那样需要符号变量。
页:
[1]