Description:
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 )
To Search:
File list (Check if you may need any files):
example4_4.m
Download