clc
close?all
clear?all
%%?-初始化參數
fs=2;???????????????????????????????????????????????????%采樣頻率
ns=1;
a=2.3;??????????????????????
tr=0.4*ns;?????????????????????
tf=2*ns;?
m=160;????????????????????????????????????????????????????%序列延遲長度
nline=4096;?????????????????????????????????????????????%序列點長（需要為偶數）
noise_n=14;?????????????????????????????????????????????%單頻噪聲數量
fir_n=100;??????????????????????????????????????????????%fir濾波器點數
dt=1/fs;
t=-nline/2*dt-m*dt:dt:nline/2*dt-dt;????????????????????%離散時間點
t_temp=-nline/2*dt-dt-m*dt:dt:nline/2*dt-dt;????????????%離散時間點
nt=length（t）;
%%?-sn?vn?yn信號序列產生
%-------------sn---------------
s=1./（exp（-a*t_temp/tr）+exp（a*t_temp/tf））;
sn=5*（s（2:end）-s（1:end-1））;?????????????????????%有用信號
%-------------yn---------------
freq=[0.611.82.134.85.25.76.16.46.777.89.3]‘*fs/20;
ampli=repmat（[0.5110.50.10.30.51110.50.31.50.5]‘1nt）;
%?rand（‘seed‘1954）;
phase=repmat（2*pi*rand（noise_n1）1nt）;
yn=sum（ampli.*cos（2*pi*freq*t+phase））;??????????%單頻噪聲信號
%-------------vn---------------
vn=wgn（1n