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信號相加的形式調(diào)制的簡單高斯信號的混合信號
figure（1）
%plot（tabs（x））;
%title（‘初始信號?‘）;

u=?awgn（x10‘measured‘）;?%u是信號x加入信噪比為10db的噪聲后的信號，加入前預(yù)估信號的功率（強(qiáng)度）
plot（tabs（u））


r=0.01;a=[0:r:2];
G=zeros（length（a）length（u））;?????????%建立一個(gè)201x301的矩陣
h=0;
for?n=1:length（a）
????T=frft（ua（n））;???????????????????%T為一個(gè)列向量
G（n:）=abs（T（:））;?????????????????????%G（n:）為一個(gè)行