QPSKmatlab代码
- 格式：docx
- 大小：15.68 KB
- 文档页数：3
QPSKmatlab代码
QPSKCode = zeros(size(ConCode));
for i = 1 : Length
L = 2*i-1;
R = 2*i;
if ConCode(L:R) == [1 1]
QPSKCode(L:R) = [1 1];
elseif ConCode(L:R) == [0 1]
QPSKCode(L:R) = [-1 1];
elseif ConCode(L:R) == [0 0]
QPSKCode(L:R) = [-1 -1];
else
QPSKCode(L:R) = [1 -1];
end
end
%双极性不归零码
%时间单位 us 频率单位 MHz 码元速率 Mb/s
close all %关闭所有的窗口
clear all %清除所有的变量 k=14;
N=2^k; %总采样点数
L=256; %每码元采样点数
M=N/L; %码元数
Rb=2; %码元速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域分辨率dt(时域取样间隔)
df=1.0/(N*dt); %频域分辨率df，满足df*dt=1/N
T=N*dt; %计算时间宽度T
Bs=N*df/2; %系统带宽
t=linspace(-T/2,T/2,N); %产生时间取样点矩阵
f=linspace(-Bs,Bs,N); %产生频域取样点矩阵
EEp=zeros(size(f)); %定义全零距阵，将来作为求功率的累加器
for ii=1:20
n=1;%赋初值
while rem(n,2)~=0;%保证仿真时有偶数个1
nrz=round(rand(1,M));%随机产生单极性归零码
n=length(find(nrz==1));
end
po0=find(nrz==0);%找出nrz中0的位置
nrz(po0)=-1;%将0赋为-1
copy0=ones(1,L);%定义复制的次数L
temp0=nrz(copy0,:); %将nrz第一行复制L次，生成L*M的矩阵
dnrz=reshape(temp0,1,L*M);%将temp0重排成1*L*M的矩阵
Z=t2f(dnrz,dt);%傅式变换
PP=Z.*conj(Z)/T;%求功率谱
EEp=(EEp*(ii-1)+PP)/ii;%求平均功率谱
end
%开窗口1
figure(1)
set(1,'position',[10,50,300,200])%设置窗口位置和大小
plot(t,dnrz)
title('双极性不零码')
axis([0,10,-1.5,1.5])%设置坐标范围
figure(2)
set(2,'position',[350,50,300,200])%设置窗口位置和大小
bb=30+10*log10(EEp+eps);
plot(f,bb)
grid on%显示网格
axis([-10,+10,-50,50])%设置坐标范围
