通信原理课程设计---2FSK数字调制系统仿真和分析

《通信原理课程设计---2FSK数字调制系统仿真和分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、X丄-*•垃n倒丿久&孑课程设计课程设计名称：—专业班级：学生姓名：学号：一指导教师：课程设计时间：1需求分析二进制频移键控(2FSK)数字调制系统:1、主要功能：对信号编码形成的0、1序列通过两种不同频率的波表达出来。经过信道有噪声的加入，在这里我们模拟为高斯口噪声，在接收端收到信号后通过分离、滤波、判决，进而还原信号。2、主要内容：对二进制数字信源进行数字调制(2FSK),画出信号波形及功率谱，并分析其性能。2概要设计整个设计包描三个部分：1、信源部分：在这里要求2FSK信号的形成，根据公式s

2、=cos(2*pi*(F+m*f).*t)(m=0或1)生成信源信号,设置F二20HZ,m=a(ceil(t+0.0005)),f=100HZ,t二0：0.0005：9.9995o生成的s为一维数列，有两万组数据。2、信道部分：在信道部分有噪声的加入，这里用如下函数n=0.01*randn(l,20000);z二s+n;此处n为噪声函数，randn函数用于生成随机分布的一维数列，z为加入了噪声的合成信号。3、信宿部分：在这一部分，要求接收到的信号z经过滤波函数filter后，再过判决，得到原信号。

3、3运行环境硬件环境：CPU、内存、硬盘、CD-ROM驱动器和鼠标软件环境：Windows98/NT/2000>WindowXP或Win74开发工具和编程语言MATLAB开发实验箱MATLAB汇编语言5详细设计第一部分：生成原信号clearall；clc；警设置时间参数、主频、偏频、抽样频率t=0:0.0005:9.9995;F=20;f=100;FS=l:500;%产生二进制信号和模拟信号a=randint(1,10,2)；m=a(ceil(t+0.0005))；s=cos(2*pi*(F+m*f

4、).*t);皆求信号的频谱和功率谱S=fft(s)；pss=S.*conj(S)/512；forj=1:2000Pss(j)=pss(10*j);end労画出二进制信号、模拟信号、功率谱的图像figure(1)；subplot(311)；plot(tzm)；axis([0,10,0,1.2])；title('~进制倍源波形J;xlabel('时间')；ylabel(1幅度■);subplot(312);plot(t,s);titleC模拟信源波形J；xlabel('时间')；ylabel(1幅度■

5、);subplot(313);plot(FSzPss(l:500));title(，原信号功率谱，);xlabel「频率T；ylabel(-幅度■)；第二部分：经过信道加入噪声咎生成噪声信号与合成信号n=0•01*randn(1,length(t));z=s+n；咎画出噪声信号图形、合成信号图形、原信号图形figure(2)；subplot(311)；plot(tzn)；title('噪声波形');xlabel('时间')；ylabel('幅度J；subplot(312)；plot(tzz)；ti

6、tle('合成信号波形');xlabel('时间')；ylabel('幅度T;subplot(313);plot(tzs);title('原信号波形');xlabel('时间')；ylabel('幅度J；咎对合成信号做快速傅立叶变换务并倆出其功率谱及原信号功率谱Z=fft(z);pzz=Z.*conj(Z)/512；forj=l：2000Pzz(j)=pss(10*j)；endfigure(3);subplot(211)；plot(FS,Pss(1:500));title('原信号功率谱T;xla

7、bel('频率')；ylabel('幅度J；subplot(212);plot(FS,Pzz(1:500));titlef合成信号功率谱J;xlabel(

8、频率J；ylabel('幅度J；第三部分：对接收到的信号滤波、判决，进而得到原信号咎求滤波器参数、并将信号通过滤波器耸分出两个不同频率的信号bl=firl(101,[15/1000,25/1000]);b2=firl(101,[115/1000,125/1000]);hl=filter(bl/l/z)；h2=filter(b2,1,z)；号作出

9、两个不同频率信号及原信号的图形figure(4)；subplot(311)；plot(tzhl)；title(，经通过频率为2OHZ的滤波器示的波形T;xlabel('时间')；ylabel('幅度■)；subplot(312)；plot(t,h2)；title(，经通过频率为120HZ的滤波器后的波形T;xlabel('时间')；ylabel('幅度J；subplot(313)；plot(t,s)；title('原信号波形');xlabel('时间')；ylabel('幅度J；咎