clear?all
N=256;???????????????????????????????%觀測樣本個數??????????????????
noise=sqrt（0.01）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
?????????????????????????????????????%噪聲信號
noise2=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
noise3=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
noise4=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
f1=0.1;??????????????????????????????%信號歸一化頻率
f2=0.25;
f3=0.27;
A1=0.1;????????????????????????????????%信號幅度???
A2=0.1;
A3=0.1;
snr1=-15;
snr?=power（10snr1/10）;%db轉化十進制
pn?=?（std（noise））.^2;
A=sqrt（pn*snr）;
signal1=A1*exp（1i*2*pi*f1*（0:N-1））;??%有用信號
signal2=A1*exp（1i*2*pi*f2*（0:N-1））;
signal3=A1*exp（1i*2*pi*f3*（0:N-1））;
for?j=1:30
????noise1=（0.5*i）*sqrt（20）*（（randn（1N）+1i*randn（1N））/sqrt（2