一、实验目的
1、掌握用matlab 中高斯白噪声信道的产生方法。 4、掌握白噪声的消除方法。
二、实验步骤
1、产生并分析理想低通高斯白噪声的时频特性,绘制其时频波形。
2、产生并分析理想带通高斯白噪声的时频特性,绘制其时频波形。
3、任意输入一个信号如()cos(2*15)s t t π=,绘制其时频域波形。
4、产生一个正态分布的高斯白噪声()n t ,绘制()()()r t s t n t =+的时频域波形。
5、将()()()r t s t n t =+通过带通滤波器,绘制通过带通滤波器后的时频域波形。
老婆饼机6、计算高斯白噪声信道的信噪比为-10dB 到10dB ,每2dB 递进,发送端信号为0、1时的误码率情况。蓝牙门禁
三、实验内容
(1) 理想低通高斯白噪声的产生和分析
程序:
打印头校准dt=0.01; %时间间隔
纺织机轴承
N=1024; %傅立叶变换点数
et=N*dt-dt; %结束时间
t=bt:dt:et; %时间域
TT=et-bt;
df=1/TT; %频率间隔
nt=-TT/2:dt:TT/2;
Tf=N*df;
f=-Tf/2+df: df :Tf/2;
fh=10;
ts=2*fh*sinc(2*fh*nt);
figure(1);
subplot(2,1,1); plot(nt,ts);
title('低通理想白噪声时域');
xlabel('时间(s)');
axis([-1 1 min(ts) max(ts)]);
grid on;
tf=fftshift(fft(ts));
subplot(2,1,2);
plot(f,abs(tf));
title('低通理想白噪声频域');
xlabel('频率(Hz)');
axis([-30 30 min(abs(tf)) max(abs(tf))]);
grid on;
(2)理想带通高斯白噪声的产生和分析程序:
fs=20;
B=10;
ts=2*B*sinc(B*nt).*cos(2*pi*fs*nt);
figure(2);
subplot(2,1,1);
plot(nt,ts);
title('带通理想白噪声时域');
xlabel('时间(s)');
axis([-1 1 min(ts) max(ts)]);
grid on;
tf=fftshift(fft(ts));
subplot(2,1,2);
plot(f,abs(tf));
title('带通理想白噪声频域');
xlabel('频率(Hz)');
axis([-30 30 min(abs(tf)) max(abs(tf))]); grid on;
(3) 信号()cos(2*15)s t t π=,绘制其时频域波形
程序:
bt=0; %开始时间
dt=0.01; %时间间隔
N=4096; %傅立叶变换点数
et=bt+N*dt-dt; %结束时间
t=bt:dt:et; %时间域
TT=et-bt; %总的时间
nt=-TT/2:dt:TT/2;
y = sin(2*pi*10*t); %待分析的信号
df=1/TT; %频率间隔
Tf=N*df %分析的频宽 f=-Tf/2+df:df:Tf/2; %频率域
figure(1);
subplot(2,1,1);
plot(t,y);
title('信号');
xlabel('时间(s)');
axis([0 2 -1.1 1.1]);
grid on;
Y =fftshift(fft(y,N));
subplot(2,1,2);
plot(f,abs(Y));
title('频谱');
xlabel('频率(Hz)');
axis([-30 30 min(abs(Y)) max(abs(Y))]);
grid on;ktv点歌台
(4)()()()
=+的时频域波形
r t s t n t
程序:
bt=0; %开始时间
dt=0.01; %时间间隔
N=4096; %傅立叶变换点数
et=bt+N*dt-dt; %结束时间
t=bt:dt:et; %时间域
隧道喷浆
TT=et-bt; %总的时间
nt=-TT/2:dt:TT/2;
y=sin(2*pi*10*t)+0.5*randn(1,length(t)); %待分析的信号
df=1/TT; %频率间隔
Tf=N*df %分析的频宽
f=-Tf/2+df:df:Tf/2; %频率域
subplot(2,1,1);
plot(t,y);
title('加噪信号');
xlabel('时间(s)');
axis([0 2 min(y) max(y)]);
grid on;
Y =fft(y,N);
subplot(2,1,2);
plot(f,abs(fftshift(Y)));
title('加噪频谱');
xlabel('频率(Hz)');
axis([-30 30 min(abs(Y)) max(abs(Y))]);
grid on;
(5)()()()
=+通过带通滤波器后的时频域波形
r t s t n t
程序:
fs=10;
B=2;
pf=zeros(1,length(f));
pf(length(f)/2+round((fs-B)/df):length(f)/2+round((fs+B)/df))=1;
pf(length(f)/2-round((fs+B)/df):length(f)/2-round((fs-B)/df))=1;
Y=Y.*fftshift(pf);
y=ifft(Y);
subplot(2,1,1);