通信原理实验项目名称:QPSK的调制解调一、实验任务任意输入长度为64比特的二进制信息,采用QPSK系统传输。码元速率为1Bps,载波频率为10Hz,采样频率为40Hz,利用Matlab画出:(1)调制后的信号波形;(2)经信道传输后的信号波形(假设加性高斯白噪声,其功率为信号功率1/10);(3)(3)任意解调方法解调后的信号波形。二、流程图开始为变量赋初值四相调制产生QPSK画出QPSK波形图输入64比特的二进制信息加入高斯白噪声画出加入噪声之后的波形输入相干解调结束三、完整程序Fd=1;%码元速率Fc=10;%载波频率Fs=40;%采样频率N=Fs/Fd;df=10;x=[110110];%任意输入64比特的二进制信息M=2;%进制数SNRpBit=10;%加性高斯白噪声,其功率为信号功率的1/10,即信噪比为10SNR=SNRpBit/log2(M);%转换为码元速率seed=[1234554321];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.51.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.51.5]);title('加入高斯白噪声后的已调信号')%带输出波形的相干M元频移键控解调subplot(212);stem([0:numPlot-1],x(1:numPlot),'bx');holdon;stem([0:numPlot-1],z1(1:numPlot),'ro');holdoff;axis([0numPlot-0.51.5]);title('相干解调后的信号')四、波形