说明：  for i=1 :Nt*Nr for l=1:1 h1=0; h2=0; for k=l:s-1 sita(k)=2*pi*rand; h1 = h1+sqrt(2)/sqrt(s-1/2)*sin(pi*k/(s-1))*cos(2*pi*fd*cos(pi*k/(2*s-1))*t+sita(k)); h2 = h2+sqrt(2)/sqrt(s-1/2)*cos(pi*k/(s-1))*cos(2*pi*fd*cos(pi*k/(2*s-1))*t+sita(k)); end sita(s)=rand; h1=h1+1/(sqrt(2)*sqrt(s-1/2))*cos(2*pi*fd*t+sita(s)); h2=h2+l/(sqrt(2)*sqrt(s-1/2))*cos(2*pi*fd*t+sita(s)); h(i,1)=h1 +j*h2; end end