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

1、多天线时变信道下基于卡尔曼滤波跟踪的数字图像通信系统仿真 140个时隙的信道跟踪值和实际的信道值2.接收到的图像效果（信噪比SNR=15dB） 源图像 接收图像3程序如下主程序：%图象信号产生picture=10*ones(64,64);for p=17:48 for q=17:48 picture(p,q)=0; endendsig_tx=picture(1:end);%(2,1,7)卷积编码bin_bits=hex2bi(sig_tx);tx_encondedbits=convencoder(bin_bits);%bit交织interleaver_depth=length(tx_encon

2、dedbits);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

3、);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); d

4、ata_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,le

5、ngth(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_es

6、t_result=h2_est_result h2_est(1); end h1_est(row+1)=a*h1_est(row); h1_est(row+2)=a2*h1_est(row); h2_est(row+1)=a*h2_est(row); h2_est(row+2)=a2*h2_est(row); %三次卡尔曼滤波进行信道跟踪 for m=1:3 H=h1_est(row+1) h2_est(row+1);conj(h2_est(row+2) conj(h1_est(row+2); if row+1<16 y=data_rx(row+1,volumn);conj(data_r

7、x(row+2,volumn); elseif volumn+1<=length(data_rx)/16 y=data_rx(row+1,volumn);conj(data_rx(1,volumn+1); else y=data_rx(row+1,volumn);0; end r=H'*y; s_dec2=r; s1_dec=s_dec2; s2_dec=-conj(s_dec2(2);conj(s_dec2(1); h_est1=y./s1_dec; h_est2=y./s2_dec; h1_est(row+1)=h_est1(1);h2_est(row+1)=h_est1(2

8、); h1_est(row+2)=h_est2(1);h2_est(row+2)=h_est2(2); end h1_est_result=h1_est_result h1_est(row+1); h1_est_result=h1_est_result h1_est(row+2); h2_est_result=h2_est_result h2_est(row+1); h2_est_result=h2_est_result h2_est(row+2); s_dec_result=s_dec_result conj(s1_dec'); end h1_est_result=h1_est_re

9、sult(1:end-1); h2_est_result=h2_est_result(1:end-1); s_dec_result=s_dec_result(1:end-1);end%16QAM解调rx_demapped=demapping(s_dec_result);%解交织bits_to_deinter=hex2bi(rx_demapped);rx_deinterleaveredbits=deinterleaver(bits_to_deinter,interleaver_depth);%卷积码解码rx_decodedbits=convdecoder(rx_deinterleaveredbi

10、ts);hex_data=bi2hex(rx_decodedbits);sig_rx=reshape(hex_data,64,64);%显示信道跟踪值与信道实际值figure(1)plot(real(h1_est_result(1:40),'k-');hold on;plot(real(h1(1:40);figure(2)plot(real(h1_est_result(1:40),'k-');hold on;plot(real(h2(1:40);figure(3)plot(imag(h1_est_result(1:40),'k-');hold on;plot(imag(h1(1:40);figure(4)plot(imag(h1_est_result(1:40),'k-');hold on;plot(imag(h2(1:40);%显示源图像与接收图像figure(5)image(picture)figure(6)image(sig_