bpsk.m

自己编写的MATLAB对调制解调的仿真

%*******************modulation****************
D=[1 0 1 1 0 0 1 0];%产生随机数
a=[];
N=length(D);%基带信号符号数目
t=0:0.01:7.99;
y=cos(2*pi*t);
d0=D;d1=2*D-1;%单极性信号变双极性1 -1 1 1 -1 -1 1 -1
D0=[];
x=ones(1,100);%100个1
for k=1:N,
        D0=[D0,d1(k)*x];
end%开关信号,为便于观察
z=y.*D0;%bpsk的表达式
%*****************demodulation***************
recos=z.*y;
meg=find(recos>0);
ch(meg)=1;
meg=find(recos<0);
ch(meg)=-1;
demo=(ch+1)/2;
%***************drawing**********************
subplot(2,1,1);plot(z);grid ;title('BPSK modulation');axis([0,800,-2,2]);
subplot(2,1,2);plot(demo);grid ;title('BPSK demodulation');axis([0,800,-2,2]);

xlabel('1 0 1 1 0 0 1 0');