单枝节阻抗匹配的Matlab实现

%单枝节阻抗匹配设计
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