基于MATLAB的电力电子技术仿真分析

1、PAGE PAGE 44电力电子技术课程设计报告题 目： 基于MATLAB的电力电子技术 仿真分析 院 （系）： 机电与自动化学院 专业班级： 电气工程及其自动化 学生姓名： 学 号： 指导教师： 2014年1月13日至2014年1月17日华中科技大学武昌分校电力电子技术课程设计任务书一、设计（调查报告/论文）题目基于MATLAB的电力电子技术仿真分析二、设计（调查报告/论文）主要内容1.晶闸管的仿真模型及以单相半波整流器为例，说明晶闸管元件应用系统的建模与仿真方法。2.晶闸管三相桥式整流带电阻性负载时系统的建模与仿真。3.绝缘栅双极型晶体管元件的仿真模型及一个由IGBT元件组成的Boost变

2、换器的建模与仿真。4.相位控制的晶闸管单相交流调压器带电阻性负载时系统的建模与仿真。三、原始资料MATLAB仿真软件四、要求的设计（调查/论文）成果编写详细的设计说明书（附上本次设计心得体会） 说明书中完成相应系统模型的建模、参数设置及仿真调试，写出设计报告。1.晶闸管的仿真模型、参数设定方法、以单相半波整流器为例说明晶闸管元件应用系统的建模与仿真方法，记录相应波形。2.晶闸管三相桥式整流带电阻性负载时系统的建模过程与仿真调试，记录波形。3.绝缘栅双极型晶体管元件的仿真模型、参数设定方法、及由IGBT元件组成的Boost变换器的建模与仿真，记录波形。4.相位控制的晶闸管单相交流调压器带电阻性负

3、载时系统的建模与仿真，记录波形。以上4个类型的仿真过程中皆需包含具体电路形式，工作过程分析，在MATLAB中的建模，各组成环节的参数设置过程，仿真波形，波形分析。五、进程安排1. 下达设计任务书，讲解设计要求、进度安排、指导时间、注意事项等，提供参考资料。（0.5天）2. 学习并熟练MATLAB Simulink/Power System 工具箱等相关内容（1.5天）3. 典型电力电子器件的仿真模型建模及仿真实例（0.5天）；4. 典型电力电子变换器的应用仿真。（1天）5. 撰写课程设计报告。（0.5天）6. 答辩。（1天）六、主要参考资料1 王兆安，刘进军.电力电子技术（第五版）.北京：机械

4、工业出版社，2010.2 周渊深. 电力电子技术与MATLAB仿真.北京：中国电力出版社, 2005.3 林飞,杜欣. 电力电子应用技术的MATLAB仿真.北京：中国电力出版社,2009.4 洪乃刚. 电力电子、电机控制系统的建模和仿真. 北京：机械工业出版社,2010.指导教师（签名）：20 年 月 日课程设计成绩评定表成绩评定项 目比例得 分平时成绩（百分制记分）30%业务考核成绩（百分制记分）70%总评成绩（百分制记分）100%评定等级优 良 中 及格 不及格指导教师（签名）：20 年 月 日目录1课程设计目的12课程设计主要内容13课程设计题目描述与要求13.1课程设计题目描述13.2

5、课程设计要求24 各电路的建模与仿真24.1单相半波可控整流器24.2晶闸管三相桥式整流电路94.3Boost变换器164.4相位控制的晶闸管单相交流调压器 195课程设计总结221 课程设计目的通过电力电子计术的课程设计达到以下几个目的：（1）培养学生文献检索的能力，特别是如何利用Internet检索需要的文献资料；（2）培养学生综合分析问题、发现问题和解决问题的能力；（3）培养学生运用知识的能力和工程设计的能力；（4）提高学生课程设计报告撰写水平；（5）提高学生通过实验测试、研究分析和完善设计的水平。2 课程设计主要内容（1）晶闸管的仿真模型及以单相半波整流器为例，说明晶闸管元件应用系统的

