資源描述:
《任意發(fā)送天線,任意接收天線stbc-mimo-ofdmmatlab程序》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、任意發(fā)送天線,任意接收天線的STBC-MIMO-OFDM的MATLAB程序%%STBC-MIMO-OFDM.mclearall;i=sqrt(-1);IFFT_bin_length=512;%傅立葉變換抽樣點數(shù)目carrier_count=100;%子載波數(shù)目symbols_per_carrier=66;%符號數(shù)/載波cp_length=10;%循環(huán)前綴長度addprefix_length=IFFT_bin_length+cp_length;M_psk=4;bits_per_symbol=log2(M
2、_psk);%位數(shù)/符號%O=[12;-2+j1+j];%AlamoutiScheme[x1x2;-x2*x1*]二天線發(fā)送矩陣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*]三天線發(fā)送矩陣co_time=size(O,1);Nt=size(O,2);%發(fā)射天線數(shù)目Nr=2;%接收天線數(shù)目num_X=1;forcc_ro=1:co_timeforcc_co=1:Ntnum_X=max(num_X,abs
3、(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
4、(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,co
5、n_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_sym
6、bols=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%開始仿真clc
7、disp('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
8、=round(rand(baseband_out_length,bits_per_symbol));%生成隨機數(shù)用于仿真de_data=bi2de(baseband_out);%二進制向十進制轉(zhuǎn)換data_buf=pskmod(de_data,M_psk,0);%PSK調(diào)制carrier_matrix=reshape(data_buf,carrier_count,symbols_per_carrier);fortt=1:Nt:symbols_per_ca