幅度调制信号的的解调性能仿真与分析文档Word下载.docx
- 文档编号:21610177
- 上传时间:2023-01-31
- 格式:DOCX
- 页数:17
- 大小:184.85KB
幅度调制信号的的解调性能仿真与分析文档Word下载.docx
《幅度调制信号的的解调性能仿真与分析文档Word下载.docx》由会员分享,可在线阅读,更多相关《幅度调制信号的的解调性能仿真与分析文档Word下载.docx(17页珍藏版)》请在冰豆网上搜索。
去掉,即可输出抑制载波双边带信号,简称双边带信号(DSB)。
DSB调制器模型如图1所示。
图1DSB调制器模型
其中,设正弦载波为
式中,
为载波幅度;
为载波角频率;
为初始相位(假定
为0)。
调制过程是一个频谱搬移的过程,它是将低频信号的频谱搬移到载频位置。
而解调是将位于载频的信号频谱再搬回来,并且不失真地恢复出原始基带信号。
双边带解调通常采用相干解调的方式,它使用一个同步解调器,即由相乘器和低通滤波器组成。
在解调过程中,输入信号和噪声可以分别单独解调。
相干解调的原理框图如图2所示:
图2相干解调器的数学模型
信号传输信道为高斯白噪声信道,其功率为
DSB信号调制过程分析
假定调制信号
的平均值为0,与载波相乘,即可形成DSB信号,其时域表达式为
的平均值为0。
DSB的频谱为
DSB信号的包络不再与调制信号的变化规律一致,因而不能采用简单的包络检波来恢复调制信号,需采用相干解调(同步检波)。
另外,在调制信号
的过零点处,高频载波相位有180°
的突变。
除了不再含有载频分量离散谱外,DSB信号的频谱与AM信号的频谱完全相同,仍由上下对称的两个边带组成。
所以DSB信号的带宽与AM信号的带宽相同,也为基带信号带宽的两倍,即
为调制信号的最高频率。
(3)SSB信号仿真分析设计概要
信号的产生
由题意可知,未调信号的频率f=1Hz,功率P=1W,载波频率10Hz。
设采用时间为0.001S,频率分辨率为0.1。
由于正弦信号的功率与幅值有以下关系:
,可以求出未调信号幅值。
所以未调信号表达为:
m=Am*cos(2*pi*ft*t)。
信号的调制
由于SSB是通过滤波法实现。
通过公式
实现DSB信号,并通过傅立叶变换得其频谱,然后去除上边频分量得到下边频分量LSSB,再通过傅立叶反变换即可产生携带下边频的单边带调幅信号即u信号。
在MATLAB中fftseq函数可以实现傅立叶变换,iff函数可以实现傅立叶反变换。
信号的解调
单边带信号的时域表达式为:
,将已调信号u与同频同相的载波
相乘后可以得到含源信号的表达试
,通过截至频率合理的低通滤波器就能将源信号恢复出来。
在MATLAB中,低通滤波器可以floor函数实现。
3运行环境
硬件:
计算机
软件:
MATLAB R2010b
4开发工具和编程语言
开发工具:
计算机,MATLAB R2010b
编程语言:
matlab
5详细设计
AM信号及调制解调原程序
clc
clearall
closeall
Fs=100000;
%samplefrequency采样频率
Fc=1000;
%carrierfrequency载波频率
Ts=0.1;
%sampletime抽样时间
Ns=Fs*Ts;
%lengthofsignal信号总长度
t=0:
1/Fs:
(Ns-1)/Fs;
A0=2;
forsnr=-5:
5
g1=2*cos(2*pi*50*t);
%g1=3*cos(50*pi*t).*cos(50*pi*t);
%g2=0;
g2=cos(2*pi*20*t);
%g2=3*sin(pi*100*(t-0.05))./(pi*100*(t-0.05));
g=g1+g2;
%原始调制信号
figure
(1)
subplot(2,1,1)
plot(g)
x=ammod(g,Fc,Fs,0,A0);
%已调AM信号
subplot(2,1,2)
plot(x)
x1=hilbert(x);
x1=awgn(x,snr);
%按照信噪比加入高斯白噪声
y=x1(:
);
t=(0:
1/Fs:
(size(y,1)-1)/Fs)'
;
t=t(:
ones(1,size(y,2)));
z=y.*cos(2*pi*Fc*t);
[num,den]=butter(5,Fc*2/Fs);
%设计低通滤波器
fori=1:
size(y,2)
z(:
i)=filtfilt(num,den,z(:
i))*2;
end
xx=z-A0;
figure
(2)
holdon
plot(xx,'
g'
)
holdoff
end
DSB信号及调制解调原程序
%DSB信号调制程序
clc;
%清除窗口中的图形
ts=0.01;
%定义变量区间步长
t0=2;
%定义变量区间终止值
t=-t0+0.0001:
ts:
t0;
%定义变量区间
fc=10;
%给出相干载波的频率
A=2;
%定义输入信号幅度
fa=1;
%定义调制信号频率
mt=A*cos(2*pi*fa.*t);
%输入调制信号表达式
ct=cos(2*pi*fc.*t);
psnt=mt.*cos(2*pi*fc.*t);
%输出调制信号表达式
subplot(3,1,1);
%划分画图区间
plot(t,mt,'
%画出输入信号波形
title('
输入信号波形'
xlabel('
Variablet'
ylabel('
Variablemt'
subplot(3,1,2);
plot(t,ct,'
b'
输入载波波形'
Variablect'
subplot(3,1,3);
plot(1:
length(psnt),psnt,'
r'
%length用于长度匹配
已调信号波形'
%画出已调信号波形
Variablepsnt'
%DSB信号加入噪声后解调程序
clf;
xzb=2;
%输入小信躁比(dB)
snr=10.^(xzb/10);
[h,l]=size(mt);
%求调制信号的维数
fangcha=A*A./(2*snr);
%由信躁比求方差
nit=sqrt(fangcha).*randn(h,l);
%产生小信噪比高斯白躁声
psmt=mt.*cos(2*pi*fc.*t);
psnt=psmt+nit;
%输出叠加小信噪比已调信号波形
xzb=20;
%输入大信躁比(dB)
snr1=10.^(xzb/10);
fangcha1=A*A./(2*snr1);
%由信躁比求方差
nit1=sqrt(fangcha1).*randn(h,l);
%产生大信噪比高斯白躁声
psnt1=psmt+nit1;
%输出已调信号波形
subplot(2,2,1);
plot(t,nit,'
小信噪比高斯白躁声'
Variablenit'
subplot(2,2,2);
plot(t,psnt,'
叠加小信噪比已调信号波形'
subplot(2,2,3);
plot(t,nit1,'
%length用于长度匹配
大信噪比高斯白躁声'
%画出输入信号与噪声叠加波形
subplot(2,2,4);
plot(t,psnt1,'
k'
叠加大信噪比已调信号波形'
%画出输出信号波形
Variablepsmt'
SSB信号调制解调原程序
t0=1;
ts=0.001;
fs=1/ts;
df=0.3;
t=[-t0/2:
t0/2];
m=sqrt
(2)*cos(2*pi*t);
c=cos(2*pi*fc.*t);
%定义载波同相分量
b=sin(2*pi*fc.*t);
v=m.*c+imag(hilbert(m)).*b;
u=m.*c-imag(hilbert(m)).*b;
[M,m,dfl]=fftseq(m,ts,df);
M=M/fs;
[U,u,dfl]=fftseq(u,ts,df);
U=U/fs;
[V,v,dfl]=fftseq(v,ts,df);
V=V/fs;
f=[0:
dfl:
dfl*(length(m)-1)]-fs/2;
holdon;
figure
(1)
subplot(1,2,1)
plot(t,m(1:
length(t)));
axis([-1,1,-2,2]);
时间'
未调信号'
subplot(1,2,2);
plot(f,abs(fftshift(M)))
频率'
未调信号的频谱'
figure
(2);
plot(t,c(1:
axis([-0.1,0.1,-2,2])
载波'
figure(3);
subplot(2,2,1)
plot(t,u(1:
axis([-0.2,0.2,-1.5,1.5]);
上边带已调信号'
)
subplot(2,2,2)
plot(t,v(1:
length(t)))
下边带已调信号'
subplot(2,2,3)
plot(f,abs(fftshift(U)))
上边带已调信号的频谱'
plot(f,abs(fftshift(V)))
下边带已调信号的频谱'
6调试分析
我在该课程设计中曾遇到过一些问题,具体如下:
其一是频谱显示问题。
起初我得出的频谱显示波形并非是频率为0点居中,而是0点在最左端。
这样不利于我观察频谱特性,为了能让频谱图居中显示,我对程序设置进行分析,发现关键在于频率矢量的设置上。
起初所设置的值f=(0:
length(U)-1)*fs/length(U),不能将图形显示在0Hz两边,经过改进为f=(0:
length(U)-1)*fs/length(U)-fs/2,使得频谱居中显示成功。
其中-fs/2就是将图形向右搬移半个长度,这样0点就能居中了。
其二是滤波器截止频率fl的设置问题。
开始我将截止频率设置为fl=50Hz,但仿真结果显示解调信号的幅度明显小于基带信号。
这说明基带信号没有完全被还原,有部分损失。
经分析我认为是滤波器的截止频率太低,将基带信号的成分也滤去了一些。
然而滤波器的截止频率也不宜过大,这样会使它无法滤除混频信号中的2
高频分量。
对于其合适的参数我进行了多次测试,力求在可以滤去2
高频分量的情况下使解调信号损失最小。
最终我发现将fl设为700Hz较为合理,解决了此问题。
7测试结果
1.AM信号及解调信号图
为了验证算法的有效性,本文进行了大量的仿真实验。
信号载波频率为1kHz,采样频率为100kHz;
调制信号为双音信号,表达式为:
,调制仿真了3种情况下AM信号,即满调幅情况下的,欠调幅情况下的以及过调幅情况下的已调的AM信号。
同时在满调幅情况下的AM信号进行了相干解调仿真,信号持续时间为0.1秒,即每次采集10000点进行处理。
加入噪声为高斯白噪声,信噪比从-5dB到5dB,步进为1dB。
图4.1,图4.2和图4.3分别给出了调制信号波形以及在满调幅情况下的,欠调幅情况下的以及过调幅情况下的已调的AM信号。
图4.4给出了信噪比为-5dB时的过调幅情况下的AM信号的相干解调结果。
图4.5,图4.6和图4.7分别给出了满调幅情况下的AM信号在信噪比为-5dB,0dB和5dB情况下的相干解调结果。
图4.1满调幅情况下的调制信号及已调AM信号
4.2调制信号与欠调幅情况下的已调AM信号
图4.3调制信号与过调幅情况下的已调AM信号
图4.4信噪比为-5dB时过调幅AM信号的相干解调结果
图4.5信噪比为-5dB时满调幅AM信号相干解调结果
图4.6信噪比为0dB时满调幅AM信号的相干解调结果
图4.7信噪比为5dB时满调幅AM信号相干解调结果
从图4.1到图4.3可以看到在满调幅与欠调幅情况下信号的包络没有发生失真,而在过调幅情况下信号的包络发生了失真,因此,对于过调幅AM信号来说,它并不适合用包络检波的方法进行解调。
但是利用相干解调的方法可以将过调幅AM信号正确解调出来。
图4.4给出了信噪比为-5dB时过调幅AM信号的相干解调结果。
从结果可以看到,利用相干解调方法对过调幅AM信号进行解调,除了噪声引起的失真以外,可以正确解调。
从图4.5到图4.7可以看到AM信号的解调结果随着信噪比的变化而变化。
信噪比越大,解调信号越接近于原始的调制信号,也就是说,噪声越大对信号的解调结果的影响越大,噪声越小,对信号的解调结果影响越小。
2.DSB信号及解调信号图
可以清晰地看出,加大噪声后,解调信号的波形杂乱无章,起伏远大于加小噪声时的波形。
造成此现象的原因是当信噪比较小时,噪声的功率在解调信号中所占比重较大,所以会造成杂波较多的情况;
而信噪比很大时,噪声的功率在解调信号中所占比重就很小了,噪声部分造成的杂乱波形相对就不是很明显,甚至可以忽略。
综上所述,叠加噪声会造成解调信号的失真,信噪比越小,失真程度越大。
所以当信噪比低于一定大小时,会给解调信号带来严重的失真,导致接收端无法正确地接收有用信号。
所以在解调的实际应用中,应该尽量减少噪声的产生。
3.SSB信号及解调信号图
图1:
未调信号及其频谱
图2:
载波信号
图3:
已调信号及其频谱
仿真分析
从图1中可以看出,未调信号频率为1Hz。
从频谱角度上看,在0处出现冲击,这是由正弦信号的频率特性确定的。
由抑制载波双边带调幅调制出来的信号是以正弦信号为包络的不等值正弦信号组成的波形,经过滤除上边带后得到图3的携带下边带信号的已调信号,及经过滤除下边带后得到图3的吓呆上边带信号的已调信号。
调制过程就是将信号的频谱进行搬移,将其搬运到载波附近的频率点上。
下图是分别对应上、下边带信号的频谱图。
与图1中的频谱图进行比较能清晰地看出,频域信号已经被搬移到载波的附件。
图2为载波信号的波形。
从图4功率普密度图中看到,已调信号的功率主要集中在10Hz附件,主要是因为信号频率已经被搬运到载波附件的原因。
经过滤波器后,信号的大致形状已经被恢复,但由于调制和解调的过程中信号产生了相移,所以与未调信号相比,相位发送了较大的变化。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 幅度 调制 信号 解调 性能 仿真 分析 文档