：滤波 fs=400 采样频率,必须大于两倍基带信号最高频率 ts=1/fs 采样时间间隔 T=5 时间窗大小 t=0:ts:T -hf1=zeros(1,N) hf1(f<=B1&f>=-B1)=1 y1=hf1.*fftshift(y) figure subplot(211) plot(f,abs(fftshift(y))/N,f,hf1, r-- ) axis([-fs/2,fs/2,0,1.2]) title( 原信号

例4：在sinc(t)信号中叠加噪声，分析对应频谱 fs=100 采样频率,必须大于两倍基带信号最高频率 ts=1/fs 采样时间间隔 T=2 时间窗大小 - x=x+0.1*randn(1,N) x=awgn(x,-10, measured ) figure plot(t,x) title( 时域信号图 ) xlabel( t /s ) y=fft(x,N) figure if mod(N,2)~=0