%双枝节阻抗匹配设计
syms s1 l1 s2 l2 s3 l3 s4 l4 %声明变量s1和s2分别对应两个方程，l1和l2分别是它们的解
wavelen=input('请输入波长值：');
phase=2*pi/wavelen; %相移常数
Z0=input('请输入特性阻抗值：');
a=input('请输入负载阻抗的实部：');
b=input('请输入负载阻抗的虚部：');
d=input('请输入两枝节间的距离（此距离为波长的倍数）：');
Y0=1/Z0;
Zl=a+i*b; %负载阻抗
Yl=1/Zl;
Gl=real(Yl);
Bl=imag(Yl);
zl=Zl/Z0; %归一化负载阻抗
yl=1/zl; %归一化导纳
t=tan(phase*d);
B1=-Bl+(Y0+sqrt((1+t*t)*Gl*Y0-Gl*Gl*t*t))/t; %第一个短截线的电纳
s1=tan(phase*l1)-B1/Y0; %关于l1的方程
e=solve(s1,l1); %解方程
l1=vpa(e,6); %取6位有效值
if B1<0; %通过B1判断l1的正负，如果l1小于0，则为其加上二分之一的波长，后面同理
l1=l1+wavelen/2;
end
l1=vpa(l1,6)
B2=(Y0*sqrt((1+t*t)*Gl*Y0-Gl*Gl*t*t)+Gl*Y0)/(Gl*t); %第二个短截线的电纳
s2=tan(phase*l2)-B2/Y0; %关于l2的方程
e=solve(s2,l2); %解方程
l2=vpa(e,6); %取6位有效值
if B2<0;
l2=l2+wavelen/2;
end
l2=vpa(l2,6)
B3=-Bl+(Y0-sqrt((1+t*t)*Gl*Y0-Gl*Gl*t*t))/t; %第一个短截线的电纳的另一个解
s3=tan(phase*l3)-B3/Y0; %关于l1另一个解的方程
e=solve(s3,l3); %解方程
l3=vpa(e,6); %取6位有效值
if B3<0;
l3=l3+wavelen/2;
end
l1f=vpa(l3,6)
B4=(-Y0*sqrt((1+t*t)*Gl*Y0-Gl*Gl*t*t)+Gl*Y0)/(Gl*t); %第二个短截线的电纳的另一个解
s4=tan(phase*l4)-B4/Y0; %关于l2另一个解的方程
e=solve(s4,l4); %解方程
l4=vpa(e,6); %取6位有效值
if B4<0;
l4=l4+wavelen/2;
end
l2f=vpa(l4,6)
