xinhao1.m

来自「OFDM系统中时域和频域波束形成算法」· M 代码 · 共 37 行

function[s,d_pilot,preamble]=xinhao1
M=4;
N=64;
d1=zeros(1,4);
for i=1:4
temp=rand;
if(temp<0.5),
d1(i)=1;
else
d1(i)=-1;
end
end;
d_pilot=zeros(4,M);%提取希望用户的导频
for i=1:M
for j=1:4
d_pilot(j,i)=d1((j-1)*16+1+(i-1)*M);
end
end
%d2,d3,d4,d5产生同上
s1=zeros(1,P);%进行IFFT变换
a=zeros(1,N);
aa=zeros(1,N);
for i=1:M
for k=1:N,
a(k)=d1((i-1)*N+k);
end
aa=ifft(a,N);
for k=1:N,
s1((i-1)*N+k)=aa(k);
end
end
preamble=zeros(1,Q_preamble);%提取期望用户的训练序列
for i=1:Q_preamble
preamble(i)=s1(i);
end    %s2,s3,s4,s5产生同上
s=[s1;s2;s3;s4;s5];