CSI信道状态信息处理入门1（附源代码）

一、1组CSI数据的SNR

csi_trace = read_bf_file('源数据');%源数据为你用CSI Tools工具采集到的原始数据

for l=1:1 %取50个数据包的数据
csia=get_scaled_csi(csi_trace{l});
for i=1:1 %1个发射天线
for j=1:2 %2个接收天线
for k=1:30 %30个子载波数据
B(i,j,k)=csia(i,j,k);
end
end
end
plot(db(abs(squeeze(B).')))
hold on
end

legend('RX Antenna A', 'RX Antenna B', 'Location', 'SouthEast' );
xlabel('Subcarrier index');
ylabel('SNR [dB]')
hold off;

二、时域柱状图

实际情况时域柱状图
csi_trace = read_bf_file('源数据');
for l=1:50 %取50个数据包的数据
csia=get_scaled_csi(csi_trace{l});
for k=1:30 %30个子载波数据
B(1,1,k)=csia(1,1,k);
end
csi_one=squeeze(B).';

for ki=1:30 %30个子载波数据
csi_amp(ki,l)=csi_one(ki); %第l个数据包的第k个子载波的相位值 【子载波，数据包】
end
end

csi_ifft=ifft(csi_amp(:,1));

T_amp=abs(csi_ifft);

bar(T_amp);

set(gca,'XTick',[0 10 20 30]);

set(gca,'xticklabel',{'0','0.5','1','1.5'});

xlabel('Delay (ms)');

ylabel('Amplitude(dB)');

三、未经过解卷绕的相位图

csi_trace = read_bf_file('源数据');
for l=1:50 %取50个数据包的数据
csia=get_scaled_csi(csi_trace{l});
for k=1:30 %30个子载波数据
B(1,1,k)=csia(1,1,k);
csi_one=squeeze(B).';
csi_phase=angle(csi_one);%angle求复数矩阵相位角的弧度值，取值-pi到pi
csi_amplitude=abs(csi_one);%angle求复数矩阵的绝对值
phase(k,l)=csi_phase(k); %第l个数据包的第k个子载波的相位值 【子载波，数据包】
amplitude(k,l)=csi_amplitude(k);
end
end
plot(180*phase/pi,'.','Markersize',25);

xlabel('Subcarrier index');
ylabel('Phase (deg)')

四、经过解卷绕后的相位图

csi_trace = read_bf_file('源数据');
for l=1:50 %取50个数据包的数据
csia=get_scaled_csi(csi_trace{l});
for k=1:30 %30个子载波数据
B(1,1,k)=csia(1,1,k);
csi_one=squeeze(B).';
csi_phase=angle(csi_one);%angle求复数矩阵相位角的弧度值，取值-pi到pi
csi_amplitude=abs(csi_one);%angle求复数矩阵的绝对值
phase(k,l)=csi_phase(k); %第l个数据包的第k个子载波的相位值 【子载波，数据包】
amplitude(k,l)=csi_amplitude(k);
end
plot((180*unwrap(phase(:,l))/pi),'b');
%plot(unwrap(phase(:,l)),'b');
hold on;
end

xlabel('Subcarrier index');
ylabel('Phase (deg)')
%ylabel('Phase (rad)')