qam.m

来自「基于卷积信号的MIMO系统盲信号估计」· M 代码 · 共 38 行

function q=QAM(N,L,seed1,seed2)

% ***********************************************
% FUNCTION q = QAM(N,L);
% -----------------
% Generates L-QAM complex sequences Xr(n)+jXi(n) 
% where Xr(n) and Xi(n) are two i.i.d. random 
% sequences independent of each other and each 
% takes the values +/- {1,3,5,...,(L-1)}
% ***********************************************

%clf;
rand('seed',seed1);
xr=rand(1,N); 
rand('seed',seed2);
xi=rand(1,N);
v=[0:1/L:1];
th=[-(L-1):2:L-1];
for i=1:L
  f=find(xr>=v(i) & xr<v(i+1));
  xr(f)=th(i)*ones(length(f),1);
  k=find(xi>=v(i) & xi<v(i+1));
  xi(k)=th(i)*ones(length(k),1);
end

x=[-L:0.5:L];
%plot([0 0],[-L L],'w');
%hold on;
%plot(x,zeros(1,length(x)),'w');
clear j;
%plot(xr+xi*j,'b*');
%axis([-L L -L L]);
q=xr+xi*j;