基于Matlab的ASK仿真

基于MATLAB的2ASK调制

根据数字带通调制原理，运用MATLAB进行编程，仿真调制过程，记录并分析仿真结果。

形成设计报告。

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

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

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

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;

signal=source(n,N);

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

1.随机信号的产生

2.ASK信号调制

3.信道加噪声后信号

4.接收信号的解调

5.解调出的二进制信号

3.4、误码率分析1.ASK误码率分析

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))

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('载波信号'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%调制

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('解调后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决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('抽样判决后波形')