MATLAB仿真的ASKPSKFSK性能比较和分析报告Word格式文档下载.docx

1、2 / 26 1.2 2FSK 2FSK可以看做是2个不同频率的2ASK的叠加，其调制与解调方法与2ASK差不多，主要频率F1和F2，不同的组合产生所要求的2FSK调制信号。cosA1? 公式如下：k1s（t）?FSK2?t，当a?0cosA? k21.3 2PSK 2PSK以载波的相位变化为基准，载波的相位随数字基带序列信号的1或者0而改变，通常用已经调制完的载波的0或者表示数据1或者0，每种相位与之一3 / 26 一对应。二数字调制技术的仿真实现 本课程设计需要借助MATLAB的M文件编程功能，对2ASK.2PSK.2FSK进行调制与解调的设计，并绘制出调制与解调后的波形，误码率的情况分析

2、，软件仿真可在已有平台上实现。1.2ASK代码主函数 allclose allclear n=16; fc=1000000; bitRate=1000000; N=50; %noise=ti; noise=10; 生成二进制代码signal=source（n,N）; % 调制后信号transmittedSignal=askModu（signal,bitRate,fc,N）;% 加噪声signal1=gussian（transmittedSignal,noise）;% configueSignal=demoASK（signal1,bitRate,fc,n,N）;source代码 sendSign

3、al=source（n,N）function4 / 26 sendSignal=randint（1,n） bit=; i=1:length（sendSignal） for sendSignal（i）=0 if bit1=zeros（1,N）; else bit1=ones（1,N）; end bit=bit,bit1; figure（1） ;on plot（1:length（bit）,bit）,title（transmitting of binary）,grid axis（0,N*length（sendSignal）,-2,2）;askModu代码 function transmittedSi

4、gnal=askModu（signal,bitRate,fc,N）%signal为输入信号，bitrate为bit速率，fc调制信号频率，N %signal=0 0 1 0 1 1 0 1; % bitRate=1000000; % fc=1000000;% N=32; t=linspace（0,1/bitRate,N）; c=sin（2*pi*t*fc）; transmittedSignal=;length（signal） for transmittedSignal=transmittedSignal,signal（i）*c; 画调制图figure（2） %Modulation plot（1

5、:length（transmittedSignal）,transmittedSignal）;title（ ;）;grid onof ASK 画频谱实部figure（3）% m=0:length（transmittedSignal）-1; F=fft（transmittedSignal）; ）;plot（m,abs（real（F）,title（ASK_frequency-domain analysis real ;grid on 画频谱虚部 %figure（4） %plot（m,imag（F）;title（ASK_frequency-domain analysis imag %grid on;

6、endCheckRatePe代码 5 / 26 PeWrong=CheckRatePe（signal1,signal2,s）function rights=0; wrongs=0; ki=1:s-2for （signal1（ki）=signal2（ki） if rights=rights+1; else wrongs=wrongs+1; end PeWrong=wrongs/（wrongs+rights）;demoASK代码 bitstream=demoASK（receivedSignal,bitRate,fc,n,N） function num load signal1=receivedSi

7、gnal; %? signal2=abs（signal1）; %LPF,2 signal3=filter（num1,1,signal2）;3 IN=fix（length（num1）/2）; bitstream=; LL=fc/bitRate*N; i=IN+LL/2;D?while （i=0.5; i=i+LL; figure（6） 接收波形 subplot（3,1,1）; % plot（1:length（signal1）,signal1）;Wave of receiving onterminal（including noise）grid 接收整流后波形 subplot（3,1,2）;% ;g

8、rid oncommutate plot（1:length（signal2）,signal2）;Wave of 包络检波波形 subplot（3,1,3）;length（signal3）,signal3）;Wave of LPF bit=;length（bitstream） for bitstream（i）=0if bit1=zeros（1,N）; else6 / 26 bit1=ones（1,N）; 解调后的二进制波形 figure（7）% ;）,grid on plot（bit）,title（binary of receiving terminal axis（0,N*length（bits

9、tream）,-2.5,2.5）;gussian代码 %加高斯白噪声 signal=gussian（transmittedSignal,noise）function signal=sqrt（2）*transmittedSignal; signal=awgn（signal,noise）; figure（5） plot（1:length（signal）,signal）;onWave including noise）,grid title（ end/ fsk主函数代码/ 二进制代码长度n=16;% 1f1=18000000;%频率 2频率f2=6000000;% 速率bitRate=1000000;

10、%bit 码元宽度N=50;% %noise=ti; 家性噪声大小noise=10;% 产生二进制代码signal=source（n,N）;% 调制transmittedSignal=fskModu（signal,bitRate,f1,f2,N）;% 解调configueSignal=demoFSK（signal1,bitRate,f1,f2,N）;% source代码%二进制信号产生函数 sendSignal=source（n,N）function sendSignal=randint（1,n） bit=;length（sendSignal） for sendSignal（i）=0 if7

