資源簡介
function y=pqmod(M,N,flag)
%本程序完成基帶信號星座圖映射
%參數(shù)說明
% y 生成的基帶調(diào)制行向量
% M 調(diào)制階數(shù)
% N 碼元個數(shù)
% flag flag=1 PSK調(diào)制
% flag=2 QAM調(diào)制
% flag=3 OQPSK調(diào)制
%函數(shù)體
x=randint(N,1,M);
if flag==1
y=pskmod(x,M);
% scatterplot(y);
elseif flag==2
y=qammod(x,M);
% scatterplot(y)
elseif flag==3
M=4;
y=oqpskmod(x);
y=y(2:end-1);
% scatterplot(y)
end;
%本程序完成基帶信號星座圖映射
%參數(shù)說明
% y 生成的基帶調(diào)制行向量
% M 調(diào)制階數(shù)
% N 碼元個數(shù)
% flag flag=1 PSK調(diào)制
% flag=2 QAM調(diào)制
% flag=3 OQPSK調(diào)制
%函數(shù)體
x=randint(N,1,M);
if flag==1
y=pskmod(x,M);
% scatterplot(y);
elseif flag==2
y=qammod(x,M);
% scatterplot(y)
elseif flag==3
M=4;
y=oqpskmod(x);
y=y(2:end-1);
% scatterplot(y)
end;

代碼片段和文件信息
%QPSK信號的四階循環(huán)累量特征?
%生成QPSK信號?
clear?all;?
clc?
M=16;?
fs=7200;?
fc=1800;?
R=600;?
N=200;?
x=pqmod(MN2)‘;?
scatterplot(x);?
y=MPSKQAMmod(xNfsRfc);?
y=awgn(y15‘measured‘);?
%?y=real(y);?
yft=abs(fft(y4096));?
figure(1)?
plot(fs*(0:2047)/4096yft(1:2048))?
title(‘頻譜‘);?
%計算四階矩?
alpha=(fc-R/2):0.5:(fc+R/2);?
S4=y.^4;?
S2=y.^2;?
m=size(y);?
k=(1:m(2))‘;?
n=0;?
for?i=1:(2*R+1)?
????e1(:i)=exp(-j*8*pi*alpha(i)*k/fs);?
????su1(:i)=S4‘.*e1(:i);?
????C1(i)=mean(su1(:i));?
????e2(:i)=exp(-j*4*pi*alpha(i)*k/fs);?
????su2(:i)=S2‘.*e2(:i);?
????C2(i)=mean(su2(:i));?
????n=n+1?
end?
C=C1-3*(C2.^2);?
figure(3)?
plot(alphaabs(C));?
title(‘8PSKSNR=20fc=32000‘)
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????566??2011-07-13?21:36??MPSK信號基于高階循環(huán)譜估計載波頻率matlab的代碼\MPSKQAMmod.m
?????文件????????730??2011-07-13?21:37??MPSK信號基于高階循環(huán)譜估計載波頻率matlab的代碼\cylic_cumulate.m
?????文件????????485??2011-07-13?21:37??MPSK信號基于高階循環(huán)譜估計載波頻率matlab的代碼\pqmod.m
?????目錄??????????0??2011-07-13?21:36??MPSK信號基于高階循環(huán)譜估計載波頻率matlab的代碼
-----------?---------??----------?-----??----
?????????????????1781????????????????????4
評論
共有 條評論