%%STBC-MIMO-OFDM.mclearall;i=sqrt(-1);IFFT_bin_length=512;%傅立叶变换抽样点数目carrier_count=100;%子载波数目symbols_per_carrier=66;%符号数/载波cp_length=10;%循环前缀长度addprefix_length=IFFT_bin_length+cp_length;M_psk=4;bits_per_symbol=log2(M_psk);%位数/符号%O=[12;-2+j1+j];%AlamoutiScheme[x1x2;-x2*x1*]二天线发送矩阵O=[1-2-3;2+j1+j0;3+j01+j;0-3+j2+j];%for[x1-x2-x3;x2*x1*0;x3*0x1*;0-x3*x2*]三天线发送矩阵co_time=size(O,1);Nt=size(O,2);%发射天线数目Nr=2;%接收天线数目num_X=1;forcc_ro=1:co_timeforcc_co=1:Ntnum_X=max(num_X,abs(real(O(cc_ro,cc_co))));endendco_x=zeros(num_X,1);forcon_ro=1:co_time%计算delta,epsilon,etaandconjmatricesforcon_co=1:Nt%用于确定矩阵“O”中元素的位置，符号以及共轭情况ifabs(real(O(con_ro,con_co)))~=0delta(con_ro,abs(real(O(con_ro,con_co))))=sign(real(O(con_ro,con_co)));epsilon(con_ro,abs(real(O(con_ro,con_co))))=con_co;co_x(abs(real(O(con_ro,con_co))),1)=co_x(abs(real(O(con_ro,con_co))),1)+1;eta(abs(real(O(con_ro,con_co))),co_x(abs(real(O(con_ro,con_co))),1))=con_ro;coj_mt(con_ro,abs(real(O(con_ro,con_co))))=imag(O(con_ro,con_co));endendendeta=eta.';eta=sort(eta);eta=eta.';carriers=(1:carrier_count)+(floor(IFFT_bin_length/4)-floor(carrier_count/2));%坐标：(1to100)+14=(15:114)conjugate_carriers=IFFT_bin_length-carriers+2;%坐标：256-(15:114)+1=257-(15:114)=(242:143)tx_training_symbols=training_symbol(Nt,carrier_count);baseband_out_length=carrier_count*symbols_per_carrier;snr_min=3;%最小信噪比snr_max=15;%最大信噪比graph_inf_bit=zeros(snr_max-snr_min+1,2,Nr);%绘图信息存储矩阵graph_inf_sym=zeros(snr_max-snr_min+1,2,Nr);forSNR=snr_min:snr_max%开始仿真clcdisp('WaituntilSNR=');disp(snr_max);SNRn_err_sym=zeros(1,Nr);n_err_bit=zeros(1,Nr);Perr_sym=zeros(1,Nr);Perr_bit=zeros(1,Nr);re_met_sym_buf=zeros(carrier_count,symbols_per_carrier,Nr);re_met_bit=zeros(baseband_out_length,bits_per_symbol,Nr);baseband_out=round(rand(baseband_out_length,bits_per_symbol));%生成随机数用于仿真de_data=bi2de(baseband_out);%二进制向十进制转换data_buf=pskmod(de_data,M_psk,0);%PSK调制carrier_matrix=reshape(data_buf,carrier_count,symbols_per_carrier);fortt=1:Nt:symbols_per_carrier%取数为空时编码做准备，此处每次取每个子载波上连续的两个数data=[];forii=1:Nttx_buf_buf=carrier_matrix(:,tt+ii-1);data=[data;tx_buf_buf];endXX=zeros(co_time*carrier_count,Nt);forcon_r=1:co_time%进行空时编码forcon_c=1:Ntifabs(real(O(con_r,con_c)))~=0ifimag(O(con_r,con_c))==0XX((con_r-1)*carrier_count+1:con_r*carrier_count,con_c)=data((abs(real(O(con_r,con_c)))-1)*carrier_count+1:abs(real(O(con_r,con_c)))...*carrier_count,1)*sign(real(O(con_r,con_c)));elseXX((con_r-1)*carrier_count+1:con_r*carrier_count,con_c)=conj(data((abs(real(O(con_r,con_c)))-1)*carrier_count+1:abs(real(O(con_r,con_c)))...*carrier_count,1))*sign(real(O(con_r,con_c)));endendendend%空时编码结束XX=[tx_training_symbols;XX];%添加训练序列rx_buf=zeros(1,addprefix_length*(co_time+1),Nr);forrev=1:Nrforii=1:Nttx_buf=reshape(XX(:,ii),carrier_count,co_time+1);IFFT_tx_buf=zeros(IFFT_bin_length,co_time+1);IFFT_tx_buf(carriers,:)=tx_buf(1:carrier_count,:);IFFT_tx_buf(conjugate_carriers,:)=conj(tx_buf(1:carrier_count,:));time_matrix=ifft(IFFT_tx_buf);time_matrix=[time_matrix((IFFT_bin_length-cp_length+1):IFFT_bin_length,:);time_matrix];tx=time_matrix(:)';%------------------------------------------------------------------------%d=randint(1,4,[1,7]);%4多经信道模拟%...
