function[A]=smithch（argzKmla）;
%用法一：A=smithch（z1）是已知?dú)w一化特性阻抗z1
%用法二：A=smithch（[VSWR?zmin?Zc]?是已知駐波比、電壓波節(jié)點(diǎn)的位置zmin（單位是波長(zhǎng)）和特性阻抗Zc
if?nargin<4;?la=1/8;?end;????%1/8波長(zhǎng)
if?nargin<3;?Km=0;?end;
if?nargin<1;?arg=1;?end;
if?length（arg）==3;
????z1=arg（1）;?zmin=arg（2）;?Zc=arg（3）;?Km=0;?ph=（1-zmin*4）*pi;?????%arg（）是取幅角
????gm=（z1-1）/（z1+1）*exp（-j*ph）;?z=（1+gm）/（1-gm）;?Kz=1;???%gm是反射系數(shù)
else;
????z1=arg（1）;?Zc=1;?Kz=0;
end;
%..................................................................繪制圓圖
t=（0:360）*pi/180;?x=cos（t）;?y=sin（t）;?r=[0.99?1.02];
plot（xy‘k‘）;?hold?on;?whitebg（‘w‘）;?axis?equal;?axis?off;???%繪制最外的等圓，whitebg為用色風(fēng)格設(shè)置
H=plot（[-1.02?1.02][0?0]）;?set（H‘color‘[0?0.5?0]）;
for?i=1:6;???????????????????????????????????????????????????%將外圓等分為24份
????u=r*cos（pi*i/12）;?v=r*sin（pi*i/12）;
????plot（uv‘k‘u-v‘k‘-uv‘k‘-u-v‘k‘）;
end;
text（-1.060‘0‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...?????%給主要的4個(gè)分點(diǎn)做標(biāo)注
????‘color‘[0?0?0.5]‘fonta‘‘italic‘）;
text（-0.061.04‘0.125‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...
????‘color‘[0?0?0.5]‘fonta‘‘italic‘）;
text（1.030‘0.25‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...
????‘color‘[0?0?0.5]‘fonta‘‘italic‘）;
text（-0.06-1.04‘0.375‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...
????‘color‘[0?0?0.5]‘fonta‘‘italic‘）;
i=150:170;?xi=x（i）*1.1;?yi=y（i）*1.1;
H=plot（[-1?xi][0.5?yi]‘k‘[-1?xi][-0.5?-yi]‘k‘）;???%繪制向源和向負(fù)載箭頭
text（-10.78‘向源‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘9）;???%箭頭類型標(biāo)注
text（-1-0.7‘向負(fù)載‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘9）;
plot（0.5*x+0.50.5*y‘b‘）;???????%繪制r=1的圓，因?yàn)橥鈭A半徑為1，所以這里用0.5做半徑
text（0.020.33‘1‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘b‘）;
if?Km<1;
????i=1:10:361;?xi=x（i）;?yi=y（i）;
????gm=1/3;?H=plot（gm*xigm*yi‘:‘）;?set（H‘color‘[0?0.5?0]）;????%繪制模值為1/3和2/3的反射系數(shù)圓
????gm=2/3;?H=plot（gm*xigm*yi‘:‘）;?set（H‘color‘[0?0.5?0]）;
text（-0.230.19‘1/3‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...
?????‘color‘[0?0.5?0]）;
text（-0.580.48‘2/3‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10?...
?????‘color‘[0?0.5?0]）;
a=1/（1+2）;?plot（a*x+1-aa*y‘b‘）;???%r=2，繪制電阻圓圖
a=1/（1+0.5）;?plot（a*x+1-aa*y‘b‘）;??%r=0.5，繪制電阻圓圖
text（0.360.033‘2‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘b‘）;
text（-0.420.033‘0.5‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘b‘）;
text（-0.970.033‘0‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘b‘）;
xx=x+1;?yy=y+1;?I=find（（xx.*xx+yy.*yy）<=1.001）;???%用find函數(shù)返回所有小于等于1.001的圓圖
plot（xx（I）yy（I）‘r‘xx（I）-yy（I）‘m‘）;???%x=-11，繪制電抗圓圖
a=1/2;?xx=a*x+1;?yy=a*y+a;?I=find（（xx.*xx+yy.*yy）<=1.001）;
plot（xx（I）yy（I）‘r‘xx（I）-yy（I）‘m‘）;?%x=-22，電抗圓圖
a=1/0.5;?xx=a*x+1;?yy=a*y+a;?I=find（（xx.*xx+yy.*yy）<=1.001）;
plot（xx（I）yy（I）‘r‘xx（I）-yy（I）‘m‘）;????%x=-0.5?0.5，電抗圓圖
text（0.530.79‘2‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘r‘）;
text（-0.040.95‘1‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘r‘）;
text（-0.630.72‘1/2‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘r‘）;
text（0.52-0.79‘-2‘‘fontname‘‘a(chǎn)rial‘‘fontsize‘10‘color‘‘m‘）;
text（-0.05-0.95‘-1‘‘fontnam

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????9091??2009-04-09?20:37??smithch.m