asech.m

来自「计算动力学系统的分岔图」· M 代码 · 共 18 行

function x = asech(x)
% asech for adtayl objects.


[m,n,p] = size(x.tc);

for i=1:m*n    
    x1.tc= reshape(x.tc(i,:),[1, 1,p]); 
    x1=class(struct('tc',x1.tc),'adtayl');
    asech0 = asech(x1.tc(1,1));
    xx=-x1*(x1+1)*((1-x1)/(1+x1))^0.5;
    dx1 = x1.tc(:,2:end).*(1:p-1);       
    y = filter( 1,xx.tc,dx1)./(1:p-1);
    xx1.tc(i,:)=[asech0 y]; 
   
 end

x.tc = reshape(xx1.tc, m,n,p);