clear all;
ly1=[];ly2=[];
C=[];
x=0.1;y=0.1;
w=eye(2,2);
r=0.11;
for p=-0.1:0.001:0.15
N=1000; % NUMBER OF ITERATIONS
sl1=0; sl2=0;
for i=1:N
xprev=x;
yprev=y;
x=1/(xprev^2+0.1)-p*xprev+r*xprev*yprev;
y=1/(yprev^2+0.1)-p*yprev+r*xprev*yprev;
jac=[-2*x/(x^2+0.1)-p+r*y r*x
r*y -2*y/(y^2+0.1)-p+r*x]; %%% JACOBIAN
F=w*jac;
[w,r]=qr(F);
sl1=sl1+log(abs(diag(r)));
l1=sl1/N;
end
if N==1000
ly1=[ly1;l1(1)];
ly2=[ly2;l1(2)];
C=[C;p];
%fprintf(1,'l1=%f\n',ly1);
end
end
subplot(2,1,1)
hold on
p=-0.1:0.001:0.15;
line(p,0,'Color','k','LineWidth',4)
plot(C,ly1,'k',C,ly2,'k')
grid
ChaChing哥,有时间帮我看看这个程序问题在哪?就是找不出问题啊,谢谢老兄