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:）為一個行