基于matlab-simulink的buck变换器的研究与设计电力电子课程设计

1、电力电子课程设计2012.5BUCK变换器的研究与设计基于MATLAB/SIMULINK指导老师:王武09级电信学院目录 TOC o 1-3 h z u HYPERLINK l _Toc327288250 课程设计任务书 课程设计任务书学生姓名： 路长鑫 专业班级： 09级电气 3班 指导教师： 王武 工作单位： 电气工程及其自动化 题目：BUCK变换器的研究与设计初始条件输入电压：2030V，输出电压：0-15V，输出负载电流：0.11A，工作频率：30KHz，采用降压斩波主电路。要求完成的主要任务1. 直流供电电源设计。2. 降压斩波主电路设计（包括电路结构形式，全控型器件的选择）并讨论主

2、电路的工作原理。3. 脉宽调制电路（如SG3525集成PWM控制器）及驱动电路设计。4. 分析PWM控制原理及波形。5. 提供电路图纸至少一张。课程设计说明书应严格按统一格式打印，资料齐全，坚决杜绝抄袭，雷同现象。应画出单元电路图和整体电路原理图，给出系统参数计算过程，图纸、元器件符号及文字符号符合国家标准。时间安排2012.3.262012.4.8 收集资料，确定设计方案2012.4.92012.5.27 系统设计2012.5.282012.6.10 撰写课程设计论文及提交指导教师签名： 年 月 日 成员分工成员分工情况如下表： 分工及完成情况成员名单分工完成情况李雪灵、张萌直流供电电源设计

3、完成李海龙 、杨浩琼、董卓奇降压斩波主电路设计完成郭道民、刘晓东脉宽调制电路完成路长鑫硬件设计及仿真完成前 言概述直流斩波电路（DC Chopper）的功能是将直流电变为另一固定电压或可调电压的直流电，也称为直接直流-直流变换器（DC/DC Converter）。直流斩波电路一般是指直接将直流电变为另一直流电的情况，不包括直流-交流-直流的情况。习惯上，DC-DC变换器包括以上两种情况。直流斩波电路的种类较多，包括6种基本斩波电路：降压斩波电路，升压斩波电路，升降压斩波电路，Cuk斩波电路，Sepic斩波电路和Zeta斩波电路，其中前两种是最基本的电路。一方面，这两种电路应用最为广泛，另一方面

4、，理解了这两种电路可为理解其他的电路打下基础。利用不同的基本斩波电路进行组合，可构成复合斩波电路，如电流可逆斩波电路、桥式可逆斩波电路等。利用相同结构的基本斩波电路进行组合，可构成多相多重斩波电路。直流斩波电路广泛应用于直流传动和开关电源领域，是电力电子领域的热点。全控型器件选择绝缘栅双极晶体管（IGBT）综合了GTR和电力MOSFET的优点，具有良好的特性。目前已取代了原来GTR和一部分电力MOSFET的市场，应用领域迅速扩展，成为中小功率电力电子设备的主导器件。所以，此课程设计选题为：设计使用全控型器件为IGBT的降压斩波电路。主要讨论电源电路、降压斩波主电路、控制电路、驱动电路和保护电路

5、的原理与设计。研究的意义从八十年代末起，工程师们为了缩小DC/DC变换器的体积，提高功率密度，首先从大幅度提高开关电源的工作频率做起，但这种努力结果是大幅度缩小了体积，却降低了效率。发热增多，体积缩小，难过高温关。因为当时MOSFET的开关速度还不够快，大幅提高频率使MOSFET的开关损耗驱动损耗大幅度增加。工程师们开始研究各种避开开关损耗的软开关技术。虽然技术模式百花齐放，然而从工程实用角度仅有两项是开发成功且一直延续到现在。一项是VICOR公司的有源箝位ZVS软开关技术；另一项就是九十年代初诞生的全桥移相ZVS软开关技术。有源箝位技术历经三代，且都申报了专利。第一代系美国VICOR公司的有