6、建模与仿真方法。（2）晶闸管三相桥式整流带电阻性负载时系统的建模与仿真。（3）绝缘栅双极型晶体管元件的仿真模型及一个由IGBT元件组成的Boost变换器的建模与仿真。（4）相位控制的晶闸管单相交流调压器带电阻性负载时系统的建模与仿真。3 课程设计题目描述与要求3.1课程设计题目描述本次课程设计包含了六个内容的建模与仿真：1晶闸管的仿真模型及以单相半波整流器为例，说明晶闸管元件应用系统的建模与仿真方法；2晶闸管三相桥式整流系统的建模与仿真；3. 可关断晶闸管的仿真模型及以可关断晶闸管元件组成的Buck变换器为例的仿真过程；4绝缘栅双极型晶体管元件的仿真模型及一个由IGBT元件组成的Boost变换

7、器的建模与仿真；5相位控制的晶闸管单相交流调压器系统的建模与仿真；6晶闸管三相半波有源逆变器的建模与仿真。这六个内容基本包含了电力变换的四大类，从中能比较全面的掌握电力电子MATLAB仿真的方法。此仿真实验主要涉及到以下四个方面，而基于MATLAB的电力电子技术仿真则是一下几个内容很好的结合。电力电子器件：电力电子器件是一系列固态高电压、大电流的电子器件，被控对象的设备功率很大。按可控性可分为三类：不控器件（二极管）、半控器件（晶闸管）、全控器件（GTR、GTO、IGBT、MOSFET等）。电力电子技术应用：该技术广泛应用于多种形式的电源、电力拖动控制、电网电能质量技术提高以及大功率电能传输。

8、MATLAB仿真：MATLAB程序设计语言是美国Math Works公司在20世纪80年代中期推出的搞性能数值计算软件，2005年8月该公司就推出MATLAB7.1版，现已成为线性代数、自动控制理论、数理统计、数字信号分析与处理、动态系统仿真等各种课程的基本数学工具。电力电子技术MATLAB实践：电力电子技术中有关电能的变换与控制过程，有各种电路原理的分析与研究、大量的计算、电能变换的波形测量、绘制与分析等，都离不开MATLAB。首先，它的运算功能强大，应用于交流电的可控整流、直流电的有源逆变与无源逆变中存在的整流输出的平均值、有效值、与电路功率计算、控制角、导通角计算。其次，MATLAB的S

9、impowerSystems实体图形化仿真模型系统，把代表晶闸管、触发器、电阻、电容、电源、电压表等实物的特有符号连接成一个整流装置电路或是一个系统，更简单方便，节省设计制作时间和成本等。并且，交流技术讨论的电能转换与控制，需要对各种电压与电流波形进行测量、绘制与分析，MATLAB提供了功能强大且方便使用的图形函数，特别适合完成这项任务。3.2课程设计要求编写详细的设计说明书,说明书中完成相应系统模型的建模、参数设置及仿真调试，写出设计报告。（1）晶闸管的仿真模型、参数设定方法、以单相半波整流器为例说明晶闸管元件应用系统的建模与仿真方法，记录相应波形。（1）晶闸管三相桥式整流带电阻性负载时系统

10、的建模过程与仿真调试，记录波形。（1）绝缘栅双极型晶体管元件的仿真模型、参数设定方法、及由IGBT元件组成的Boost变换器的建模与仿真，记录波形。（1）相位控制的晶闸管单相交流调压器带电阻性负载时系统的建模与仿真，记录波形。以上4个类型的仿真过程中皆需包含具体电路形式，工作过程分析，在MATLAB中的建模，各组成环节的参数设置过程，仿真波形，波形分析。4各电路的建模与仿真4.1单相半波可控整流器4.1.1晶闸管的仿真 晶闸管模型晶闸管是一种门极信号触发导通的半导体器件。晶闸管有两个输入端和两个输出端，第一个输入与输出是阳极媏（a）与阴极端（k），第二个输入（g）是门极控制信号端如图4-1，当

11、勾选“Show measurement port”项时便显示第二个输出端（m）如图4-2，这是晶闸管检测输出向量Iak Uak端，可连接仪表检测流经晶闸管的电流（Iak）与晶闸管的正向压降（Uak），晶闸管组件的符号和仿真模型图如图所示。 图 4-1 图 4-2 晶闸管参数及其设置在模型结构图中，当鼠标双击模型时，则弹出晶闸管参数对话框，如下图4-3所示图 4-3“Resistance Ron(Ohms)”：晶闸管导通电阻Ron（）。“Inductance Lon（H）”：晶闸管元件内电感Lon（H）。电感参数与电阻参数不能同时设 为0“Forward voltage Vf（V）”：晶闸管元件

