该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
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);%计算发送信号,且未经过衰