通信原理的MATLAB实验-QPSK的调制解调报告.doc

通信原理实验
项目名称：QPSK的调制解调
一、实验任务
任意输入长度为64比特的二进制信息，采用QPSK系统传输。码元速率为1Bps，载波频率为10Hz，采样频率为40 Hz， 利用Matlab画出：
（1）调制后的信号波形；
（2）经信道传输后的信号波形（假设加性高斯白噪声，其功率为信号功率1/10）；（3）（3）任意解调方法解调后的信号波形。

二、流程图

开始
为变量赋初值
四相调制产生QPSK
画出QPSK波形图
输入64比特的二进制信息
加入高斯白噪声
画出加入噪声之后的波形
输入相干解调
结束

三、完整程序
Fd=1; %码元速率
Fc=10; %载波频率
Fs=40; %采样频率
N=Fs/Fd;
df=10;
x=[ 1 1 0 1 1 0];%任意输入64比特的二进制信息
M=2; %进制数
SNRpBit=10;%加性高斯白噪声，其功率为信号功率的1/10，即信噪比为10
SNR=SNRpBit/log2(M); %转换为码元速率
seed=[12345 54321];
numPlot=length(x);
figure(1)%画出输入二进制序列
subplot(211);
stem([0:numPlot-1],x(1:numPlot),'bx');
title('输入波形’)
%调制
y=dmod(x,Fc,Fd,Fs,'fsk',M,df);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
subplot(212);%画出调制后的信号
plot(t,y(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('调制后的信号')
%在已调信号中加入高斯白噪声
randn('state',seed(2));
y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured',[],'dB');
%相干解调
figure(2)
subplot(211);
plot(t,y(1:length(t)),'b-');%画出经过信道的实际信号
axis([min(t) max(t) -1.5 1.5]);
title('加入高斯白噪声后的已调信号')
%带输出波形的相干M元频移键控解调
subplot(212);
stem([0:numPlot-1],x(1:numPlot),'bx');
hold on;
stem([0:numPlot-1],z1(1:numPlot),'ro');
hold off;
axis([0 numPlot -0.5 1.5]);
title('相干解调后的信号')
四、波形