基于Matlab的ASK仿真
基于Ｍａｔｌab的ASK仿真
————————————————————————————————作者: ————————————————————————————————日期:
?
基于ＭATLAB的2ASＫ调制
根据数字带通调制原理，运用MATLAB进行编程，仿真调制过程，记录并分析仿真结果。
形成设计报告。
课程设计需要运用MAＴLAB编程实现2ＡＳK,2FSK,２ＰSK，2DPSK调制解调过程,并且输出其调制后的波形,画出频谱、功率谱密度图，并比较各种调制的误码率情况,讨论其调制效果。
二进制数字调制技术原理
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。通常使用键控法来实现数字调制，比如对载波的振幅、频率和相位进行键控。
2ASK
２ＡSK信号的产生方法通常有两种:模拟调制和键控法。解调有相干解调和非相干解调。Ｐ=1时f(t）=Ａｃｏswｔ;p=0时ｆ(t)=0;其功率谱密度是基带信号功率谱的线性搬移
2.2数字调制技术的仿真实现?
MATLAＢ是一种功能强大的科学计算和工程仿真软件,它的交互式集成界面能够帮助用户快速的完成数值分析、数字信号处理、仿真建模、和优化等功能。本课程设计需要运用ＭAＴLAB编程实现2ASK,2FＳK,2ＰSK，２DＰＳK 调制解调过程，并且输出其调制后的波形,画出频谱、功率谱密度图,并比较各种调制的误码率情况，讨论其调制效果。
3.程序与调试图形
3.1、ＡSＫ调制解调
３.1．1ASK程序:
３.3.2ＡSK图形：
closｅａll
clｅar all
n=1６;
ｆｃ＝1000000;％ｆc>＝bitRaｔｅfc/bitRate为每个二进制包含ｓin 周期个数
biｔRate=10００000;
Ｎ=５０;
%ｎｏiｓｅ=ti;
ｎoise=10;
siｇnal＝source(n，N)；
tranｓmittedSignaｌ＝asｋMoｄu（sｉgｎal,bitRatｅ,fc,N);
signal1＝gｕsｓian(transmittｅｄＳiｇnal,noｉｓe);
configｕeSｉgnａl=ｄemｏAＳＫ（ｓigｎal1,biｔRａte,fc，n，N);
1.随机信号的产生
?
2.ASK信号调制
?
3.信道加噪声后信号
4.接收信号的解调
?
5.解调出的二进制信号
3．4、误码率分析1.AＳK误码率分析
ＡSＫ
ｃlear ａlｌ
ｃloｓｅall
i=10;
j=5000；
ｔ=linsｐace（0,5，ｊ);%取0，０.0０1，０.00２，。。。5 共50０１个点
fc=１０；
ｆm=i/５；
B=2＊fc; %%％%%%%%%%%%%%%%％％％％%%％%%％%%％%%%%%％％%%产生基带信号
a＝rｏund(rａnd(1,ｉ)）;％随机序列
figure（2)
ｐlot(rａnd（1，ｉ))
sｔ=t;
for n=1：10
if ａ（n)<1;
ｆor m=j/i*（ｎ-1)+1:j/i*n %一个信号取50０个点
st（m）=0;
ｅnd
eｌse
fｏｒm=j/i*(ｎ-１）+1:j/i＊ｎ
ｓｔ(m)=1;
enｄ
enｄ
end
fiｇｕre(1);
subplot(42１）;
pｌｏt(ｔ,st）;
ａxｉｓ（[0,5，-1,２]);
tｉｔｌe('基带信号＇)；％%%％％%%%％%%%%%％%%%%％%%％%%%%%％%％%％%%%%%％%%％%％％载波
s1＝cos（2*ｐi*fc*t)；
ｓubｐlｏt（422)；
plot（s1);
ｔitle（'载波信号'); %％%％％%％%%%％%%%%%%％%%%%%％%％%％%%%%%%%%%% %%%%%%%调制
e_２ask=st.*s1;
subｐlot（423）;
ｐloｔ(t，ｅ_2ask);
ｔitle(＇已调信号＇）;
noiｓe＝ｒaｎd(1，j)；
ｅ＿2ask=ｅ_2aｓk+noise;%加入噪声
suｂplｏt(424);
ploｔ(t,e_２ask);
ｔitle('加入噪声的信号＇); %％%%％%%%％％%%％%%%％%%%%%%%%%%%%%%%%%％相干解调
at=e＿2aｓｋ.*coｓ(2*pi*fc*ｔ);
ａｔ=aｔ-mｅan(ａt）;
ｓubpｌot（425);
pｌoｔ(ｔ,at);
title（＇相乘后信号')；
[f,af]＝T2F(t,aｔ);%通过低通滤波器
[t，at] = ｌpf（f，af,２＊ｆｍ);
subｐｌoｔ(4２６);
pｌot(t，ａt)；
tｉtlｅ（'解调后波形'）; %%％%%%％%%％％%%%%%%%%％%%％%％%％%％％%%%%抽样判决
foｒｍ=０:i-1;
ｉｆat（1,ｍ*５00+250)＋0.5<0.５；
fｏr j=m*5００+１：(m+1)*５00;
ａｔ(1,j)=0;
ｅnd
else
ｆor j＝m*5０0＋1:（ｍ+１）*５00;
at(１,j)=1;
eｎd
eｎｄ
end
ｓｕbploｔ(427);
ploｔ（t,at);
ａxis([0,5,-１,2]);
title（'抽样判决后波形'）