计算机仿真试卷

桂林电子科技大学研究生试卷

一、填空（每空1分，共20分）。

1、仿真是指用代替进行实验和研究。

2、为使仿真的结果被证实是真实可靠的，也就是结果是可信的，仿真所遵循的基本原则是相似原理：、和。

3、按照仿真实验时间标尺τ与实际系统的时间标尺t比例来分类，将的仿真称为实时仿真，而τ/t≠1的仿真为非实时仿真。一般有介入的半实物仿真属于实时仿真，而纯计算机仿真为非实时仿真。

4、凡是系统的输入量是的连续函数，无论其输出量是连续单调函数或是单值或是多值的函数，均称这个系统是连续系统。用表示的系统可以是非线性或线性系统，而对于表示的系统，只适用于单输入－单输出的线性定常系统，所以表示有一定的局限性。

5、各种数值积分法的截断误差与积分方法的和有关。是由于的字长是有限制的，计算只能限制在有限位数内，因而产生了舍入误差。

6、对于一个稳定的微分方程，经过转换得到的差分方程的计算结果是稳定的，不同的数值积分方法的计算稳定性相同。

7、按环节离散化就是将系统分成若干个环节，在每个环节的入口处加一个虚拟的，并立即跟一个信号重构过程（加），以便使信号恢复为形式，这时系统实际上已成为一个系统。当采样周期足够小时，这个采样系统就近似等价于原系统。

二、简述系统仿真的特点（8分）

三、线性定常连续系统的数学模型有哪几种形式，分别叙述：

1、每种形式数学模型的模型表示；（3分）

2、各种数学模型所表示的系统特征是什么？（6分）

3、模型之间的相互转化关系是什么？（3分）

四、在采样控制系统的数字仿真中，

1、画出典型的采样控制系统图；（4分）

2、采样控制系统工作的特点；（5分）

3、采样控制系统仿真前应处理的问题；（5分）

T之间的关系。（6分）

4、如何确定计算步长T和采样控制系统自己固定的采样周期

s

五、应用MATLAB中的模型表示：

1、分别阐述函数ss()、tf()和zpk()的功能；（5分）

2、写出sys=tf({-5;[1 -5 6]},{[1 -1];[1 1 0]})的运行结果；（5分）

3、写出sys1=ss(sys)的运行结果；（5分）

4、写出sys=zpk({[];[2 3]},{1;[0 -1]},[-5;1]) 的运行结果；（5分）

六、用图形表示连续调制波形)

t

y 。

sin(t

)

9

sin(

1、写出应用四种不同采样频率的MATLAB实现的程序；（10分）

2、调试、写出正确的运行结果，并对各结果作出比较分析。（10分）

一、填空

1、仿真是指用模型（物理模型或数学模型）代替实际系统进行实验和研究。

2、为使仿真的结果被证实是真实可靠的，也就是结果是可信的，仿真所遵循的基本原则是相似原理：几何相似、环境相似和性能相似。

3、按照仿真实验时间标尺τ与实际系统的时间标尺t比例来分类，将τ/t=1的仿真称为实时仿真，而τ/t≠1的仿真为非实时仿真。一般有实物介入的半实物仿真属于实时仿真，而纯计算机仿真为非实时仿真。

4、凡是系统的输入量是时间t的连续函数，无论其输出量是连续单调函数或是单值或是多值的函数，均称这个系统是连续系统。用微分方程表示的系统可以是非线性或线性系统，而对于传递函数表示的系统，只适用于单输入－单输出的线性定常系统，所以表示有一定的局限性。

5、各种数值积分法的截断误差与积分方法的阶次和步长有关。是由于计算机的字长是有限制的，计算只能限制在有限位数内，因而产生了舍入误差。

6、对于一个稳定的微分方程，经过转换得到的差分方程的计算结果不一定是稳定的，不同的数值积分方法的计算稳定性也不相同。

7、按环节离散化就是将系统分成若干个环节，在每个环节的入口处加一个虚拟的采样开关，并立即跟一个信号重构过程（加保持器），以便使信号恢复为连续形式，这时系统实际上已成为一个采样系统。当采样周期足够小时，这个采样系统就近似等价于原连续系统。

二、简述系统仿真的特点

1仿真的实验性质

仿真基本上是一种通过实验来研究系统的综合实验技术，具有一般实验的性质。不论是利用仿真技术进行系统分析还是设计，都必须通过一系统实验来完成。所建立的仿真模型应具有实验的性质，即模型与原型的功能及参数之间应有相似性和对应性2数字仿真与解析法

解析法是应用数学推导、演绎去求解数学模型的方法，而数字仿真是通过计算机在数学模型上进行一系列实验来研究问题的方法。

解析法求解问题得出对问题的通解，而数字仿真的每一次运行只能给出在特定条件的数值解（特解）。用解析法所能求解的问题是有限的，而仿真适应范围广。

3数字仿真与数值求解。数字仿真是指建立系统的数学模型并在计算机上运行和分析的整个过程；数值求解则是仿真所采用的在计算机上求解数学模型的方法。

三、四（略）

五、应用MATLAB中的模型表示：

1、分别阐述函数ss()、tf()和zpk()的功能；（5分）

ss()的功能：生成或将其它模型转换为状态空间模型

tf()的功能：生成或将其它模型转换为传递函数模型

zpk()的功能：生成或将其它模型转换为零极点增益模型

2、写出sys=tf({-5;[1 -5 6]},{[1 -1];[1 1 0]})的运行结果；（5分）

Transfer function from input to output...

-5

#1: -----

s - 1

s^2 - 5 s + 6

#2: -------------

s^2 + s

3、写出sys1=ss(sys)的运行结果；（5分）

a =

x1 x2 x3

x1 1 0 0

x2 0 -1 0

x3 0 1 0

b =

u1

x1 2

x2 4

x3 0

c =

x1 x2 x3

y1 -2.5 0 0

y2 0 -1.5 1.5

d =

u1

y1 0

y2 1

Continuous-time model.

4、写出sys=zpk({[];[2 3]},{1;[0 -1]},[-5;1]) 的运行结果；（5分）

Zero/pole/gain from input to output...

-5

#1: -----

(s-1)

(s-2) (s-3)

#2: -----------

s (s+1)

六、用图形表示连续调制波形)

