基于MATLAB的2FSK数字通信系统设计.doc

青 岛 农 业 大 学 理学与信息科学学院 通 信 原 理 课 程 设 计 设 计 题 目 学生专业班级 学生姓名（学号） 指 导 教 师 完 成 时 间 实 习（设计）地点 2011年月日  四、实现方法 通过在Matlab环境下利用相关程序对二进制信号进行抽样产生二进制随机序列，并在此时确定抽样间隔、抽样个数，进行显示。然后对二进制信号进行2FSK调制，利用Matlab函数库中awgn函数在信道中加入高斯白噪声，并对此时的信号进行显示，接着进行相干解调，输出显示解调后序列，对比信号调制前后的变化，统计误码率，并分析调制优劣。以下为详细过程： 调整模块 使用键控法产生2FSK信号，即由电子开关在两个独立的频率源之间转换形成，这样的相邻码元之间相位不一定连续。 产生2FSK调制信号原理图： Fc=10; %载频 Fs=40; %系统采样频率 Fd=1; %码速率 N=Fs/Fd; df=10; numSymb=25; %进行仿真的信息代码个数 M=2; %进制数 SNRpBit=60; %信噪比 SNR=SNRpBit/log2(M); seed=[12345 54321]; numPlot=25; 如上图，使用产Bernoulli Binary Generator模块生二进制随机比特流，通过键控法进行2FSK调制。 figure(1) %产生25个二进制随机码 x=randsrc(numSymb,1,[0:M-1]);%产生25个二进制随机码 subplot(2,2,1) stem([0:numPlot-1],x(1:numPlot),'bx'); title('二进制随机序列') xlabel('Time'); ylabel('Amplitude'); 2、 传输模块 （1）2FSK调制阶段 在二进制数字调制中，若正弦载波的频率随二进制基带信号在f1和f2两个频率点间变化，则产生二进制移频键控信号（2FSK信号）。二进制移频键控信号的时间波形如下图所示： 图中波形g可分解为波形e和波形f，即二进制移频键控信号可以看成是两个不同载波的二进制振幅键控信号的叠加。 若二进制基带信号的1符号对应于载波频率f1，0符号对应于载波频率f2，则二进制移频键控信号的时域表达式为： %调制 y=dmod(x,Fc,Fd,Fs,'fsk',M,df); numModPlot=numPlot*Fs; t=[0:numModPlot-1]./Fs; subplot(2,2,2) plot(t,y(1:length(t)),'b-'); axis([min(t) max(t) -1.5 1.5]); title('调制后的信号') xlabel('Time'); ylabel('Amplitude'); （2）在已调信号中加入高斯白噪声 在二进制数字调制系统抗噪声性能分析中，假设信道特性是恒参信道，在信号的频带范围内其具有理想矩形的传输特性(可取传输系数为K)。 噪声为等效加性高斯白噪声， 其均值为零，方差为σ2。 加入高斯噪声信道的总图如下 其中，Scope图为原随机二进制比特流图形，Scope3图为高频率f1信号调制后的波形，Scope2图为低频率f2信号调制后的波形，Scope14图为经过2FSK调制后未加高斯白噪声前的波形图，Scope13为经过加有高斯白噪声信道后的波形图，Scope5图为经过解调抽样判决后的波形信号。 %在已调信号中加入高斯白噪声 randn('state',seed(2)); y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured',[],'dB'); %在已调信号中加入高斯白噪声 subplot(2,2,4) plot(t,y(1:length(t)),'b-'); %画出经过信道的实际信号 axis([min(t) max(t) -1.5 1.5]); title('加入高斯白噪声后的已调信号') xlabel('Time'); ylabel('Amplitude'); 3、 解调模块 （1）解调输出 本设计使用相干解调。其解调原理是将2FSK信号分解为上下两路2ASK信号分别进行解调，然后进行判决。抽样判决是直接比较两路信号抽样值的大小，判决规则应与调制规则相呼应，调制时若规定“1”符号对应载波频率f1，则接收时上支路的样值较大，应判为“1”；反之则判为“0”。 将加入的高斯白噪声滤除，并抽样判决出解调