【導(dǎo)讀】統(tǒng)的性能提出了更高的要求。選擇性衰落等方面具備優(yōu)越的性能，是未來通信領(lǐng)域中的關(guān)鍵技術(shù)?；贛ATLAB軟件對所建立的MIMO系統(tǒng)的信道容量進(jìn)。MIMO-OFDM系統(tǒng)進(jìn)行了完善并采用MATLAB對其性能進(jìn)行了仿真，結(jié)果顯示，

　　

【正文】 )。 %Calculate the recieved data with AWGN Rx=sqrt(64)*HD+sqrt(N0(n)/2)*NOISE。 %weights calculation % MMSE 均衡， Q*R % 均衡矩陣 Q=inv(H*H39。+N0*I)*H，這里先求出 inv(H*H39。+ N0*I) for i=1:N_sym Hfi=H(:,i)。 Ri=Hfi*Hfi39。+Rnn。 % H*H39。+N0*I invRi=inv(Ri)。 invP1((i1)*L+1:i*L,(i1)*L+1:i*L)=invRi。 end。 南京林業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)（論文） 26 %estimate the transmitted signal estmodDt=mod_d。 for k=1:K Hfk=H((k1)*L+1:k*L,:)。 Wk=invP1*reshape(Hfk,L*N_sym,1)。 Wk=reshape(Wk,L,N_sym)。 %multiply received signal with weights fdefk=sum(conj(Wk).*Rx((k1)*L+1:k*L,:), 1)。 fdetk=fdefk。 fdet(k,:)=fdetk。 end %get the estimated modulated data %解調(diào) for i=1:N_sym ei=transpose(fdet(:,i))。 estmodDti=(sign(real(ei))+1i*sign(imag(ei)))/sqrt(2*Es)。 estmodDt(i,:)=estmodDti。 end。 %calculate the average BER error=transpose(reshape(round((mod_destmodDt)/sqrt(2*Es)),N_sym,K))。 BERk(:,n)=BERk(:,n)+sum(abs(error).^2,2)。 end。 end。 %calculate the average BER BERk=BERk./(ones(K,1)*count*N_sym*log2M)。 BER=mean(BERk,1)。 %plot the figure figure semilogy(SNR_dB,BER,39。bo39。)。 hold on。 grid on。 legend(39。OFDM MMSE 1x139。)。 ylabel(39。Average BER39。)。 xlabel(39。SNR(dB)39。)。 任意收發(fā)天線下的 基于 STBC 的 MIMOOFDM 系統(tǒng) clear all。 i=sqrt(1)。 IFFT_bin_length=512。 南京林業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)（論文） 27 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=[1 2。2+j 1+j]。 %Alamouti Scheme [x1 x2。x2* x1*] O=[1 2 3。2+j 1+j 0。3+j 0 1+j。0 3+j 2+j]。 %for [x1 x2 x3。x2* x1* 0。x3* 0 x1*。0 x3* x2*] co_time=size(O,1)。 Nt=size(O,2)。 Nr=2。 num_X=1。 for cc_ro=1:co_time for cc_co=1:Nt num_X=max(num_X,abs(real(O(cc_ro,cc_co))))。 end end co_x=zeros(num_X,1)。 for con_ro=1:co_time for con_co=1:Nt if abs(real(O(con_ro,con_co)))~=0 delta(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))。 end end end eta=eta.39。 eta=sort(eta)。 eta=eta.39。 南京林業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)（論文） 28 carriers = (1: carrier_count) + (floor(IFFT_bin_length/4) floor(carrier_count/2))。 conjugate_carriers=IFFT_bin_lengthcarriers+2。 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_maxsnr_min+1,2,Nr)。 graph_inf_sym=zeros(snr_maxsnr_min+1,2,Nr)。 for SNR=snr_min:snr_max clc disp(39。Wait until SNR=39。)。disp(snr_max)。 SNR n_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)。 carrier_matrix=reshape(data_buf,carrier_count,symbols_per_carrier)。 for tt=1:Nt:symbols_per_carrier data=[]。 for ii=1:Nt tx_buf_buf=carrier_matrix(:,tt+ii1)。 data=[data。tx_buf_buf]。 end XX=zeros(co_time*carrier_count,Nt)。 for con_r=1:co_time for con_c=1:Nt if abs(real(O(con_r,con_c)))~=0 if imag(O(con_r,con_c))==0