基于Matlab的ASK仿真

基于MATLAB的2ASK调制

根据数字带通调制原理，运用MATLAB进行编程，仿真调制过程，

记录并分析仿真结果。

形成设计报告。

课程设计需要运用MATLAB编程实现2ASK,2FSK,2PSK，2DPSK调制解调过程，并且输出其调制后的波形，画出频谱、功率谱密度图，并比较各种调制的误码率情况，讨论其调制效果。

二进制数字调制技术原理

数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。通常使用键控法来实现数字调制，比如对载波的振幅、频率和相位进行键控。

2ASK

2ASK信号的产生方法通常有两种：模拟调制和键控法。解调有相干解调和非相干解调。P=1时f(t)=Acoswt;p=0时f(t)=0;其功率谱密度是基带信号功率谱的线性搬移

专业整理

word格式文档

2.2 数字调制技术的仿真实现

MATLAB是一种功能强大的科学计算和工程仿真软

件，它的交互式集成界面能够帮助用户快速的完成数值分析、数字信号处理、仿真建模、和优化等功能。本课程设计需要运用MATLAB编程实现2ASK,2FSK,2PSK，2DPSK调制解调过程，并且输出其调制后的波形，画出频谱、功率谱密度图，并比较各种调制的误码率情况，讨论其调制效果。

3．程序与调试图形 3.1、ASK调制解调 3.1.1 ASK程序： 3.3.2 ASK图形：

close all clear all n=16;

fc=1000000; %fc>=bitRate fc/bitRate为每个二进制包含sin周期个数 bitRate=1000000; N=50; %noise=ti; noise=10;

专业整理

word格式文档

signal=source(n,N);

transmittedSignal=askModu(signal,bitRate,fc,N); signal1=gussian(transmittedSignal,noise); configueSignal=demoASK(signal1,bitRate,fc,n,N);

1. 随机信号的产生

2. ASK信号调制

专业整理

word格式文档

3. 信道加噪声后信号

4. 接收信号的解调

专业整理

word格式文档

5. 解调出的二进制信号

3.4、误码率分析

1. ASK误码率分析

专业整理

word格式文档

ASK clear all close all i=10; j=5000;

t=linspace(0,5,j);%取0，0.001，0.002，。。。5 共5001个点 fc=10; fm=i/5; B=2*fc;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号 a=round(rand(1,i));%随机序列 figure(2) plot(rand(1,i))

专业整理

word格式文档

st=t; for n=1:10 if a(n)<1;

for m=j/i*(n-1)+1:j/i*n %一个信号取500个点 st(m)=0; end else

for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end figure(1); subplot(421); plot(t,st); axis([0,5,-1,2]); title('基带信号');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波 s1=cos(2*pi*fc*t); subplot(422); plot(s1);

title('载波信号');

专业整理

word格式文档

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制 e_2ask=st.*s1; subplot(423); plot(t,e_2ask); title('已调信号'); noise =rand(1,j);

e_2ask=e_2ask+noise;%加入噪声 subplot(424); plot(t,e_2ask);

title('加入噪声的信号');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调 at=e_2ask.*cos(2*pi*fc*t); at=at-mean(at); subplot(425); plot(t,at);

title('相乘后信号');

[f,af] = T2F(t,at);%通过低通滤波器 [t,at] = lpf(f,af,2*fm); subplot(426); plot(t,at);

title('解调后波形');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决

专业整理

word格式文档

for m=0:i-1;

if at(1,m*500+250)+0.5<0.5; for j=m*500+1:(m+1)*500; at(1,j)=0; end else

for j=m*500+1:(m+1)*500; at(1,j)=1; end end end

subplot(427); plot(t,at); axis([0,5,-1,2]); title('抽样判决后波形')

专业整理