信号分析处理课程设计-基于matlab的模拟信号单边带幅度调制(ssb)与解调分析(精品

信号分析处理课程设计-基于matlab的模拟信号单边带幅度调制(ssb)与解调分析(精品内容摘要：

()( wjwH  的滤波器称为希尔伯特滤波器。 )(wH 传递函数的模和相位特性如图 24 所示。 从图 24 可见，希尔伯特滤波器是 一个宽带 90o 移相网络，是正交变换网络。 )( wHw00 w)( w 图 24 希尔伯特滤波器的传递函数 武汉理工大学《信号分析处理》课程设计说明书 4 4 解调原理 相干解调 相干解调也叫同步检波。 解调与调制的实质一样，均是频谱搬移。 调制是把基带信号的谱搬到了载频位置，这一过程可以通过一个相乘器与载波相乘来实现。 解调则是调制的反过程，即把在载波位置的已调信号的谱搬回到原始基带位置，因此同样可以用相乘器与载波相乘来实现。 相干解调器的一般模型如图 25 所示。 )(tSm )(tSpLPF )(tS d twtccco s)(  图 25 相干解调 器的一般模型 2 级单边带调制解调 在实际应用中，有时需要将信号调制到较高频率的载波上进行传输，但一般设备很难一次性调制成功，所以这时需要将信号分两级调制。 第 1 级调制时，信号先与频率相对较低 1 级载波相乘，再通过带通滤波器滤除下边带频谱得单边带调制信号；第 2 级调制时， 2 级载波频率相对第 1 级较高，再通过带通滤波器后输出高频单边带信号。 解调时， 1 级解调时先与 2 级载波相乘，再通过低通滤波器， 2 级解调时与 1 级载波相乘再通过低通滤波器，最后得解调信号。 3. SSB 调制与解调的 MATLAB 程序实现 函数的使用 正弦 波函数 cos()： 调用格式为 x = a*cos(b*2*pi*t)； abs(x)：纯量的绝对值或向量的长度 ,abs 函数返回一个数的绝对值。 用法为： result = abs(number) 其中 result 是 number 参数的绝对值 ； filter ：一维数字滤波 filter(fb,fa,mo)，这里 fb,fa 分别为滤波器的上下限截止频率，而 mo 为滤波器的输入信号。 其他简单函数： cla 清除当前坐标轴 clc 清除命令窗口显示 clf 清除 当前图形窗口 demo 运行 MATLAB 演示程序 grid 给图形加网格线 gtext 在鼠标指定的位置加文字说明 hold 当前图形保护模式 length 查询向量的维数 linspace 构造线性分布的向量 武汉理工大学《信号分析处理》课程设计说明书 5 5 logspace 构造等对数分布的向量 pi 圆周率 π plot 线性坐标图形绘制 subplot 将图形窗口分成若干个区域 title 给图形加标题 xlabel 给图形加 x 坐标说明 ylabel 给图形加 y 坐标说明 MATLAB 程序实现 这里调制信号为正弦信号， 调制信号幅度 =载波幅度 ，程序如下： Fs=100000。 %抽样频率 t=[0:1/Fs:]。 %抽样间隔 Fc=30000。 %载波频率 a=。 m=a*cos(300*2*pi*t)。 %调制信号 X=fft(m)。 X=abs(X(1:length(X)/2+1))。 %调制信号频谱 frqX=(0:length(X)1)*Fs/length(X)/2 sm = modulate(m,Fc,Fs,39。 amssb39。 )。 %对信号进行调制 Y=fft(sm)。 Y=abs(Y(1:length(Y)/2+1))。 %已调信号频谱 frqY=(0:length(Y)1)*Fs/length(Y)/2 k1=awgn(sm,1)。 %对已调信号加噪声 k2=awgn(sm,5)。 k3=awgn(sm,10)。 k4=awgn(sm,15)。 k5=awgn(sm,20)。 Y1=fft(k1)。 Y1=abs(Y1(1:length(Y1)/2+1))。 %加入噪声后已调信号频谱 frqY1=(0:length(Y1)1)*Fs/length(Y1)/2。 sn=demod(sm,Fc,Fs,39。 amssb39。 )。 %无噪声 ssb信号解调 sn1=demod(k1,Fc,Fs,39。 amssb39。 )。 %加噪声 ssb信号解调 sn2=demod(k2,Fc,Fs,39。 amssb39。 )。 sn3=demod(k3,Fc,Fs,39。 amssb39。 )。 sn4=demod(k4,Fc,Fs,39。 amssb39。 )。 sn5=demod(k5,Fc,Fs,39。 amssb39。 )。 dy1=k1sm。 武汉理工大学《信号分析处理》课程设计说明书 6 6 snr1=var(sm)/var(dy1)。 dy2=sn1sn。 snr2=var(sn)/var(dy2)。 dy12=k2sm。 snr12=var(sm)/var(dy12)。 dy22=sn2sn。 snr22=var(sn)/var(dy22)。 dy13=k3sm。 snr13=var(sm)/var(dy13)。 dy23=sn3sn。 snr23=var(sn)/var(dy23)。 dy14=k4sm。 snr14=var(sm)/var(dy14)。 dy24=sn4sn。 snr24=var(sn)/var(dy24)。 dy15=k5sm。 snr15=var(sm)/var(dy15)。 dy25=sn5sn。 snr25=var(sn)/var(dy25)。 in=[snr1,snr12,snr13,snr14,snr15]。 out=[snr2,snr22,snr23,snr24,snr25]。 Z=fft(sn)。 Z。