矢量网络分析仪（VNA）测量S参数，利用MATLAB处理Active S11

矢量网络分析仪（VNA）测量S参数，利用MATLAB处理Active S11

测试S参数的注意点：
1.要将所有的连接器都进行校准，最好的状态是校准完之后直接装天线，可以减小甚至消除因为连接器带来的相位误差；
2.每测量一次要进行一次校准，为什么测量后有偏差？
3.开始校准时要进行频率的选择和步进的确定；
4.特别分批次测量时，如果选择步进不一致，数据处理有麻烦；
5.线长对应相位的相位周期数计算公式，应该是的吧
用matlab处理：
1.当储存的数据分别在两个CSV后缀文件里，想把他们合并在同一文件里面。例如实部虚部被分开储存了，例如幅度相位被分开了。
处理想法是，将第二文件的第二列加到第一个文件的第三列，或者将第一个文件的第三列替换成第二个文件的第二列；
2.要处理幅度和相位改成实部和虚部。直接带入公式很简单，注意角度输入要转化成弧度带入；下面这段文字主要是针对幅度相位的，顺带求出了实部虚部，并且以多种形式储存了。

clear all;
close all;

X1=importdata('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_MAG.csv');
X2=importdata('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_P.csv');

res=zeros(201,3);%这个分号使命令行窗口不出现0矩阵，简洁
for i=1:201
   res(i,1)=X1.data(i,1);
end
for i=1:201
   res(i,2)=X1.data(i,2);
end
for i=1:201
   res(i,3)=X2.data(i,2);
end
 
Sum=zeros(201,2);
for i=1:201
 temp1=res(i,3)*pi/180;
 Sum(i,1)=res(i,2)*cos(temp1);
end
for i=1:201
 temp2=res(i,3)*pi/180;
 Sum(i,2)=res(i,2)*sin(temp2);
end
xlswrite('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_1',res);
xlswrite('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_2',Sum);
%或者整合成一个文件；导出为XLS文件。
clear all;
close all;
 
X1=importdata('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_MAG.csv');
X2=importdata('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_P.csv');
 
res=zeros(201,5)
for i=1:201
   res(i,1)=X1.data(i,1);
end
for i=1:201
   res(i,2)=X1.data(i,2);
end
for i=1:201
   res(i,3)=X2.data(i,2);
end
for i=1:201
 temp1=res(i,3)*pi/180;
 res(i,4)=res(i,2)*cos(temp1);
end
for i=1:201
 temp2=res(i,3)*pi/180;
 res(i,5)=res(i,2)*sin(temp2);
end
xlswrite('C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12_P2',res);
保存成.csv格式.xls格式.xlsx格式都可以。
t=table(res);
writetable(t,'C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S12.csv');
t=table(res);
writetable(t,'C:\Users\Administrator\Desktop\2021.1.15\2021.1.15\N5_S122.xlsx');

3将数据部分处理完成后就是对S参数的整合过程，并将整合的S参数用图形绘制出来。

clear all;
close all;
 
X1=importdata('C:\Users\Administrator\Desktop\New folder 2021.1.13\1_X1_S12_MAG.csv');
X2=importdata('C:\Users\Administrator\Desktop\New folder 2021.1.13\1_X2_S12_MAG.csv');
X3=importdata('C:\Users\Administrator\Desktop\New folder 2021.1.13\1_X3_S12_MAG.csv');
X4=importdata('C:\Users\Administrator\Desktop\New folder 2021.1.13\1_X4_S12_MAG.csv');
S11=importdata('C:\Users\Administrator\Desktop\New folder 2021.1.13\1_S11_MAG.csv');
 
Sum=zeros(1,201);
for i=1:201
temp=S11(i,4)+X1(i,4)+X2(i,4)+X3(i,4)+X4(i,4)+1j*(S11(i,5)+X1(i,5)+X2(i,5)+X3(i,5)+X4(i,5));
Sum(1,i)=20*log10(abs(temp));
end
 
freq=zeros(1,201);
for i=1:201
   freq(1,i)=S11(i,1)/1000000000;
end
 
plot(freq,Sum,':','LineWidth',3);
grid on;
xlabel('fre/GHz');
ylabel('S11/dB');

成品