%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%設置參數%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all
close?all
clc
pf=0.1;????????????????????????????????%虛警概率為0.1
fs=1000;???????????????????????????????%采樣頻率為1000
N=1024;????????????????????????????????%數據點數為1024

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%求信號的頻譜%%%%%%%%%%%%%%%%%%%%%%%%%%
haoduan=zeros（110）;???????????????????
haoduan（2）=1;??????????????????????????%原始信號的第2段有信號
haoduan（4）=1;??????????????????????????%原始信號的第4段有信號
haoduan（7）=1;??????????????????????????%原始信號的第7段有信號
haoduan（9）=1;??????????????????????????%原始信號的第9段有信號

x=0;
m=0;
fc=180:180:500;????????????????????????%分配載頻
for?i=1:2
????n=0:N-1;
????t=n/fs;????????????????????????????%時間序列
????m=2*sin（2*pi*fc（i）*t）;?????????????????
????x=x+m;?????????????????????????????%無噪聲信號???
end
y=x+2*randn（1length（t））;??????????????%含有噪聲的信號
Y=fft（yN）;????????????????????????????%對含有噪聲的信號進行快速Fourier變換
mag=abs（Y）;????????????????????????????%得到振幅值
f=n*fs/N;??????????????????????????????%頻率序列
figure（1）