12、的正向管压降Vf（V）。“Initial current Ic（A）”：初始电流Ic（A）。“Snubber resistance Rs（ohms）”：缓冲电阻Rs（）。“Snubber capacitance Cs（F）”：缓冲电容Cs（F）。可对Rs与Cs设置不同的数值以改变或者取消吸收电路。“Show measurement port”为设置是否显示检测端（m）。需要说明的是，含有晶闸管模型的电路仿真时，最好采用特定的算法Ode23tb与Oder15s，而当电路进行离散化处理时，晶闸管的内电感量应设为0。4.1.2单相半波可控整流电路的仿真电路图及工作原理 单相半波可控整流电路（阻感负载

13、）图 如上图所示，当晶闸管VT处于断态时，电路中电流Id=0，负载上的电压为0，U2全部加在VT两端，在触发角处，触发VT使其导通，U2加于负载两端，由于电感L的存在使电流id不能突变，id从0开始增加同时L的感应电动势试图阻止id增加，这时交流电源一方面供给电阻R消耗的能量，一方面供给电感L吸收的电磁能量，到U2由正变负的过零点处处id已经处于减小的过程中，但尚未降到零，因此VT仍处于导通状态，当id减小至零，VT关断并承受反向压降，电感L延迟了VT的关断时刻使Ud波形出现负的部分。（2）建立仿真模型根据原理图用matalb软件画出正确的仿真电路图，整体模型如下图所示单相半波晶闸管可控整流电

14、路（阻感负载）的仿真模型仿真参数：选择ode23tb算法，将相对误差设置为1e-3，开始仿真时间设置为0，停止仿真时间设置为0.12，如下图所示 模型参数简介与设置 交流电压源提取路径： SimulinkSimPoweSystenElectricalAC Voltage Source “Peak amplitude”:正弦电压峰值Um，单位V， “Phase”：正弦电压初相角,单位 度， “Frequency”：正弦电压频率f，单位 Hz， “Sample time”：采样时间 ，单位 s， 本实验参数设置为频率50Hz，电压幅值220V，其他为默认设置，如下图所示。 晶闸管提取路径：Simu

15、linkSimPoweSystenPower Electronics Thyristor设置“Snubber resistance Rs（ohms）”缓冲电阻Rs=500，“Snubber capacitance Cs（F）”：缓冲电容Cs为无穷大inf ，其他为默认设置，如下图所示。 RLC元件提取路径：SimulinkSimPoweSystenElements Series RLC Branch设置“Resistance (Ohms)”电阻R=1，“Inductance Lon（H）”电感L=5e-3H，“capacitance（F）”电容为无穷大inf，“measurements”测量选

