大规模MIMO下LS信道估计算法程序

%{ 本程序仿真多小区单用户情况下的LS信道估计算法 %} close all; clc; clear all; %% 一些参数设置 Cells=3; %小区数 Users=1; %每小区的用户数 Base_antenna=64; %基站天线数 Modutle_order=4; %qam调制阶数（4qam==qpsk） Pilot_length=16; %导频长度 Pr=1; %用户的发射功率 Large_scale_fading=[1 0 0; 0.8 0.1 0.1; 0.5 0.25 0.25]; %大尺度衰落系数beta,第一列为目标小区的 %Large_scale_fading=[0.8 0.1 0.1]; %大尺度衰落系数beta,第一列为目标小区的 Type=size(Large_scale_fading,1); SNR=-5:5:30; %信噪比的范围 Length_SNR=length(SNR); %信噪比的长度 Itor_number=400; %每一个信噪比下的循环迭代次数 MSE_LS=zeros(Type,Length_SNR); %LS信道估计误差 for itor_type=1:Type for itor_SNR=1:Length_SNR snr=SNR(itor_SNR); % 此时的信噪比 % 显示仿真进度 fprintf('共 %d 组,程序运行第 %d 组下信噪比为 %d \n',Type,itor_type,snr); Mse=0; for itor_Itor_number=1:Itor_number %产生上行信道 [H_F_up]=Gena_channal(Cells,Users,Base_antenna,Large_scale_fading(itor_type,:)); % 发送过程 产生导频数据，调制 [pilot]=Transmit(Users,Pilot_length,Modutle_order); % 过信道 所有用户的数据都过目标用户的信道 [Data_after_pass_chan]=pass_channel(Cells,Users,H_F_up,pilot,Base_antenna,Pilot_length); %加噪声 [Nosize_signal]=Add_Nosize(Data_after_pass_chan,snr); % LS信道估计 只估计小区1用户1的信道 [H_LS]=Channel_estimate(Nosize_signal,pilot(1,:),Pilot_length); % 计算信道估计误差 [H_LS_MSE]=Channel_Estimate_Mse(H_F_up{1,1},H_LS); Mse=Mse+H_LS_MSE; end MSE_LS(itor_type,itor_SNR)=Mse/Itor_number; end end MSE_LS figure(1); semilogy(SNR,MSE_LS(1,:),'-ks');hold on; semilogy(SNR,MSE_LS(2,:),'-k<');hold on; semilogy(SNR,MSE_LS(3,:),'-k*');hold on; legend('LS，大尺度衰落系数1,0,0','LS，大尺度衰落系数0.8,0.1,0.1','LS，大尺度衰落系数0.5,0.25,0.25'); title('LS信道估计的性能仿真'); xlabel('信噪比/SNR'); ylabel('均方误差/MSE'); set(gca, 'XTick', [-5 0 5 10 15 20 25 30]); grid on;