多天线时变信道下基于卡尔曼滤波跟踪的数字图像通信系.docVIP

多天线时变信道下基于卡尔曼滤波跟踪的数字图像通信系统仿真 1．40个时隙的信道跟踪值和实际的信道值 2.接收到的图像效果（信噪比SNR=15dB） 源图像 接收图像 3．程序如下 = 1 \* GB3 ①主程序： %图象信号产生 picture=10*ones(64,64); for p=17:48 for q=17:48 picture(p,q)=0; end end sig_tx=picture(1:end); %(2,1,7)卷积编码 bin_bits=hex2bi(sig_tx); tx_encondedbits=convencoder(bin_bits); %bit交织 interleaver_depth=length(tx_encondedbits); tx_interleaveredbits=interleaver(tx_encondedbits,interleaver_depth); data_to_map=bi2hex(tx_interleaveredbits); %16QAM调制 tx_mappeddata=mapping(data_to_map); data_tx=tx_mappeddata; len=length(data_tx); %初始化信道 h1=[];h2=[]; v=sqrt(0.002/2)*(randn(1,len)+rand(1,len)*i); a=besselj(0,2*pi*0.025)*exp(j*2*pi*0.2); h1(1)=-0.15-0.85*i+(1-a)*v(1); h2(1)=-0.90+0.95*i+(1-a)*v(1); for order=2:len; h1(order)=a*h1(order-1)+(1-a)*v(order); h2(order)=a*h2(order-1)+(1-a)*v(order); end %计算信噪比 sig_pow=std(data_tx); noise_pow=std(v); SNR=10*log10(sig_pow/noise_pow) %空时编码信号通过时变信道 data_rx=[]; for t1=0:length(data_tx)-1 n1=fix(t1/2); data_rx(2*n1+1)=h1(2*n1+1)*data_tx(2*n1+1)+h2(2*n1+1)*data_tx(2*(n1+1))+v(2*n1+1); data_rx(2*(n1+1))=-h1(2*(n1+1))*conj(data_tx(2*(n1+1)))+h2(2*(n1+1))*conj(data_tx(2*n1+1))+v(2*(n1+1)); end %-------------------------------------- %信号接受、信道跟踪、信号恢复 %-------------------------------------- data_rx=reshape(data_rx,16,length(data_rx)/16); h1=reshape(h1,16,length(h1)/16); h2=reshape(h2,16,length(h2)/16); h1_est=[]; h2_est=[]; h1_est_result=[]; h2_est_result=[]; s_dec_result=[]; for volumn=1:size(data_rx,2) for row=1:2:16 %信道训练 if row==1 s_dec1=data_tx((volumn-1)*16+1); s_dec_result=[s_dec_result s_dec1]; h1_est(1)=h1(row,volumn); h2_est(1)=h2(row,volumn); h1_est_result=[h1_est_result h1_est(1)]; h2_est_result=[h2_est_result h2_est(1)]; end h1_est(row+1)=a*h1_est(row); h1_est(row+2)=a^2*h1_est(row); h2_est(row+1)=a*h2_est(row);