16、None。如下图所示 脉冲信号发生器提取路径：SimulinkSimlinkSourcePulse Generator“Amplitude”：脉冲幅值，“Period(secs)”：周期（秒），“Pulse Width(% of Period”：脉冲宽度（周期的百分数），“Phase delay(secs)”：相位延迟（秒）。 振幅A=3V，周期T=0.02，占空比10%，时相延迟（1/50）x（/360）s，如下图所示， 为移相控制角示波器设置Number of axes 为5，显示5段波形，分别为脉冲电压Ug，晶闸管两端电压UVT，负载电流id，负载电压ud，电源电压U2，如下图所示。 仿

17、真结果设置触发脉冲分别为0、30、60、90、120。其产生的相应波形分别如图所示。在波形图中第一列为脉冲电压Ug波形，第二列为晶闸管两端电压UVT波形，第三列为负载电流id波形，第四列为负载电压ud波形，第五列为电源电压U2 波形。 阻感负载 触发角=0阻感负载 触发角=30阻感负载 触发角=60阻感负载 触发角=90阻感负载 触发角=1204.2晶闸管三相桥式整流电路电路图及工作原理 以=0为例，6个晶闸管 的导通顺序为VT1-VT2-VT3-VT4-VT5-VT6，触发脉冲为宽脉冲宽度大于60，保证了每个时刻均有两个晶闸管导通，当VT1-VT2导通时桥臂输出电压为Uac，然后VT2-VT

18、3导通输出电压为Ubc，VT3-VT4导通输出电压为Uba，VT4-VT5导通输出电压为Uca，VT5-VT6 导通输出电压为Ucb， VT6-VT1导通输出电压为Uab。因此输出整流电压Ud波形为线电压在正半周的包络线。三相桥式整流电路原理图建立仿真模型根据原理图用matalb软件画出正确的仿真电路图，整体模型如下图所示 三相桥式全控整流系统（电阻负载）的仿真模型仿真参数：选择ode23tb算法，将相对误差设置为1e-3，开始仿真时间设置为0，停止仿真时间设置为0.05，模型参数简介与设置交流电压源三相交流电源通过三个频率为50Hz、幅值为220V、相位两两相差120，A相的设置如右图所示，

19、另外两相设置为B相相位滞后A相120，Phase设置为-120，C相相位超前A相120，Phase设置为120，测量“measurements”三相都要选Voltage，以便使用万用表测量电压。 通用桥输入端A,B,C为三相交流的相电压输入端子，输入端g为触发脉冲输入端子，+,-为整流器输出正负极端子。“Number of bridge arms”：通用整流桥臂的相数，“Snubber resistance Rs（ohms）”：缓冲电阻Rs（），“Snubber capacitance Cs（F）”：缓冲电容Cs（F），“Power Electronic device”：电力电子器件的种类，默

20、认晶闸管“Ron(Ohms)”：器件内电阻（），“Lon(H)”：器件内电感（H），“Forward voltage Vf(V)”：整流桥门槛电压(伏）.在测量“Measurements”选“All voltages and currents”（全部电压和电流）以便测量桥臂内晶闸管的电压和电流，其他参数为默认值.如下图所示。 常量三相桥式全控整流系统仿真模型要使用两个常量模块，一个提供触发角的值，一个设置为0连接同步6脉冲触发器的使能端Block，使其能正常工作。如下图所示。同步6脉冲触发器输入端alpha_deg为移相控制角给定信号，用常量模块constant输入控制角，输入端AB,BC,C

21、A是同步线电压输入端，输入端Block是触发器的使能端，当此端置0时，则输出脉冲，输出端Pulse是触发脉冲的输出，它是一个6维向量，即6个触发脉冲 “Frequency of synchronisation voltages(Hz)”：同步电压频率（Hz),“Pulse width(degrees)”：脉冲宽度（度）频率设置为50Hz，脉冲用宽脉冲设置为80，如下图所示。万用表三相桥式全控整流系统仿真模型使用了两个万用表，其中一个万用表的参数如下图所示，选中Isw1和Usw1，点击【】移入右侧的对话框中，分别测量iVT1，uVT1 ，另一个万用表选择Usrc：A，Usrc：B，Usrc：C，

22、分别测量A,B,C三相电压。 示波器三相桥式全控整流系统仿真模型使用了两三个示波器，最主要的一个设置Number of axes 为4，显示4段波形，分别为负载电压ud，负载电流id，脉冲信号电压Ug，A,B,C三相电压，与万用表连接的示波器，设置Number of axes 为2，显示2段波形，分别为晶闸管VT1的电压和电流，另一个示波器设置Number of axes 为3，显示3段波形，分别为A,B,C三相的电流。电压电流测量由于同步6脉冲触发器的AB,BC,CA端为同步线电压输入端，而三相电源提供的是相电压所以要通过三个电压表进行转换，其他电流电压测量无需设置直接使用RLCR=10，L

23、=0H,C=inf（无穷大） 仿真结果设置触发脉冲=0，负载电压ud，负载电流id，脉冲信号电压Ug，A,B,C三相电压晶闸管VT1的电压和电流通过A,B,C三相的电流设置触发脉冲分别为30、60、90、110，产生的相应波形分别如图所示，第一列为负载电压ud波形，第二列负载电流id波形，第三列脉冲信号第四列电压Ug波形，A,B,C三相电压波形。电阻负载 触发角=30电阻负载 触发角=60电阻负载 触发角=90电阻负载 触发角=1104.3 Boost变换器4.3.1绝缘栅双极型晶体管（IGBT）的仿真绝缘栅双极型晶体管模型绝缘栅双极型晶体管（IGBT）是一种栅极信号触发导通的全控型器件。晶闸

24、管模型有两个输入端和两个输出端，第一个输入与输出是集电极（C）与发射极（E），第二个输入（g）是栅极控制信号端如图，当勾选“Show measurement port”项时便显示第二个输出端（m）如图，这是IGBT检测输出向量Iak Uak端，可连接仪表检测流经IGBT的电流（Iak）与正向压降（Uak），IGBT组件的符号和仿真模型图如下图所示。 绝缘栅双极型晶体管参数及其设置在模型结构图中，当鼠标双击模型时，则弹出晶闸管参数对话框，如下图所示由图可知，IGBT的参数设置与普通晶闸管的参数设置几乎完全相同，另有2个参数类似GTO参数设置。“Current 10% fall time Tf(s

25、)”：电流下降时间Tf。“Current tail time Tt(s)”：电流拖尾时间Tt。对于绝缘栅双极型晶体管（IGBT）模型的电路仿真时，同样宜采用Ode23tb与Oder15s算法。4.3.2 Boost变换器的仿真 电路图及工作原理首先假设电路中电感L的值很大，电容C值 也很大。当IGBT处于通态时，电源E向电感L充电，充电电流基本恒定为I1，同时电容C上的电压向负载R供电。因C值很大，基本保持输出电压u0为恒值，记为U0 。设IGBT处于通态的时间为ton，此阶段电感L上积蓄的能量为EI1ton。当IGBT处于断态时E和L共同向电容C充电并向负载R提供能量。设IGBT处于断态的时

26、间为toff，则在此期间电感L释放的能量为（U0 -E）I1toff。当电路工作于稳态时，一个周期T中电感L上积蓄的能量与释放的能量相等 EI1ton=（U0 -E）I1toff 化简为 U0=T*E/toff 输出电压高于电源电压 升压斩波电路（电阻负载）原理图 建立仿真模型 仿真参数：选择ode23tb算法，将相对误差设置为1e-3 开始仿真时间设置为0，停止仿真时间设置为0.004。模型参数简介与设置直流电压源 设置A=100V，“measurements”选None（不测量电压），如下图所示。 脉冲信号发生器设置振幅A=3V，周期T=0.0001，占空比40%，时相延迟0s，如下图所示

27、。绝缘栅双极型晶体管（IGBT） 勾选“Show measurement port”项，其他为默认设置，如下图所示。二极管勾选“Show measurement port”项，其他为默认设置，如下图所示。RLC元件 R=50,L=0H，C=3e-6F，如下图所示。示波器设置Number of axes 为5，显示5段波形，如下图所示。 仿真结果产生的相应波形分别如图所示。在波形图中第一列为通过电感的电流i1波形，第二列为负载电流i0波形，第三列为负载电压u0波形，第四列为IGBT电流iV波形，第五列为脉冲电压Ug波形。4.4相位控制的晶闸管单相交流调压器电路图及工作原理在交流电源U1的正半周和

28、负半周，分别对VT1和VT2的触发延迟角进行控制，使得输出电压波形为正弦电压的一部分，从而实现调节输出电压的目的，负载阻抗角=arctan(L/R),负载电压相位滞后于晶闸管输出电压相位，把=0的时刻定在电源电压过零的时刻，显然阻感负载下稳态时的移相范围为-。 建立仿真模型根据原理图用matalb软件画出正确的仿真电路图，整体模型如下图所示。仿真参数：选择ode23tb算法，将相对误差设置为1e-3,开始仿真时间设置为0，停止仿真时间设置为0.12。模型参数设置交流电压源 参数设置为频率50Hz，电压幅值220V,“measurements”测量选“Voltage”，其他为默认设置，如下图所示

29、。 脉冲信号发生器振幅A=1.1V，周期T=0.02，占空比0.001，时相延迟（1/50）x（/360）其他为默认设置，为移相控制角两个脉冲信号发生器相位相差180,如下图所示。晶闸管 不勾选“Show measurement port”其他均为默认设置。RLC元件R=1，H=0，C=inf。示波器设置Number of axes 为5，显示5段波形，如下图所示 仿真结果设置Pulse1触发脉冲角1分别为60、90，Pulse2触发脉冲角2则应对应为240、270，产生的相应波形分别如下图所示，第一列为交流电源U1波形，第二列负载电流i0波形，第三列电压u0波形，第四列晶闸管两端电压uVT波

30、形。电阻负载 触发角=60电阻负载 触发角=905.课程设计总结我们这次的课题是应用Matlab软件对一些典型的电路进行仿真，刚拿到题目时有些无从下手，毕竟从来没有学过Matlab的使用，后来在查阅辅导书后我基本明白了仿真方法，感觉并不太难，于是以为做起来应该没什么问题，但在实际动手的过程中我才发现，想起来容易做起来难，理论到实践之间还有着巨大的努力空间。在这次的课程设计中，通过在Matlab软件对一些典型的电路进行仿真，我极大地锻炼了对于Matlab的应用能力。在实际操作过程中我由于不熟练、粗心以及对理论知识的掌握不到位，遇到了形形色色的问题，为了使电路正常工作，必须灵活运用原理找出解决方法

31、。在同学的帮助下，我一一解决了这些问题。有很多课堂上并没有讲授的知识应用中需要用到，我则通过去图书馆借阅资料、在网上搜索等形式自己吸收补充。这整个过程很好地巩固了我们学过的专业知识，使我对电路系统结构也有了更进一步的了解和认识，同时对Matlab等系列知识都有了一定的了解，也使我们把理论与实践从真正意义上相结合了起来,考验了我们借助互联网络搜集、查阅相关文献资料、和组织材料的综合能力。这次的课程设计使我受益良多，最后我谨向在这个过程中给予我帮助的老师以及各位同学们致以最诚挚的谢意。附录资料：matlab绘图指令大全绘图指令1 二维曲线图1.1 绘制折线图plot指令图例Y=1,3,6,5,9,

32、0,2;plot(Y);X=0: pi/10: pi*2;Y=sin(X);plot(X,Y);X=0: pi/10: pi*2;Y1=sin(X);Y2=cos(X);Plot(X,Y1,X,Y2);调整坐标范围：axisaxis(0,300,0,2)1.2 绘制自定义函数DrawCircle.mfunction DrawCircle(Point,Radius) Hold on t=0: pi/10: 2*pi; x=Point(1)+ Radius*cos(t); y=Point(2)+ Radius*sin(t); plot(x,y);DrawCircle(10 10,1)DrawCir

33、cle(20 10,2)DrawCircle(10 20,3)1.3 绘制符号函数显函数ezplot(sin(x),0,2*pi)隐函数ezplot(x2+y2-10,-5,5,-6,6)参数方程ezplot(cos(t)3,sin(t)3,0,2*pi)1.4 绘制自定义函数function y=myf1(x) y=sqrt(100-x2);fplot(myf1,-15 15)fplot(sin(x) cos(x) myf1(x),-15 15)1.5 图形修饰 设置颜色 y m c r g b w k 设置线型 - : -. - 设置标记 . o x + * 指令图例Y=1,3,6,5,9

34、,0,2;plot(Y, r-+);X=0: pi/10: pi*2;Y=sin(X);plot(X,Y, b-.);X=0: pi/10: pi*2;Y1=sin(X); Y2=cos(X);plot(X,Y1,r+-, X,Y2,b-*); 在指定坐标处，书写文字：text(3.5, 0.6, 曲线比较);x=1.6*pi, 1.6*pi; y=-0.3, 0.8;s=曲线cos; 曲线sin; text(x,y,s);1.6 更多类型的二维图指令图例bar直方图X=0:pi/10:2*pi;Y=sin(X);bar(X,Y);polar极坐标图T=0: pi/10: 4*pi;R=T;p

35、olar(T, R);误差棒棒图X=0:pi/10:2*pi;Y=sin(X);e=0.2*rand(size(X);errorbar(X,Y,e);火柴杆图X=0:pi/10:2*pi; Y=sin(X);stem(X,Y);stairs楼梯图X=0:pi/10:2*pi; Y=sin(X);stairs(X,Y);多边形填色图X=1,2,3,4,5; Y=3,5,2,1,6;fill(X,Y,r);hold on; % 保持图形plot(X,Y,o)1.7 数值函数的二维图 可用于绘图，更可用于采样取点。 fplot(0.5*cos(x),-pi,pi) % 绘图X,Y = fplot(0

36、.5*cos(x),-pi,pi); % 返回点坐标fplot(cos(x),-pi,pi,r-+); % 观察点的位置控制采样点的密度fplot(cos(x),-pi,pi,r-+,0.05);fplot(cos(x),-pi,pi,r-+,0.1); 可绘制系统函数，也可绘制自定义函数的图形。2 三维曲线图2.1 三维曲线plot3指令图例X=0: 0.1: 8*pi;Y=sin(X);Z=cos(X); plot3(X,Y,Z,r);X=0: 0.1: 8*pi;Y=sin(X);Z1=cos(X);Z2=2*cos(X); plot3(X,Y,Z1,r, X,Y,Z2,b);2.2 三

37、维面填色fill3指令图例X1=2,2,1;Y1=0,2,1; Z1=0,0,1;fill3(X1,Y1,Z1,r);hold on;X2=1,0,0;Y2=1,2,0;Z2=1,0,0;fill3(X2,Y2,Z2,r);X3=0,2,1;Y3=2,2,1;Z3=0,0,1;fill3(X3,Y3,Z3,b);text(1,1,1,1,1,1);3 曲面图形3.1 网格点坐标的表示x=1:2:7y=2:2:6X,Y = meshgrid(x,y)X = 1 3 5 7 1 3 5 7 1 3 5 7Y = 2 2 2 2 4 4 4 4 6 6 6 63.2 三维网格mesh、meshc、m

38、eshz 用途：数据场的观察分析命令图例随机数据的网格Z=rand(5,5);mesh(Z);% 设置颜色colormap(1,0,0);自定义函数的网格x=-4: 1: 4;y=-5: 1: 5;X,Y=meshgrid(x,y);Z=X.2+Y.2;mesh(X,Y,Z); 消影开关：hidden on / hidden off 利用peaks(50)作为模拟数据矩阵;命令图例 带等高线的网格Z=peaks(50);meshc(Z);Z=peaks(50);meshc(Z);colormap(1,0,0);带基准面的网格Z=peaks(50);meshz(Z);剪孔Z=peaks(50);

39、Z(30:45,15:30)=NaN*ones(16,16);meshc(Z);3.3 着色表面图surf、surfc命令图例表面着色的网格Z=peaks(50);surf(Z);自定义函数的着色网格x=-2: 0.1: 2;y=-2: 0.1: 2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);surfc(X,Y,Z);3.4 二元函数的伪彩色图pcolor 用途：污染浓度场的观察分析。命令图例Z=peaks(50); pcolor(Z);colorbar(hor);colorbar(vec);3.5 等高线contour不仅可用于绘图，更可以用以求截面数据。命令图例以

40、矩阵下标为x、y分量的等高线Z=peaks(50);C=contour(Z);colormap(1,0,0);C：保存了全部等高线上的点坐标。均分n条等高线，并标注之Z=peaks(50);n=5;C=contour(Z,n);colormap(1,0,0);clabel(C);在指定高度绘制等高线Z=peaks(50);V=-10: 2: 10;C=contour(Z, V);colormap(1,0,0);clabel(C);完整图形数据的等高线x=-2: 0.1: 2;y=-2: 0.1: 2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);n=5;C=contour

