LS.zip_LS 均衡_LS均衡_OFDM均衡_OFDM的均衡

clear all;

close all;

fprintf('OFDM 信道估计仿真');

carrier_count=64;%-----------载波数目

num_symbol=50;%--------------OFDM 符号个数

Guard=8;%--------------------循环前缀

pilot_Inter=8;%--------------导频间隔

modulation_mode=16;%---------调制方式

SNR=[0:2:20];%-------------信噪比取值

MSE2=zeros(length(SNR),NumLoop);

%%%%%%%%%%%%%%%%%%%%%%% 主 程 序 循

环%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for c1=1:length(SNR)

fprintf(' 仿真信噪比为%f',SNR(c1));

for num1=1:NumLoop

%---------------产生发送的随机序列——————————————-

BitsLen=carrier_count*num_symbol;

BitsTx=randint(1,BitsLen);

%---------------符号调制---------------------------------------

pilot_symbols(1,i)=pilot_symbols(1,i)-1;

else

pilot_symbols(1,i)=pilot_symbols(1,i);

end

end

pilot_symbols=pilot_symbols';

%----------------计算导频和数据数目----------------------------

num_pilot=ceil(num_symbol/pilot_Inter);

if rem(num_symbol,pilot_Inter)==0

num_pilot=num_pilot+1;

pilot_Indx=zeros(1,num_pilot);

Data_Indx=zeros(1,num_pilot*(pilot_Inter+1));

pilot_Indx(1,num_pilot)=num_data;

for j=0:num_pilot

Data_Indx(1,(1+j*pilot_Inter):(j+1)*pilot_Inter)=(2+j*(pilot_Inter+1)):((j+1)*(pilot_Inter+1));

end

%----------------导频插入-------------------------------------

piloted_ofdm_syms(:,pilot_Indx)=repmat(pilot_symbols,1,num_pilot);

%----------------IFFT 变换———————————————————

time_signal=sqrt(carrier_count)*ifft(piloted_ofdm_syms);

%----------------加循环前缀------------------------------------

add_cyclic_signal=[time_signal((carrier_count-Guard+1:carrier_count),:);time_signal];

Tx_data_trans=reshape(add_cyclic_signal,1,(carrier_count+Guard)*num_data);

%----------------信道处理--------------------------------------

Tx_signal_power=sum(abs(Tx_data_trans(:)).^2)/length(Tx_data_trans(:));

noise_var=Tx_signal_power/(10^(SNR(c1)/10));

Rx_data=awgn(Tx_data_trans,SNR(c1),'measured');

Rx_pilot=Rx_carriers(:,pilot_Indx);

Rx_fre_data=Rx_carriers(:,Data_Indx);

pilot_patt=repmat(pilot_symbols,1,num_pilot);

pilot_esti=Rx_pilot./pilot_patt;

%----------------LS 估计的线性插值———————————————

int_len=pilot_Indx;

len=1:num_data;

for ii=1:carrier_count