t

sin(t

y 。

9

sin(

)

1、写出应用四种不同采样频率的MATLAB实现的程序；（10分）

t1=(0:11)/11*pi; % <1> y1=sin(t1).*sin(9*t1);

t2=(0:100)/100*pi; % <3> y2=sin(t2).*sin(9*t2);

subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图 (1)') subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图 (2)') subplot(2,2,3),plot(t1,y1,t1,y1,'r.')

axis([0,pi,-1,1]),title('子图 (3)')

subplot(2,2,4),plot(t2,y2)

axis([0,pi,-1,1]),title('子图 (4)')

2、调试、写出正确的运行结果，并对各结果作出比较分析。（10分）

计算机仿真课程设计报告

、 北京理工大学珠海学院 课程设计任务书 2010 ～2011 学年第 2学期 学生姓名：林泽佳专业班级：08自动化1班指导教师：钟秋海工作部门：信息学院一、课程设计题目 ： 《控制系统建模、分析、设计和仿真》 本课程设计共列出10个同等难度的设计题目，编号为：[0号题]、[1号题]、[2号题]、[3号题]、[4号题]、[5号题]、[6号题]、[7号题]、[8号题]、[9号题]。 学生必须选择与学号尾数相同的题目完成课程设计。例如，学号为8xxxxxxxxx2的学生必须选做[2号题]。 二、课程设计内容 （一）《控制系统建模、分析、设计和仿真》课题设计内容|

！ " [2 有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 （二）《控制系统建模、分析、设计和仿真》课题设计要求及评分标准【共100分】 , 1、求被控对象传递函数G(s)的MATLAB描述。（2分） 2、求被控对象脉冲传递函数G(z)。（4分） 3、转换G(z)为零极点增益模型并按z-1形式排列。（2分） 4、确定误差脉冲传递函数Ge(z)形式,满足单位加速度信号输入时闭环稳态误差为零和实际 闭环系统稳定的要求。（6分） 5、确定闭环脉冲传递函数Gc(z)形式，满足控制器Dy(z)可实现、最少拍和实际闭环系统稳 定的要求。（8分）

