（全新）微波期中matlab仿真smith圆图.docVIP

微波工程期中课程设计 班 级： 2010 姓 名： 学 号： 序 号： 02 日 期： 2013年 月 EMAIL ：  目录 一：软件界面 2 二：功能介绍 2 1. 输入模块 2 2. 归一化阻抗和导纳求解 2 3. 反射系数和驻波比求解 2 4. 单支节双支节匹配 3 5. 演示部分 3 6. 清除 6 三：操作步骤 6 四：软件原理 6 A. 公式 6 B. 代码原理 8 1：输入部分 8 2：计算归一化阻抗和反射系数 9 3：单支节匹配 11 4：双支节匹配 15 五：总结 16 六：致谢 17  一：界面 ： 输入模块 此处输入负载的实部和虚部以及特征阻抗。 和导纳求解 按下计算按钮后此处显示归一化阻抗和导纳 和驻波比求解 、 此处显示反射系数和驻波比的计算结果 单支节双支节匹配 此处显示单支节和双支节匹配的计算结果（双支节需要输入d1和d2） 演示部分 归一化演示 单支节匹配演示 双支节演示 此键将左侧图重置，并清除所有数据。 操作步骤 输入负载的实部和虚部，写成不能空着。 计算按钮，得到归一化阻抗导纳及驻波比和反射系数 单支节匹配得到单支节匹配结果 distance1和distance2点击双支节匹配可得到 四：软件原理 公式 1: 求反射系数的公式。在onvertZ2AbsL等函数中用到 求反射系数的模，在[AbsL，h]函数中用。 与驻波比的关系。 求输入阻抗，在双支节匹配中，由y1旋转到y1‘用到，用于求第一个直接匹配后的导纳。 4： 用于计算按钮中，求驻波比和 5： 单支节匹配中用于求distance和length 函数中用到。 输入阻抗虚部不为用到。 6： 主要用于绘制圆的函数， 代码原理输入部分 输入部分类似都是用str2double函数读取edit类型框中的输入数值并进行检错，如果不是数字则报错。Handle.R1和X1以及Zc分别存入负载和特性阻抗的值。 2：计算归一化阻抗和反射系数 利用公式Zl=Z??∕0，得到归一化阻抗的实部ZIR虚部ZIX。（此处只能用handle来储存值，不用随意定义变量储存） 为了接下来计算方便（不使用handle计算，打字麻烦）读取刚才的结果值给Z_Re,Z_Im，这么做是因为不能直接将handle赋值给变量 上面这么做无法通过编 此处计算导纳值，complex为复数函数。 ConvertZ2AbsL_Ph 为计算|??|和??的函数，因为后面多次用到，所以独立成了一个函数，具体内容如下。 这个函数通过阻抗的实部和虚部返回反射系数和相角。然后用set将结果返回到GU框中，利用公式??=（1+??）/（1???）和K=1/??得到驻波比，同时返回给GUI的edit8和edit11 最后 这个Draw函数完成绘图，细节 Smithchart为的函数，每次赋予从绘。 Draw_Ratio完成绘制等反射系数Theta从pi一周x_zINy_zIN为以??的单位上。P定义了 绘制实数圆与以上大致相同。 绘制虚数部分，根据书上公式，虚数部分是以1,1/Z_Im为圆心，1/Z_Im为半径的一部分圆。所以r=1/abs（Z_Im），x_g和y_g为以r为半径的圆加上C_X和C_Y的偏移量，然后用find找到x_g，y_g在半径为1的园内的点。最后用plot函数画出。 3：单支节匹配 同计算，首先读取负载和特性阻抗。并赋值给Z_Re和Z_Im用于后面的计算 根据书上公式，首先判断导纳的虚部是否为0， if Y_Im == 0 如果为=（??/2±/（1-gL） L1=atan(sqrt(Y_Re)/(1-Y_Re))/2/pi; L2=-atan(sqrt(Y_Re)/(1-Y_Re))/2/pi; 如果不为则=（??/2±/（1-??） L1=atan(sqrt(P)/(1-P))/2/pi; L2=atan(-sqrt(P)/(1-P))/2/pi; 上面用来计算distance。 如果Y_Im为0，则 d= （/4）arccos（（gl-1）/（gl+1））。 由于matlab中acos只返回一个值，所以需要加上 theta2=2*pi-theta1; 当Y_Im不为0时， 令 theta1=acos((K-1)