计算机仿真试题及答案

南京林业大学研究生试卷

课程号 33351 课程名称 计算机仿真 20 16 ~20 17 学年第 1 学期 一、（20分）某系统的传递函数为

)

4)(3)(2)(14(1

43)(2++++++=

S S S S S S S G 试用级联法、串联法和并联法中的两种画出系统模拟结构图，并列出相应的状态方程和输出方程。自拟仿真，比较两者结果。

答：选用级联法和并联法写出传递函数，公式如下：

)

4(1)3(43)2()

31()4

1

()1(/6/4

122/4113/4371/41

/1/43)4)(3)(2)(14(143)(4324

322

+?

+?++?++=++++++=

++++++=s s s s s s s s s s s s s s s s s s s s G

分别画出级联法、并联法的系统模拟结构图及其状态方程和输出方程，图式如下：

432342312432114

143'''64

122

4113437'x x x y x x x x x x x x x x u x ++=

===----

=

4

4

343

2132

121

144

3

'335

43'243

'41

'x y x x x x x x u x x x u x x u x =-=--+=-+=-=

使用龙格—库塔法程序进行仿真，令u=1，系统在t=0-2s 之间（仿真步长0.1s ）则输出y

截图如下：

比较后发现两者仿真结果接近。

二、（20分）已知系统的状态方程及输出方程 1123102123

20312330123

320220220

22x x x x x x

x x x x x x x x x y x x x =+-=??=-+=??

=+-+=??=+-?

拟用四阶龙格-库塔法求系统在t=0-2s 之间的输出值（仿真步长取0.01s 或0.1s ），试编写仿真程序，对仿真结果截图。（单号C++，双号VB ）

答：用c++进行四阶龙格-库塔法，系统在t=0-2s间的（仿真步长0.01）编程：#include

#include

using namespace std;

float x1;

float x2;

float x3;

float x4;

float t;

float f1(float x1,float x2,float x3,float x4,float t)

{

float f1;

return f1=3*x1+2*x2-x3;

}

float f2(float x1,float x2,float x3,float x4,float t)

{

float f2;

return f2=2*x1-x2+2*x3;

}

float f3(float x1,float x2,float x3,float x4,float t)

{

float f3;

return f3=x1+2*x2-x3+2;

}

float f4(float x1,float x2,float x3,float x4,float t)

{

float f4;

return f4=0;

}

int main()

{

float k11,k12,k13,k14,k21,k22,k23,k24,k31,k32,k33,k34,k41,k42,k43,k44;

float t=0,h=0.01;

double y;

float x1=0,x2=0,x3=0,x4=0;

for(int i=1;i<200;i++)

{

t=t+h;

k11 = h * f1(x1, x2, x3, x4, t);

k21 = h * f2(x1, x2, x3, x4, t);

k31 = h * f3(x1, x2, x3, x4, t);

k41 = h * f4(x1, x2, x3, x4, t);

k12 = h * f1(x1 + k11 / 2, x2 + k21 / 2, x3 + k31 / 2, x4 + k41 / 2, t + h / 2);

k22 = h * f2(x1 + k11 / 2, x2 + k21 / 2, x3 + k31 / 2, x4 + k41 / 2, t + h / 2);

k32 = h * f3(x1 + k11 / 2, x2 + k21 / 2, x3 + k31 / 2, x4 + k41 / 2, t + h / 2);

k42 = h * f4(x1 + k11 / 2, x2 + k21 / 2, x3 + k31 / 2, x4 + k41 / 2, t + h / 2);

k13 = h * f1(x1 + k12 / 2, x2 + k22 / 2, x3 + k32 / 2, x4 + k42 / 2, t + h / 2);

k23 = h * f2(x1 + k12 / 2, x2 + k22 / 2, x3 + k32 / 2, x4 + k42 / 2, t + h / 2);

k33 = h * f3(x1 + k12 / 2, x2 + k22 / 2, x3 + k32 / 2, x4 + k42 / 2, t + h / 2);

k43 = h * f4(x1 + k12 / 2, x2 + k22 / 2, x3 + k32 / 2, x4 + k42 / 2, t + h / 2);

k14 = h * f1(x1 + k13, x2 + k23, x3 + k33, x4 + k43, t + h);

k24 = h * f2(x1 + k13, x2 + k23, x3 + k33, x4 + k43, t + h);

k34 = h * f3(x1 + k13, x2 + k23, x3 + k33, x4 + k43, t + h);

k44 = h * f4(x1 + k13, x2 + k23, x3 + k33, x4 + k43, t + h);

x1 = x1 + (k11 + 2 * k12 + 2 * k13 + k14) / 6;

x2 = x2 + (k21 + 2 * k22 + 2 * k23 + k24) / 6;

x3 = x3 + (k31 + 2 * k32 + 2 * k33 + k34) / 6;

}

y=x1+2*x2-2*x3;

cout<<"y="<

return 0;

}

仿真结果截图如下：

三、（20分）已知单位反馈系统的开环传递函数为()

()

()()3

4

1

50

+

+

+

=

s

s

s

s

s G

试用Z 域离散相似法和时域离散相似法进行仿真，编写仿真程序，对仿真结果截图。（输出初态为0，输入为3，仿真步长为0.02S,100-200个周期）

（单号VB ，双号C++）

答：先将函数()()()()34150+++=

s s s s s G 化成为：50(1)1+150

()(4)(3)s+4s+3

s s G s s s s s +==??

++ 画出系统结构图：

使用z 域离散相似法计算得：

3

/))*3exp(1(*50*)()*3exp(*)()1(4/))*4exp(3(*)()1()*4exp(*)()1(*3)()1(t n x t n y n y t n w n w t n x n x t

n w n w --+-=+-+-++-=++=+ 仿真程序如下：

Private Sub Command1_Click() Dim n As Single Dim t As Single

Dim w(200) As Single Dim x(200) As Single Dim z(200) As Single Dim y(200) As Single Dim u(200) As Single t = 0.02 y(0) = 0

For n = 0 To 199

w(n + 1) = w(n) + t * 3

x(n + 1) = x(n) * Exp(-4 * t) + w(n + 1) - w(n) * (3 + Exp(-4 * t)) / 4 y(n + 1) = y(n) * Exp(-3 * t) + x(n) * 50 * (1 - Exp(-3 * t)) / 3 Next

Text1.Text = y(200) End Sub

仿真结果截图如下：

使用时域相似法计算得：

3

/))*3exp(1(*)(*50)(*)*3exp()1())*4exp(*25.025.0(*)3)(()*4exp(*)()1(3

*)()1(T n x n y T n y T n w T n x n x T n w n w --+-=+--++-=++=+ 仿真程序如下：

Private Sub Command1_Click() Dim w(2000) As Single Dim x(2000) As Single Dim y(2000) As Single Dim T As Single Dim n As Single

T = 0.02

y(0) = 0

For n = 0 To 199

w(n + 1) = w(n) + T * 3

x(n + 1) = x(n) * Exp(-4 * T) + (w(n) + 3) * (0.25 - 0.25 * Exp(-4 * T))

y(n + 1) = Exp(-3 * T) * y(n) + 50 * x(n) * (1 - Exp(-3 * T)) / 3

Next

Text1.Text = y(200)

End Sub

仿真结果截图如下：

四、（20分）投掷硬币实验，至少投掷多少次，才能使正面朝上的频率在

0.4-0.6之间的概率不小于0.9 ？试编写仿真程序，对仿真结果截图。（仿真8000-12000次）

答：（仿真12000次）程序如下：Private Sub Command1_Click()

Dim a As Single

Dim i As Integer

Dim j As Integer

Dim t As Integer

Dim s As Single

For j = 1 To 12000

n = 0

For i = 1 To 25

r = Rnd()

If r >= 0.5 Then

n = n + 1

End If

Next

If 0.4 <= n / 25 <= 0.6 Then

m = m + 1

End If

If m / 25 >= 0.9 Then

s = s + 1

End If

Next

Text1.Text = s

End Sub

仿真结果截图如下：

五、（20分）某自选市场，平均每分钟有6位顾客到来，每位顾客购物时间服从（4.2，7.2）分钟均布（不包括交费时间）。有十个收银台，每位顾客服务时间服从（1.2，2）分钟均布。试用GPSS仿真其排队情况。要求画出仿真模型的程序块框图，列出程序清单,对仿真结果截图。(仿真4000-5000次)

答：程序清单及程序块框图如下：

SIMULATE

1 STORAGE 10

GENERATE 1/6,0

ADVANCE 5.7,1.5

QUEUE 1

ENTER 1

DEPART 1

MARK

ADVANCE 1.6,.4

LEAVE 1

TABULATE 1

TERMINATE 1

1 TABLE IA,5,5,10

START 4000

END

计算机仿真课程设计报告

、 北京理工大学珠海学院 课程设计任务书 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：完成程序仿真调试和图形仿真调试；完成课程设计说明书；课程设计答 辩总结。 [ 四、基本要求

计算机仿真技术与CAD习题答案

第0章绪论 0-1 什么是仿真？它所遵循的基本原则是什么？ 答： 仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。 它所遵循的基本原则是相似原理。 0-2 仿真的分类有几种？为什么？ 答： 依据相似原理来分：物理仿真、数学仿真和混合仿真。 物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。 数学仿真：就是应用数学相似原理，构成数学模型在计算机上进行研究。它由软硬件仿真环境、动画、图形显示、输出打印设备等组成。 混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿真系统，这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理仿真。 0-3 比较物理仿真和数学仿真的优缺点。 答： 在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较简单，主要是受控系统的建模、控制方式的确立和计算机编程。数学仿真实验所需的时间比物理仿真大大缩短，实验数据的处理也比物理仿真简单的多。 与数学仿真相比，物理仿真总是有实物介入，效果直观逼真，精度高，可信度高，具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性不强。 0-4 简述计算机仿真的过程。 答： 第一步：根据仿真目的确定仿真方案 根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。 第二步：建立系统的数学模型 对于简单的系统，可以通过某些基本定律来建立数学模型。而对于复杂的系统，则必须利用实验方法通过系统辩识技术来建立数学模型。数学模型是系统仿真的依据，所以，数学模型的准确性是十分重要。

通信对抗原理大作业题目

通信对抗原理仿真大作业题目 基本要求：仿真大作业分组完成，每个组3~5人，至少选择4个题目，并且在每一类中至少选择一个题目。利用MATLAB完成计算机仿真，并且撰写仿真实验报告。大作业完成情况将作为评价平时成绩的依据。 第一类：测频方法仿真 1.FFT法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以上， 基于FFT法进行载波频率测量。画出信号的时域、频域波形，给出FFT法测量的结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比与测量误差的关系。 2.互相关法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以上， 基于互相关法进行载波频率测量。画出信号的时域、频域波形，给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比与测量误差的关系。 3.相位差分法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以 上，基于相位差分法法进行载波频率测量。画出信号的时域、频域波形，给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比与测量误差的关系。 第二类：测向方法仿真 4.相位干涉仪测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上， 基于相位干涉仪测向方法，对不同方向到达的通信信号进行测向。画出信号的时域、频域波形，给出到达方向测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比、到达角与测量误差的关系。 5.到达时差测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上，基 于到达时差测向方法，对不同方向到达的通信信号进行测向。画出信号的时域、频域波形，观察相关函数，给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比、到达角与测量误差的关系。6.多普勒测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上，基于 多普勒测向方法，对不同方向到达的通信信号进行测向。画出信号的时域、频域波形，给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比、到达角与测量误差的关系。 7.沃森-瓦特测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上，基 于沃森-瓦特测向方法，对不同方向到达的通信信号进行测向。画出信号的时域、频域波形，给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线，分析信噪比、到达角与测量误差的关系。 第三类：信号处理技术仿真 8.信号带宽和幅度测量方法仿真。仿真模拟通信信号或者数字通信信号三种以 上，基于FFT法进行信号带宽、信号相对幅度测量。画出信号的时域、频域

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

计算机仿真技术的发展概述及认识 摘要：随着经济的发展和社会的进步，计算机技术高速发展，使人类社会进入了信息时代，计算机作为后期新秀渗入到人们生活中的每一个领域，给人们的生活带来了前所未有的变化。作为新兴的技术，计算机技术在人类研究的各个领域起到了只管至关重要的作用，帮助人类解决了许多技术难题。在科研领域，计算机技术与仿真技术相结合，形成了计算机仿真技术，作为人们科学研究的一种新型方法，被人们应用到各个领域，用来解决人们用纯数学方法或者现实实验无法解决的问题，对科研领域技术成果的形成有着积极地促进作用。 本文在计算机仿真技术的理论思想基础上，分析了计算机仿真技术产生的基本原因，也就是人们用计算机模拟解决问题的优点所在，讨论了模拟、仿真、实验、计算机仿真之间的联系和区别，介绍了计算机仿真技术的发展历程，并查阅相关资料介绍了计算机仿真技术在不同领域的应用，分析并预测了计算机仿真的未来发展趋势。经过查阅大量数据资料并加以分析对比，这对于初步认识计算机仿真技术具有重要意义。 关键词：计算机仿真；模拟；仿真技术；发展 一、引言 计算机仿真技术是以多种学科和理论为基础，以计算机及其相应的软件为工具，通过虚拟试验的方法来分析和解决问题的一门综合性技术。计算机仿真（模拟）早期称为蒙特卡罗方法，是一门利用随机数实验求解随机问题的方法。其原理可追溯到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所示，主要由三个模块组成，即前台多媒体培训仿真界面、信息响应和计算机后台分析计算模块。后台分析计算是计算机仿真培训系统的核心，它包括计算机仿真模型、仿真模块、仿真结果分析、资源库和多媒体计算机网络。用户在前台操作，通过响应模块将用户的操作信息传递给后台，在后台通过仿真模型的判断和运算，做出决策，并将决策信息由后台经过信息处理与响应模块传递到用户操作前台多媒体培训仿真界面，实现交互操作训练。

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 数字滤波器。 本次课程设计加深理解和巩固通信原理、数字信号处理课上所学的有关基本概念、基本理论和基本方法，并锻炼分析问题和解决问题的能力。

计算机仿真期末大作业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 } }

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

五、（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);

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)

计算机仿真课程设计

附件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号题] 控制系统建模、分析、设计和仿真

控制系统数字仿真大作业.

《控制系统数字仿真》课程 大作业 姓名： 学号： 班级： 日期： 同组人员：

目录 一、引言 (2) 二、设计方法 (2) 1、系统数学模型 (2) 2、系统性能指标 (4) 2.1 绘制系统阶跃响应曲线、根轨迹图、频率特性 (4) 2.2 稳定性分析 (6) 2.3 性能指标分析 (6) 3、控制器设计 (6) 三、深入探讨 (9) 1、比例-微分控制器(PD) (9) 2、比例-积分控制（PI） (12) 3、比例-微分-积分控制器（PID） (14) 四、设计总结 (17) 五、心得体会 (18) 六、参考文献 (18)

一、引言 MATLAB语言是当今国际控制界最为流行的控制系统计算机辅助设计语言，它的出现为控制系统的计算机辅助分析和设计带来了全新的手段。其中图形交互式的模型输入计算机仿真环境SIMULINK，为MATLAB应用的进一步推广起到了积极的推动作用。现在，MATLAB语言已经风靡全世界，成为控制系统CAD领域最普及、也是最受欢迎的软件环境。 随着计算机技术的发展和应用，自动控制理论和技术在宇航、机器人控制、导弹制导及核动力等高新技术领域中的应用也愈来愈深入广泛。不仅如此，自动控制技术的应用范围现在已发展到生物、医学、环境、经济管理和其它许多社会领域中，成为现代社会生活中不可或缺的一部分。随着时代进步和人们生活水平的提高，在人类探知未来，认识和改造自然，建设高度文明和发达社会的活动中，控制理论和技术必将进一步发挥更加重要的作用。作为一个自动化专业的学生，了解和掌握自动控制的有关知识是十分必要的。 利用MATLAB软件及其SIMULINK仿真工具来实现对自动控制系统建模、分析与设计、仿真，能够直观、快速地分析系统的动态性能和稳态性能，并且能够灵活的改变系统的结构和参数，通过快速、直观的仿真达到系统的优化设计，以满足特定的设计指标。 二、设计方法 1、系统数学模型 美国卡耐尔基-梅隆大学机器人研究所开发研制了一套用于星际探索的系统，其目标机器人是一个六足步行机器人，如图(a)所示。该机器人单足控制系统结构图如图(b)所示。 要求： （1）建立系统数学模型； （2）绘制系统阶跃响应曲线、根轨迹图、频率特性； （3）分析系统的稳定性，及性能指标； （4）设计控制器Gc(s)，使系统指标满足：ts<10s,ess=0,，超调量小于5%。

MATLAB计算机仿真设计

《计算机仿真技术》 课程设计 姓名： 学号： 班级： 1 专业： 学院： 2016年12月24日

目录 一、设计目的 (1) 二、设计任务 (1) 三、具体要求 (1) 四、设计原理概述 (1) 五、设计内容 (2) 六、设计方案及分析 (2) 1、观察原系统性能指标 (2) 2、手动计算设计 (6) 3、校正方案确定 (8) 七、课程设计总结 (14)

模拟随动控制系统的串联校正设计 一、设计目的 1、通过课程设计熟悉频域法分析系统的方法原理。 2、通过课程设计掌握滞后-超前校正作用与原理。 3、通过在实际电路中校正设计的运用，理解系统校正在实际中的意义。 二、设计任务 控制系统为单位负反馈系统，开环传递函数为) 1025.0)(11.0()(G ++=s s s K s ,设计校正装置，使系统满足下列性能指标：开环增益100K ≥；超调量30%p σ<； 调节时间ts<0.5s 。 三、具体要求 1、使用MATLAB 进行系统仿真分析与设计，并给出系统校正前后的 MATLAB 仿真结果，同时使用Simulink 仿真验证； 2、使用EDA 工具EWB 搭建系统的模拟实现电路，分别演示并验证校正前 和校正后的效果。 四、设计原理概述 校正方式的选择：按照校正装置在系统中的链接方式，控制系统校正方式分 为串联校正、反馈校正、前馈校正和复合校正4种。串联校正是最常用的一种校 正方式，这种方式经济，且设计简单，易于实现，在实际应用中多采用这种校正 方式。串联校正方式是校正器与受控对象进行串联链接的。本设计按照要求将采 用串联校正方式进行校正。 校正方法的选择：根据控制系统的性能指标表达方式可以进行校正方法的确 定。本设计要求以频域指标的形式给出，因此采用基于Bode 图的频域法进行校 正。 几种串联校正简述：串联校正可分为串联超前校正、串联滞后校正和滞后- 超前校正等。 超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提

通信原理大作业

通信原理大作业 1、说明 在通信原理课程中，介绍了通信系统的基本理论，主要包括信道、基带传输、调制 / 解调方法等。为了进一步提高和改善学生对课程基本内容的掌握，进行课程作业方法的改革的试点，设立计算机仿真大作业。成绩将计入平时成绩。 2、要求 参加的同学3~5人一组，选择1?2个题目，协作和共同完成计算机编程和仿真，写出计算机仿真报告。推荐的计算机仿真环境为MATLAB也可以 选择其它环境。 3、大作业选题 (1) 信道噪声特性仿真产生信道高斯白噪声，设计信道带通滤波器对高斯白噪 声进行滤波， 得到窄带高斯噪声。对信道带通滤波器的输入输出的噪声的时域、频域特性进行统计和分析，画出其时域和频域的图形。 (2) 基带传输特性仿真利用理想低通滤波器作为信道，产生基带信号，仿真验证奈氏第一准则的给出的关系。改变低通滤波器的特性，再次进行仿真，验证存在码间干扰时的基带系统输出，画出眼图进行观察。加入信道噪声后再观 察眼图。 (3) 2ASK言号传输仿真 按照2ASK产生模型和解调模型分别产生2ASK言号和高斯白噪声，经过信道传

输后进行解调。对调制解调过程中的波形进行时域和频域观察，并且对解调结果进行误码率测量。2ASK信号的解调可以选用包络解调或者相干解调法。(4) 2FSK信号传输仿真 按照2FSK产生模型和解调模型分别产生2FSK信号和高斯白噪声，经过信道传输后进行解调。对调制解调过程中的波形进行时域和频域观察，并且对解调结果进行误码率测量。2FSK信号的解调可以选用包络解调或者相干解调法。(5) 2PSK信号传输仿真 按照2PSK产生模型和解调模型分别产生2PSK言号和高斯白噪声，经过信道传输后进行解调。对调制解调过程中的波形进行时域和频域观察，并且对解调结果进行误码率测量。2PSK信号的解调选用相干解调法。 ⑹2DPSK言号传输仿真 按照2DPSK产生模型和解调模型分别产生2DPSK言号和高斯白噪声，经过信道传输后进行解调。对调制解调过程中的波形进行时域和频域观察，并且对解调结果进行误码率测量。2DPSK信号的解调可以选用非相干解调或者相干解调法。 (7) 模拟信号的数字传输 产生模拟语音信号，进行PCM编码过程的计算机仿真。仿真发送端采样、 量化编码的过程、仿真接收端恢复语音信号的过程。按照有或者无信道噪 声两种情况分别进行仿真。

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

、数值计算，编程完成以下各题(共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].

计算机仿真课程设计

计算机仿真课程设计 Prepared on 22 November 2020

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

学生必须选择与学号尾数相同的题目完成课程设计。例如，学号为8xxxxxxxxx2的学生必须选做[2号题]。 二、课程设计内容 （一）《控制系统建模、分析、设计和仿真》课题设计内容 [0 [1号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用一阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [2号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [3号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用一阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [4号题] 控制系统建模、分析、设计和仿真 设连续被控对象的实测传递函数为： 用零阶保持器离散化，采样周期取秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 [5号题] 控制系统建模、分析、设计和仿真

大作业题目

控制系统仿真大作业 1、曲线拟合的Matlab实现和优化度检验 通过一个实际的例子，介绍最小二乘曲线拟合法的基本原理，对最小二乘曲线拟合法的Matlab实现方法进行研究，并给出曲线拟合Matlab实现的源程序。论述了Matlab软件在做曲线拟合时的用法，并进行曲线的拟合和相应的图像。 2、基于Matlab的液位串级控制系统 运用组态王和Matlab混合编程的方法设计了一个双容(两个水箱串联)液位串级在线控制系统，由组态王编制人机交互界面，用Matlab完成控制算法，二者通过DDE进行实时数据交换；采用串级控制策略,减小二次干扰的影响，验证其方法的有效性。 3、基于Matlab的变压器差动保护闭环仿真研究 应用Matlab建立了微机保护仿真系统,并对不同原理的变压器差动保护进行了仿真和比较.仿真系统采用积木式结构,根据微机保护的实现原理构建模块,实现保护的闭环仿真,对保护的动作过程进行分析. 4、基于MATLAB/SIMULINK的交流电机调速系统建模与仿真 根据直接转矩控制原理，利用MATLAB/SIMULINK软件构造了一个交流电机调速系统，该系统能够很好地模拟真实系统，实现高效的调速系统设计。仿真结果验证该方法的有效性。 5、基于MCGS和MATLAB的薄膜厚度控制系统仿真 以MCGS组态软件和MATLAB为平台,设计和仿真了一个薄膜厚度控制系统.MCGS完成硬件接口的设置、数据的实时采集、人机对话、以动画的方式显示控制系统的运行情况,MATLAB完成PID参数的自动整定,并利用动态数据交换(DDE)技术建立两者间的通讯.并分析其仿真结果。 6、Matlab在动态电路分析中的应用 用Matlab计算动态电路,可得到解析解和波形图.一阶电路先计算3要素,后合成解

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

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

电力电子技术综合了电子电路、电机拖动、计算机控制等多学科知识，是一门实践性和应用性很强的课程。由于电力电子器件自身的开关非线性，给电力电子电路的分析带来了一定的复杂性和困难，一般常用波形分析的方法来研究。仿真技术为电力电子电路的分析提供了崭新的方法。 我们在电力电子技术课程的教学中引入了仿真，对于加深学生对这门课程的理解起到了良好的作用。掌握了仿真的方法，学生的想法可以通过仿真来验证，对培养学生的创新能力很有意义，并且可以调动学生的积极性。实验实训是本课程的重要组成部分，学校的实验实训条件毕竟是有限的，也受到学时的限制。而仿真实训不受时间、空间和物质条件的限制，学生可以在课外自行上机。仿真在促进教学改革、加强学生能力培养方面起到了积极的推动作用。 【关键字】电力电子，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的功率都还较小但这代表了电力电子技术发展的一个重要方向 利用电力电子器件实现工业规模电能变换的技术，有时也称为功率电子技术。一般情况下，它是将一种形式的工业电能转换成另一种形式的工业电能。例如，将交流电能变换成直流电能或将直流电能变换成交流电能；将工频电源变换为设备所需频率的电源；在正常交流电源中断时，用逆变器（见电力变流器）将蓄电池的直流电能变换成工频交流电能。应用电力电子技术还能实现非电能与电能之间的转换。例如，利用太阳电池将太阳辐射能转换成电能。与电子技术不同，电力电子技术变换的电能是作为能源而不是作为信息传感的载体。因此人们关注的是所能转换的电功率。 电力电子技术是建立在电子学、电工原理和自动控制三大学科上的新兴学科。因它本身是大功率的电技术，又大多是为应用强电的工业服务的，故常将它归属于电工类。电力电子技术的内容主要包括电力电子器件、电力电子电路和电力电子装置及其系统。电力电子器件以半导体为基本材料，最常用的材料为单晶硅；它的理论基础为半导体物理学；它的工艺技术为半导体器件工艺。近代新型电力电子器件中大量应用了微电子学的技术。电力电子电路吸收了电子学的理论基础，根据器件的特点和电能转换的要