6、源箝位ZVS技术，其专利已经于2002年2月到期。VICOR公司利用该技术，配合磁元件，将DC/DC的工作频率提高到1MHZ，功率密度接近200W/in3，然而其转换效率却始终没有超过90%，主要原因在于MOSFET的损耗不仅有开关损耗，还有导通损耗和驱动损耗。特别是驱动损耗随工作频率的上升也大幅度增加，而且因1MHZ频率之下不易采用同步整流技术，其效率是无法再提高的。因此，其转换效率始终没有突破90%大关。为了降低第一代有源箝位技术的成本，IPD公司申报了第二代有源箝位技术专利。它采用P沟MOSFET在变压器二次侧用于forward电路拓朴的有源箝位。这使产品成本减低很多。但这种方法形成的M

7、OSFET的零电压开关（ZVS）边界条件较窄，在全工作条件范围内效率的提升不如第一代有源箝位技术，而且PMOS工作频率也不理想。为了让磁能在磁芯复位时不白白消耗掉，一位美籍华人工程师于2001年申请了第三代有源箝位技术专利，并获准。其特点是在第二代有源箝位的基础上将磁芯复位时释放出的能量转送至负载。所以实现了更高的转换效率。它共有三个电路方案：其中一个方案可以采用N沟MOSFET。因而工作频率较高，采用该技术可以将ZVS软开关、同步整流技术、磁能转换都结合在一起，因而它实现了高达92%的效率及250W/in3以上的功率密度。MATLABMATLAB是矩阵实验室（Matrix Laborator

8、y）的简称，是美国MathWorks公司出品的商业HYPERLINK /view/1658637.htm t _blank数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，SIMULINK是MATLAB软件的扩展，它是实现动态系统建模和仿真的一个软件包，本课程设计的仿真即需要在SIMULINK中来完成电路的仿真与计算。通过系统建模和仿真，掌握和运用MATLAB/SIMULINK工具分析系统的基本方法。致谢诚挚的感谢老师的细心指导，诚挚的感谢作出努力的每一位组员，使得本学期电力电子课程设计的顺利完成，使我更进一步的明白了团队力量的内涵，为以后的工作和学习奠定

9、了基础，树立了标杆！第1章 设计思路与框图设计思路本课程设计主要应用了MATLAB 软件及其组件之一SIMULINK进行系统的设计与仿真。系统主要包括：直流稳压电源部分、BUCK降压斩波主电路部分、PWM控制部分和负载。直流稳压电源部分采用SIMULINK中的交流电源模块，模拟工频50Hz的220V交流电源，经过变压器降压，把220V电压变为30V左右的交流电，为主电路提供电源。BUCK降压斩波主电路部分拖动带反电动势的电阻负载，模拟现实中一般的负载，若实际负载中没有反电动势，只需令其为零即可。PWM控制部分为主电路部分提供脉冲信号，控制全控器件IGBT的导通和关断，实现整个系统的运行。在SI

10、MULINK中完成各个功能模块的绘制后，即可进行仿真和调试，用SIMULINK提供的示波器观察波形，进行相应的电压和电流等的计算，最后进行总结，完成整个BUCK变换器的研究与设计。系统框图系统框图如图1.1所示。图1.1 BUCK变换器系统结构总框图第2章 直流稳压电源设计电源设计原理设计电路采用的是电容滤波的单项不可控整流电路，小功率稳压电源由电源变压器、整流电路、滤波电路和稳压电路四个部分组成，其原理框图如图2.1,波形图如图2.2所示。 图2.1 直流稳压电源原理框图图2.2 直流稳压波形图 电路的工作原理及其波形分析a) 电路 b) 波形图2.3 电容滤波的单相桥式不可控整流电路及其工

11、作波形基本工作过程在u2正半周过零点至wt=0期间，因u2 50V,额定工作电流ID=1； 滤波电容：容量470F,耐压50V 第3章 降压斩波电路课程设计降压斩波电路主电路图图3.1降压斩波电路的原理图及波形a）电路图 b）电流连续时的波形c）电流断续时的波形电路分析直流斩波电路(即斩波器)的功能是将直流电变换为另一固定电压或可调电压的直流电。本课程设计使用一个全控型器件VIGBT，若采用晶闸管，需设置使晶闸管关断的辅助电路。另外设置了续流二极管VD，在V关断时给负载中电感电流提供通道。该电路主要用于电子电路的供电电源，也可拖动直流电动机或带蓄电池负载等，后两种情况下负载中均会出现反电动势，

