mysignal.m

wigner分布的计算

function signal = mysignal(varargin)
%signal
%高斯包络的线性调频信号LFM
n = varargin{1};
t1 =5;
afa = 0.001;
beta = 1;
w0 = 8;
t0 = 0;
if nargin >1
    afa = varargin{2};
elseif nargin >2
    beta = varargin{3};
elseif nargin > 3
    w0 = varargin{4};
elseif nargin > 4
    t0 = varargin{5};
elseif  nargin > 5
    t1 = varargin{5};
end      
t =[-n:2:n]/n*t1;
j = sqrt(-1);
signal = (afa/pi)^(1/4)*exp(-afa*(t-t0).^2/2+j*beta*(t-t0).^2/2+j*w0*(t-t0));
% signal = (afa/pi)^(1/4)*exp(-afa*(t-t0).^2/2+j*beta*(t-t0).^2/2+j*w0*(t-t0)).*(t-t0);