二进制数字调制与解调系统的设计DOCWord格式文档下载.docx
《二进制数字调制与解调系统的设计DOCWord格式文档下载.docx》由会员分享，可在线阅读，更多相关《二进制数字调制与解调系统的设计DOCWord格式文档下载.docx（15页珍藏版）》请在冰豆网上搜索。
非相干解调
2ASK信号非相干解调过程的时间波形
振幅键控是利用载波的幅度变化来传递数字信息，而频率和初始相位保持不变。
在2ASK中：
S2ask=m（t）*cos（2*pi*f*t），
其中m（t）为数字信号，后者为载波。
载波在二进制基带信号控制下通断变化，所以又叫通-断键控（OOK）。
2ASK的产生方法有两种：
模拟调制和键控法而解调也有两中基本方式：
非相干解调（包络检波）和相干解调（同步检测法）DS2ask=s（t）*cos（2*pi*f*t）
=0.5*m（t）+0.5*m（t）*cos（2*wc*t）乘以相干载波后，只要滤去高频部分就可以了本次仿真使用相干解调方式：
2ask信号带通滤波器与与载波相乘低通滤波器抽样判决输出以下就是matlab的仿真结果极其频谱图（省去了带通filter）可以看到解调后的信号与信源有一定的延时。
通过观察频谱图，用放大镜可以清楚的看到，2ask实现了频谱的搬移，将基带信号
搬移到了fc=150hz的频率上，而且若只计频谱的主瓣则有：
B2ask=2fs，fs=1/Ts
其中Ts为一个码元宽度
即：
2ask信号的传输带宽是码元传输速率的2倍
Matlab程序实现
clc;
clearall;
closeall;
%信源
a=randint（1,15,2）;
t=0:
0.001:
0.999;
m=a（ceil（15*t+0.01））;
subplot（511）
plot（t,m）;
axis（[01.2-0.21.2]）;
title（'
信源'
）;
%载波
f=150;
carry=cos（2*pi*f*t）;
%2ASK调制
st=m.*carry;
subplot（512）;
plot（t,st）
axis（[01.2-1.21.2]）
2ASK信号'
）
%加高斯噪声
nst=awgn（st,70）;
%解调部分nst=nst.*carry;
subplot（513）plot（t,nst）axis（[01.2-0.21.2]）;
title（'
乘以相干载波后的信号'
）%低通滤波器设计wp=2*pi*2*f*0.5;
ws=2*pi*2*f*0.9;
Rp=2;
As=45;
[N,wc]=buttord（wp,ws,Rp,As,'
s'
[B,A]=butter（N,wc,'
%低通滤波
h=tf（B,A）;
%转换为传输函数dst=lsim（h,nst,t）;
subplot（514）plot（t,dst）axis（[01.2-0.21.2]）;
经过低通滤波器后的信号'
%判决器
k=0.25;
pdst=1*（dst>
0.25）;
subplot（515）plot（t,pdst）axis（[01.2-0.21.2]）;
经过抽样判决后的信号'
）%频谱观察
%调制信号频谱T=t（end）;
df=1/T;
N=length（st）;
f=（-N/2:
N/2-1）*df;
sf=fftshift（abs（fft（st）））;
figure
（2）subplot（411）plot（f,sf）title（'
调制信号频谱'
）%信源频谱mf=fftshift（abs（fft（m）））;
subplot（412）plot（f,mf）title（'
信源频谱'
%乘以相干载波后的频谱
mmf=fftshift（abs（fft（nst）））;
subplot（413）
plot（f,mmf）
乘以相干载波后的频谱'
%经过低通滤波后的频谱dmf=fftshift（abs（fft（dst）））;
subplot（414）
plot（f,dmf）
经过低通滤波后的频谱'
二、2FSK调制解调
频移键控是利用载波的频率来传递数字信号，在2FSK中，载波的频率随着二进制基带信号在f1和f2两个频率点间变化，频移键控是利用载波的频移变化来传递数字信息的。
在2FSK中，载波的频率随基带信号在f1和f2两个频率点间变化。
故其表达式为：
Acos（1tn）
e2FSK（t）{Acos（t）
Acos（2tn）
典型波形如下图所示。
由图可见。
2FSK信号可以看作两个不同载频的ASK信号的叠加。
因此2FSK信号的时域表达式又可以写成：
s2FSK（t）[ang（tnTs）]cos（1tn）[ang（tnTs）]cos（2tn）
nn
cos（w1t+θn）
t
s1（t）cos（w1t+θn）
s2（t）cos（w2t+φn）
2FSK信号
2FSK数字系统的调制原理
2FSK调制就是使用两个不同的频率的载波信号来传输一个二进制信息序列。
可以用二进制“1”来对应于载频f1，而“0”用来对应于另一相载频w2的已调波形，而这个可以用受矩形脉冲序列控制的开关电路对两个不同的独立的频率源w1、f2进行选择通。
如下原理
图：
2FSK的调制原理图
2FSK的解调方式
2FSK的解调方式有两种:
相干解调方式和非相干解调方式.下面我们将详细的介绍
1非相干解调
经过调制后的2FSK数字信号通过两个频率不同的带通滤波器f1、f2滤出不需要的信号，
然后再将这两种经过滤波的信号分别通过包络检波器检波，最后将两种信号同时输入到抽样判决器同时外加抽样脉冲，最后解调出来的信号就是调制前的输入信号。
其原理图如下图所示：
带通滤波器F1包络检波器
输出
抽样脉冲抽样判决器输出
带通滤波器F2包络检波器
非相干方式原理图
2相干解调
根据已调信号由两个载波f1、f2调制而成，则先用两个分别对f1、f2带通的滤波器对已调信号进行滤波，然后再分别将滤波后的信号与相应的载波f1、f2相乘进行相干解调，再分别低通滤波、用抽样信号进行抽样判决器即可。
原理图如下：
cos2πf1t
cos2πf2t
相干方式原理图
Fc=150;
%载频
Fs=40;
%系统采样频率
Fd=1;
%码速率
N=Fs/Fd;
df=10;
numSymb=25;
%进行仿真的信息代码个数
M=2;
%进制数
SNRpBit=60;
%信噪比
SNR=SNRpBit/log2（M）;
%60
seed=[1234554321];
numPlot=15;
x=randsrc（numSymb,1,[0:
M-1]）;
%产生25个二进制随机码
figure
（1）
stem（[0:
numPlot-1],x（1:
numPlot）,'
bx'
%显示15个码元，杆图，从x的前十五个
随机数中选取
二进制随机序列'
xlabel（'
Time'
ylabel（'
Amplitude'
%调制
y=dmod（x,Fc,Fd,Fs,'
fsk'
M,df）;
%数字带通调制
numModPlot=numPlot*Fs;
%15*40
t=[0:
numModPlot-1]./Fs;
%数组除法（仿真时间）
figure
（2）
plot（t,y（1:
length（t））,'
b-'
axis（[min（t）max（t）-1.51.5]）;
调制后的信号'
）xlabel（'
ylabel（'
在已调信号中
%在已调信号中加入高斯白噪声randn（'
state'
seed
（2））;
%生成-2到+2之间的随机数矩阵y=awgn（y,SNR-10*log10（0.5）-10*log10（N）,'
measured'
[],'
dB'
%
加入高斯白噪声figure（3）plot（t,y（1:
%画出经过信道的实际信号axis（[min（t）max（t）-1.51.5]）;
加入高斯白噪声后的已调信号'
%相干解调figure（4）z1=ddemod（y,Fc,Fd,Fs,'
fsk/eye'
相干解调后的信号的眼图'
）%带输出波形的相干M元频移键控解调figure（5）stem（[0:
holdon;
numPlot-1],z1（1:
ro'
holdoff;
axis（[0numPlot-0.51.5]）;
相干解调后的信号原序列比较'
）legend（'
原输入二进制随机序列'
'
相干解调后的信号'
%非相干解调figure（6）z2=ddemod（y,Fc,Fd,Fs,'
fsk/eye/noncoh'
非相干解调后的信号的眼图'
）%带输出波形的非相干M元频移键控解调figure（7）
holdon;
numPlot-1],z2（1:
非相干解调后的信号'
legend（'
%误码率统计
[errorSymratioSym]=symerr（x,z1）;
figure（8）
simbasebandex（[0:
1:
5]）;
相干解调后误码率统计'
[errorSymratioSym]=symerr（x,z2）;
figure（9）
非相干解调后误码率统计'
%滤除高斯白噪声
Delay=3;
R=0.5;
PropD=0;
%滞后3s
[yf,tf]=rcosine（Fd,Fs,'
fir'
R,Delay）;
%升余弦函数[yo2,to2]=rcosflt（y,Fd,Fs,'
filter'
yf）;
%加入高斯白噪声后的已调信号和经过升余弦滤波器后的已调信号t=[0:
figure（10）
r-'
plot（to2,yo2,'
%滤出带外噪声
axis（[030-1.51.5]）;
经过升余弦滤波器后的已调信号'
）title（'
升余弦滤波前后波形比较'
eyediagram（yo2,N）;
%眼图
加入高斯白噪声后的已调信号的眼图'
仿真结果
三、2psk信号调制解调
2psk信号的调制不能采用包络检测的方法，只能进行相干解调，其原理框图如下：
不考虑噪声时，带通滤波器输出可以表示为y（t）=cos（wct+Φn）
式中Φn为2psk信号某一码元的初相。
Φn=0时，代表数字“0”，Φn=π时，代表数字“1”。
与同步载波COSwct相乘后，输出为
Z（t）=COS（wct+Φn）COSwct=1/2cosΦn+1/2cos（2wct+Φn）经过低通滤波器滤除高频分量，得解调输出为
根据发送端产生2psk信号时Φn代表数字信息1或0的规定，以及接收端x（t）与Φn的关系特性，抽样判决器的判决准则为
其中，x为x（t）在抽样时刻的值。
2psk信号相干解调的过程实际上就是输入已调信号与本地载波信号进行极性比较的过程，故常称为极性比较解调。
clear;
fs=8e5;
%抽样频率
fm=20e3;
%基带频率n=2*（6*fs/fm）;
final=（1/fs）*（n-1）;
fc=2e5;
%载波频率t=0:
1/fs:
（final）;
Fn=fs/2;
%耐奎斯特频率
%用正弦波产生方波twopi_fc_t=2*pi*fm*t;
A=1;
phi=0;
x=A*cos（twopi_fc_t+phi）;
%方波
am=1;
x（x>
0）=am;
x（x<
0）=-1;
figure
（1）subplot（321）;
plot（t,x）;
axis（[02e-4-22]）;
基带信号'
gridoncar=sin（2*pi*fc*t）;
%载波
ask=x.*car;
%载波调制
subplot（322）;
plot（t,ask）;
axis（[0200e-6-22]）;
PSK信号'
gridon;
%=====================================================vn=0.1;
noise=vn*（randn（size（t）））;
%产生噪音
subplot（323）;
plot（t,noise）;
噪音信号'
axis（[0.2e-3-11]）;
askn=（ask+noise）;
%调制后加噪subplot（324）;
plot（t,askn）;
加噪后信号'
%带通滤波
fBW=40e3;
f=[0:
3e3:
4e5];
w=2*pi*f/fs;
z=exp（w*j）;
BW=2*pi*fBW/fs;
a=.8547;
%BW=2（1-a）/sqrt（a）p=（j^2*a^2）;
gain=.135;
Hz=gain*（z+1）.*（z-1）./（z.^2-（p））;
subplot（325）;
plot（f,abs（Hz））;
带通滤波器'
gridon;
Hz（Hz==0）=10^（8）;
%avoidlog（0）
subplot（326）;
plot（f,20*log10（abs（Hz）））;
Receiver-3dBFilterResponse'
axis（[1e53e5-31]）;
%滤波器系数
a=[100.7305];
%[10p]
b=[0.1350-0.135];
%gain*[10-1]faskn=filter（b,a,askn）;
figure
（2）subplot（321）;
plot（t,faskn）;
axis（[0100e-6-22]）;
通过带通滤波后输出'
cm=faskn.*car;
%解调subplot（322）;
plot（t,cm）;
通过相乘器后输出'
%低通滤波器p=0.72;
gain1=0.14;
%gain=（1-p）/2Hz1=gain1*（z+1）./（z-（p））;
Hz1（Hz1==0）=10^（-8）;
%avoidlog（0）plot（f,20*log10（abs（Hz1）））;
LPF-3dBresponse'
axis（[05e4-31]）;
%滤波器系数a1=[1-0.72];
%（z-（p））b1=[0.140.14];
%gain*[11]so=filter（b1,a1,cm）;
so=so*10;
%addgainso=so-mean（so）;
%removesDCcomponentsubplot（324）;
plot（t,so）;
axis（[08e-4-3.53.5]）;
通过低通滤波器后输出'
%比较器High=2.5;
Low=-2.5;
vt=0;
%设立比较标准error=0;
len1=length（so）;
forii=1:
len1ifso（ii）>
=vt
Vs（ii）=High;
else
Vs（ii）=Low;
endendVo=Vs;
subplot（325）;
plot（t,Vo）,title（'
解调后输出信号'
）,
axis（[02e-4-55]）gridon;
时间（s）'
）,ylabel（'
幅度（V）'
升级会员 