6、根据4、5、列写方程组，求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 （12分） 7、求针对单位加速度信号输入的最少拍有波纹控制器Dy(z)并说明Dy(z)的可实现性。 （3分） ！ 8、用程序仿真方法分析加速度信号输入时闭环系统动态性能和稳态性能。（7分） 9、用图形仿真方法(Simulink)分析单位加速度信号输入时闭环系统动态性能和稳态性能。 （8分） 10、确定误差脉冲传递函数Ge(z)形式,满足单位速度信号输入时闭环稳态误差为零和实际 闭环系统稳定的要求。（6分） 11、确定闭环脉冲传递函数Gc(z)形式，满足控制器Dw(z)可实现、无波纹、最少拍和实际 闭环系统稳定的要求。（8分） 12、根据10、11、列写方程组，求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 （12分） 13、求针对单位速度信号输入的最少拍无波纹控制器Dw(z)并说明Dw(z)的可实现性。（3分） 14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。（7分） 15、用图形仿真方法(Simulink)分析单位速度信号输入时闭环系统动态性能和稳态性能。 & （8分） 16、根据8、9、14、15、的分析，说明有波纹和无波纹的差别和物理意义。（4分） 三、进度安排 6月13至6月14：下达课程设计任务书；复习控制理论和计算机仿真知识，收集资料、熟悉仿真工具；确定设计方案和步骤。 6月14至6月16：编程练习，程序设计；仿真调试，图形仿真参数整定；总结整理设计、 仿真结果，撰写课程设计说明书。 6月16至6月17：完成程序仿真调试和图形仿真调试；完成课程设计说明书；课程设计答 辩总结。 [ 四、基本要求

计算机仿真技术的发展概述及认识

计算机仿真技术的发展概述及认识 摘要：随着经济的发展和社会的进步，计算机技术高速发展，使人类社会进入了信息时代，计算机作为后期新秀渗入到人们生活中的每一个领域，给人们的生活带来了前所未有的变化。作为新兴的技术，计算机技术在人类研究的各个领域起到了只管至关重要的作用，帮助人类解决了许多技术难题。在科研领域，计算机技术与仿真技术相结合，形成了计算机仿真技术，作为人们科学研究的一种新型方法，被人们应用到各个领域，用来解决人们用纯数学方法或者现实实验无法解决的问题，对科研领域技术成果的形成有着积极地促进作用。 本文在计算机仿真技术的理论思想基础上，分析了计算机仿真技术产生的基本原因，也就是人们用计算机模拟解决问题的优点所在，讨论了模拟、仿真、实验、计算机仿真之间的联系和区别，介绍了计算机仿真技术的发展历程，并查阅相关资料介绍了计算机仿真技术在不同领域的应用，分析并预测了计算机仿真的未来发展趋势。经过查阅大量数据资料并加以分析对比，这对于初步认识计算机仿真技术具有重要意义。 关键词：计算机仿真；模拟；仿真技术；发展 一、引言 计算机仿真技术是以多种学科和理论为基础，以计算机及其相应的软件为工具，通过虚拟试验的方法来分析和解决问题的一门综合性技术。计算机仿真（模拟）早期称为蒙特卡罗方法，是一门利用随机数实验求解随机问题的方法。其原理可追溯到1773年法国自然学家G.L.L.Buffon为估计圆周率值所进行的物理实验。根据仿真过程中所采用计算机类型的不同，计算机仿真大致经历了模拟机仿真、模拟－数字混合机仿真和数字机仿真三个大的阶段。20世纪50年代计算机仿真主要采用模拟机；60年代后串行处理数字机逐渐应用到仿真之中，但难以满足航天、化工等大规模复杂系统对仿真时限的要求；到了70年代模拟－数字混合机曾一度应用于飞行仿真、卫星仿真和核反应堆仿真等众多高技术研究领域；80年代后由于并行处理技术的发展，数字机才最终成为计算机仿真的主流。现在，计算机仿真技术已经在机械制造、航空航天、交通运输、船舶工程、经济管理、工程建设、军事模拟以及医疗卫生等领域得到了广泛的应用。 二、基本概念 模拟：（Simulation）应用模型和计算机开展地理过程数值和非数值分析。不是去求系统方程的解析解，而是从系统某初始状态出发，去计算短暂时间之后接着发生的状态，再以此为初始状态不断的重复，就能展示系统的行为模式。模拟是对真实事物或者过程的虚拟。模拟要表现出选定的物理系统或抽象系统的关键特性。模拟的关键问题包括有效信息的获取、关键特性和表现的选定、近似简化和假设的应用，以及模拟的重现度和有效性。可以认为仿真是一种重现系统外在表现的特殊的模拟。 仿真：（Emulation）利用模型复现实际系统中发生的本质过程，并通过对系统模型的实验来研究存在的或设计中的系统，又称模拟。即使用项目模型将特定于某一具体层次的不确定性转化为它们对目标的影响，该影响是在项目仿真项目

计算机仿真在电工特种作业培训中的应用(新版)

( 安全技术 ) 单位：_________________________ 姓名：_________________________ 日期：_________________________ 精品文档 / Word文档 / 文字可改 计算机仿真在电工特种作业培 训中的应用(新版) Technical safety means that the pursuit of technology should also include ensuring that people make mistakes

计算机仿真在电工特种作业培训中的应用 (新版) 传统的培训方法大多是在实际系统中进行的。随着工业和科学技术的发展,系统的规模日益庞大、复杂,系统的造价日益昂贵,训练时因操作不当引起的破坏和危险性也大大增加。因此,提高系统运行的安全性事关重大，关键在于人员培训。为了解决培训问题,需要有这样的系统,它能模拟实际系统的工作状况及环境,又可避免采用实际系统时所带来的危险性及高昂代价,这就是培训仿真系统。根据特种作业事故的基本特点，人为因素是事故发生的主要原因。要保证电工特种作业安全,必须控制电工作业人员的不安全行为；要控制电工作业人员的不安全行为,必须提高他们的安全技术素质和安全操作技能，这主要靠培训来完成。目前，电工特种作业培训主要以授课演示加实物模拟平台训练的方法。由于电的危害性，在实物模拟

平台上进行训练也存在许多危险；加之实物模拟平台数量有限，受时间、地点等限制，不能供人员反复使用。这些因素直接影响培训效率和质量。这就要求我们设计出更为有效且具有重复性好、安全性高、经济性好、效率高等优点的电工特种作业仿真培训平台。 计算机仿真培训原理 仿真是以相似原理、系统技术、信息技术及应用领域有关专业技术为基础，以计算机系统与相关物理效应设备及仿真器为工具，利用模型对系统进行研究的多学科综合性技术。 计算机仿真培训原理图。计算机仿真培训系统的组成如图1所示，主要由三个模块组成，即前台多媒体培训仿真界面、信息响应和计算机后台分析计算模块。后台分析计算是计算机仿真培训系统的核心，它包括计算机仿真模型、仿真模块、仿真结果分析、资源库和多媒体计算机网络。用户在前台操作，通过响应模块将用户的操作信息传递给后台，在后台通过仿真模型的判断和运算，做出决策，并将决策信息由后台经过信息处理与响应模块传递到用户操作前台多媒体培训仿真界面，实现交互操作训练。

计算机仿真实验

计算机仿真实验报告 专业：电气工程及其自动化班级：09电牵一班学号：22 姓名：饶坚指导老师：叶满园实验日期：2012年4月30日 一、实验名称 三相桥式SPWM逆变电路仿真 二、目的及要求 1.了解并掌握三相逆变电路的工作原理； 2.进一步熟悉MA TLAB中对Simulink的使用及构建模块； 3.掌握SPWM原理及构建调制电路模块； 4.复习在Figure中显示图形的程序编写和对图形的修改。 三、实验原理与步骤、电路图 1、实验原理图

2、电路原理（采用双极性控制方式） U、V和W三相的PWM控制通常公用一个三角波载波Uc，三相的调制信号Uru、Urv和Urw依次相差120°。 电路工作过程（U相为例）：当Uru>Uc时，上桥臂V1导通，下桥臂V4关断，则U相相对于直流电源假想中点N’的输出电压Uun’=Ud/2。当Uru

对电路模型进行封装如下图示：

其中Subsystem1为主电路，Subsystem2为负载，Subsystem3为检测电路，Subsystem4为输入信号，Subsystem5为调制电路，Scope 为示波器，Repeating Sequence为三角载波。 各子系统电路分别如下所示： Subsystem1 Subsystem2 Subsystem3

计算机仿真与建模实验报告

中南大学 计算机仿真与建模 实验报告 题目：理发店的服务过程仿真 姓名：XXXX 班级：计科XXXX班 学号：0909XXXX 日期：2013XXXX

理发店的服务过程仿真 1 实验案例 (2) 1.1 案例：理发店系统研究 (2) 1.1.1 问题分析 (3) 1.1.2 模型假设 (3) 1.1.3 变量说明 (3) 1.1.4 模型建立 (3) 1.1.5 系统模拟 (4) 1.1.6 计算机模拟算法设计 (5) 1.1.7 计算机模拟程序 (6) 1实验案例 1.1 案例：理发店模拟 一个理发店有两位服务员A和B顾客随机地到达该理发店,每分钟有一个顾客到达和没有顾客到达的概率均是1/2 , 其中60%的顾客理发仅用5分钟,另外40%的顾客用8分钟. 试对前10分钟的情况进行仿真。 （“排队论”，“系统模拟”，“离散系统模拟”，“事件调度法”）

1.1.1 问题分析 理发店系统包含诸多随机因素，为了对其进行评判就是要研究其运行效率， 从理发店自身利益来说，要看服务员工作负荷是否合理，是否需要增加员工等考 虑。从顾客角度讲，还要看顾客的等待时间，顾客的等待队长，如等待时间过长 或者等待的人过多，则顾客会离开。理发店系统是一个典型的排队系统，可以用 排队论有关知识来研究。 1.1.2 模型假设 1． 60％的顾客只需剪发，40％的顾客既要剪发，又要洗发； 2． 每个服务员剪发需要的时间均为5分钟，既剪发又洗发则花8分钟； 3． 顾客的到达间隔时间服从指数分布； 4． 服务中服务员不休息。 1.1.3 变量说明 u ：剪发时间（单位：分钟），u=5m ； v: 既剪发又理发花的时间（单位：分钟），v=8m ； T ： 顾客到达的间隔时间，是随机变量，服从参数为λ的指数分布，（单位： 分钟） T 0：顾客到达的平均间隔时间（单位：秒），T 0＝λ 1； 1.1.4 模型建立 由于该系统包含诸多随机因素，很难给出解析的结果，因此可以借助计算机 模拟对该系统进行模拟。 考虑一般理发店的工作模式，一般是上午9：00开始营业，晚上10：00左 右结束，且一般是连续工作的，因此一般营业时间为13小时左右。 这里以每天运行12小时为例，进行模拟。 这里假定顾客到达的平均间隔时间T 0服从均值3分钟的指数分布， 则有 3小时到达人数约为603 603=?人， 6小时到达人数约为1203 606=?人， 10小时到达人数约为2003 6010=?人， 这里模拟顾客到达数为60人的情况。 （如何选择模拟的总人数或模拟总时间）

Simulink系统仿真课程设计

《信息系统仿真课程设计》 课程设计报告 题目信息系统课程设计仿真 院（系）: 信息科学与技术工程学院 专业班级：通信工程1003 学生姓名： 学号: 指导教师：吴莉朱忠敏 2012年1 月14 日至2012年1 月25 日 华朴中科技大学武昌分校制 信息系统仿真课程设计任务书

20 年月日 目录 摘要 (5)

一、Simulink 仿真设计 (6) 1.1 低通抽样定理 (6) 1.2 抽样量化编码 (9) 二、MATLA仿真设计 (12) 2.1 、自编程序实现动态卷积 (12) 2.1.1 编程分析 (12) 2.1.2 自编matlab 程序： (13) 2.1.3 仿真图形 (13) 2.1.4 仿真结果分析 (15) 2.2 用双线性变换法设计IIR 数字滤波器 (15) 2.2.1 双线性变换法的基本知识 (15) 2.2.2 采用双线性变换法设计一个巴特沃斯数字低通滤波器 (16) 2.2.3 自编matlab 程序 (16) 2.2.4 仿真波形 (17) 2.2.5 仿真结果分析 (17) 三、总结 (19) 四、参考文献 (19) 五、课程设计成绩 (20) 摘要 Matlab 是一种广泛应用于工程设计及数值分析领域的高级仿真平台。它功能

强大、简单易学、编程效率高，目前已发展成为由MATLAB 语言、MATLAB 工作环境、MATLAB 图形处理系统、MATLAB 数学函数库和MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的系统。本次课程设计主要包括MATLAB 和SIMULINKL 两个部分。首先利用SIMULINKL 实现了连续信号的采样及重构，通过改变抽样频率来实现过采样、等采样、欠采样三种情况来验证低通抽样定理，绘出原始信号、采样信号、重构信号的时域波形图。然后利用SIMULINKL 实现抽样量化编码，首先用一连续信号通过一个抽样量化编码器按照A 律13折线进量化行，观察其产生的量化误差，其次利用折线近似的PCM 编码器对一连续信号进行编码。最后利用MATLAB 进行仿真设计，通过编程，在编程环境中对程序进行调试，实现动态卷积以及双线性变换法设计IIR 数字滤波器。 本次课程设计加深理解和巩固通信原理、数字信号处理课上所学的有关基本概念、基本理论和基本方法，并锻炼分析问题和解决问题的能力。

武汉理工大学-计算机仿真实验作业答案

五、（10分）已知系统的传递函数为6 168682)(232+++++=s s s s s s G 。 语言建立系统传递函数模型，并求： ⑴ 该系统的单位阶跃响应；（2分） ⑵ 输入函数为u(t)时的响应；（3分） （u(t)正弦信号，周期2秒，仿真时间8秒，采样周期0.1）； (3) 输入函数为u(t)时的响应；（3分） （u(t)方波输入信号，周期10秒，仿真时间20秒，采样周期0.05） (4) 绘出系统的波德图（Bode ）。（2分） 解答： num=[2 8 6]; den=[1 8 16 6]; sys=tf(num,den); t=0:0.1:8; y1=step(sys,t); u=sin(t*pi); y2=lsim(sys,u,t); subplot(2,2,1);plot(t,y1); grid; title('阶跃响应曲线'); xlabel('响应时间'); ylabel('响应值'); hold on; subplot(2,2,2);plot(t,y2); grid on ; title('对sin(t)的响应曲线'); xlabel('响应时间'); ylabel('响应值'); t=0:0.05:20 u=square(pi/5*t) y3=lsim(sys,u,t); subplot(2,2,3);plot(t,y3) grid on ; title('对方波信号的响应曲线'); xlabel('响应时间'); ylabel('响应值'); subplot(2,2,4);bode(sys);

grid ;title('bode 图'); 运行结果： 六、（10分）设二阶动力学系统的传递函数如下，假设将无阻尼固有频 率固定为ωn ＝1 rad/s ，将阻尼比的值分别设置成ζ＝0，0.1，0.2，0.3，…， MATLAB 语言编程，分析在这些阻尼比ζ的取值下该系统的阶跃响应。 2222)(n n n s s s G ω?ωω++= 解答：wn=1; kesi=[0:0.1:1,2,3,4,5]; figure('color',[1 1 1]); hold on for i=kesi num=wn.^2 den=[1,2*i*wn,wn.^2]; step(num,den);

计算机仿真实验-基于Simulink的简单电力系统仿真

实验七 基于Simulink 的简单电力系统仿真实验 一． 实验目的 1) 熟悉Simulink 的工作环境及SimPowerSystems 功能模块库； 2) 掌握Simulink 的的powergui 模块的应用； 3) 掌握发电机的工作原理及稳态电力系统的计算方法； 4）掌握开关电源的工作原理及其工作特点； 5）掌握PID 控制对系统输出特性的影响。 二．实验内容与要求 单机无穷大电力系统如图7-1所示。平衡节点电压0 44030 V V =∠? 。负荷功率10L P kW =。线路参数：电阻1l R =Ω；电感0.01l L H =。发电机额定参数：额定功率100n P kW =；额定电压440 3 n V V =；额定励磁电流 70 fn i A =；额定频率50n f Hz =。发电机定子侧参数：0.26s R =Ω， 1 1.14 L mH =，13.7 md L mH =，11 mq L mH =。发电机转子侧参数：0.13f R =Ω，1 2.1 fd L mH =。发电机阻尼绕组参数：0.0224kd R =Ω， 1 1.4 kd L mH =，10.02kq R =Ω，11 1 kq L mH =。发电机转动惯量和极对数分别 为224.9 J kgm =和2p =。发电机输出功率050 e P kW =时，系统运行达到稳态状态。在发电机输出电磁功率分别为170 e P kW =和2100 e P kW =时，分析发电机、平衡节点电源和负载的电流、电磁功率变化曲线，以及发电机转速和功率角的变化曲线。

