资源描述：
《练习46 绘制smith图和波形图》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、练习46绘制Smith图和波形图知识背景从本练习开始，我们来学习MATLAB在电学方面的应用。在现代社会，电工与电子技术得到了广泛的应用。从收音机到电子计算机，从医疗诊断器械到飞机、火箭，我们随处可以感觉到电工与电子科学的存在。电工与电子是个统称。具体说起来，基本电路、过渡过程、变压器、电动机、继电器、接触器等属于电工范畴；而半导体二极管、三极管、各种放大电路、数字集成电路、晶闸管应用电路等属于电子范畴。MATLAB丰富的语句函数和绘图功能以及工具箱为我们更好的研究电工和电子领域的问题提供了可能。主要内容【
2、本练习讲述知识点】本练习作为铺垫，将要介绍电工与电子常见图形在MATLAB下的绘制。我们蒋哦使用一些循环语句和嵌套结构来完成特定功能，复习MATLAB基本绘图命令（包括坐标选取、坐标轴标注和图线标注）以及傅立叶分析常用函数fft的用法。练习过程（1）我们首先来学习绘制Smith图凡是从事高频工程领域研究的人，经常会遇到Smith图。所谓Smith图，是一种反映阻抗与复反射系数之间关系的图形。Smith图能够清楚的表现出阻抗与复反射系数的练习，具有很强的直观性。因为这些特点，Smith图得到了广泛应用。我们设
3、复反射系数为Y，阻抗为z（z＝I*x），则有：Y＝（z－1）/（z＋1）＝u＋iv绘制Smith图时我们以u为横坐标，v为纵坐标。设定两个参数r和x，并令其为常量，则我们可以得到一个圆，这就是Smith图。那么所得圆的方程是什么呢？由阻抗和反射系数关系式，我们得到两个实数方程：对于第一个式子，r为常数，是x的轨迹；对于第一个式子，x为常数，是r的轨迹。如果我们相画r＝0，0.2，0.5，1，1.5，2，4时的Smith图，应该在命令区中输入：plot([0,0],[-1.1,1.1]),holdon,xla
4、bel(‘u’),plot([0,0],[-1.1,1.1]),holdon,xlabel(‘v’),axisequal,axis([-1.1,1.1,-1.1,1.1]),gridtr=2*pi*(0:0.005:1);forr=[0,0.2,0.5,1,1.5,2,4]rr=1/(1+r);cr=1-rr;plot(cr+rr*cos(tr),rr*sin(tr))endforx=[0,0.2,0.5,1,1.5,2,4]rx=1/x;cx=rx;tx=2*atan(x)*(0:0.005:1);plo
5、t(1-rx*sin(tx),cx-rx*cos(tx))plot(1-rx*sin(tx),-cx+rx*cos(tx))end按回车后，我们得到图46－1所示的图形。在图中没有标注，读者可以利用gtext命令将每条线所代表的意义标注出来。上面我们利用将复数方程化为两个实数方程的方法，得到了Smith图。那么，还有没有别的方法得到Smith图呢？答案是肯定的。MATLAB有着优秀的复数计算功能，我们就用前面给出的复数方程，还是令r和x中一者为常数。我们来看看具体是怎么实现的。在命令区内输入：p=logsp
6、ace(-2,2);x=[-p,p];q=[0,0.05,0.10,0.20,0.50,1,2,5,10];forr=qz=(r+i*x-1)./(r+i*x+1);plot(z),holdon,axisequalendr=p;forx=[-q.q]z=(r+i*x-1)./(r+i*x+1);plot(z),holdon,axisequalendgrid,holdoff得到的图形，如图46－2所示。复数法中最为关键的语句是z=(r+I*x－1)/(r+I*x+1)，这一句使得z与x建立了确定的复数域里的函
7、数关系。（1）绘制波形图在电子科学中，经常要对波形进行调制，从而得到调制波。一般的调制波方程为：使振幅增量△A，频率增量△ω和相位增量△φ分别随信号变化，就形成调幅、调频和调相。ω0称为载频，信号频率一般低于载频，出于方便的考虑，在这里，我们取信号频率为0.1ω0。我们先来绘制调幅波形。在命令区内输入：clf,formatcompact,t=0:0.001:1;A0=10,w0=100,phi0=0;x=A0.*sin(w0*t+phi0);subplot(3,1,1),plot(t,x),holdxlab
8、el(‘t’),ylabel(‘载频信号’)dA=5*sin(0.1*w0*t),dw=0;dphi=0;y=(10+dA).*sin((w0+dw).*t+dphi);subplot(3,1,2),plot(t,y),ylabel(‘调制波形’)subplot(3,1,3),Y=fft(y);plot(abs(Y)),gridaxis([0,50,0,4000]),ylabel(‘信号频谱’)得到的结果如图46