12、如图中Em所示。IGBT简介IGBT的等效电路如图3.2所示。由图可知,若在IGBT的栅极和发射极之间加上驱动正电压,则MOSFET导通,这样PNP晶体管的集电极与基极之间成低阻状态而使得晶体管导通;若IGBT的栅极和发射极之间电压为0V,则MOSFET截止,切断PNP晶体管基极电流的供给,使得晶体管截止。图3.2 IGBT的等效电路图由此可知,IGBT的安全可靠与否主要由以下因素决定：IGBT栅极与发射极之间的电压;IGBT集电极与发射极之间的电压;流过IGBT集电极发射极的电流;IGBT的结温。如果IGBT栅极与发射极之间的电压,即驱动电压过低,则IGBT不能稳定正常地工作,如果过高超过栅

13、极发射极之间的耐压则IGBT可能永久性损坏;同样,如果加在IGBT集电极与发射极允许的电压超过集电极发射极之间的耐压,流过IGBT集电极发射极的电流超过集电极发射极允许的最大电流,IGBT的结温超过其结温的允许值,IGBT都可能会永久性损坏。工作原理 t=0时刻驱动V导通，电源E向负载供电，负载电压uo=E，负载电流io按指数曲线上升。 t=t1时控制V关断，二极管VD续流，负载电压uo近似为零，负载电流呈指数曲线下降，通常串接较大电感L使负载电流连续且脉动小。基本的数量关系如下:电流连续时,负载电压的平均值为:式中，ton为V处于通态的时间，toff为V处于断态的时间，T为开关周期，a为导通

14、占空比，简称占空比或导通比。 负载电流平均值为: 电流断续时,负载电压uo平均值会被抬高，一般不希望出现电流断续的情况。控制方式根据对输出电压平均值进行调制的方式不同,斩波电路有三种控制方式(时间比控制方式):脉冲宽度调制（PWM）：T不变，改变ton。（定频调宽控制模式）频率调制：ton不变，改变T。 （定宽调频控制模式）混合型：ton和T都可调，改变占空比 （调宽调频混合控制模式）对降压斩波电路进行解析 基于分时段线性电路这一思想，按V处于通态和处于断态两个过程来分析，初始条件分电流连续和断续。电流连续时 (1) （2） 用泰勒级数近似，可得： （3）平波电抗器L为无穷大，此时负载电流最大

15、值、最小值均等于平均值。（3）式所示的关系还可从能量传递关系简单地推得，一个周期中，忽略电路中的损耗，则电源提供的能量与负载消耗的能量相等，即： 则：假设电源电流平均值为I1，则有： 其值小于等于负载电流Io，由上式得： 即输出功率等于输入功率，可将降压斩波器看作直流降压变压器。电流断续时有I10=0，且t=ton+tx时，i2=0，可以得出： 电流断续时，txElectrical Sources中选择“AC Voltage Sources”交流电源模块，用于模拟实际的交流电。Simulink中模块如图5.1所示，双击模块，打开参数设置对话框，这里，参数设置如图5.2所示。反电动势在SimPo

16、werSystemElectrical Sources选择“DC Voltage Sources”直流电原模块，用于模拟负载为直流电动机或蓄电池时出现的反电势，若负载无反电势时只需令其为零即可，模块如图5.3所示。参数设置：Amplitude（V）：1 图5.1 交流电压模块 图5.3 直流电模块 图5.2 交流电参数设置 图5.5 变压器参数设置变压器在SimPowerSystemElements选择“Linear Transformer”变压器模块，用于将工频50Hz的220V较高电压的交流电，变为适合整理电路用较低电压的交流电。Linear Transformer模块如图5.4所示，参数

