微波通信技术课程设计报告-Smith圆图的软件设计

1、微波实验课设姓名:壬變2008010430班级:通信0803班学号：微波通信技术课程设计报告smith圆图的软件设计简要说明smith圆图是p. h. smith于1939年在贝尔实验室发明的，它主要用 于计算微波网络的阻抗、导纳及网络阻抗匹配设计，还可用于设计微波 元器件。smith圆图软件不仅适用于微波工程设计，亦可用于电磁场、 微波技术及天线与电波传播等课程相关内容的教学，该软件通过形象的 演示可以深刻理解圆图的应用。微波网络的正弦稳态分析含有复数计算，运算十分繁琐和耗时。在 计算机运算速度和内存不够发达以前，图解分析法得到长足发展，其中 多年来应用最广的是smith圆图。在计算微波传输

2、线输入阻抗、导纳及 阻抗匹配等问题时，它不仅能避开繁琐的公式及复数运算，使工程设计 中相关计算简单便捷，而u图解过程物理概念清晰，所得结果直观形象。 随着计算机技术的飞速发展，图解法在计算精度上的固有缺陷日益显现, 因为，圆图的计算精度取决于圆图中必须有足够的圆周数，而过多的圆 周会导致图线过于密集，不便将阻抗、反射系数、电压驻波系数 （vswr ）及电长度等相关数据从图上直接读出。通过对圆图构成的基 本原理和应用问题的分析，利用现代计算机技术可以解决圆图计算精度 等方面存在的问题，为此设计的smith圆图软件既保留圆图计算直观、 便捷的大众性，又满足工程设计中相关参数的计算精度。在计算机应用

3、 日益普及的今天，该软件特别适合电磁场、微波技术与天线等领域的教 学和工程设计相关参数计算使用。2. 设计目的通过具体的软件编程和多媒体制作，进一步加深对微波通信技术的 理解和掌握，提高动手能力，提高解决实际问题的综合能力。3. 设计要求计算结果以图形和数据并行输出，处理包括复数的矩阵运算。为使程 序代码简单，执行运算速度快，计算精度高，选择matlab软件作为设 计技术平台较为理想（也可利用vc,vb等）。该软件数学表示、函数集丰 富且功能强大、良好的用户界面以及许多函数本身会绘制图形u自动选 取坐标刻度等显著优点，特别适用大量计算。利用mat lab强大的作图 功能容易画出完整的smi t

4、h圆图。整个圆图软件分为用户图形界面模块、 圆图计算模块、画图演示模块。上述:大模块又进一步分解，其中用户图 形界面模块分为：主页、主菜单；圆图计算模块分为反射系数计算、单 支节匹配计算、输入阻抗计算以及整个smith圆图；画图演示模块分为 等归一化电阻圆、等归一化电抗圆、反射系数圆等；确定阻抗值在圆图 上的位置、圆图的基本应用、求输入阻抗及其在圆图上的位置以及单支 节匹配等问题。4. 实验结果5. 总结这次课设的题目不简单，对于我来说，计算机语言学的不好，vb还只是 在高中学过一些，只会一些简单的语句，所以只有去找有没有相关的资料，然 后再研究，把代码的功能搞清楚，尽力去理解，更改了界面。并

5、且对于微波圆 图也进一步加深了印彖，对于这种简单的方法产生了兴趣，我认为，每一次课 设都能带给我一些新东四，不论最后得到了多少，也都比过去强一些。学到了 新知识同时还得到了更深的知识，感觉很有意思。极深刻理解了所学的知识乂 学到了新的知识提高不少，希望以后多锻炼得到更高的提高。附录：1.源程序(包括调试程序)priva te sub com mandl_clic k()endend subprivate sub comma nd2_click()if comman d2.caption ="显示波节”thencomman d2.caption ="显示波腹,p icture

6、l.li ne (500, 2 400)-(2500,2400),r gb(o, 0, 2 55)elsecom mand2.capt ion = ”显示波节”picturel .line (250 0, 2400)-(4500, 2400 ),rgb(0,255, 0)end ifend subprivate su b command3 _click()pi cturel.cispicturel.e nabled = f alsedim zo as doubledim rras doubledim rxas doub ledim gga s doubledi mgbas do ubledim

7、 ra d asdoubl edim radl as doubled im rad2 as doubledim taoas do ubledim pi as doublepi = 4* atn(l)rr=val(text3.text)rx=val(text4.text)zo =val(text7.text)tao=sqr(rr zo) a 2 +rxa 2)/ (rr+ z 0) a 2+r x a 2)tex tl.text ="0" +roun d(tao, 7)p icturel.ci rcle (2500,2400),t ao* 2000,rgb(0, 25 5,

8、255)pic turel.circ le (4500 -2000 /(r r/zo + 1), 2400),2000 /(rr /z0 + 1),rgb(o, 2 55, 255)if rx> 0 th enpicturel.circle (45 00, 2400 -2000 /rx / zo), 20 00 /rx/z0,rgb(0,255, 255), 3 * pi /2-2 *a tn(rx/ zo ), 3 * pi / 2end ififrx < 0 thenpi cturel.cir cle (4500,2400 - 20 00 /rx/zo), -2000 /rx/

