log.m

计算动力学系统的分岔图

function x= log(x)
%

[m,n,p1] = size(x.tc);
% for i=1:p1
%   if x.tc(:,i)==0
%     x.tc(:,i)=0;
%  end
% else
[m,n,p1] = size(x.tc);
p = p1-1;
mn = m*n;
x.tc = reshape(x.tc, mn,p1);

log0 = log(x.tc(:,1));
%log0 = log(x.tc(1,1));
for i=1:mn
  dx = x.tc(i,2:end).*(1:p);
  y = filter( 1, x.tc, dx)./(1:p);
  x.tc(i,:)=[log0 y];  
end 

x.tc = reshape(x.tc, m,n,p1);