17、设置如图5.5所示。 图5.4 变压器模块 图5.6 串联RLC模块串联阻感容模块在SimPowerSystemElements选择“Series RLC Branch”阻感容串联模块。在Simulink模块库中没有专用的电阻、电感、电容模块，它们均可以通过Series RLC Branch模块通过参数的设置来实现，模块如图5.6所示。参数设置：电阻RBranch type：RResistance（0hms）：15电感L1Branch type：LInductance（H）：1e-10电感L2Branch type：LInductance（H）：2电容CBranch type：CCapacit

18、ance（F）：470e-6电力二极管模块在SimPowerSystemPower Electronics选择“Diode”电力二极管模块，模块如图5.7所示，参数设置如图5.8所。 图5.7 电力二极管模块 图5.9 IGBT模块 图5.8电力二极管参数 图5.10 IGBT参数全控器件IGBT模块在SimPowerSystemPower Electronics选择全控型“IGBT”模块，模块如图5.9所示，参数设置如图5.10所。电压电流测量模块在SimPowerSystemMeasurements选择“Voltage Measurement”电压测量模块和“Current Measure

19、ment”电流测量模块，通过这些模块，可以方便的与示波器模块相连接来进行参数的测量。模块如图5.11所示。图5.11 电压和电流测量模块示波器模块在SimulinkSinks选择“Scope”示波器模块，用来与电压和电流测量模块配合使用，显示测量点的电压或电流波形。“Scope”示波器模块可以参数设置测量输入端的数目，也就是说可以同时进行多路的测量，既可以是电压，也可以是电流，仿真时可以通过双击示波器模块，打开显示波形的界面，该界面有很多按钮，可以进行X轴和Y轴的放大显示，方便观察测量的波形，选择“Parameters”按钮模块，打开示波器的属性设置窗口，在Number of axes中输入需

20、要的端口数目即可。模块如图5.12所示。 图5.12 示波器模块 图5.13 Powergui模块Powergui模块在SimPowerSystem选择“Powergui”模块，模块如图5.13所示。Powergui 是一个环境模块，是一个接口工具，用来采集信号，连接起电路图的各个元件，任何一个含有 SimPowerSystem 模块的模型中必须含有一个。它储存了电路模型的等效数学模型（状态空间方程）。没有它，仿真不能启动，会给出一个错误提示。模块如图5.13所示。参数设置：在参数设置对话框中，选择Simulation and configuration optionconfigure par

21、ameters选项，弹出的“Solver”选项卡后进行相关设置。Simulation type：DiscreteSample time：1e-6PWM脉冲模块在SimulinkSources选择“Pulse Generator”模块，用来模拟PWM控制电路和驱动电路，该模块通过参数的设置，可以实现任意周期，任意宽度，任意幅值的脉冲信号，模块如图5.14所示。参数设置如图5.15所示。 图5.14 脉冲产生模块 图5.15 脉冲产生模块参数设置仿真电路设计在Simulink中选择FileNewModel，即可创建以一个由工具栏和绘图区构成的“*.mdl”文件，将选择的各个模块从库中拖到新建的绘图

22、区，进行连线，即可完成电路图的绘制。电路图如图5.16所示。 图5.16 仿真电路图仿真波形设置仿真时间为1s，仿真过程中或仿真结束后，双击示波器模块，即可查看各个测量点的波形，如图5.17所示。图 5.17 各测量点的电压和电流波形仿真分析由仿真输出的波形,分析如下:输入侧为标准工频50Hz的220V交流电源，变压器变比为22:3，可知，变压器二次侧输出为30V的交流电，由整流输出电压波形可以看出。初期电压峰值约为30V，波形很快趋于稳定约为E=23V，位于20-30V之间，符合设计要求。对主电路而言，输入侧即为整流输出的E=23V直流电，由“Pulse Generator”模块参数设置可知

