该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

clc;

clear;

Nfft=16;

L=3;

N=11;

hl=[0.2944+1.6236j,-1.3362-0.6918j,0.7143+0.8580j];

fl=linspace(0.02,0.06,L);

SNR=linspace(-10,40,N);

mse_hl_1_dB=zeros(1,N);

mse_fl_1_dB=zeros(1,N);

mse_hl_2_dB=zeros(1,N);

mse_fl_2_dB=zeros(1,N);

nl=0:L-1;

delta_h=0.01;

delta_f=0.01;

repeat=10;

hl_initial=(1+j)*ones(2,L);

fl_initial=0.05*ones(2,L);

for i=1:N

hl_temp1=zeros(1,N);

fl_temp1=zeros(1,N);

hl_temp2=zeros(1,N);

fl_temp2=zeros(1,N);

for t=1:repeat

s=round(rand(1,2*Nfft));%产生随机的输入

X=QPSK_modulation(s,Nfft);

X=X/sqrt(X*X');

y=calcu_org(hl,fl,nl,X,Nfft,L);%计算发送信号，且未经过衰落和噪声干扰

Y=add_noise(y,Nfft,SNR(i));

[hl_1,fl_1]=channel_est_p1(X,nl,Y,Nfft,L,hl_initial,fl_initial,delta_h,delta_f);%信道估计算法1

[hl_2,fl_2]=channel_est_p2(X,nl,Y,Nfft,L,hl_initial,fl_initial,delta_h,delta_f);%信道估计算法2

hl_temp1(i)=hl_temp1(i)+mean(abs(hl_1-hl).^2);

fl_temp1(i)=fl_temp1(i)+mean(abs(fl_1-fl).^2);

hl_temp2(i)=hl_temp2(i)+mean(abs(hl_2-hl).^2);

fl_temp2(i)=fl_temp2(i)+mean(abs(fl_2-fl).^2);

end

hl_temp1_mse(i)=hl_temp1(i)/repeat;

fl_temp1_mse(i)=fl_temp1(i)/repeat;

hl_temp2_mse(i)=hl_temp2(i)/repeat;

fl_temp2_mse(i)=fl_temp2(i)/repeat;

end

mse_hl_1_dB=10*log10(hl_temp1_mse);

mse_fl_1_dB=10*log10(fl_temp1_mse);

mse_hl_2_dB=10*log10(hl_temp2_mse);

mse_fl_2_dB=10*log10(fl_temp2_mse);

figure;

plot(SNR,mse_hl_1_dB,'-ok',SNR,mse_hl_2_dB,'-*m');

xlabel('平均信噪比(dB)');

ylabel('hl估计值与理论值的均方误差(dB)');

legend('算法1','算法2');

grid on;

figure;

plot(SNR,mse_fl_1_dB,'-ok',SNR,mse_fl_2_dB,'-*m');

xlabel('平均信噪比(dB)');

ylabel('fl估计值与理论值的均方误差(dB)');

legend('OFDM LS','OFDM LMMSE');

grid on;

这是OFDM信道估计的主程序，哪位大神能告诉我这个程序里面用了哪些参数和做了什么中间做了什么处理么，急