|
楼主 |
发表于 2006-10-13 10:33
|
显示全部楼层
是我没有说清楚.看原题
clear
clc
a=[1;1]
load canshu.mat %导入S参数
m=2 %端口数
n=3 %频点数
for p=1:n %频率循环
S=[S11(p,2)*exp(j*S11(p,3)),S12(p,2)*exp(j*S12(p,3));S21(p,2)*exp(j*S21(p,3)),S22(p,2)*exp(j*S22(p,3))]
c=S*a
b(1:m,p)=c % 改成b(2,j)=c不行
for k=1:m
R=b(k,p)/a(k)
T(k,p)=R %T为复反射系数
v(k,p)=(1+abs(T(k,p)))/(1-abs(T(k,p))) %将反射系数转化为驻波
end
end
plot(S11(:,1),abs(v(1,:)),'-r',S11(:,1),abs(v(2,:)),'-b')
title('复反射系数与频率的函数')
xlabel('频率MHz')
ylabel('复反射系数')
text(2010,2.5,'红色--1端口 蓝色--2端口')
上面的程序仅仅是2维的,如果是8维的话,表示S太麻烦了,请问有没有更好的办法? |
|