G 发电机节点 V 负 荷 l R l L L P 图 7.1 单机无穷大系统结构图 输电线路 三．实验步骤 1. 建立系统仿真模型 同步电机模块有2个输入端子、1个输出端子和3个电气连接端子。模块的第1个输入端子（Pm）为电机的机械功率。当机械功率为正时，表示同步电机运行方式为发电机模式；当机械功率为负时，表示同步电机运行方式为电动机模式。在发电机模式下，输入可以是一个正的常数，也可以是一个函数或者是原动机模块的输出；在电动机模式下，输入通常是一个负的常数或者是函数。模块的第2个输入端子（Vf）是励磁电压，在发电机模式下可以由励磁模块提供，在电动机模式下为一个常数。 在Simulink仿真环境中打开Simulink库，找出相应的单元部件模型，构造仿真模型，三相电压源幅值为4403，频率为50Hz。按图连接好线路，设置参数，建立其仿真模型，仿真时间为5s，仿真方法为ode23tb，并对各个单元部件模型的参数进行修改，如图所示。

matlab课程设计报告书

《计算机仿真及应用》课程设计报告书 学号：08057102，08057127 班级：自动化081 姓名陈婷，万嘉

目录 一、设计思想 二、设计步骤 三、调试过程 四、结果分析 五、心得体会 六、参考文献

