用简单仿真链路产生 WiFi CSI（不依赖专用工具箱，matlab实现）

Nt = 2;                  % 发射天线数
Nr = 2;                  % 接收天线数
Nfft = 64;               % FFT点数（802.11a/n 20MHz）
Ng   = 16;               % 循环前缀长度
subc_idx = [-26:-1 1:26];% 52个有效子载波（跳过DC）
Nsc = numel(subc_idx);
Es  = 1;                 % 训练符平均功率
SNRdB = 25;              % 接收端SNR（可改）

三、构造训练OFDM符

用BPSK全已知子载波（实际LTF有特定序列，这里只需"已知且非零"即可）

Matlab 复制代码

Xfull = zeros(Nfft,1);
trainSymbols = (2*(randi([0 1],Nsc,1))-1); % BPSK +/-1
trainSymbols = sqrt(Es)*trainSymbols;      % 调整功率
Xfull(subc_idx + (subc_idx<0)*Nfft + 1) = trainSymbols;  % 映射到FFT bins
x_time = ifft(ifftshift(Xfull), Nfft);     % IFFT（使用ifftshift对齐DC）
tx_ofdm = [x_time(end-Ng+1:end); x_time];  % 加CP
Ns = length(tx_ofdm);

四、多径瑞利信道

瑞利衰落信道（Rayleigh fading channel）是一种无线电信号传播环境的统计模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，即"衰落"，并且其包络服从瑞利分布。这一信道模型能够描述由电离层和对流层反射的短波信道，以及建筑物密集的城市环境。

Matlab 复制代码

% 设定一个离散PDP
tapDelays = [0 1 3 5 7];          % 以采样点为单位的延迟
tapPowers_dB = [0 -3 -6 -9 -12];  % 功率谱
tapPowers = 10.^(tapPowers_dB/10);
L = tapDelays(end)+1;

% 为每个链路生成独立脉冲响应
h = cell(Nr,Nt);
for r = 1:Nr
    for t = 1:Nt
        h_rt = zeros(L,1);
        for l = 1:numel(tapDelays)
            amp = sqrt(tapPowers(l)/2)*(randn+1j*randn);
            h_rt(tapDelays(l)+1) = amp;
        end
        h{r,t} = h_rt;
    end
end

五、逐天线时分发射并接收

Matlab 复制代码

Y = zeros(Nr,Nt,Nfft);  % 保存每次接收的频域符
noiseVarLin = Es/Nsc * 10^(-SNRdB/10);

for t = 1:Nt
    % 当前仅开启第t根发射天线，其它静默
    txSig_t = tx_ofdm;
    % 通过信道到各接收天线并叠加噪声
    for r = 1:Nr
        yr = conv(txSig_t, h{r,t});
        yr = yr(1:Ns); % 对齐长度（忽略尾符泄露，训练符足够长可避免IS I）
        % 加噪
        yr = yr + sqrt(noiseVarLin/2)*(randn(size(yr))+1j*randn(size(yr)));
        % 去CP、FFT
        y_noCP = yr(Ng+1:Ng+Nfft);
        Yfreq = fftshift(fft(y_noCP, Nfft));
        Y(r,t,:) = Yfreq; % 保存
    end
end

六、CSI信道估计

Matlab 复制代码

Hhat = zeros(Nr,Nt,Nsc); % 只保留有效子载波
Xused = Xfull(subc_idx + (subc_idx<0)*Nfft + 1); % 已知训练符频域值

for t = 1:Nt
    for r = 1:Nr
        Yrt = squeeze(Y(r,t,:));              % Nfft x 1
        Yused = Yrt(subc_idx + (subc_idx<0)*Nfft + 1);
        Hhat(r,t,:) = Yused ./ Xused;         % LS估计
    end
end

七、幅值和相位可视化

使用abs函数获取幅值，使用angle函数获取相位。

Matlab 复制代码

k = 1:Nsc;
H11 = squeeze(Hhat(1,1,:));
figure; plot(k, 20*log10(abs(H11))); grid on; xlabel('子载波索引'); ylabel('幅度(dB)');
title('估计CSI幅度（r=1,t=1）');
figure; plot(k, angle(H11)); grid on; xlabel('子载波索引'); ylabel('相位(rad)');
title('估计CSI相位（r=1,t=1）');

八、模拟与真实对比

模拟产生的CSI幅值：

真实产生的CSI幅值：

模拟产生的CSI相位

真实产生的CSI相位：

博主的每篇博文都是用心去写的，喜欢的可以多多支持和收藏，创作不易，未经作者允许，请勿转载或者抄袭。看WiFi指纹定位，我只认准全村第二帅和数产小黑娃，如外站有抄袭，记得点点举报！

关于WiFi 指纹定位迷茫的童鞋可以参考博主以前文章：WiFi指纹定位发论文求生指南------浅论WiFi指纹定位的前世今生与毕业困境-CSDN博客