23、脉冲周期T=3.33e-5s，ton=50%T，由电流连续时,负载电压的平均值计算公式:可得U0=1/2E=11.5V。由负载电压波形可以看出，电压上升时间约为0.4s，上升到11V左右后趋于稳定，近乎直线，与计算值接近，可见输出波形较为理想，且位于0-15V之间，符合设计要求。由变压器二次侧电流波形与图2.3 b对比，波形基本符合。由整流输出电流波形可以看出，电流脉动频率很高，主要是由于IGBT全控开关管以30KHz的频率导通和关闭，使得电流也随之导通和关断，从而产生了很高的脉动频率，符合设计要求。反电动势Em=1V，负载电阻R=15，由负载电流平均值计算公式: 可得I0=(11.5-1)/

24、15=0.7V。负载电流波形可以看出，电流上升时间约为0.4s，上升到0.7A后趋于稳定，近乎直线，与计算值接近，可见输出波形较为理想，且位于0-1A之间，符合设计要求。 第6章 设计总结设计回顾通过这个学期的电力电子课程设计，初步掌握了如何在MATLAB/SIMULINK环境中进行电路图的绘制和仿真,基本学会了如何进行各个参数的设置和计算,了解了一些故障的解决办法，为以后的工作和学习奠定了一定的基础。心得体会通过合作，我们的合作意识得到加强。合作能力得到提高。上大学后，很多同学都没有过深入的交流，在设计的过程中，我们用了分工与合作的方式，每个人互责一定的部分，同时在一定的阶段共同讨论，以解决

25、分工中个人不能解决的问题，在交流中大家积极发言，和提出意见，同时我们还向别的同学请教。在此过程中，每个人都想自己的方案得到实现，积极向同学说明自己的想法。能过比较选出最好的方案。在这过程也提高了我们的表过能力。 在设计的过程中我们还得到了老师的帮助与意见。在学习的过程中，不是每一个问题都能自己解决，向老师请教或向同学讨论是一个很好的方法.参考文献【1】王兆安 黄俊.电力电子技术.北京:机械工业出版社, 2000【2】浣喜明 姚为正.电力电子技术.高等教育出版社，2000【3】莫正康电力电子技术应用（第3版）.机械工业出版社，2000【4】郑琼林 耿学文.电力电子电路精选.机械工业出版社，199

26、6【5】刘定建 朱丹霞.实用晶闸管电路大全.机械工业出版社，1996【6】刘祖润 胡俊达.毕业设计指导.机械工业出版社，1995【7】邵玉斌.Matlab/Simulink通信系统建模与仿真分析. HYPERLINK /Search/?publisher=清华大学出版社&page=1 o 搜索 清华大学出版社 出版的更多图书 t _blank 清华大学出版社，2008【8】张圣勤.MATLAB 7.0实用教程.机械工业出版社，2006附录资料：MATLAB函数和命令的用法Binocdf二项式累积分布函数语法格式Y = binocdf(X,N,P)函数功能Y = binocdf(X,N,P) 计

27、算X中每个X(i)的二项式累积分布函数，其中，N中对应的N(i)为试验数，P中对应的P(i)为每次试验成功的概率。Y, N, 和 P 的大小类型相同，可以是向量、矩阵或多维数组。输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。The values in N must all be positive integers, the values in X must lie on the interval 0,N, and the values in P must lie on the interval 0, 1.The binomial cdf for a given value x an

28、d a given pair of parameters n and p is The result, y, is the probability of observing up to x successes in n independent trials, where the probability of success in any given trial is p. The indicator function I(0,1,.,n)(i)ensures that x only adopts values of 0,1,.,n.示例若一个棒球队在一个赛季要比赛162场，每场比赛取胜的机会是

29、50-50，则该队取胜超过100 场的概率为： 1-binocdf(100,162,0.5)ans = 0.0010433相关函数 HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binofit.html binofit | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binoinv.html binoinv | HYPERLINK jar:file:/

30、C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binopdf.html binopdf | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binornd.html binornd | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binostat.html binostat

31、 | HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/cdf.html cdf附：二项式分布(binomial distribution )定义二项分布的概率密度函数为where k is the number of successes in n trials of a Bernoulli process with probability of success p.The binomial distribution is discrete, defined for integ

