respsp.m

来自「一种新的时频分析方法的matlab源程序。」· M 代码 · 共 34 行

function ra=respsp(x1,x2)

% The function RESPSP computes the Hilbert response function 
% from one input, x2,  and one output, x1, both are Hilbert spectra
% of the same dimensions.
% The basic formula is based on R=output/input; therefore,
%
%			Ra=respsp(x1,x2)
%
% The algorithm has put 0.01 at any point having zero for value.
%
% Calling sequence-
% ra=respsp(x1,x2)
%
% Input-
%	x1,x2	- 2-D matrix
% Output-
%	ra	    - 2-D matrix ( ra=x1/x2 )

a1=size(x1);
ix=a1(1);
iy=a1(2);
ra=ones(ix,iy);
ra=ra*0.01;

for i=1:ix
   for j=1:iy
      if(abs(x2(i,j))>.1 & abs(x1(i,j))>.1 )
         ra(i,j)=x1(i,j)/x2(i,j);
      end
   end
end