9、 zo ,rgb(o, 2 55, 255),pi / 2, pi / 2 2 *atn(rx/zo)elsepictu rel.line (500, 2400)-(4500, 24 00),rgb(2 55, 255, 0 )end ifrad 1 =atn(rx / (rr-1)rad2 = a tn(rx/ (r r + 1)if radl >= 0 andrx>=0 then '1,3radl =radlend ifif radl >= 0 and rx 0 then radl = radl + p iend ifif ndrx>= 0 then 

10、9;2,4r adl = pi +radlend ifif and rx thenradl =radl + 2 *piend ifr ad= radl rad2if r ad>2 * p i thenrad =rad- 2 * piend iftext2.text =round(r ad, 7)gg=rr / (rr a 2 +rxa2) gb = -r x/ (rra2 +rxa 2 )text5.tex t =round(gg, 7) text 6.text =r ound(gb, 7 )end subpr ivate sub command4_c lick()if c ommand

11、4.ca ption ="设定特征阻抗” the ncommand4.caption 二”取消阻抗设定"te xt7.visibl e = truela bell2.visi ble = truelabel9.vis ible = tru eelsecomma nd4.captio n = ”设定特征阻抗,，text7.te xt = ,1hte xt7.visibl e = falsel abel9.visi ble = fals elabell2.v isible = f alseend ifend subpri vate sub c ommand5_ci ick()p

12、ictu rel.enable d = falsep icturel.ci sdim zo as doubledim rras dou bledim rx as doubled imggas d oubledim g bas doubl edim rad a s doubledi m radl as doubledim rad2 as do ubledim ta oas doubl edim pi as doublegg = round(te xt5.text,刀gb= rou nd(text6.t ext, 7)pi = 4 *atn(l)zo = vai (text7.tex t)rr =

13、 gg / (gga 2 +gb a 2)r x = -gb/(gga 2 +gb a 2)tao = sqr(r r- zo) a2 +rxa 2 ) / (rr+zo) a 2 +rx a 2)t ext 1.text = "0" +ro und(tao, 7 )picturel.circle (25 00, 2400),tao* 200 0,rgb(0,255, 255)p icturel.ci rcle (4500 - 2000/(rr/zo +1), 2400), 2000/(rr/zo +l),rgb(0,255, 255)if rx> 0 thenpic

14、turel .circle (4 500, 2400 - 2000 /r x/ zo), 2 000 /rx/zo,rgb(o , 255, 255 ), 3 * pi /2 - 2 *atn(rx/ z 0), 3 * pi / 2end ifif rx thenpicturel.c ircle (450 0, 2400 -2000 /rx / zo), -20 00 /rx/zo,rgb(o,255, 255), pi / 2,pi /2 2 *atn(rx / zo)elsepicturel.line (500,2400)-(45 00, 2400),rgb(255,0, 0)end i

15、 fradl = at n(rx / (rr 1)rad2 =atn(rx / (rr+ 1)if radl >= 0 andrx >= 0 then '1,3ra dl = radlend ifif radl >= oand rx nradl = rad 1 + piend ifif radl rx >= 0 then *2,4radl =pi + radlend ifif radl and r xradl = radl + 2*pien d ifrad=radl rad 2if rad>2 * pi the nrad = rad - 2 * pie n

16、d iftext2 .text =ro und(rad, 7 )text8.tex t =round(l - tao)/ (1 +tao ), 7) text9 .text =ro und(l +1 ao) / (1 -tao), 7) e nd subpriv ate sub co mmand6_cli ck() pictur el.enabled = truepic turel.cise nd subpriv ate sub co mmand7_cli ck()form2.showend su bprivate s ub drivel_change()en d subpriva te su

17、b com mand8_clic k()form3.s howend subprivate su bform_mou semove(but ton as i nt eger, shift as integ er, x as s ingle, y a s single)d im pi as d oublepi=4 *atn(l)dim a as d oublepictu rel.circle (2500, 24 00), 2000,rgb© 0,128)a = 20picturel.circle (25 00, 2400)za,rgb(0,0,128)a = 500pictu rel.