32、ers k = 0, 1, 2, . n, where it is nonzero.背景The binomial distribution models the total number of successes in repeated trials from an infinite population under the following conditions:Only two outcomes are possible on each of ntrials.The probability of success for each trial is constant.All trials

33、are independent of each other.The binomial distribution is a generalization of the HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/brn2ivz-2.html Bernoulli distribution; it generalizes to the HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/

34、help.jar%21/brn2ivz-84.html multinomial distribution.参数Suppose you are collecting data from a widget manufacturing process, and you record the number of widgets within specification in each batch of100. You might be interested in the probability that an individual widget is within specification. Par

35、ameter estimation is the process of determining the parameter, p, of the binomial distribution that fits this data best in some sense.One popular criterion of goodness is to maximize the likelihood function. The likelihood has the same form as the binomial pdf above. But for the pdf, the parameters

36、(nandp) are known constants and the variable isx. The likelihood function reverses the roles of the variables. Here, the sample values (the xs) are already observed. So they are the fixed constants. The variables are the unknown parameters. MLE involves calculating the value of p that give the highe

37、st likelihood given the particular set of data.The function HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/binofit.html binofit returns the MLEs and confidence intervals for the parameters of the binomial distribution. Here is an example using random numbers from

38、 the binomial distribution with n=100 and p=0.9. r = binornd(100,0.9)r = 85 phat, pci = binofit(r,100)phat = 0.85pci = 0.76469 0.91355The MLE for parameterp is0.8800, compared to the true value of0.9. The 95% confidence interval forp goes from 0.7998 to0.9364, which includes the true value. In this

39、made-up example you know the true value ofp. In experimentation you do not.示例The following commands generate a plot of the binomial pdf for n = 10 and p = 1/2.x = 0:10;y = binopdf(x,10,0.5);plot(x,y,+) 相关内容 HYPERLINK jar:file:/C:/Program%20Files/MATLAB/R2010a/help/toolbox/stats/help.jar%21/bqt29ct.h

40、tml l bqt290a-1 Discrete Distributions附：二项式分布（网上）定义若某事件概率为p，现重复试验n次，该事件发生k次的概率为:P=C(k,n)pk(1-p)(n-k)C(k,n)表示组合数，即从n个事物中拿出k个的方法数。二项分布的概念考虑只有两种可能结果的随机试验，当成功的概率（）是恒定的，且各次试验相互独立，这种试验在统计学上称为贝努里试验（Bernoullitrial）。如果进行n次贝努里试验，取得成功次数为X（X=0,1,n）的概率可用下面的二项分布概率公式来描述：P=C(X,n) X(1-)(n-X)式中的n为独立的贝努里试验次数，为成功的概率，（1

41、-）为失败的 HYPERLINK /doc/5399415.html t _blank 概率，X为在n次贝努里试验中出现成功的次数，C(X,n)表示在n次试验中出现X的各种组合情况，在此称为二项系数（binomialcoefficient）。内容简介二项分布，伯努里分布：进行一系列试验,如果1. 在每次试验中只有两种可能的结果,而且是互相对立的；2. 每次实验是独立的，与其它各次试验结果无关；3. 结果事件发生的 HYPERLINK /doc/5399415.html t _blank 概率在整个系列试验中保持不变，则这一系列试验称为伯努力试验。在这试验中，事件发生的次数为一 HYPERLIN

42、K /doc/2273201.html t _blank 随机事件,它服从二次分布。二项分布可以用于可靠性试验。可靠性试验常常是投入n个相同的式样进行试验T小时，而只允许k个式样失败，应用二项分布可以得到通过试验的概率。一个事件必然出现，就说它100%要出现。100%=1，所以100%出现的含义就是出现的概率P=1。即必然事件的出现概率为1。若掷一枚硬币，正面向上的结果的概率为0.5。反面向上的结果的概率也是0.5。则出现正面向上事件或者反面向上事件的概率就是0.5+0.5=1，即二者必居其一。若掷两次硬币，由独立事件的概率乘法定理那么两次都是正面（反面）向上的概率是0.50.5=0.25。另

