MATLAB仿真 BPSK调制

matlab BPSK 调制与解调
1、调制
clear all;
g=[1 0 1 0 1 0 0 1];%基带信号
f=100; %载波频率
t=0:2*pi/99:2*pi;
cp=[];sp=[];
mod=[];mod1=[];bit=[];
for n=1:length(g);
if g(n)==0;
die=-ones(1,100); %Modulante
se=zeros(1,100); %
else g(n)==1;
die=ones(1,100); %Modulante
se=ones(1,100); %
end
c=sin(f*t);
cp=[cp die];
mod=[mod c];
bit=[bit se];
end
bpsk=cp.*mod;
subplot(2,1,1);plot(bit,'LineWidth',1.5);grid on; title('Binary Signal');
axis([0 100*length(g) -2.5 2.5]);
subplot(2,1,2);plot(bpsk,'LineWidth',1.5);grid on; title('ASK modulation');
axis([0 100*length(g) -2.5 2.5]);
2、调制解调加噪声
clc;
close all;
clear; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 假定：
% 2倍载波频率采样的bpsk信号
% 调制速率为在波频率的 N/2m
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
m=128;
N=512;
n=1:1:N;
N0=0.5*randn(1,N) %噪声
h0=zeros(1,N);
% 30阶低通滤波器 h0
f = [0 0.3 0.3 1]; w0 = [1 1 0 0];
b = fir2(30,f,w0);
[h,w] = freqz(b,1,N/2);
h0(1,1:N/2)=abs(h');
for i=1:N/2
h0(1,N-i+1)=h0(1,i);
end;
%%%%%%%%% 随机序列
a=rand(1,m);
for i=1:m
if(a(1,i)>0.5)
a(1,i)=1;
else
a(1,i)=-1;
end;
end;
%%% 生成BPSK信号
bpsk_m=zeros(1,N);
j=1;k=1;
for i=1:N
if(j==(N/m+1))
j=1;
k=k+1;
end; % 0.05*pi 为初始相位，可以任意改变
bpsk_m(1,i)=a(1,k)*sin(2*pi*0.5*i+0.05*pi)+a(1,k)*cos(2*pi*0.5* i+0.05*pi);
j=j+1;
end;
bpsk_m=bpsk_m+N0;% 信号加噪声，模拟过信道
% 接收处理用正交本振与信号相乘，变频
bpsk_m1=bpsk_m.*sin(2*pi*0.5*n);
bpsk_m2=bpsk_m.*cos(2*pi*0.5*n);
%滤波
tempx=fft(bpsk_m1);
tempx=tempx.*h0; %低通滤波
tempx=ifft(tempx);
real_x=real(tempx);
tempx=h0.*fft(bpsk_m2);
tempx=tempx.*h0; %低通滤波
tempx=ifft(tempx);
real_x1=real(tempx);
subplot(2,1,1);
plot(real_x1+real_x,'b');
axis([1 N -2.5 2.5]);
grid on;
hold on;
In=real_x1+real_x; % 可只取一路，这里取了两路之和for i=1:N % 滤波后整形
if(In(1,i)>0) % 判决，得到解调结果
In(1,i)=1;
else
In(1,i)=-1;
end;
end;
plot(In,'r');
an=zeros(1,m);
for i=1:m
an(1,i)=In(1,(i-1)*N/m+N/(2*m));
end;
subplot(2,1,2); % 比较误码
plot(an,'r*');hold on;
axis([1 m -2 2]);
plot(a,'b^');。