%雙枝節(jié)阻抗匹配設計
syms??s1?l1?s2?l2?s3?l3?s4?l4????????????????????%聲明變量s1和s2分別對應兩個方程，l1和l2分別是它們的解

wavelen=input（‘請輸入波長值：‘）;??????
phase=2*pi/wavelen;???????????????????%相移常數(shù)
Z0=input（‘請輸入特性阻抗值：‘）;?????
a=input（‘請輸入負載阻抗的實部：‘）;??
b=input（‘請輸入負載阻抗的虛部：‘）;?
d=input（‘請輸入兩枝節(jié)間的距離（此距離為波長的倍數(shù)）：‘）;?

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判斷l(xiāng)1的正負，如果l1小于0，則為其加上二分之一的波長，后面同理
?????l1=l1+wavelen/2;