43、外第一个是正第二个是反的出现概率也是0.50.5=0.25。同理第一个反第二个正的出现概率也是0.50.5=0.25。于是一正一反的概率是前面两个情况的和，即0.25+0.25=20.25=0.5。它们的合计值仍然是1。binopdf二项分布的概率密度函数语法格式Y = binopdf(X,N,P)函数功能Y = binopdf(X,N,P) 计算X中每个X(i)的概率密度函数，其中，N中对应的N(i)为试验数，P中对应的P(i)为每次试验成功的概率。Y, N, 和 P 的大小类型相同，可以是向量、矩阵或多维数组。输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。N中的值必须是正整数，

44、0 P 1 。对于给出的x和两个参数n和p，二项分布概率密度函数为其中 q = 1 p。 y为n次独立试验中成功x次的概率，其中，每次试验成功的概率为p。指示器函数 I(0,1,.,n)(x) 确保x 取值为 0, 1, ., n。 示例一个质量检查技术员一天能测试200块电路板。假设有 2% 的电路板有缺陷，该技术员在一天的测试中没有发现有缺陷的电路板的概率是多少？binopdf(0,200,0.02)ans =0.0176质量检查技术员一天中最大可能检测出有缺陷的电路板是多少块？defects=0:200;y = binopdf(defects,200,.02);x,i=max(y);de

45、fects(i) ans =4相关函数 HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binocdf.html binocdf | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binofit.html binofit | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/bin

46、oinv.html binoinv | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binornd.html binornd | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/binostat.html binostat | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/toolbox/stats/help.jar%21/p

47、df.html pdfdlmread将以ASCII码分隔的数值数据文件读入到矩阵中语法格式M = dlmread(filename)M = dlmread(filename, delimiter)M = dlmread(filename, delimiter, R, C)M = dlmread(filename, delimiter, range)函数功能M = dlmread(filename)读取有分隔符的ASCII数值数据文件filename，并把数据给到矩阵M中。文件名(filename)以字符串形式用单引号括起来。 dlmread 从文件的格式中推断分隔符。M = dlmread(f

48、ilename, delimiter)指定分隔符delimiter。使用 t 代表制表符 tab 作为分隔。M = dlmread(filename, delimiter, R, C)R和C指定了数据读取数据，其左上角的值位于文件中的第R行，第C列。R和C从0开始，所以R=0, C=0 指向文件中的第一个值。M = dlmread(filename, delimiter, range)读取由range = R1 C1 R2 C2指定的区域块， (R1,C1)是左上角，(R2,C2)是右下角。也可以使用电子表格表示法来指定，如range = A1.B7。备注All data in the inp

49、ut file must be numeric. dlmread does not read files that contain nonnumeric data, even if the specified rows and columns contain only numeric data.若没有指定分隔符，当从文件格式中推断分隔符时，连续的空格符当作一个分隔符对待。若指定了分隔符，则重复的分隔符将分别作为单独的分隔符对待。If you want to specify an R, C, or range input, but not a delimiter, set the delimit

50、er argument to the empty string, (two consecutive single quotes with no spaces in between, ). For example,M = dlmread(myfile.dat, , 5, 2)In this case, dlmread treats repeated white spaces as a single delimiter.Dlmread将用0填充没有边界的区域。有多行的数据文件，若以非空格分隔符结束，例如分号，则导入后会多产生全0的一列在最后。Dlmread在导入任何复数时，将其作为一个整体导入到一

51、个复数单元中。下面是有效的复数格式：i|jExample: 5.7-3.1ii|jExample: -7j嵌入了空格的复数是不正确的格式，空格将被认为是分隔符。示例例1Export a 5-by-8 test matrix M to a file, and read it with dlmread, first with no arguments other than the filename:M = gallery(integerdata, 100, 5 8, 0); dlmwrite(myfile.txt, M, delimiter, t)dlmread(myfile.txt)ans =

52、96 77 62 41 6 21 2 42 24 46 80 94 36 20 75 85 61 2 93 92 82 61 45 53 49 83 74 42 1 28 94 21 90 45 18 90 14 20 47 68Now read a portion of the matrix by specifying the row and column of the upper left corner:dlmread(myfile.txt, t, 2, 3)ans = 92 82 61 45 53 42 1 28 94 21 90 14 20 47 68This time, read a

