tanh.m

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

function x = tanh(x)
% Tanhx for adtayl objects.

if ~isreal(x.tc)
  error('X for adtayl objects must not have imaginary part')
end
 
[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');
    y= (exp(x1)-exp(-x1))/(exp(x1)+exp(-x1));
    
    xx1.tc(i,:)=y.tc;
   
 end
x.tc = reshape(xx1.tc, m,n,p);
x.tc =real(x.tc);