2FSK相干解调与非相干解调.docx
10级通信通信原理课程设计班 级:姓 名:学 号:指导教师:2013年1月6日星期日2FSK相干解调与非相干解调一.课程内容1.2FSK信号的产生2. 加入高斯白噪声的波形3. 上下支路FSK信号的相干解调与非相干解调二原理2FSK信号的解调原理是通过带通滤波器将2FSK信号分解为上下两路2FSK 信号后分别解调,然后进行抽样判决输出信号。三. 设计过程1. 相干解调相干解调是指利用乘法器,输入一路与载频相干的参考信号与载频相乘。根据已调信号由两个载波fl、f2调制而成,则先用两个分别对fl、f2带通 的滤波器对己调信号进行滤波,然后再分别将滤波后的信号与相应的载波fl、f2 相乘进行相干解调,再分别低通滤波、用抽样信号进行抽样判决器即可。原理图 如下:cos2n fitcos2n f2t相干方式原理图2. 非相干解调经过调制后的2FSK数字信号通过两个频率不同的带通滤波器f1、f2滤出不 需要的信号,然后再将这两种经过滤波的信号分别通过包络检波器检波,最后将两种信号同时输入到抽样判决器同时外加抽样脉冲,最后解调出来的信号就是调制前的输入信号。其原理图如下图所示:带通滤波器包络检波输入5u带通滤波器)包络检波器非相干方式原理图四. 实现效果1.相干解调仿真图咨口舍fe纹0心苞更 nEilt! Edi tIns er + Tools iLeEktop S.i ndow it&lp产生的波形0.5段 Ina-E-0.5将要通过滤波器的波形2HOI 0 12- - W-K宅都蝗Figure 2_|n| xiEile Edit Vi ew Ijiser t lools HeEktup Itindow Help q a 鼠 矿直 经过低通滤波器洲1后的波形-|a| x|Eile Edit Vi ew Ineert Tools liesktop itinilow Help小 a: h aTM o ® cFn e0.80.6他04暨0.212345678910经过低通滤波器辿2后的波形-|n| x| h e fe q n ® « 口Eile Edit 归 ew Insert Io ol s Be ski op Btimlow Help2.非相干解调仿真图-! xEile Edit View Ins er + lools LLeskt op Window Help也 q e玲T 蚁 0咨©匐 bT 口经过带通滤波器f2后的波形012345678910t五. 附录(原程序)fs=2000;%采样频率fc=900;dt=1/fs;f1=20;f2=120;%两个信号的频率a=round(rand(1,10);%随机信号g1=ag2=a;%信号反转,和g1反向g11=(ones(1,2000)'*g1;% 抽样g1a=g11(:)'g21=(ones(1,2000)'*g2;g2a=g21(:)'t=0:dt:10-dt;t1=length(t);fsk1=g1a.*cos(2*pi*f1.*t);fsk2=g2a.*cos(2*pi*f2.*t);fsk=fsk1+fsk2;%产生的信号sn=awgn(fsk,20);subplot(211);plot(t,fsk);title('产生的波形')ylabel('幅度')subplot(212);plot(t,sn);title ('将要通过滤波器的波形')ylabel('幅度的大小')xlabel('t')figure(2) %FSK 解调 b1=fir1(101,10/800 20/800);b2=fir1(101,90/800 110/800); % 设置带通参数H1=filter(b1,1,sn);H2=filter(b2,1,sn);%经过带通滤波器后的信号subplot(211); plot(t,H1); title('经过带通滤波器f1后的波形') ylabel('幅度') xlabel('t') subplot(212);plot(t,H2);title('经过带通滤波器f2后的波形')ylabel('幅度')xlabel('t')sw1=H1.*H1;sw2=H2.*H2; %经过相乘器figure(3)subplot(211);plot(t,sw1);title(经过相乘器h1后的波形')ylabel('幅度')subplot(212);plot(t,sw2);title(经过相乘器h2后的波形')ylabel(' 幅度')xlabel('t')bn=fir1(101,2/800 10/800);%经过低通滤波器figure(4)st1=filter(bn,1,sw1);st2=filter(bn,1,sw2);subplot(211);plot(t,st1);title('经过低通滤波器sw1后的波形')ylabel('幅度')subplot(212);plot(t,st2);title('经过低通滤波器sw2后的波形')ylabel('幅度')xlabel('t')%判决for i=1:length(t)if(st1(i)>=st2(i) st(i)=0;else st(i)=st2(i);endendfigure(5)st=st1+st2;subplot(211);plot(t,st);title(经过抽样判决器后的波形')ylabel('幅度')subplot(212);plot(t,sn);title('原始的波形')ylabel('幅度')xlabel('t')%非相干figure(6) %FSK 解调 b1=fir1(101,10/800 20/800);b2=fir1(101,90/800 110/800); % 设置带通参数H1=filter(b1,1,sn);H2=filter(b2,1,sn);%经过带通滤波器后的信号subplot(211); plot(t,H1); title('经过带通滤波器f1后的波形') ylabel('幅度') xlabel('t') subplot(212);plot(t,H2);title('经过带通滤波器f2后的波形')ylabel('幅度')xlabel('t')figure(7)st1=ssbdemod(H1,fc,fs,2./pi);% 包络解波 subplot(211);plot(t,st1);title('经过包络检波的波形')ylabel('幅度')xlabel('t')st2=ssbdemod(H2,fc,fs,2./pi);% 包络解波 subplot(212);plot(t,st2);title('经过包络检波的波形')ylabel('幅度')xlabel('t')%判决for i=1:length(t)if(st1(i)>=st2(i)st(i)=0;else st(i)=st2(i);endendfigure(8)st=st1+st2;subplot(211);plot(t,st);title(经过抽样判决器后的波形')ylabel('幅度')subplot(212);plot(t,sn);title('原始的波形')ylabel('幅度')xlabel('t')