53、 different part of the matrix using a range specifier:dlmread(myfile.txt, t, C1.G4)ans = 62 41 6 21 2 80 94 36 20 75 93 92 82 61 45 74 42 1 28 94例2Export matrix M to a file, and then append an additional matrix to the file that is offset one row below the first:M = magic(3);dlmwrite(myfile.txt, M*5

54、M/5, )dlmwrite(myfile.txt, M/3, -append, . roffset, 1, delimiter, )type myfile.txt40 5 30 1.6 0.2 1.215 25 35 0.6 1 1.420 45 10 0.8 1.8 0.4 2.6667 0.33333 21 1.6667 2.33331.3333 3 0.66667When dlmread imports these two matrices from the file, it pads the smaller matrix with zeros:dlmread(myfile.txt)

55、40.0000 5.0000 30.0000 1.6000 0.2000 1.2000 15.0000 25.0000 35.0000 0.6000 1.0000 1.4000 20.0000 45.0000 10.0000 0.8000 1.8000 0.4000 2.6667 0.3333 2.0000 0 0 0 1.0000 1.6667 2.3333 0 0 0 1.3333 3.0000 0.6667 0 0 0替代作为dlmread的替代，可使用导入向导。选择菜单 “File | Import Data” 激活导入向导。相关函数 HYPERLINK jar:file:/D:/Pr

56、ogram%20Files/MATLAB/help/techdoc/help.jar%21/ref/dlmwrite.html dlmwrite | HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/techdoc/help.jar%21/ref/textscan.html textscanfprintf写数据到文本文件或输出到命令窗口语法格式fprintf(fileID, format, A, .)fprintf(format, A, .)count = fprintf(.)函数功能fprintf(fileID, format, A, .)

57、 fileID为文件句柄，指定要写入的文件；format是用来控制所写数据格式的格式符；A是用来存放数据的矩阵。applies the format to all elements of array A and any additional array arguments in column order, and writes the data to a text file. fprintf uses the encoding scheme specified in the call to HYPERLINK jar:file:/D:/Program%20Files/MATLAB/help/t

58、echdoc/help.jar%21/ref/fopen.html fopen.fprintf(format, A, .) 将A中的数据按格式format输出到命令窗口。count = fprintf(.) 返回fprintf写的字节数。输入参数fileID取值为下面中的一种：从文件打开时所得到的文件句柄，为一个整数值。1 为标准输出设备（屏幕）。2 for standard error.fileID缺省时，输出到屏幕，即1为默认值。format单引号括起来的是 “转换控制字符串 ”，指定输出的格式。由下面给出的各部分的组合而成：百分号后紧跟的格式字符，如%s 。输出区域的宽度，精度，和其它选

59、项（附加格式说明符）。要原样输出的文本（普通字符）。转义字符，包括：单引号% 百分号反斜线a响铃b退格符fForm feedn换行r回车t水平制表符v垂直制表符xN十六进制数, NN八进制数, N转换字符和可选的操作符按以下顺序排列（包括空格）：下表列出可用的转换字符和限定类型字符。数据类型转换说 明Integer, signed%d 或 %i十进制整数%ld 或 %li64-bit base 10 values%hd 或 %hi16-bit base 10 valuesInteger, unsigned%u十进制整数%o八进制整数(octal)%x十六进制整数 (hexadecimal),

60、使用字母 af%X十六进制整数 (hexadecimal), 使用字母AF%lu%lo%lx or %lX64-bit values, 十、八、十六进制整数%hu%ho%hx or %hX16-bit values, 十、八、十六进制整数Floating-point number%f定点表示法%e指数表示法，如 3.141593e+00%E同 %e，但用E，如 3.141593E+00%g%e 或 %f的紧凑格式，没有尾部的0%G%E 或 %f的紧凑格式，没有尾部的0%bx or %bX%bo%bu双精度的十六、八、十进制数值如：%bx prints pi as 400921fb54442d1