基于Matlab的ASK仿真.docx
《基于Matlab的ASK仿真.docx》由会员分享，可在线阅读，更多相关《基于Matlab的ASK仿真.docx（7页珍藏版）》请在冰点文库上搜索。
基于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．1 ASK程序:
３.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;
forn=1：
10
ifａ（n）<1;
ｆorm=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ｏrj=m*5００+１：
（m+1）*５00;
ａｔ（1,j）=0;
ｅnd
else
ｆorj＝m*5０0＋1:
（ｍ+１）*５00;
at（１,j）=1;
eｎd
eｎｄ
end
ｓｕbploｔ（427）;
ploｔ（t,at）;
ａxis（[0,5,-１,2]）;
title（'抽样判决后波形'）
升级会员 