基于matlab的史密斯圆图演示仿真图

1.算法描述

   史密斯图表(Smith chart，又称史密斯圆图)是在反射系散平面上标绘有归一化输入阻抗（或导纳）等值圆族的计算图。是一款用于电机与电子工程学的图表，主要用于传输线的阻抗匹配上。该图由三个圆系构成，用以在传输线和某些波导问题中利用图解法求解，以避免繁琐的运算。一条传输线(transmission line)的电阻抗力(impedance)会随其长度而改变，要设计一套匹配(matching)的线路，需要通过不少繁复的计算程序，史密夫图表的特点便是省略一些计算程序。

    史密斯图表的基本在于以下的算式：

   当中的Γ代表其线路的反射系数(reflection coefficient)，即S参数（S-parameter）里的S11，ZL是归一负载值，即ZL / Z0。当中，ZL是线路本身的负载值，Z0是传输线的特征阻抗（本征阻抗）值，通常会使用50Ω。
   图表中的圆形线代表电阻抗力的实数值，即电阻值，中间的横线与向上和向下散出的线则代表电阻抗力的虚数值，即由电容或电感在高频下所产生的阻力，当中向上的是正数，向下的是负数。图表最中间的点(1+j0)代表一个已匹配(matched)的电阻数值(ZL)，同时其反射系数的值会是零。图表的边缘代表其反射系数的长度是1，即100%反射。在图边的数字代表反射系数的角度(0-180度)和波长(由零至半个波长)。 

   史密斯圆图： 反射系数圆图+阻抗/导纳圆图。输入阻抗不易直接测得，通过测量反射系数间接获得输入阻抗值。史密斯圆图中反射系数和输入阻抗一一对应。

史密斯圆图的应用：

    已知传输线的特性阻抗=50Ω，传输线的负载阻抗为=(50+j50)Ω，求离负载z=0.1λ处的等效阻抗。解：=/=1+j1 R=1 X=1

    找到R=1的电阻圆和X=1的电抗圆的交点A，连接OA并延长交电长度刻度圆于B，OA顺时针旋转电长度0.1λ交电长度刻度圆于C，连接OC，交c1于P（据此读出|r|和Φ），过P点的R和X就是离负载z=0.1λ处的归一化输入阻抗。

2.仿真效果预览
matlab2022a仿真结果如下：

3.MATLAB核心程序

    %case1
    shuntL=Zo/(s*w)
    seriesC=1/(w*abs(deltaZP*Zo))
    figure(2)
    smith3
    gama1=(zInterceptP-1)/(zInterceptP+1) %reflection coeffients of the interscetion to 1+js
    t=num2str(s);
    s4=[ '  Y=1 - j' t]
    plot(gama1,'rx')
    text(real(gama1),imag(gama1),s4) %plot intersection point
    s1=['Shunt L ' num2str(shuntL) '  Series C '  num2str(seriesC)] %list 2 element match on figure
    title(s1)
    plot(gama,'rx') %plot desired complex impedance
    text(real(gama),imag(gama),s3)
    text(.5,1,s5)
    %check solution
    [zTest]=shuntL_seriesC(shuntL,seriesC,w,Zo)
    xlabel(['actual impdeance=' num2str(zTest)])
    
    %*******************************************
elseif (deltaZP < 0 & imag(s)==0) % Shunt L Series L zInterceptP
    %case2
    shuntL=Zo/(w*s)
    seriesL=abs(deltaZP*Zo)/w
    figure(3)
    smith3
    gama1=(zInterceptP-1)/(zInterceptP+1)
    t=num2str(s);
    s4=[ '  Y=1 - j' t]
    plot(gama1,'rx')
    text(real(gama1),imag(gama1),s4)
    s1=['Shunt L ' num2str(shuntL) '  Series L '  num2str(seriesL)]
    title(s1)
    plot(gama,'rx')
    text(real(gama),imag(gama),s3)
    text(.5,1,s5)
    [zTest]=shuntL_seriesL(shuntL,seriesL,w,Zo)
    xlabel(['actual impdeance=' num2str(zTest)])
    %*******************************************
end
 
    if (deltaZM < 0 & imag(s)==0) % Shunt C Series L zInterceptM
    %case3
    shuntC=s/(Zo*w)
    seriesL=abs(deltaZM*Zo)/w
    figure(4)
    smith3
    gama1=(zInterceptM-1)/(zInterceptM+1)
    t=num2str(s);
    s4=[ '  Y=1 +  j' t]
    plot(gama1,'rx')
    text(real(gama1),imag(gama1),s4)
    s1=['Shunt C ' num2str(shuntC) '  Series L '  num2str(seriesL)]
    title(s1)
    plot(gama,'rx')
    text(real(gama),imag(gama),s3)
    text(.5,1,s5)
    [zTest]=shuntC_seriesL(shuntC,seriesL,w,Zo)
    xlabel(['actual impdeance=' num2str(zTest)])
   
    %*******************************************
elseif (deltaZM > 0 & imag(s)==0) % Shunt C Series C zInterceptM
    %case4
    shuntC=s/(Zo*w)
    seriesC=1/(Zo*deltaZM*w)
    figure(5)
    smith3
    gama1=(zInterceptM-1)/(zInterceptM+1)
    t=num2str(s);
    s4=[ '  Y=1 + j' t]
    plot(gama1,'rx')
    text(real(gama1),imag(gama1),s4)
    s1=['Shunt C ' num2str(shuntC) '  Series C '  num2str(seriesC)]
    title(s1)
    plot(gama,'rx')
    text(real(gama),imag(gama),s3)
    text(.5,1,s5)
    [zTest]=shuntC_seriesC(shuntC,seriesC,w,Zo)
    xlabel(['actual impdeance=' num2str(zTest)])
   
    %*******************************************
end