zuixiaoercheng.rar_最小二乘算法

%**************************MIMD最小二乘估计******************************* %**********hg估计信道 clc; clear all; %1.产生一个【16，3】矩阵,并且产生3个用户8psk随机信号 for i=1:3 Xi=randint(16,1,8); Xi=j*sin(pi*Xi/4)+cos(pi*Xi/4); X(:,i)=Xi; % 组成16行3列的数组 end %2.在接收端0产生叠加信号 h=[1 0.5 0.2]'; y=X*h; %3.产生白噪声 y1=zeros(1,13); y2=zeros(1,13); y3=zeros(1,13); M=5000; for zz=1:M for snr=0:12 Y=awgn(y,snr); %加高斯白噪声 %4.求出估计信道值 Xw=inv(X'*X)*X'; %X伪逆 hg=Xw*Y; %信道估计 %5.求均方误差 MSE=abs(hg-h).^2; y11(snr+1)=MSE(1,:); %h00 y22(snr+1)=MSE(2,:); %h10 y33(snr+1)=MSE(3,:); %h20 end for k=1:13 y1(k)=y1(k)+y11(k); %将h00信道的13个值依次放入数组中 y2(k)=y2(k)+y22(k); y3(k)=y3(k)+y33(k); end end y1=y1/M; y2=y2/M; y3=y3/M; ym=(y1+y2+y3)/3; n1=0; n2=12; nn=[n1:n2]; %x坐标 subplot(1,3,1) plot(nn,ym) grid xlabel('snr'); ylabel('MSE');