%單枝節阻抗匹配設計
syms??s1?l1?s2?l2?????????????????????%聲明變量s1和s2分別對應兩個方程，l1和l2分別是它們的解

wavelen=input（‘請輸入波長值：‘）;??????
phase=2*pi/wavelen;???????????????????%相移常數
Z0=input（‘請輸入特性阻抗值：‘）;?????
a=input（‘請輸入負載阻抗的實部：‘）;??
b=input（‘請輸入負載阻抗的虛部：‘）;?

Zl=a+i*b;??????????????????????%負載阻抗
zl=Zl/Z0;??????????????????????%歸一化負載阻抗
yl=1/zl;???????????????????????%歸一化導納

if?a==Z0???????????????%如果負載阻抗的實部等于特性阻抗則只有一個根
????t1=-b/（2*Z0）;
else?
????t1=（b+sqrt（a*（（Z0-a）*（Z0-a）+b*b）/Z0））/（a-Z0）;???%對應兩個根
????t2=（b-sqrt（a*（（Z0-a）*（Z0-a）+b*b）/Z0））/（a-Z0）;
end
if??t1>=0
????????d1=wavelen*atan（t1）/（2*pi）???????????%求d的公式，對應兩種情況
else????d1=wavelen*（atan（t1）+pi）/（2*pi）
end

if??t2>=0
????????d2=wavelen*ata