qpsk_raylaigh1.asv

qpsk在raylaigh信道下的误码率曲线

%Program 1
%qpsk_fading.m
%******************************* Preparation part ************************  
sr=1000; %设码元速率为1000
ml=2;
fc=10;    %fc是载波频率
nd=1000;  %设置传输码元个数，
br=sr*ml;  %bit率
ebn0=5;   %EB/N0=5
%****************************** data generation **************************
data1=rand(1,1000*ml)>0.5;  %产生服从0—1分布随机数
%***************************** qpsk modulation ***************************
t=-500*1/sr+1/sr:1/sr:500*1/sr;
n=length(t);
data2=data1*2-1;        %转化成－1，1的序列
[ich,qch]=qpskmod(data2,1,nd,ml);
tich=ich.*cos(2*pi*fc*t); %Qpsk调制，I signal码元乘载波
tqch=qch.*sin(2*pi*fc*t); %Q signal码元乘载波
ampitude_modulated=tich+tqch;
                               %%%提取星座图1
%**************************** Rayleigh fading *************************
[tich2,tqch2]=ray(tich,tqch); %调用Rayleigh fading的产生
n1=sqrt(tich2.*tich2+tqch2.*tqch2);  %Rayleigh fading I和Q信号的幅度均方值
tich2=tich.*n1;
tqch2=tqch.*n1;
%***************************** Attenuation Calculation******************
spow=sum(tich.*tich+tqch.*tqch)/nd;
attn=0.5*spow*sr/br*10.^(-ebn0/10);
attn=sqrt(attn);
%***************************** AGWN ************************************
[ich2,qch2]=comb(tich2,tqch2,attn); %加性高斯白噪声   
recieved_amplitude=ich2+qch2;
                             %%%提取星座图2
%**************************** qpsk demodulation ************************
rich=ich2.*cos(2*pi*fc*t);       %对I信号分量解调
rqch=qch2.*sin(2*pi*fc*t);       
ampitude_demodulated=rich+rqch;
%[demodata]=qpskdemod(ich2,qch2,1,nd,ml);
%*************************** BER ********************************
%demodata=zeros(1,ml*nd);      
demodata((1:1),(1:ml:ml*nd-1))=rich((1:1),(1:nd))>=0;  %对I分量判决
demodata((1:1),(2:ml:ml*nd))=rqch((1:1),(1:nd))>=0;    %对Q分量判决
%noe=sum(abs(data1-demodata));  %计算误码率
%*************************** draw the figure **********************
figure(2); %%%%%source图形
plot(data2);
title('Source figure');

figure(3);
plot(t,ampitude_modulated);%%%qpsk调制后的图形
title('figure after modulation');
figure(100);   %%%%星座图1
plot(tich,tqch,'.');
title('constellation before channel');

figure(4);%%%%%%%%%%% qpsk+raylaigh_fading图形
plot(t,recieved_amplitude);
title('figuer qpsk+Raylaigh _ fading');
figure(101);   %%%%星座图2
plot(ich2,qch2,'.');
title('constellation after channel');

figure(5);%%%%%%%%%%%%%% 解调信号图形
plot(t,ampitude_demodulated);
title('figure after demodulation');

figure(6);
plot(demodata);
title('figure after demodulation11');