11、/ 26 axis（0,N*length（sendSignal）,-2.5,2.5）;fskModu代码%频率调制函数 transmittedSignal=fskModu（signal,bitRate,f1,f2,N）function t=linspace（0,1/bitRate,N）; 1调制信号 c1=sin（2*pi*t*f1）;% 2调制信号 c2=sin（2*pi*t*f2）;% transmittedSignal=; 调制 i=1:length（signal）% for signal（i）=1 if transmittedSignal=transmittedSignal,c1; t

12、ransmittedSignal=transmittedSignal,c2; 画调制后波形图figure（2） % plot（1:Modulation onof FSKfigure（3） %画调制后频谱图 m=0:ongrid demoFSK代码 bitstream=demoFSK（receivedSignal,bitRate,f1,f2,N） function num load ，得到高通分量通过 signal2=filter（gaotong,1,signal1）; %HPF8 / 26 整流% signal3=abs（signal2）; 通过低通，形成包络% signal3=filter（

13、lowpass,1,signal3）; IN1=fix（length（lowpass）/2）+fix（length（gaotong）/2）; %延迟时 间 bitstream1=; 的抽样点数每个bit LL=N; % i=IN1 +LL/2; 判决%while （i bitstream1 figure（5） subplot（3,1,1）;grid onterminal（including noise） subplot（3,1,2）;grid onAfter Passing HPFtitle（ subplot（3,1,3）;onPassing LPFgrid plot（1:After 得到低频

14、分量BPF signal4=filter（daitong,1,signal1）; %通过 整流% signal5=abs（signal4）; ，形成包络LPF%通过 signal5=filter（lowpass,1,signal5）; 延迟时间% IN2=fix（length（lowpass）/2）+fix（length（daitong）/2）; bitstream2=; bit的的抽样点数 LL=N; %每个 i=IN2 +LL/2; 判决 （i bitstream2 figure（6） subplot（3,1,1）;grid terminal（including noise）on subp

15、lot（3,1,2）;onPassing After BPFlength（signal4）,signal4）;on）;length（signal5）,signal5）;After Passing LPF 判决 i=1:min（length（bitstream1）,length（bitstream2） for%9 / 26 （bitstream1（i）bitstream2（i） if bitstream（i）=1; bitstream（i）=0; bitstream 接收端波形% bit=;length（bitstream）for bitstream（i）=0if else figure（7）

16、; end CheckRatePe代码 s-2for （signal1（ki）=signal2（ki）if rights=rights+1;gussian代码 figure（4） ）;length（signal）,signal）,title（Adding Noise10 / 26 grid on end2psk主函数代码 二进制码长n=16;% 载波频率fc=1000000;% 信息频率bitRate=1000000; 码宽N=50;% 信道加性噪声大小noise=10;% 生成二进制代码signal=source（n,N）;transmittedSignal=bpskModu（signal,

17、bitRate,fc,N）;对信号进行调制并进行频谱 分析 加信道噪声signal1=gussian（transmittedSignal,noise）% 信号解调configueSignal=demoBPSK（signal1,bitRate,fc,n,N）;length（sendSignal） for sendSignal（i）=0if on）,grid plot（bit）,title（bpskModu代码 transmittedSignal=bpskModu（signal,bitRate,fc,N）function t=linspace（0,1/bitRate,N）; c1=sin（2*pi

18、*t*fc）; c2=sin（2*pi*t*fc + pi）;11 / 26 length（signal） for signal（i）=1 if transmittedSignal=transmittedSignal,c1; 画调制图figure（2） % onof BPSK 画频谱图figure（3）% m=0:BPSK_frequency-domain analysis realplot（m,abs（real（F）,title（ ;demoBPSK代码 bitstream=demoBPSK（receivedSignal,bitRate,fc,n,N） function 存储的低通滤波用的数

19、据num%读取 load num signal1=receivedSignal; signal=;nfor signal=signal,c;12 / 26 sin%乘同频同相 signal2=signal1.*signal; 3%LPF,包络检波 signal3=filter（num1,1,signal2）; 延迟时间 IN=fix（length（num1）/2）;=length（signal3） % while bitstream=bitstream,signal3（i）=0; figure（5） 画接收的包含噪声的波形 subplot（3,1,1）;% plot（1: 相干解调波形 subplot（3,1,2）;After Multipling sin onFuctiontitle（ length（bitstream） for bitstream（i）=0 if bit1=zeros（1,N）; 二进制接收信号波形 figure（6） ;grid on plot（bit）; axis（0,N*length（bitstream）,-2.5,2