MATLAB仿真BPSK调制_原创精品文档.pdfVIP

MATLAB仿真BPSK调制

matlabBPSK调制与解调

1、调制

clearall;

g=[10101001];%基带信号

f=100;%载波频率

t=0:2*pi/99:2*pi;

cp=[];sp=[];

mod=[];mod1=[];bit=[];

forn=1:length(g);

ifg(n)==0;

die=-ones(1,100);%Modulante

se=zeros(1,100);%

elseg(n)==1;

die=ones(1,100);%Modulante

se=ones(1,100);%

end

c=sin(f*t);

cp=[cpdie];

mod=[modc];

bit=[bitse];

end

bpsk=cp.*mod;

subplot(2,1,1);plot(bit,LineWidth,1.5);gridon;title(Binary

Signal);

axis([0100*length(g)-2.52.5]);

subplot(2,1,2);plot(bpsk,LineWidth,1.5);gridon;title(ASK

modulation);

axis([0100*length(g)-2.52.5]);

2、调制解调加噪声

clc;

closeall;

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=[00.30.31];w0=[1100];

b=fir2(30,f,w0);

[h,w]=freqz(b,1,N/2);

h0(1,1:N/2)=abs(h);

fori=1:N/2

h0(1,N-i+1)=h0(1,i);

end;

%%%%%%%%%随机序列

a=rand(1,m);

fori=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;

fori=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;%