行业知识
通信原理与MATLAB(十):QPSK的调制解调
2026-11-04 06:23:09
来源：互联网转载
目录
- 1.QPSK的调制原理
- 2.QPSK的解调原理
- 3.QPSK代码
- 4.结果图
- 5.特点
- 6.加星座图的QPSK代码
1.QPSK的调制原理
QPSK调制原理如下图所示，QPSK相当于两个正交的BPSK相加而成。其调制原理是将基带码元分成I、Q两路，I路是原始基带码元的奇数位置码元，Q路是原始基带码元的偶数位置码元，然后两条支路分别和对应的载波相乘实现BPSK的调制，然后将两条支路相加实现QPSK的调制。
2.QPSK的解调原理
QPSK的解调原理如下图所示，DPSK信号再分为I、Q两路和对应的载波相乘，然后经过低通滤波器后进行抽样判决，相当于作两路的BPSK解调。判决之后的I、Q路码元进行合并，I路为最终码元序列的奇数位置码元，Q路为最终码元序列的偶数位置码元，恢复出原始的码元序列。
3.QPSK代码
clear all; % 清除所有变量close all; % 关闭所有窗口clc; % 清屏%% 基本参数M=20; % 产生码元数 L=100; % 每码元复制L次,每个码元采样次数Ts=0.001; % 每个码元的宽度,即码元的持续时间Rb=1/Ts; % 码元速率1Kdt=Ts/L; % 采样间隔TotalT=M*Ts; % 总时间t=0:dt:TotalT-dt; % 时间TotalT2=(M/2)*Ts; % 总时间2t2=0:dt:TotalT2-dt; % 时间2Fs=1/dt; % 采样间隔的倒数即采样频率%% 产生双极性波形wave=randi([0,1],1,M); % 产生二进制随机码,M为码元个数wave=2*wave-1; % 单极性变双极性fz=ones(1,L); % 定义复制的次数L,L为每码元的采样点数x1=wave(fz,:); % 将原来wave的第一行复制L次，称为L*M的矩阵jidai=reshape(x1,1,L*M); % 产生双极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵%% I、Q路码元% I路码元是基带码元的奇数位置码元，Q路码元是基带码元的偶数位置码元I=[]; Q=[];for i=1:M if mod(i, 2)~=0 I = [I, wave(i)]; else Q = [Q, wave(i)]; endendx2 = I(fz,:); % 将原来I的第一行复制L次，称为L*(M/2)的矩阵I_lu = reshape(x2,1,L*(M/2));% 将刚得到的L*(M/2)矩阵，按列重新排列形成1*(L*(M/2))的矩阵x3 = Q(fz,:); % 将原来Q的第一行复制L次，称为L*(M/2)的矩阵Q_lu = reshape(x3,1,L*(M/2));% 将刚得到的L*(M/2)矩阵，按列重新排列形成1*(L*(M/2))的矩阵figure(1); % 绘制第1幅图subplot(311); % 窗口分割成3*1的，当前是第1个子图 plot(t,jidai,'LineWidth',2);% 绘制基带码元波形，线宽为2title('基带信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT,-1.1,1.1]) % 坐标范围限制subplot(312); % 窗口分割成3*1的，当前是第2个子图 plot(t2,I_lu,'LineWidth',2);% 绘制基带码元波形，线宽为2title('I路信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-1.1,1.1]) % 坐标范围限制subplot(313); % 窗口分割成3*1的，当前是第3个子图 plot(t2,Q_lu,'LineWidth',2);% 绘制基带码元波形，线宽为2title('Q路信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-1.1,1.1]) % 坐标范围限制%% QPSK调制fc=2000; % 载波频率2kHz zb1=cos(2*pi*fc*t2); % 载波1psk1=I_lu.*zb1; % PSK1的调制 zb2=sin(2*pi*fc*t2); % 载波2psk2=Q_lu.*zb2; % PSK2的调制 qpsk=psk1+psk2; % QPSK的实现figure(2); % 绘制第2幅图subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(t2,qpsk,'LineWidth',2);% 绘制基带码元波形，线宽为2title('QPSK信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-2.1,2.1]) % 坐标范围限制%% 信号经过高斯白噪声信道tz=awgn(qpsk,20); % 信号qpsk中加入白噪声，信噪比为SNR=20dBsubplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(t2,tz,'LineWidth',2); % 绘制2ASK信号加入白噪声的波形axis([0,TotalT2,-2.5,2.5]); % 坐标范围设置title('通过高斯白噪声信道后的信号');% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 解调部分figure(3);tz1=tz.*zb1; % 相干解调，乘以相干载波subplot(211) % 窗口分割成2*1的，当前是第1个子图 plot(t2,tz1,'LineWidth',2) % 绘制I路乘以相干载波后的信号axis([0,TotalT2,-2.5,2.5]); % 设置坐标范围title("I路乘以相干载波后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签tz2=tz.*zb2; % 相干解调，乘以相干载波subplot(212) % 窗口分割成2*1的，当前是第2个子图 plot(t2,tz2,'LineWidth',2) % 绘制Q路乘以相干载波后的信号axis([0,TotalT2,-2.5,2.5]); % 设置坐标范围title("Q路乘以相干载波后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 加噪信号经过滤波器% 低通滤波器设计fp=2*Rb; % 低通滤波器截止频率，乘以2是因为下面要将模拟频率转换成数字频率wp=Rb/(Fs/2)b=fir1(30, fp/Fs, boxcar(31));% 生成fir滤波器系统函数中分子多项式的系数% fir1函数三个参数分别是阶数，数字截止频率，滤波器类型% 这里是生成了30阶(31个抽头系数)的矩形窗滤波器[h,w]=freqz(b, 1,512); % 生成fir滤波器的频率响应% freqz函数的三个参数分别是滤波器系统函数的分子多项式的系数，分母多项式的系数(fir滤波器分母系数为1)和采样点数(默认)512lvbo1=fftfilt(b,tz1); % 对信号进行滤波，tz1是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数lvbo2=fftfilt(b,tz2); % 对信号进行滤波，tz2是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数figure(4); % 绘制第4幅图 subplot(311); % 窗口分割成3*1的，当前是第1个子图 plot(w/pi*Fs/2,20*log(abs(h)),'LineWidth',2); % 绘制滤波器的幅频响应title('低通滤波器的频谱'); % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度/dB'); % y轴标签subplot(312) % 窗口分割成3*1的，当前是第2个子图 plot(t2,lvbo1,'LineWidth',2); % 绘制I路经过低通滤波器后的信号axis([0,TotalT2,-2.1,2.1]); % 设置坐标范围title("I路经过低通滤波器后的信号");% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签subplot(313) % 窗口分割成3*1的，当前是第3个子图 plot(t2,lvbo2,'LineWidth',2); % 绘制Q路经过低通滤波器后的信号axis([0,TotalT2,-2.1,2.1]); % 设置坐标范围title("Q路经过低通滤波器后的信号");% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 抽样判决figure(5);k=0; % 设置抽样限值pdst1=1*(lvbo1>0); % 滤波后的向量的每个元素和0进行比较，大于0为1，否则为0pdst2=1*(lvbo2>0); % 滤波后的向量的每个元素和0进行比较，大于0为1，否则为0subplot(311) % 窗口分割成3*1的，当前是第1个子图 plot(t2,pdst1,'LineWidth',2)% 画出经过抽样判决后的信号axis([0,TotalT2,-0.1,1.1]); % 设置坐标范用title("I路经过抽样判决后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签subplot(312) % 窗口分割成3*1的，当前是第2个子图 plot(t2,pdst2,'LineWidth',2)% 画出经过抽样判决后的信号axis([0,TotalT2,-0.1,1.1]); % 设置坐标范用title("Q路经过抽样判决后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% I、Q路合并I_zong = [];Q_zong = [];% 取码元的中间位置上的值进行判决for j=L/2:L:(L*M/2) if pdst1(j)>0 I_zong=[I_zong,1]; else I_zong=[I_zong,-1]; endend% 取码元的中间位置上的值进行判决for k=L/2:L:(L*M/2) if pdst2(k)>0 Q_zong=[Q_zong,1]; else Q_zong=[Q_zong,-1]; endendcode = [];% 将I路码元为最终输出的奇数位置码元，将Q路码元为最终输出的偶数位置码元for n=1:M if mod(n, 2)~=0 code = [code, I_zong((n+1)/2)]; else code = [code, Q_zong(n/2)]; endendx4=code(fz,:); % 将原来code的第一行复制L次，称为L*M的矩阵dout=reshape(x4,1,L*M); % 产生单极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵subplot(313); % 窗口分割成3*1的，当前是第3个子图 plot(t,dout,'LineWidth',2);% 绘制基带码元波形，线宽为2title('I、Q路合并信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT,-1.1,1.1]) % 坐标范围限制%% 绘制频谱%% 信源频谱figure(6) % 绘制第6幅图T=t(end); % 时间df=1/T; % 频谱分辨率N=length(jidai); % 采样长度f=(-N/2:N/2-1)*df; % 频率范围mf=fftshift(abs(fft(jidai)));%对信源信号采用快速傅里叶变换并移到矩阵中心subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f,mf,'LineWidth',2); % 绘制信源频谱波形title("基带信号频谱"); % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% QPSK信号频谱T2=t2(end); % 时间2df2=1/T2; % 频谱分辨率2N2=length(qpsk); % 采样长度2f2=(-N2/2:N2/2-1)*df2; % 频率范围sf=fftshift(abs(fft(qpsk)));% 对QPSK信号采用快速傅里叶变换并将0-fs频谱移动到-fs/2-fs/2subplot(212) % 窗口分割成2*1的，当前是第2个子图 plot(f2,sf,'LineWidth',2) % 绘制QPSK调制信号频谱title("QPSK信号频谱") % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf]) % 坐标范围限制%% I路乘以相干载波后的频谱mmf=fftshift(abs(fft(tz1)));% 对相干载波信号采用快速傅里叶变换并移到矩阵中心figure(7) % 绘制第7幅图subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f2,mmf,'LineWidth',2) % 画出乘以相干载波后的频谱title("I路乘以相干载波后的频谱")xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% Q路乘以相干载波后的频谱mmf2=fftshift(abs(fft(tz2)));% 对相干载波信号采用快速傅里叶变换并移到矩阵中心subplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(f2,mmf2,'LineWidth',2) % 画出乘以相干载波后的频谱title("Q路乘以相干载波后的频谱")xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% 经过低通滤波后的频谱figure(8);dmf=fftshift(abs(fft(lvbo1)));%对低通滤波信号采用快速傅里叶变换并移到矩阵中心subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f2,dmf,'LineWidth',2) % 画出经过低通滤波后的频谱title("I路经过低通滤波后的频谱");xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签dmf2=fftshift(abs(fft(lvbo2)));%对低通滤波信号采用快速傅里叶变换并移到矩阵中心subplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(f2,dmf2,'LineWidth',2) % 画出经过低通滤波后的频谱title("Q路经过低通滤波后的频谱");xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签4.结果图
结果图中QPSK信号是经过信道，加了高斯白噪声的。
如果不想加噪声，把下面这行代码注释即可。
tz=awgn(qpsk,20); % 信号qpsk中加入白噪声，信噪比为SNR=20dB
5.特点
随着进制数增加，抗噪声性能下降，但是传输速率提高，牺牲通信index.php/tags-13.html" >系统的可靠性，提高有效性。
6.加星座图的QPSK代码
clear all; % 清除所有变量close all; % 关闭所有窗口clc; % 清屏%% 基本参数M=20; % 产生码元数 L=100; % 每码元复制L次,每个码元采样次数Ts=0.001; % 每个码元的宽度,即码元的持续时间Rb=1/Ts; % 码元速率1Kdt=Ts/L; % 采样间隔TotalT=M*Ts; % 总时间t=0:dt:TotalT-dt; % 时间TotalT2=(M/2)*Ts; % 总时间2t2=0:dt:TotalT2-dt; % 时间2Fs=1/dt; % 采样间隔的倒数即采样频率%% 产生双极性波形wave=randi([0,1],1,M); % 产生二进制随机码,M为码元个数wave=2*wave-1; % 单极性变双极性fz=ones(1,L); % 定义复制的次数L,L为每码元的采样点数x1=wave(fz,:); % 将原来wave的第一行复制L次，称为L*M的矩阵jidai=reshape(x1,1,L*M); % 产生双极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵%% I、Q路码元% I路码元是基带码元的奇数位置码元，Q路码元是基带码元的偶数位置码元I=[]; Q=[];for i=1:M if mod(i, 2)~=0 I = [I, wave(i)]; else Q = [Q, wave(i)]; endendx2 = I(fz,:); % 将原来I的第一行复制L次，称为L*(M/2)的矩阵I_lu = reshape(x2,1,L*(M/2));% 将刚得到的L*(M/2)矩阵，按列重新排列形成1*(L*(M/2))的矩阵x3 = Q(fz,:); % 将原来Q的第一行复制L次，称为L*(M/2)的矩阵Q_lu = reshape(x3,1,L*(M/2));% 将刚得到的L*(M/2)矩阵，按列重新排列形成1*(L*(M/2))的矩阵figure(1); % 绘制第1幅图subplot(311); % 窗口分割成3*1的，当前是第1个子图 plot(t,jidai,'LineWidth',2);% 绘制基带码元波形，线宽为2title('基带信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT,-1.1,1.1]) % 坐标范围限制subplot(312); % 窗口分割成3*1的，当前是第2个子图 plot(t2,I_lu,'LineWidth',2);% 绘制基带码元波形，线宽为2title('I路信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-1.1,1.1]) % 坐标范围限制subplot(313); % 窗口分割成3*1的，当前是第3个子图 plot(t2,Q_lu,'LineWidth',2);% 绘制基带码元波形，线宽为2title('Q路信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-1.1,1.1]) % 坐标范围限制%% QPSK调制fc=2000; % 载波频率2kHz zb1=cos(2*pi*fc*t2); % 载波1psk1=I_lu.*zb1; % PSK1的调制 zb2=sin(2*pi*fc*t2); % 载波2psk2=Q_lu.*zb2; % PSK2的调制 qpsk=psk1+psk2; % QPSK的实现figure(2); % 绘制第2幅图subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(t2,qpsk,'LineWidth',2);% 绘制基带码元波形，线宽为2title('QPSK信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT2,-2.1,2.1]) % 坐标范围限制%% 信号经过高斯白噪声信道tz=awgn(qpsk,20); % 信号qpsk中加入白噪声，信噪比为SNR=20dBsubplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(t2,tz,'LineWidth',2); % 绘制2ASK信号加入白噪声的波形axis([0,TotalT2,-2.5,2.5]); % 坐标范围设置title('通过高斯白噪声信道后的信号');% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 解调部分figure(3);tz1=tz.*zb1; % 相干解调，乘以相干载波subplot(211) % 窗口分割成2*1的，当前是第1个子图 plot(t2,tz1,'LineWidth',2) % 绘制I路乘以相干载波后的信号axis([0,TotalT2,-2.5,2.5]); % 设置坐标范围title("I路乘以相干载波后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签tz2=tz.*zb2; % 相干解调，乘以相干载波subplot(212) % 窗口分割成2*1的，当前是第2个子图 plot(t2,tz2,'LineWidth',2) % 绘制Q路乘以相干载波后的信号axis([0,TotalT2,-2.5,2.5]); % 设置坐标范围title("Q路乘以相干载波后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 加噪信号经过滤波器% 低通滤波器设计fp=2*Rb; % 低通滤波器截止频率，乘以2是因为下面要将模拟频率转换成数字频率wp=Rb/(Fs/2)b=fir1(30, fp/Fs, boxcar(31));% 生成fir滤波器系统函数中分子多项式的系数% fir1函数三个参数分别是阶数，数字截止频率，滤波器类型% 这里是生成了30阶(31个抽头系数)的矩形窗滤波器[h,w]=freqz(b, 1,512); % 生成fir滤波器的频率响应% freqz函数的三个参数分别是滤波器系统函数的分子多项式的系数，分母多项式的系数(fir滤波器分母系数为1)和采样点数(默认)512lvbo1=fftfilt(b,tz1); % 对信号进行滤波，tz1是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数lvbo2=fftfilt(b,tz2); % 对信号进行滤波，tz2是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数figure(4); % 绘制第4幅图 subplot(311); % 窗口分割成3*1的，当前是第1个子图 plot(w/pi*Fs/2,20*log(abs(h)),'LineWidth',2); % 绘制滤波器的幅频响应title('低通滤波器的频谱'); % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度/dB'); % y轴标签subplot(312) % 窗口分割成3*1的，当前是第2个子图 plot(t2,lvbo1,'LineWidth',2); % 绘制I路经过低通滤波器后的信号axis([0,TotalT2,-2.1,2.1]); % 设置坐标范围title("I路经过低通滤波器后的信号");% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签subplot(313) % 窗口分割成3*1的，当前是第3个子图 plot(t2,lvbo2,'LineWidth',2); % 绘制Q路经过低通滤波器后的信号axis([0,TotalT2,-2.1,2.1]); % 设置坐标范围title("Q路经过低通滤波器后的信号");% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% 抽样判决figure(5);k=0; % 设置抽样限值pdst1=1*(lvbo1>0); % 滤波后的向量的每个元素和0进行比较，大于0为1，否则为0pdst2=1*(lvbo2>0); % 滤波后的向量的每个元素和0进行比较，大于0为1，否则为0subplot(311) % 窗口分割成3*1的，当前是第1个子图 plot(t2,pdst1,'LineWidth',2)% 画出经过抽样判决后的信号axis([0,TotalT2,-0.1,1.1]); % 设置坐标范用title("I路经过抽样判决后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签subplot(312) % 窗口分割成3*1的，当前是第2个子图 plot(t2,pdst2,'LineWidth',2)% 画出经过抽样判决后的信号axis([0,TotalT2,-0.1,1.1]); % 设置坐标范用title("Q路经过抽样判决后的信号")% 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签%% I、Q路合并I_zong = [];Q_zong = [];% 取码元的中间位置上的值进行判决for j=L/2:L:(L*M/2) if pdst1(j)>0 I_zong=[I_zong,1]; else I_zong=[I_zong,-1]; endend% 取码元的中间位置上的值进行判决for k=L/2:L:(L*M/2) if pdst2(k)>0 Q_zong=[Q_zong,1]; else Q_zong=[Q_zong,-1]; endendcode = [];% 将I路码元为最终输出的奇数位置码元，将Q路码元为最终输出的偶数位置码元for n=1:M if mod(n, 2)~=0 code = [code, I_zong((n+1)/2)]; else code = [code, Q_zong(n/2)]; endendx4=code(fz,:); % 将原来code的第一行复制L次，称为L*M的矩阵dout=reshape(x4,1,L*M); % 产生单极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵subplot(313); % 窗口分割成3*1的，当前是第3个子图 plot(t,dout,'LineWidth',2);% 绘制基带码元波形，线宽为2title('I、Q路合并信号波形'); % 标题xlabel('时间/s'); % x轴标签ylabel('幅度'); % y轴标签axis([0,TotalT,-1.1,1.1]) % 坐标范围限制%% 绘制频谱%% 信源频谱figure(6) % 绘制第6幅图T=t(end); % 时间df=1/T; % 频谱分辨率N=length(jidai); % 采样长度f=(-N/2:N/2-1)*df; % 频率范围mf=fftshift(abs(fft(jidai)));%对信源信号采用快速傅里叶变换并移到矩阵中心subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f,mf,'LineWidth',2); % 绘制信源频谱波形title("基带信号频谱"); % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% QPSK信号频谱T2=t2(end); % 时间2df2=1/T2; % 频谱分辨率2N2=length(qpsk); % 采样长度2f2=(-N2/2:N2/2-1)*df2; % 频率范围sf=fftshift(abs(fft(qpsk)));% 对QPSK信号采用快速傅里叶变换并将0-fs频谱移动到-fs/2-fs/2subplot(212) % 窗口分割成2*1的，当前是第2个子图 plot(f2,sf,'LineWidth',2) % 绘制QPSK调制信号频谱title("QPSK信号频谱") % 标题xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf]) % 坐标范围限制%% I路乘以相干载波后的频谱mmf=fftshift(abs(fft(tz1)));% 对相干载波信号采用快速傅里叶变换并移到矩阵中心figure(7) % 绘制第7幅图subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f2,mmf,'LineWidth',2) % 画出乘以相干载波后的频谱title("I路乘以相干载波后的频谱")xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% Q路乘以相干载波后的频谱mmf2=fftshift(abs(fft(tz2)));% 对相干载波信号采用快速傅里叶变换并移到矩阵中心subplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(f2,mmf2,'LineWidth',2) % 画出乘以相干载波后的频谱title("Q路乘以相干载波后的频谱")xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签axis([-20000,20000,-inf,inf])% 坐标范围限制%% 经过低通滤波后的频谱figure(8);dmf=fftshift(abs(fft(lvbo1)));%对低通滤波信号采用快速傅里叶变换并移到矩阵中心subplot(211); % 窗口分割成2*1的，当前是第1个子图 plot(f2,dmf,'LineWidth',2) % 画出经过低通滤波后的频谱title("I路经过低通滤波后的频谱");xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签dmf2=fftshift(abs(fft(lvbo2)));%对低通滤波信号采用快速傅里叶变换并移到矩阵中心subplot(212); % 窗口分割成2*1的，当前是第2个子图 plot(f2,dmf2,'LineWidth',2) % 画出经过低通滤波后的频谱title("Q路经过低通滤波后的频谱");xlabel('频率/Hz'); % x轴标签ylabel('幅度'); % y轴标签figure(9); % 绘制第9幅图scatter(I,Q); % 绘制星座图，如果没有4个点，是因为码元数量不够的原因% 让坐标轴放在坐标原点ax = gca;ax.XAxisLocation = 'origin';ax.YAxisLocation = 'origin';