选题一、 考虑如下图所示的电机拖动控制系统模型，该系统有双输入，给定输入)(t R 和负载输入)(t M 。 1、 编制MATLAB 程序推导出该系统的传递函数矩阵。 2、 若常系数增益为：C 1＝Ka ＝Km ＝1，Kr ＝3，C2＝0.8，Kb ＝1.5，时间常数T 1＝5， T 2＝0.5，绘制该系统的根轨迹、求出闭环零极点，分析系统的稳定性。若)(t R 和)(t M 分别为单位阶跃输入，绘制出该系统的阶跃响应图。（要求C 1，Ka ，Km ，Kr ，C2，Kb ， T 1，T 2所有参数都是可调的） 一．设计思想 题目分析： 系统为双输入单输出系统，采用分开计算，再叠加。 要求参数均为可调，而matlb 中不能计算未赋值的函数，那么我们可以把参数设置为可输入变量，运行期间根据要求赋值。 设计思路： 使用append 命令连接系统框图。 选择‘参数=input('inputanumber:')’实现参数可调。 采用的方案： 将结构框图每条支路稍作简化，建立各条支路连接关系构造函数，运行得出相应的传递函数。 在得出传递函数的基础上，使用相应的指令求出系统闭环零极点、画出其根轨迹。 通过判断极点是否在左半平面来编程判断其系统是否稳定。 二．设计步骤 （1）将各模块的通路排序编号

（2）使用append命令实现各模块未连接的系统矩阵 （3）指定连接关系 （4）使用connect命令构造整个系统的模型 三．调试过程 出现问题分析及解决办法： 在调试过程出现很多平时不注意且不易寻找的问题，例如输入的逗号和分号在系统运行时不支持中文格式，这时需要将其全部换成英文格式，此类的程序错误需要细心。 在实现参数可调时初始是将其设为常量，再将其赋值进行系统运行，这样参数可调性差，后用‘参数=input('inputanumber:')’实现。 最后是在建立通路连接关系时需要细心。 四．结果分析 源代码： Syms C1 C2 Ka Kr Km Kb T1 T2 C1=input('inputanumber:') C2=input('inputanumber:') Ka=input('inputanumber:') Kr=input('inputanumber:') Km=input('inputanumber:') Kb=input('inputanumber:') T1=input('inputanumber:') T2=input('inputanumber:') G1=tf(C1,[0 1]); G2=tf(Ka*Kr,[0 1]); G3=tf(Km,[T1 1]); G4=tf(1,[T2 1]); G5=tf(1,[1 0]); G6=tf(-C2,1); G7=tf(-Kb,1); G8=tf(-1,1); Sys=append(G1,G2,G3,G4,G5,G6,G7,G8) Q=[1 0 0;2 1 6;3 2 7;4 3 8;5 4 0;6 5 0;7 4 0;8 0 0;]; INPUTS1=1; OUTPUTS=5; Ga=connect(Sys,Q,INPUTS1,OUTPUTS) INPUTS2=8; OUTPUTS=5; Gb=connect(Sys,Q,INPUTS2,OUTPUTS) rlocus(Ga)

计算机仿真期末大作业Mersenne Twister随机数发生器及随机性测试

Mersenne Twister随机数发生器及随机性测试 一、实验目的 用MATLAB实现Mersenne Twister随机数发生器,并对其随机性进行测试。二、实验原理 伪随机数的产生，首先是选取种子，然后是在此种子基础上根据具体的生成算法计算得到一个伪随机数，然后利用此伪随机数再根据生成算法递归计算出下二个伪随机数，直到将所有不重复出现的伪随机数全部计算出来。这个伪随机数序列就是以后要用到的伪随机数序列。上面的计算过程可以一次性计算完毕，也可以使用一次递归计算一次，每次生成的伪随机数就是这个伪随机数序列中的一个，不过不管怎么样，只要确定了种子，确定了生成算法，这个序列就是确定的了。所谓种子，就是一个对伪随机数计算的初始值。 Mersenne Twister算法是一种随机数产生方法，它是移位寄存器法的变种。该算法的原理：Mersenne Twister算法是利用线性反馈移位寄存器(LFSR)产生随机数的，LFSR的反馈函数是寄存器中某些位的简单异或，这些位也称之为抽头序列。一个n位的LFSR能够在重复之前产生2^n-1位长的伪随机序列。只有具有一定抽头序列的LFSR才能通过所有2^n-1个内部状态，产生2^n - 1位长的伪随机序列，这个输出的序列就称之为m序列。为了使LFSR成为最大周期的LFSR，由抽头序列加上常数1形成的多项式必须是本原多项式。一个n阶本原多项式是不可约多项式，它能整除x^(2*n-1)+1而不能整除x^d+1，其中d能整除2^n-1。例如(32,7,5,3,2,1,0)是指本原多项式x^32+x^7+x^5+x^3+x^2+x+1，把它转化为最大周期LFSR就是在LFSR小邓第32，7，5，2，1位抽头。利用上述两种方法产生周期为m的伪随机序列后，只需要将产生的伪随机序列除以序列的周期，就可以得到(0，1)上均匀分布的伪随机序列了。 伪代码如下： // 建立624位随机序列数组 int[0..623] MT int index = 0 //初始化随机序列数组 function initializeGenerator(int seed) { MT[0] := seed for i from 1 to 623 { MT[i] := last 32 bits of(1812433253 * (MT[i-1] xor(right shift by 30 bits(MT[i-1]))) + i) // 0x6c078965 } }

《计算机仿真技术》试题(含完整答案)

、数值计算，编程完成以下各题(共20分，每小题5 分) 1、脉冲宽度为d，周期为T的矩形脉冲的傅里叶级数如下式描述: d［i.^= sin(^d/T)cos(^：n.) T n」n rd /T 当n =150，d..「T =1；4，- 1/2 ：：：.：：： 1/2，绘制出函数f(.)的图形。 解： syms n t; f=((si n(n *pi/4))/( n*pi/4))*cos(2*pi* n*t); s=symsum(f, n,1,150); y=(1+2*s)/4; x=-0.5:0.01:0.5; Y=subs(y,'t',x); plot(x,Y) 2 0 05x2 5 ■ 5 2、画出函数f (x)二(sin 5x) e .- 5x cos1.5x 1.5x 5.5 x 在区间［3, 5］的图形，求出该函数在区间［3, 5］中的最小值点X min和函数的最小值f min . 解：程序如下 x=3:0.05:5; y=(si n(5*x).A2).*exp(0.05*x.A2)-5*(x.A5).*cos(1.5*x)+1.5*abs(x+5.5)+x.A2.5; mix_where=fi nd(y==mi n(y)); xmin=x(mix_where); hold on; plot(x,y); plot(xmi n,min (y),'go','li newidth',5); str=strcat('(' ,nu m2str(xmi n),',' ,nu m2str(mi n(y)),')'); text(xmi n,min (y),str);

Ylabel('f(x)') 经过运行后得到的图像截图如下： 运行后的最小值点X min =4.6 , f m in = -8337.8625 3、画出函数f (x) = cos2x「e^'x — 2.5 X在口，3]区间的图形, 解该非线 并用编程求性方程 f (x) = 0的一个根，设初始点为X o = 2 . 解： x=1:0.02:3; x0=2; y=@(x)(cos(x).A2).*exp(-0.3*x)-2.5*abs(x); fplot(y,[1,3]); Xlabel('x') Ylabel('f(x)') X仁fzero('(cos(x).A2).*exp(-0.3*x)-2.5*abs(x)',x0) 运行后求得该方程的一个根为z=0.3256 。 4、已知非线性方程组如下，编程求方程组的解，设初始点为[1 0.5 -1].

计算机仿真实训实验报告实验1-4

实验一 熟悉MATLAB 工作环境 16电气5班 周树楠 20160500529 一、实验目的 1.熟悉启动和退出MATLAB 软件的方法。 2.熟悉MATLAB 软件的运行环境。 3.熟悉MATLAB 的基本操作。 二、实验设备及条件 计算机一台（带有MATLAB6.0以上的软件境）。 三、实验内容 1.练习下面指令： cd,clear,dir,path,help,who,whos,save,load 。 2.建立自己的工作目录MYBIN 和MYDATA ，并将它们分别加到搜索路径的前面或者后面。 3.求23)]47(*212[÷-+的算术运算结果。 4.M 文件的建立，建立M 文件，求出下列表达式的值： ?? ????-+=++=+= 545.0212),1ln(21 185sin 2222 1i x x x z e z o 其中

5.利用MATLAB的帮助功能分别查询inv、plot、max、round函数的功能和用法。 四、运行环境介绍及注意事项 1.运行环境介绍 打开Matlab软件运行环境有图1-1所示的界面

图1-1 MATLAB的用户界面 操作界面主要的介绍如下： 指令窗（ Command Window ），在该窗可键入各种送给 MATLAB 运作的指令、函数、表达式，并显示除图形外的所以运算结果。 历史指令窗（ Command History ），该窗记录已经运行过的指令、函数、表达式；允许用户对它们进行选择复制、重运行，以及产生 M 文件。 工作空间浏览器（ Workspace Browser ），该窗口罗列出 MATLAB 工作空间中所有的变量名、大小、字节数；并且在该窗中，可对变量进行观察、编辑、提取和保存。 其它还有当前目录浏览器（ Current Directory Browser ）、 M 文件编辑 / 调试器（Editor/Debugger ）以及帮助导航/ 浏览器（Help Navigator/Browser ）等，但通常不随操作界面的出现而启动。 利用 File 菜单可方便对文件或窗口进行管理。其中 File | New 的各子菜单， M-file （ M 文件）、 Figure （图形窗口）、或 Model （ Simulink 编辑界面）分别可创建对应文件或模块。 Edit 菜单允许用户和 Windows 的剪切板交互信息。 2.在指令窗操作时应特别注意以下几点 1）所有输入的指令、公式或数值必须按下回车键以后才能执行。例如： >>(10*19+2/4-34)/2*3 （回车） ans= 234.7500 2）所有的指令、变量名称都要区分字母的大小写。 3）%作为MATLAB注释的开始标志，以后的文字不影响计算的过程。 4）应该指定输出变量名称，否则MATLAB会将运算结果直接存入默认的输出变量名ans。 5）MATLAB可以将计算结果以不同的精确度的数字格式显示，可以直接在指令视窗键入不同的数字显示格式指令。例如：>>format short (这是默认的) 6）MATLAB利用了↑↓二个游标键可以将所输过的指令叫回来重复使用。按下↑则前一次输入的指令重新出现，之后再按Enter键，即再执行前一次的指令。

计算机仿真实验报告7

山东工商学院计算机仿真及应用实验报告 实验七 MATLAB的基本应用(二)及Simulink仿真 （验证性实验） 学院： 专业班级： 实验时间： 学号： 姓名：

一、实验目的 1、掌握连续信号的仿真和傅里叶分析方法 2、掌握连续系统的分析方法（时域分析法，拉氏变换法和傅里叶分析法）； 3、掌握离散信号的仿真和分析运算方法 4、掌握离散系统的分析方法（时域分析法）； 5、掌握符号运算方法； 6、掌握Simulink仿真工具； 二、实验原理 1、连续信号的仿真和分析法，参考教材第6.1节，重点： 单位冲激信号的仿真方法；单位阶跃信号的仿真方法；复指数信号的仿真方法 2、连续系统的分析方法，参考教材第6.1节，重点： 例6.2，LTI系统的零输入响应的求解方法； 例6.3，LTI系统的冲激响应的求解方法 例6.5，LTI系统的零状态响应的求解方法 例6.6，系统中有重极点时的计算 3、系统的频域分析方法，参考教材第6.2节，重点： 例6.7，方波分解为多次正弦波之和 例6.8：全波整流电压的频谱 例6.10：调幅信号通过带通滤波器 例6.12：用傅里叶变换计算滤波器的响应和输出 4、离散信号的仿真和分析法，参考教材第6.3节，7.1节，重点： 单位脉冲序列impseq，单位阶跃序列stepseq 例7.1：序列的相加和相乘 例7.2：序列的合成与截取 例7.3：序列的移位和周期延拓运算 三、实验内容（包括内容，程序，结果） 以自我编程练习实验为主，熟悉各种方法和设计，结合课堂讲授，实验练习程序代码。 1、根据教材第6.1节的内容，练习连续信号和系统的时域分析和拉氏变换方法。 q602 clear,clc a=input('输入分母系数向量a=[a1,a2,...]= '); n=length(a)-1; Y0=input('输入初始条件向量Y0=[y0,Dy0,D2y0,...]= '); p=roots(a);V=rot90(vander(p));c=V\Y0'; dt=input('dt= ');tf=input('tf= '); t=0:dt:tf;y=zeros(1,length(t));

计算机仿真课程设计

附件1： 北京理工大学珠海学院 《计算机仿真》课程设计说明书题目: 控制系统建模、分析、设计和仿真 学院：信息学院 专业班级： 学号： 学生姓名： 指导教师： 2012年6 月16 日 附件2： 北京理工大学珠海学院 课程设计任务书 2011 ～2012 学年第2学期 学生姓名：专业班级： 指导教师：工作部门：信息学院 一、课程设计题目 《控制系统建模、分析、设计和仿真》 本课程设计共列出10个同等难度的设计题目，编号为：[0号题]、[1号题]、[2号题]、[3号题]、[4号题]、[5号题]、[6号题]、[7号题]、[8号题]、[9号题]。 学生必须选择与学号尾数相同的题目完成课程设计。例如，学号为8xxxxxxxxx2的学生必须选做[2号题]。

[0号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [1号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用一阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [2号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [3号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用一阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [4号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [5号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用一阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹 控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [6号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹 控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [7号题] 控制系统建模、分析、设计和仿真

电力电子系统的计算机仿真

《电力电子系统的计算机仿真》题目：方波逆变电路的计算机仿真

电力电子技术综合了电子电路、电机拖动、计算机控制等多学科知识，是一门实践性和应用性很强的课程。由于电力电子器件自身的开关非线性，给电力电子电路的分析带来了一定的复杂性和困难，一般常用波形分析的方法来研究。仿真技术为电力电子电路的分析提供了崭新的方法。 我们在电力电子技术课程的教学中引入了仿真，对于加深学生对这门课程的理解起到了良好的作用。掌握了仿真的方法，学生的想法可以通过仿真来验证，对培养学生的创新能力很有意义，并且可以调动学生的积极性。实验实训是本课程的重要组成部分，学校的实验实训条件毕竟是有限的，也受到学时的限制。而仿真实训不受时间、空间和物质条件的限制，学生可以在课外自行上机。仿真在促进教学改革、加强学生能力培养方面起到了积极的推动作用。 【关键字】电力电子，MATLAB，仿真。

第一章电力电子与MATLAB软件的介绍 一、电力电子概况 二、MATLAB软件介绍 第二章电力电子器件介绍 一、电力二极管特性介绍 二、晶闸管特性介绍 三、IGBT特性介绍 第三章主电路工作原理 一、单相桥式逆变电路 二、三相桥式逆变电路 三、PWM控制基本原理 第四章仿真模型的建立 一、单极性SPWM触发脉冲波形的产生 二、双极性SPWM触发脉冲波形的产生 三、单极性SPWM方式下的单相桥式逆变电路 四、双极性SPWM方式下的单相桥式逆变电路第五章仿真结果分析 第六章心得体会 第七章参考文献

第一章电力电子与MATLAB软件的 介绍 一、电力电子概况 电力电子技术是一门新兴的应用于电力领域的电子技术，就是使用电力电子器件（如晶闸管，GTO，IGBT等）对电能进行变换和控制的技术。电力电子技术所变换的“电力”功率可大到数百MW甚至GW，也可以小到数W甚至1W以下，和以信息处理为主的信息电子技术不同电力电子技术主要用于电力变换。 电力电子技术分为电力电子器件制造技术和交流技术（整流，逆变，斩波，变频，变相等）两个分支。 一般认为，电力电子技术的诞生是以1957年美国通用电气公司研制出的第一个晶闸管为标志的，电力电子技术的概念和基础就是由于晶闸管和晶闸管变流技术的发展而确立的。此前就已经有用于电力变换的电子技术，所以晶闸管出现前的时期可称为电力电子技术的史前或黎明时期。70年代后期以门极可关断晶闸管（GTO），电力双极型晶体管（BJT），电力场效应管（Power-MOSFET）为代表的全控型器件全速发展（全控型器件的特点是通过对门极既栅极或基极的控制既可以使其开通又可以使其关断），使电力电子技术的面貌焕然一新进入了新的发展阶段。80年代后期，以绝缘栅极双极型晶体管（IGBT 可看作MOSFET和BJT的复合）为代表的复合型器件集驱动功率小，开关速度快，通态压降小，在流能力大于一身，性能优越使之成为现代电力电子技术的主导器件。为了使电力电子装置的结构紧凑，体积减小，常常把若干个电力电子器件及必要的辅助器件做成模块的形式，后来又把驱动，控制，保护电路和功率器件集成在一起，构成功率集成电路（PIC）。目前PIC的功率都还较小但这代表了电力电子技术发展的一个重要方向 利用电力电子器件实现工业规模电能变换的技术，有时也称为功率电子技术。一般情况下，它是将一种形式的工业电能转换成另一种形式的工业电能。例如，将交流电能变换成直流电能或将直流电能变换成交流电能；将工频电源变换为设备所需频率的电源；在正常交流电源中断时，用逆变器（见电力变流器）将蓄电池的直流电能变换成工频交流电能。应用电力电子技术还能实现非电能与电能之间的转换。例如，利用太阳电池将太阳辐射能转换成电能。与电子技术不同，电力电子技术变换的电能是作为能源而不是作为信息传感的载体。因此人们关注的是所能转换的电功率。 电力电子技术是建立在电子学、电工原理和自动控制三大学科上的新兴学科。因它本身是大功率的电技术，又大多是为应用强电的工业服务的，故常将它归属于电工类。电力电子技术的内容主要包括电力电子器件、电力电子电路和电力电子装置及其系统。电力电子器件以半导体为基本材料，最常用的材料为单晶硅；它的理论基础为半导体物理学；它的工艺技术为半导体器件工艺。近代新型电力电子器件中大量应用了微电子学的技术。电力电子电路吸收了电子学的理论基础，根据器件的特点和电能转换的要

电路计算机仿真实验报告

电路计算机仿真分析 实验报告

实验一直流电路工作点分析和直流扫描分析 一、实验目的 1、学习使用Pspice软件，熟悉它的工作流程，即绘制电路图、元件类别的选择及其参数的赋值、分析类型的建立及其参数的设置、Probe窗口的设置和分析的运行过程等。 2、学习使用Pspice进行直流工作点分析和直流扫描分析的操作步骤。 二、原理与说明 对于电阻电路，可以用直观法（支路电流法、节点电压法、回路电流法）列写电路方程，求解电路中各个电压和电流。PSPICE软件是采用节点电压法对电路进行分析的。 使用PSPICE软件进行电路的计算机辅助分析时，首先在capture环境下编辑电路，用PSPICE 的元件符号库绘制电路图并进行编辑、存盘。然后调用分析模块、选择分析类型，就可以“自 动”进行电路分析了。需要强调的是，PSPICE软件是采用节点电压法“自动”列写节点电 压方程的，因此，在绘制电路图时，一定要有参考节点（即接地点）。此外，一个元件为一 条“支路”（branch）,要注意支路（也就是元件）的参考方向。对于二端元件的参考方向定 义为正端子指向负端子。 三、示例实验 应用PSPICE求解图1-1所示电路个节点电压和各支路电流。 图1-1 直流电路分析电路图

4.000V R2 1 2.000A 0V Idc2 4Adc 4.000A 6.000V R1 1 4.000A Idc1 2Adc 2.000A R3 3 2.000A 图1-2 仿真结果 四、选做实验 1、实验电路图 （1）直流工作点分析，即求各节点电压和各元件电压和电流。 （2）直流扫描分析，即当电压源Us1的电压在0-12V之间变化时，求负载电阻R L中电流I RL随电压源Us1的变化曲线。 R4 3 Is3 2Adc 0Vs2 10Vdc RL 1 Is1 1Adc Is2 1Adc R1 4 I Is5 3Adc R2 2 12Vdc IPRINT Vs3 5Vdc Vs4 7Vdc 图1-3 选做实验电路图 2、仿真结果