实验1LSB

说明：

1) 原理图上部分为解调部分，下部分为调制部分（即将双边带信号通过一个满足在载波频

率互补对称的低通系统）。sm(t)sLSB(t)

2) 带通滤波器的设计：考虑下边带信号的带宽为2KHz，且在小于20KHz的部分。带通滤

波器的带宽与此相符即可。

3) 低通滤波器设计：考虑到要恢复出的信号带宽2KHz，只要设计出带宽略大于此的低通

滤波器即可。

2.实验结果分析

图1： 下边带（LSB）信号及其频谱

图2：信号在不同过程的时域图

图3：信号在不同过程的频域图

3.实验源程序 %LSB调制解调 clc clear

%基带信号

T=4e-3; fs=1e6; dt=1/fs; t=0:dt:T; f1=2e3;

s=cos(2*pi*f1*t);

%载波信号 fc=20e3;

s1=cos(2*pi*fc*t);

%调制信号

n=hilbert(s);

x=s.*s1+imag(n).*sin(2*pi*fc*t); %高斯白噪声（SNR=5dB） SNR=5;

y=awgn(x,SNR);

%带通滤波 fsamp=1e6;

fcuts=[(20e3-4e3) (20e3-3e3) (20e3+1e3) (20e3+4e3)]; mags=[0 1 0];

devs=[0.05 0.01 0.05];

[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp); hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');

st_p=fftfilt(hh,y); %相干解调

so_step=st_p.*s1;

%**********低通滤波 fsamp=1e6;

fcuts=[3e3 1e4]; mags=[1 0];

devs=[0.01 0.05];

[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp); hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale'); %figure %freqz(hh)

so=fftfilt(hh,so_step);

figure

subplot(2,1,1) plot(t,x)

title('时域LSB波形') subplot(2,1,2) Y=fft(x);

f=(0:(length(Y)-1))*fs/(length(Y))-fs/2; plot(f,fftshift(abs(Y)))

title('LSB信号频域图') figure hold on

subplot(4,1,1) Y=fft(y);

f=(0:(length(Y)-1))*fs/(length(Y))-fs/2; plot(f,fftshift(abs(Y)))

title('加噪后的信号频域图') subplot(4,1,2) Y=fft(st_p);

f=(0:(length(Y)-1))*fs/(length(Y))-fs/2; plot(f,fftshift(abs(Y)))

title('加噪信经BPF后频域图') subplot(4,1,3) Y=fft(so_step);

f=(0:(length(Y)-1))*fs/(length(Y))-fs/2; plot(f,fftshift(abs(Y)))

title('加噪信号与载波相乘后频域图') subplot(4,1,4) Y=fft(so);

f=(0:(length(Y)-1))*fs/(length(Y))-fs/2; plot(f,fftshift(abs(Y))) title('加噪信号经LPF后频域图') figure

subplot(4,1,1) plot(t,y)

title('加噪后的信号时域图') subplot(4,1,2) plot(t,st_p)

title('加噪信经BPF后时域图') subplot(4,1,3) plot(t,so_step)

title('加噪信号与载波相乘后时域图') subplot(4,1,4)

plot(t,so)

title('加噪信号经LPF后时域图')

hold off