example4_4

代码说明：

   例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 N=N-1 end f=linspace(0,fs/2,N/2) plot(f,abs(y(1:N/2))*2/N) title( 叠加噪声的sinc(t)信号频谱 ) xlabel( f /Hz ) figure y0=y(1:N/2) P=y0.*conj(y0)/N P=10*log10(P) plot(f,P) title( 功率谱 ) xlabel( f /Hz ) )