【数字信号调制】基于matlab二进制数字振幅调制与解调（2ASK）【含Matlab源码 996期】.zip

clear all close all i=5;%5个码元 j=5000; t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即分成5000份 fc=2;%载波频率 fm=i/4;%码元速率 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号 x=(rand(1,i))%rand函数产生在0-1之间随机数，共1-10个 a=round(x);%随机序列，round取最接近小数的整数 st=t; for n=1:i if a(n)<1; for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end figure(1); subplot(221); plot(t,st); axis([0,5,-0.2,1.2]); title('基带信号'); %载波 s1=cos(2*pi*fc*t); subplot(222); plot(t,s1); axis([0,5,-1,1]); title('载波信号'); %调制 e_2ask=st.*s1; subplot(223); plot(t,e_2ask); axis([0,5,-1,1]); title('已调信号'); %相干解调 at=e_2ask.*cos(2*pi*fc*t); at=at-mean(at);%因为是单极性波形，还有直流分量，应去掉 subplot(223); [f,af]= T2F(t,at);%通过低通滤波器 [t,at]= lpf(f,af,2*fm); %抽样判决 for m=0:i-1; if at(1,m*1000+500)+0.5<0.5; for j=m*1000+1:(m+1)*1000; at(1,j)=0; end else for j=m*1000+1:(m+1)*1000; at(1,j)=1; end end end subplot(224); plot(t,at); axis([0,5,-0.2,1.2]); title('相干解调后波形')