18、circle (2500, 24 00), a,rg b(0, 0,12 8)a = 1200picturel.c ircle (250 0, 2400),a,rgb(0,0,128)a =300pictur el.circle (4500 - a,2400), a,rgb(0, 0,128)a = 700picturel .circle (4 500 - a, 2 400), a,r gb(0, 0,1 28)a = 130 0picturel.circle (45 00 a, 24 00), a,rg b(0, 0,12 8)a = 1800picturel.c ircle (450 0

19、a, 240 0), a,rgb (0, 0,128 )a = 500pi cturel.cir cle (4500,2400 a), a,rgb(0,0,128),3 * pi /2 - 2 *at n(2000 / a ), 3 * pi /2a = 120 0picturel.circle (45 00, 2400 -a), a,rg b(0, 0,12 8), 3 * pi / 2 2 *atn(2000 / a), 3 *pi /2a =3500pictur el.circle (4500, 240 0 - a), a,rgb(0, 0,128), 3 *pi/2 -2 *atn(2

20、0 00 /a), 3 * pi/2a = 9000pic turel.circ le (4500,2400 a),a,rgb(0,0,128),3 * pi/2 2 *atn (2000/a), 3* pi /2a = -500picturel.c ircle (450 0, 2400 -a), -a,rg b(0, 0,12 8), pi / 2 , pi / 2 -2 *atn(2 000 / a)a = -1200pic turel.circ le (4500,2400 - a),-a,rgb(0z 0,128),pi / 2, p i / 2-2 *atn(2000 / a)a =

21、-3500pictur el.circle (4500, 240 0 - a), -a ,rgb(0, 0,128), pi / 2, pi /2 2 *a tn(2000 /a)a = -900 0picturel.circle (45 00, 2400 -a), -a,r gb(0, 0,1 28), pi /2, pi/2 -2 *atn(2000 / a)e nd subpriv ate sub pi cturel_cli ck()pictur el.enabled = falseen d subpriva te sub pic turel_mous emove(butt on as

22、inte ger, shift as intege r, x as si ngle, y as single)di m a as dou bledim zo as doubled im xx as d oubledim y yas doubl edim pi as doublepi = 4 *atn(l)dim tao as doubled imtaouas doubledim taovas d oubledim r ad asdoub ledim rra s doubledi mrxas do ubledim gg as doubledim gbas doubledim midi as do

23、 ubledim mi d2 as doub lezo = vai (text7.tex t)xx = cdb l(x)yy = c dbl(y)pict urel.cispi cturel.cir cle (2500,2400), 20 00,rgb(0,0,128)a = 20pictur el.circle (2500, 240 0), a,rgb (0, 0,128 )a = 500pi cturel.cir cle (2500,2400), a,rgb© 0,128)a = 1 200picture l.circle (2500, 2400 ), a,rgb(0, 0, 1

24、28)a = 300pic turel.circ le (4500 -a, 2400),a,rgb(0q 128)a = 700pictu rel.circle (4500 a , 2400), a ,rgb(0, 0,128)a =1300pictur el.circle (4500 - a,2400), a,rgb(0, 0,128)a = 1 800picture l.circle (4500 - a,2400), a,rgb© 0,128)a = 50 0picturel.circle (45 00, 2400 -a), a,rg b(0, 0,12 8), 3 * pi /

25、 2-2 *atn(2000 / a), 3 *pi /2a =1200pictur el.circle (4500, 240 0 a), a,rgb(0, 0,128), 3 *pi/2 -2 *atn(20 00 / a)z 3 * pi /2a = 3500pic turel.circ le (4500,2400 a)a,rgb(0q 128),3 * pi/2 2 *atn (2000/a), 3* pi/2 a = 9000picturel.c ircle (450 0, 2400 -a), a,rgb (0, 0, 128 ), 3 * pi / 2-2 *atn(2000 /a)

26、, 3 * p i /2a = -500picture l.circle (4500, 2400 - a), -a,rgb(0, 0,128), pi / 2, pi /2-2 *at n(2000 / a ) a = -1200picturel.c ircle (450 0, 2400 -a), -a,rg b(0, 0,12 8), pi/2, pi/2 -2 *atn(2 000 / a) a = -3500pic turel.circ le (4500,2400 - a),-a,rgb(0,0,128),pi / 2, p i / 2-2 *atn(2000 / a) a = -900

27、0pictur el.circle (4500, 240 0 a), -a ,rgb(0, 0,128), pi / 乙 pi /2 2 *a tn(2000 /a) midi = s qr(xx - 2 500) a 2 +(yy- 240 0) a 2)if 00 then '查找功能实现tao= mid 1 / 2000te xtl.text ="0" + rou nd(tao, 7)'反射系数计算pi cturel.cir cle (2500,2400), 20 00 *tao,rgb(255, 0,0)if xx enrad =atn(240 0 -yy) /(xx - 2500 )elserad = pi * (24 00 -yy) /abs(yy-2400) / 2e nd ifif xx >= 2500a nd yy 00 thenrad = radend ifif x thenrad=pi +raden d ifif xx >= 2500an d yy> 240 0 thenrad =rad+ 2 * piend ifif rad othenrad =rad+ 2 * p