pilot.m

OFDM系统中时域和频域波束形成算法

function[pro]=Pilot(SNRindB,fd)
M=4;
N=64;
[s,d_pilot]=xinhao1;
U=xindao1(s,SNRindB,fd);
j=sqrt(-1);   %产生4*64的傅立叶变换矩阵
F4_64=zeros(4,N);
for i=1:4
for k=1:64
F4_64(i,k)=exp((-j)*2*pi*(k-1)*(i-1)*16/N);
end
end
w=[0,0,0,0]';
y1=zeros(1,N);
jieshou=zeros(1,N);
panjue=zeros(1,Q);
for l=1:M
V=U(1:4,((l-1)*N+1):(l*N));
y1=w'*V;
Y1=zeros(1,N);   %FFT变换
Y1=fft(y1,N);
for i=1:N      %判决
if(real(Y1(i))>=0)
jieshou(i)=1;
else
jieshou(i)=-1;
end
end
panjue(1,((l-1)*N+1):l*N)=jieshou;
pilot=zeros(4,1);          %提取解调后期望信号中的pilot
for i=1:4
pilot(i,1)=jieshou((i-1)*16+1);
end
d_pilot1=F4_64'*d_pilot(1:4,l);
pilot1=F4_64'*pilot;
w=w+2*buchang*V*(d_pilot1-pilot1);
end
%误码率计算同上