使用MATLAB进行CSI数据处理

前言：Intel5300网卡获取CSI数据的教程强烈推荐该博文
安装抓取WIFI CSI的工具

• CSI幅值数据处理
  可以使用MATLAB自带函数abs实现。

% 读取自己的数据集
close all;clear;clc
load("datasets\test_1011\bow_ant1.mat","bow.mat")
bow_amp = abs(bow);

• CSI相位数据处理

感谢@weixin_44408931 的提醒，使用MATLAB中的unwrap函数需要注意自己的数据格式。若子载波是按行，可将unwrap(YourData,pi,2)；若按列，则设置为unwrap(YourData,pi,1)。具体请参考

https://ww2.mathworks.cn/help/matlab/ref/unwrap.html

function [CSI_phase_recovery] = GetPhase(R)
% phase
phase_matrix = angle(R);
no_of_subcarriers = 30;
m = [-28,-26,-24,-22,-20,-18,-16,-14,-12,-10,-8,-6,-4,-2,-1,1,3,5,7,9,11,13,15,17,19,21,23,25,27,28]';
CSI_phase_recovery = zeros(size(phase_matrix));
CSI_phase_temp = zeros(no_of_subcarriers,1);
CSI_phase_unwarp_temp = zeros(no_of_subcarriers,1);
for line_index = 1:(size(phase_matrix,1)/no_of_subcarriers)
    for row_index = 1:size(phase_matrix,2)
        CSI_phase_temp = phase_matrix(no_of_subcarriers*(line_index-1)+1:no_of_subcarriers*line_index,row_index);
        CSI_phase_unwarp_temp = unwrap(CSI_phase_temp,pi,1);
        b = mean(CSI_phase_unwarp_temp);
        k = (CSI_phase_unwarp_temp(no_of_subcarriers)-CSI_phase_unwarp_temp(1))/(m(no_of_subcarriers)-m(1));
        H = zeros(no_of_subcarriers,1);
        for i = 1:no_of_subcarriers
            H(i) = CSI_phase_unwarp_temp(i)-k*m(i)-b;
        end
        CSI_phase_recovery(no_of_subcarriers*(line_index-1)+1:no_of_subcarriers*line_index,row_index) = H;
    end
end

end

上述代码仅针对Wi-Fi 20MHz频宽，如果使用更高的频宽和更多的子载波数，可以参考

https://blog.csdn.net/a_beatiful_knife/article/details/119247331

Reference:

PhaseFi: Phase Fingerprinting for Indoor Localization with a Deep Learning Approach