qpsk.m

matlab仿真程序

clc;
clear all
N=16;
g=0;
m=1;
k=1;
I=[];
Q=[];
fc=1;
t=0:0.01:0.99;
r1=cos(2*pi*fc*t);
stem(r1)
r2=sin(2*pi*fc*t);
x=round(rand([1,N]));
for i=1:16
    if(x(i)==1)
        x(i)=1;
    else x(i)=-1;        
    end   
end
for i=1:16
    if(mod(i,2)==0)
       Q(m)=x(i);
        m=m+1;    
    else
         I(k)=x(i);
         k=k+1;
    end
end 
for i=1:8
    h=1;
    for j=1+g:100+g
       I2(j)=I(i).*r1(h);    
       Q2(j)=Q(i).*r2(h);
       h=h+1;
    end
    g=g+100;
end
g=0;
d=1;
for d=1:16;
if(mod(d,2)~=0)
    for j=1+g:100+g
        y(j)=I2(j-100*floor(d/2));
    end
    g=g+100;
   
else
    for j=1+g:100+g
        y(j)=Q2(j-100*floor(d/2));
    end
    g=g+100;
    
end
end
figure(1)
stem(x)
figure(2)
plot(y)