二维傅里叶变换二维峰值检索MATLAB代码

fs=30;N=151;
t=-（N-1）/fs:1/fs:（N-1）/fs;%????????
x=exp（-pi*（t-4）.^2）+exp（-j*pi*t.^2）;  %以chirp信号相加的形式调制的简单高斯信号的混合信号
figure（1）
%plot（tabs（x））;
%title（‘初始信号 ‘）;

u= awgn（x10‘measured‘）; %u是信号x加入信噪比为10db的噪声后的信号，加入前预估信号的功率（强度）
plot（tabs（u））


r=0.01;a=[0:r:2];
G=zeros（length（a）length（u））;         %建立一个201x301的矩阵
h=0;
for n=1:length（a）
    T=frft（ua（n））;                   %T为一个列向量
G（n:）=abs（T（:））;                     %G（n:）为一个行