st2tr.m

可以运行的m4k20功率分配图

function [st_out1,st_out2]=st2tr(input,fftlen)
% 空时编码程序
% 对输入序列按帧进行space-time 编码

input_len=length(input);                             %输入信号长度
frame_num=input_len/fftlen;                          %帧数

st_out1=zeros(1,input_len);
st_out2=zeros(1,input_len);

for n=0:(frame_num/2-1)                              %按帧进行space-time 编码
    st_out1(2*n*fftlen+1:(2*n+1)*fftlen)=input(2*n*fftlen+1:(2*n+1)*fftlen);
    st_out1((2*n+1)*fftlen+1:2*(n+1)*fftlen)=-conj(input((2*n+1)*fftlen+1:2*(n+1)*fftlen));
    st_out2(2*n*fftlen+1:(2*n+1)*fftlen)=input((2*n+1)*fftlen+1:2*(n+1)*fftlen);
    st_out2((2*n+1)*fftlen+1:2*(n+1)*fftlen)=conj(input(2*n*fftlen+1:(2*n+1)*fftlen));    
end

%st_out1=st_out1/sqrt(2);                             %每根天线上的功率都是原来的1/2
%st_out2=st_out2/sqrt(2);