41、(X,Y,Z,n);三维等高线x=-2:0.1:2;y=-2:0.1:2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);n=10;C=contour3(X,Y,Z,n); 3.6 矢量场图quiver用于挖掘数据变化趋势。命令图例构造起伏跌宕的曲面x=-2:0.2:2;y=-1:0.2:1;X,Y=meshgrid(x,y);Z=X.*exp(-X.2-Y.2);mesh(X,Y,Z);xlabel(X轴);ylabel(Y轴);colormap(1,0,0);计算曲面的梯度px,py= gradient(Z,0.2,0.2);绘制矢量场图quiver(x,y,px,py

42、);3.7 视角控制view视点控制方式及效果：view(1 1 1)view(2 1 1)view(3 1 1)view(1 1 1)view(1 2 1)view(1 3 1)view(1 1 1)view(1 1 2)view(1 1 3)方位角、仰角控制方式及效果：缺省为(-37.5,30)。view(-37.5,30)view(-17.5,30)view(-5.5,30)view(-37.5,30)view(-37.5,45)view(-37.5,60)3.8 多视区控制subplotsubplot(2,1,1); mesh(X,Y,Z);subplot(2,1,2); quiver(x,y,px,py);3.9 制作、播放动画x,y,z=peaks(30); surf(x,y,z)% 制作动画m=moviein(5); % 5帧画面的动画变量for i=1:5 view(1 1 i) % 不断调整视点 m(:,i)=getframe; % 将当前画面作为帧保存到序列中end% 播放动画movi