%單枝節(jié)阻抗匹配設(shè)計(jì)
syms??s1?l1?s2?l2?????????????????????%聲明變量s1和s2分別對應(yīng)兩個(gè)方程，l1和l2分別是它們的解

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

Zl=a+i*b;??????????????????????%負(fù)載阻抗
zl=Zl/Z0;??????????????????????%歸一化負(fù)載阻抗
yl=1/zl;???????????????????????%歸一化導(dǎo)納

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

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