计算radon-wigner变换

计算radon-wigner变换的Matlab程序

clear;
clc;

fp=24;%构造函数
a=-10;
dt=0.002;
TT=0.5;
t=0:dt:TT;
s=exp（a*t）.*sin（2*pi*fp*t）;


nsmp=length（t）;
N=512;
for k=nsmp+1:N
    s（k）=0.0;
end

xf=fft（sN）;

for i=1:N
    if i==1
        xf（i）=xf（i）;
    elseif i>1&&i        xf（i）=2*xf（i）;
    else
        xf（i）=0;
    end
end
i=1:N;

% 
% x=ifft（xfN）;


% xr=real（x）;
% xi=imag（x）;
% aa=abs（x）;

% ts=i*dt;

% p=-511:1:512;
% [PI]=meshgrid（pi）;
% ans=P.*（I.^2）;
% XF=meshgrid（xf）;
% G=XF.*exp（0.5j.*ans）;
% g=ifft（GN）;
% D=sqrt（1+P.^2）.*（abs（g））.^2;
% 
PP=500;
p=-PP:1:PP;
np=length（p）;
G=zeros（Nnp）;
Z=xf;
[PI]=meshgrid（pi）;
for w=1:N
     for p=1:PP
         G（wPP+p）=Z（w）.*exp（0.5j.*（p）.*w^2）;
     end
end
 for w=1:N
     for p=1:PP+1
         G（wp）=Z（w）.*exp（0.5j.*（p-PP-1）.*w^2）;
     end
 end
g=ifft（GN）;
M=zeros（12*PP+1）;
D=zeros（12*PP+1）;
for p=1:PP
    M（p+PP）=sqrt（1+p.^2）;
end
for p=1:PP+1
    M（p）=sqrt（1+（p-PP-1）.^2）;
end
for w=1:N
     for p=1:PP
         D（wPP+p）=M（p）.*（abs（g（wPP+p）））.^2;
     end
end
 for w=1:N
     for p=1:PP+1
         D（wp）=M（p）.*（abs（g（wp）））.^2;
     end
 end


% D=M.*（abs（g））.^2;

% angH=angle（x）;
% dxdt_angH=diff（angH）/dt;
% dxdt_angH（512）=0;
% 
% figure
% plot（ixf）;axis（[-100600-3030]）;

figure
plot（iabs（xf））;
axis（[0300050]）;
% 
% title（‘频域分布X（k）‘）;grid on;
% figure
% plot（tsxr‘-‘tsxi‘:‘tsaa‘--‘）;axis（[00.5-11]）;title（‘解析函数及其包络‘）;grid on;
% figure
% plot（tsangH）;axis（[00.5-55]）;title（‘瞬时相位‘）;grid on;
% figure
% plot（tsdxdt_angH）;axis（[00.5-4000800]）;title（‘瞬时频率‘）;grid on;

figure
surf（IPD）;colormap（jet）;shading interp;



% figure
% subplot（311）;plot（tsxr‘-‘tsxi‘:‘tsaa‘--‘）;
% axis（[01.1-11]）;
% subplot（312）;plot（tsangH）;
% % axis（[01.1-55]）;
% subplot（313）;plot（tsdxdt_angH）;
% % axis（[01.1-40001000]）;

% for k=1:nsmp
%     bb（k）=b（k）;
%     yr（k）=real（x（k））;
%     yi（k）=imag（x（k））;
%     tts（k）=k*dt;
%     aaa（k）=abs（x（k））;
%     aangH（k）=angH（k）;
%     adxdt_angH（k）=dxdt_angH（k）;
% end
% figure
% 
% subplot（311）;plot（ttsyr‘-‘ttsyi‘:‘ttsaaa‘--‘）; axis（[00.6-11]）;
% subplot（312）;plot（ttsaangH）;axis（[00.6-55]）;
% subplot（313）;plot（ttsadxdt_angH）;axis（[00.6-4000500]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2355  2013-05-17 20:09  rwtransform.m