双枝节阻抗匹配的Matlab实现

%双枝节阻抗匹配设计
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（s1l1）;          %解方程
l1=vpa（e6）;              %取6位有效值
if B1<0;         %通过B1判断l1的正负，如果l1小于0，则为其加上二分之一的波长，后面同理
     l1=l1+wavelen/2;