信道估计( channel estimation）

clear; clc; y=abs('abcdefghijklmn');%将字符串转换成ASCII码 y1=dec2bin(y,8); %将ASCII码转换成8位二进制数，‘char’ % y2=y1-48; %将二进制代码转换为矩阵，即将字符串转为数据,'double'，或者用str2num y2=y1+0-'0'; %char to double,char,其实也是ascII码，就给char数字加0，就将char变为double矩阵了，但是ascII码值，再减去‘0’就行了 x1=reshape(y2,1,[]);%将y2转为一维矩阵，上述两步为重中之重，只有转换为矩阵才能进行下面的运算 a=length(x1); b=a/2; c=ceil(b); if c~=b error('出错，请重新输入'); end L=3; t=0:20; for n=1:2:a if x1(n)==0 && x1(n+1)==0; d((n+1)/2)=-1/sqrt(2); d1((n+1)/2)=-1/sqrt(2); elseif x1(n)==0 && x1(n+1)==1; d((n+1)/2)=-1/sqrt(2); d1((n+1)/2)=1/sqrt(2); elseif x1(n)==1 && x1(n+1)==0; d((n+1)/2)=1/sqrt(2); d1((n+1)/2)=-1/sqrt(2); elseif x1(n)==1 && x1(n+1)==1; d((n+1)/2)=1/sqrt(2); d1((n+1)/2)=1/sqrt(2); end end s=(d+d1*1i).'; for SNR=0:20 sigma=L/exp(SNR/10); for u=1:100 % h1=normrnd(0,1,1,L); % h2=normrnd(0,1,1,L); % h=((h1+h2*1i)/sqrt(2)).'; h=wgn(L,1,0,'complex'); x=comsystem(s,L,h,SNR); h_estmvu=mvu(x,s,L); h_estmmse=mmse(x,s,L,sigma); var_mvu1=abs(h-h_estmvu).^2/3; var_mmse1=abs(h-h_estmmse).^2/3; var_mvu2(u)=sum(var_mvu1); var_mmse2(u)=sum(var_mmse1); end errormvu=sum(var_mvu2); errormmse=sum(var_mmse2); var_minmvu(SNR+1) = errormvu/100; var_minmmse(SNR+1)=errormmse/100; end semilogy(t,var_minmvu,'r-o'); %用圆圈标志点位置 hold on semilogy(t,var_minmmse,'s-'); xlabel('SNRdB'); ylabel('误差'); legend('MVU估计','MMSE估计');