自己编写的调试通过的radon_wigner变换.txt

自己编写的调试通过的radon变换＋wigner变换求得峰值点坐标

%radon transform
clear all;
%
 N=800;
 n=1:N;
 fs=200;
 t=n/fs;
x1=exp(j*2*pi*(5*t+0.5*5*t.^2));
x2=exp(j*2*pi*(5*t+0.5*15*t.^2));
x=x1+x2;

%****************************************RWT
tfr=tfrwv(x);
%hough transmate
figure(2);
htl(abs(tfr));
[wh,rho,theta]=htl(abs(tfr));
% colormap([0,0,0]);
% mesh(wh);
% xlabel('极半径');
% ylabel('角度');
%找出峰值点的坐标
b=max(max(wh)); 
[u,a]=find(wh=b);