Matlab对于2XSK信号的产生处理

【实例简介】本实例对通信原理中2ASK、2PSK、2FSK以及2DPSK信号都做了类似的处理，描述了信号调制、加噪声、到最后解调比较的全过程
【实例截图】
【核心代码】a=randint(1,20);
figure(1)
m=a(ceil(t 0.01));
plot(t,m)
title('产生随机八比特二进制比特序列');
axis([0,8,-1.5,1.5]);
%*************调制************%
fc=4000;          %载波频率
fs=80000;         %采样频率
ts=0:1/fs:(800-1)/fs;
ts1=0:1/fs:(100-1)/fs;
tzxh1=cos(2*pi*fc*ts);
%**************2ASK调制************%
ask=m.*tzxh1;
figure(2)
plot(t,ask)
title('2ASK调制波形');
axis([0,8,-1.5,1.5]);
%**************叠加噪声************%
e_2ask=awgn(ask,10);
figure(3)
plot(t,e_2ask)
title('2ASK调制信号叠加噪声波形');
%**************2ASK相干解调************%
[b11,a11]=ellip(5,0.5,60,[2000,6000]*2/80000);%带通椭圆滤波器设计
e_ask1=filter(b11,a11,e_2ask);%相乘器
e_ask2=2*e_ask1.*tzxh1*2;%低通滤波器
[b12,a12]=ellip(5,0.5,60,1000*2/80000);       %低通滤波器设计
ask_xgjt=filter(b12,a12,e_ask2);       %相干解调后，抽样判决前的结果
figure(4)
subplot(2,1,1)
plot(t,ask_xgjt)
title('2ASK调制信号相干解调后通过低通滤波器，抽样判决前的信号');
axis([0,8,-1.5,1.5]);
%**************2ASK非相干解调************%  
e_ask1_f=abs(e_ask1);%包络检波
ask_fxgjt=filter(b12,a12,e_ask1_f);       %抽样判决前的结果
subplot(2,1,2)
plot(t,ask_fxgjt)
title('2ASK调制信号非相干解调后通过低通滤波器，抽样判决前的信号');
axis([0,8,-1.5,1.5]);
%**************2ASK的相干解调法的抽样判决结果与原数据比较*************%
for i=0:7
    if ask_xgjt(1,(i 1)*100)>0.5
        ask_hyjt(i*100 1:i*100 100)=ones(1,100);
    else
        ask_hyjt(i*100 1:i*100 100)=zeros(1,100); 
    end
end   
figure(5)
subplot(2,1,1);
plot(t,ask_hyjt)
title('2ASK调制信号相干解调，抽样判决后的信号');
axis([0,8,-1.5,1.5]);
%**************2ASK的非相干解调法的抽样判决结果与原数据比较*************%
for i=0:7
    if ask_fxgjt((i 1)*100)>0.5
        ask_fhyjt(i*100 1:i*100 100)=ones(1,100);
    else       
        ask_fhyjt(i*100 1:i*100 100)=zeros(1,100);  
    end
end   
subplot(2,1,2);
plot(t,ask_fhyjt)
title('2ASK调制信号非相干解调，抽样判决后的信号');
axis([0,8,-1.5,1.5]);