MATLAB系统建模与仿真实验

系统建模与仿真实验报告

报告一： 产生10中独立分布的随机数，并检验其一、二阶距的性质。

1、[0,1]区间的均匀分布

采用乘同余法产生均匀分布在（0,1）之间的随机数。乘同余法的递推公式为：

1(mod )n n x x M λ+=

一般情况下，323a λ=±，a 为整数，M 于计算机的字长有关，2m M =,m 为16或32，x 的初值为(1)21b x =+。利用MATLAB 实现，代码如下： function u=undistribution(a,b,m);%乘同余法 lam=8*a-3; M=pow2(m); x(1)=pow2(b)+1; for i=2:10000; y=lam*x(i-1); x(i)=mod(y,M); end u=x/M; end

调用函数，并检验产生随机数的数字特性。 y=undistribution(3,2,32); hist(y,50); E=mean(y); D=var(y);

title('0-1均匀分布直方图');

text(0,-20,strcat('均值为',num2str(E)));

text(0.77,-20,strcat('均值为',num2str(D)));

00.10.20.30.40.50.60.7

0.80.910

50

100

150

200

250

0-1均匀分布直方图

均值为0.50227

方差为0.082933

2、标准正态分布

高斯分布的概率密度函数：

22

()2()x u P x σ-=

；

首先利用前面产生均匀分布随机数的方法生成两组均匀分布的随机数u1,u2

；利用公式：2)Z u π=，Z 服从高斯分布。MATLAB 实现代码如下：

u1=undistribution(3,2,32); u2=undistribution(2,3,32);

z=sqrt(-2*log(u1)).*cos(2*pi*u2); hist(z,100); E=mean(z); D=var(z);

title('标准正态分布直方图')

text(-6,-40,strcat('均值为',num2str(E))); text(6,-40,strcat('方差为',num2str(D)));

-6-4-2024

68050

100150200250300350400

450标准正态分布直方图

均值为0.0081833

方差为1.0025

3、指数分布

指数分布的概率密度函数如下：

,0()0,x e x P x λλ-?≥=??

其它 ；

首先利用前面产生均匀分布随机数的方法生成一组均匀分布的随机数u ；则数列1

ln y u λ

=-

，为均值为λ，方差为2λ的指数分布随机数列。MATLAB 实现

代码如下：

u=undistribution(3,2,32); lam=8;

y=-log(u)/lam; hist(y,150);

E=mean(y); D=var(y);

title('指数分布直方图')

text(0,-130,strcat('均值为',num2str(E))); text(2.3,-130,strcat('方差为',num2str(D)));

00.51 1.52 2.53

200

400

600

800

1000

1200

1400

指数分布直方图

均值为0.12474方差为0.017187

4、广义指数分布

广义指数分布的概率密度函数如下：

0()())p x I x s =-+，

设有两组高斯分布的数列x 、y ，

另1x x =其中s 为信噪比，则221z x y =+服从广义指数分布。MATLAB 实现代码入下： u1=undistribution(3,2,32); u2=undistribution(2,3,32);

x=sqrt(-2*log(u1)).*cos(2*pi*u2); y=sqrt(-2*log(u1)).*sin(2*pi*u2); s=2;

x1=x+sqrt(2*s);

z=x1.^2+y.^2; hist(z,100) E=mean(z); D=var(z);

title('广义指数分布直方图');

text(0,-60,strcat('均值为',num2str(E))); text(63,-60,strcat('方差为',num2str(D)));

010

2030405060

70800

100200300400500600700

800广义指数分布直方图

均值为6.0285

方差为21.1718

5、瑞利分布

瑞利分布的概率密度函数为：

2

22exp(),0

()20,0x x x f x x σσ

?-≥?=??

； 瑞利分布数组可以利用高斯分布数列来产生，设有高斯分布的数列x 、y ，

则z =服从瑞利分布。利用MATLAB 实现代码如下： u1=undistribution(3,2,32);

u2=undistribution(2,3,32);

x=sqrt(-2*log(u1)).*cos(2*pi*u2); y=sqrt(-2*log(u1)).*sin(2*pi*u2); z=sqrt(x.^2+y.^2); hist(z,50) E=mean(z); D=var(z);

title('瑞利分布直方图');

text(0,-65,strcat('均值为',num2str(E))); text(5.5,-65,strcat('方差为',num2str(D)));

012345

670

100200300400500600700

800瑞利分布直方图

均值为1.2496

方差为0.43442

6、广义瑞利分布

广义瑞利分布的概率密度函数如下：

220

222

+()exp -(),02a a p εεεεεσσσ=≥（）I ；

同样，可以利用产生瑞利分布的方法来生产广义瑞利分布数列。设有高斯分

布数列x、y，则z=a为常数。MATLAB 实现代码如下：

u1=undistribution(3,2,32);

u2=undistribution(2,3,32);

x=sqrt(-2*log(u1)).*cos(2*pi*u2);

y=sqrt(-2*log(u1)).*sin(2*pi*u2);

a=2;

x1=x+a;

z=sqrt(x1.^2+y.^2);

hist(z,100)

E=mean(z);

D=var(z);

title('广义瑞利分布直方图');

text(0,-30,strcat('均值为',num2str(E)));

text(7.4,-30,strcat('方差为',num2str(D)));

01234567

890

50100150200250300350

400广义瑞利分布直方图

均值为2.2767

方差为0.84506

7、韦布尔分布

韦布尔分布的概率密度函数如下：

1

()()exp()a a n n x x x x a p x b b b

---=

-； 其中0,0,0x a b ≥>>，

则n z x =+其中i u 是服从N(0,1)的正态分布。MATLAB 实现代码如下： u1=undistribution(3,2,32); u2=undistribution(2,3,32);

x=sqrt(-2*log(u1)).*cos(2*pi*u2); E=mean(x);D=var(x);

y=(x-E)/sqrt(D);%将高斯数列x 转为为服从N(0,1)的高斯分布数列y E1=mean(y);D1=var(y); xn=3;a=3;b=2;

z=xn+b*(-log(y)).^(1/a); hist(z,100); E=mean(z);

D=var(z);

title('?¤2???·?2??±·?í?');

text(3,-60,strcat('?ù?μ?a',num2str(E))); text(6.6,-60,strcat('·?2??a',num2str(D)));

3 3.5

4 4.5

5 5.5

6 6.5

77.50

100200300400500600700

800韦布尔分布直方图

均值为5.0848-0.85012i

方差为1.0795

8、对数正态分布

对数正态分布概率密度函数如下：

2

ln(/)

())2x u p x σ=

-

。 要生成具有对数正态分布的随机数列，只需对高斯分布的数列进行指数变换。MATLAB 实现代码如下： u1=undistribution(3,2,32); u2=undistribution(3,2,16);

x=sqrt(-2*log(u1)).*cos(2*pi*u2); z=exp(x); hist(z,2000);

axis([0,15,0,3000]) E=mean(z); D=var(z);

title('对数正态分布直方图')

text(0,-300,strcat('均值为',num2str(E))); text(12,-300,strcat('方差为',num2str(D)));

51015

500

1000

1500

2000

2500

3000

对数正态分布直方图

均值为1.772

方差为59.456

9、Swerling Ⅱ型分布：

设有两组高斯分布的数列x 、y ，则22()2

r

z x y =+就是Swerling Ⅱ型分布随

机数，其中r 可以取全1数列。MATLAB 实现代码如下： u1=undistribution(3,2,32); u2=undistribution(2,3,32);

x=sqrt(-2*log(u1)).*cos(2*pi*u2); y=sqrt(-2*log(u1)).*sin(2*pi*u2); r=ones(1,10000);

z=(r/2).*(x.^2+y.^2); hist(z,100)

E=mean(z); D=var(z);

title('SwerlingII 分布直方图');

text(0,-150,strcat('均值为',num2str(E))); text(20,-150,strcat('方差为',num2str(D)));

0510

15

20250

20040060080010001200140016001800

2000SwerlingII 分布直方图

均值为0.99789

方差为1.0999

10、χ2

分布

2χ分布的概率密度函数如下：

12

21

()exp()22()

2

n n x p x x n

p +=-

首先生成n 组服从标准正态分布的随机数组i x ，则2n

j i i

x ε=∑服从2χ分布。

MATLAB 实现代码如下：

u1=undistribution(3,2,32); u2=undistribution(2,3,32);

u3=undistribution(5,4,32); u4=undistribution(4,5,32);

g1=sqrt(-2*log(u1)).*cos(2*pi*u2); g2=sqrt(-2*log(u1)).*sin(2*pi*u2); g3=sqrt(-2*log(u3)).*cos(2*pi*u4); g4=sqrt(-2*log(u3)).*sin(2*pi*u4);

x(1,1:10000)=g1;x(2,1:10000)=g2;x(3,1:10000)=g3;x(4,1:10000)=g4;

for i=1:10000

z(i)=x(1,i)^2+x(2,i)^2+x(3,i)^2+x(4,i)^2; end

hist(z,100) E=mean(z); D=var(z);

title('x^2分布直方图');

text(0,-130,strcat('均值为',num2str(E))); text(65,-130,strcat('方差为',num2str(D)));

010

2030405060

70800

500

1000

1500

x 2分布直方图

均值为3.9773

方差为9.2147

现代信号处理Matlab仿真——例611

例6.11 利用卡尔曼滤波估计一个未知常数 题目： 设已知一个未知常数x 的噪声观测集合，已知噪声v(n)的均值为零， 方差为 ，v(n)与x 不相关，试用卡尔曼滤波估计该常数 题目分析： 回忆Kalman 递推估计公式 由于已知x 为一常数，即不随时间n 变化，因此可以得到： 状态方程： x(n)=x(n-1) 观测方程： y(n)=x(n)+v(n) 得到A(n)=1,C(n)=1, , 将A(n)=1,代入迭代公式 得到：P(n|n-1)=P(n-1|n-1) 用P(n-1)来表示P(n|n-1)和P(n-1|n-1),这是卡尔曼增益表达式变为 从而 2v σ1??(|1)(1)(1|1)(|1)(1)(1|1)(1)()()(|1)()[()(|1)()()]???(|)(|1)()[()()(|1)](|)[()()](|1)H w H H v x n n A n x n n P n n A n P n n A n Q n K n P n n C n C n P n n C n Q n x n n x n n K n y n C n x n n P n n I K n C n P n n --=----=----+=--+=-+--=--2()v v Q n σ=()0w Q n =(|1)(1)(1|1)(1)()H w P n n A n P n n A n Q n -=----+21 ()(|1)[(|1)]v K n P n n P n n σ-=--+22(1)()[1()](1)(1)v v P n P n K n P n P n σσ-=--=-+

《工程系统建模》实验报告

《工程系统建模与仿真》实验报告

实验一扭摆法测定物体的转动惯量 实验名称 扭摆法测定物体的转动惯量 同组成员 学号姓名 XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX 三、实验器材 1）转动惯量测试仪 2）数字式电子台秤 3）游标卡尺 4）扭摆及几种有规则的待测转动惯量的物体：金属载物圆盘、塑料圆柱体、木球、验证转动惯量平行轴定理用的金属细杆，杆上有两块可以自由移 动的金属滑块。 四、实验原理 转动惯量的测量，一般都是使刚体以一定形式运动，通过 表征这种运动特征的物理量与转动惯量的关系，进行转换测 量。本实验使物体作扭转摆动，动周期及其它参数的测定计算 出物体的转动惯量。 扭摆的构造如图1-1所示，在垂直轴1上装有一根薄片状 的螺旋弹簧以产生恢复力矩。在轴的上方可以装上各种待测物 体以降低摩擦力矩。3为水平仪，用来调整系统平衡。 将物体在水平面内转过一定角度9 后，在弹簧的恢复力矩作用下物体就 开始绕垂直轴作周期往返扭转运动。 根据虎克定律，弹簧受扭转而产生的恢复力矩M与所转 过的角度9成正比，即：M=-K9 ⑴ 上式中，K为弹簧的扭转常数。 由转动定律M= I卩得：萨M/I ⑵ 令J2=K/I，忽略轴承的 摩擦阻力矩，由式⑴、⑵得： d2K 2 dt2I 上述方程表示扭摆运动具有角简谐振动的特性，角加速度 与角位移成正比，且方向相反。此方程的解为：9=Acos（3+?）。 式中，A为谐振动的角振幅，?为初相位角，①为角速度，此谐振动的周期由于摆 2,用垂直轴与支座间装有轴 承，

系统建模与仿真习题2

系统建模与仿真习题二 1. 考虑如图所示的典型反馈控制系统框图 (1)假设各个子传递函数模型为 66.031.05 .02)(232++-+=s s s s s G ，s s s G c 610)(+=，2 1)(+=s s H 分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法求该系统的传递函数模型。 (2) 假设系统的受控对象模型为s e s s s G 23 )1(12 )(-+=，控制器模型为 s s s G c 32)(+=，并假设系统是单位负反馈，分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法能求出该系统的传递函数模型？如果不能，请近似该模型。 2. 假定系统为： )(0001)(111000100001024269)(t u t x t x ????? ???????+????????????----= [])(2110)(t x t y = 请检查该系统是否为最小实现，如果不是最小实现，请从传递函数的角度解释该模型为何不是最小实现，并求其最小实现。 3. 双输入双输出系统的状态方程：

)(20201000)()(20224264)(75.025.075.125 .1125.15.025.025.025.125.425.25.025.1525.2)(t x t y t u t x t x ??????=????? ???????+????????????------------= （1）试将该模型输入到MATLAB 空间，并求出该模型相应的传递函数矩阵。 （2）将该状态空间模型转化为零极点增益模型，确定该系统是否为最小实现模型。如果不是，请将该模型的传递函数实现最小实现。 （3）若选择采样周期为s T 1.0=，求出离散后的状态方程模型和传递函数模型。 （4）对离散的状态空间模型进行连续变化，测试一下能否变回到原来的系统。 4. 假设系统的传递函数模型为： 222 )(2+++=s s s s G 系统状态的初始值为?? ????-21，假设系统的输入为t e t u 2)(-=。 （1）将该传递函数模型转化为状态空间模型。 （2）利用公式 ?--+=t t t A t t A d Bu e t x e t x 0 0)()()()(0)(τττ求解],0[t 的状态以及系统输出的解析解。 （3）根据上述的解析解作出s ]10,0[时间区间的状态以及系统输出曲线。 （4）采用lsim 函数方法直接作出s ]10,0[时间区间的状态以及系统输出曲线，并与（3）的结果作比较。 5. 已知矩阵 ???? ??????----=212332110A （1）取1:1.0:0=t ，利用expm(At)函数绘制求A 的状态转移矩阵，看运行的速度如何？ （2）采用以下程序绘制A 的状态转移矩阵的曲线，看运行的速度如何？ clc;clear; A=[0 1 -1;-2 -3 3;2 1 -2]; t=0:0.1:2; Nt=length(t);

Matlab仿真实例-卫星轨迹

卫星轨迹 一．问题提出 设卫星在空中运行的运动方程为： 其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示，通过仿真，研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面（r=6400KM ，θ=0）分别以v=8KM/s,v=10KM/s,v=12KM/s 发射时，卫星绕地球运行的轨迹。 二．问题分析 1．卫星运动方程一个二阶微分方程组，应用Matlab 的常微分方程求解命令ode45求解时，首先需要将二阶微分方程组转换成一阶微分方程组。若设，则有： 2．建立极坐标如上图所示，初值分别为：卫星径向初始位置，即地球半径：y(1,1)=6400；卫星初始角度位置：y(2,1)=0；卫星初始径向线速度：y(3,1)=0；卫星初始周向角速度：y(4,1)=v/6400。 3．将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解，可得到一定时间间隔的卫星的径向坐标值y(1)向量；周向角度坐标值y(2)向量；径向线速度y(3)向量；周向角速度y(4)向量。 4．通过以上步骤所求得的是极坐标下的解，若需要在直角坐标系下绘制卫星的运动轨迹，还需要进行坐标变换，将径向坐标值y(1)向量；周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。 5．卫星发射速度速度的不同将导致卫星的运动轨迹不同，实验将绘制卫星分别以v=8KM/s ，v=10KM/s ，v=12KM/s 的初速度发射的运动轨迹。 三．Matlab 程序及注释 1．主程序 v=input('请输入卫星发射速度单位Km/s ：\nv=');%卫星发射速度输入。 axis([-264007000-1000042400]);%定制图形输出坐标范围。 %为了直观表达卫星轨迹，以下语句将绘制三维地球。 [x1,y1,z1]=sphere(15);%绘制单位球。 x1=x1*6400;y1=y1*6400;???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ θθ2)(2 22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2() 3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)] 2(sin[)1(Y )]2(cos[)1(X y y y y

基于MATLAB的变压器仿真 与分析

于MATLAB_Simulink的牵引变压器建模与仿真 基于MATLAB/Simulink的牵引变压器建模与仿真徐（西安铁路局安康供电段新陕西汉中 723000）摘要：针对多种牵引变压器接线方式，建立数学模型，基于Matlab/Simulink仿真软件，建立牵引变压器的仿真模型，并验证数学模型和仿真模型的一致性。利用所建立仿真模型对不同接线形式牵引变压器在不同条件下对公用电网产生的谐波和负序影响进行仿真试验，对研究各种类型的牵引变压器特性在我国电气化铁路的应用提供条件。关键词：牵引变压器；数学模型；仿真模型；Matlab/Simulink 中图分类号：U223.6 文献标识码：A 文章编号：1671－7597（2011）0610061－03 牵引变压器按其特性可分为平衡接线和不平衡接线。其中不平衡接线有单相接线、Vv接线和YNd11接线；平衡接线是试图实现三相两相对称变换而提出的，主要代表方式有Scott，Leblanc、Kubler、Wood-bridge、阻抗匹配接线等。本次主要总结了常用牵引变压器的特点并建立数学模型，包括每种牵引变压器的原理结构、原次边电气量关系等，基于Matlab/Simulink软件建立牵引变压器仿真模型，并对牵引变压器在不同条件下的负序、谐波特性的进行了研究. 1 牵引变压器数学模型研究 1.1 YNd11接线 YNd11变压器接线原理如下图所示，如果忽略激磁电流及其漏阻抗压降，二次侧绕组ac相与一次侧绕组A相同相，cb相与C相同相。由于变压器一次侧绕组A，B，C相与电力系统的相序一致，A相滞后C相，对应的二次侧ac也滞后cb相[2]。其中Z为牵引端口对应变压器漏抗，和β相的端口电压。 1.2 Vv接线 Vv接线牵引变压器接线原理如图2所示。为二次侧空载相即α相图2 Vv接线牵引变压器设Vv接线变压器一次侧、二次侧绕组匝数分别为可得电流输入输出关系[3]：和，电压输入输出关系如下：图1 YNd11接线牵引变压器设YNd11接线变压器一次侧、二次侧绕组匝数分别为和假设变压器原边中性点接地，可以得出一次侧三相电流。，其中为牵引端口对应变压器漏抗，为二次侧空载相即α相和β相的端口电压。 1.3 Scott接线 Scott接线变压器（又称T形接法变压器）属于能完成三相-两相变换的平衡变压器，Scott接线牵引变压器接线原理如图3所示。图3 Scott牵引变压器接线原理图 1 61 设一次侧绕组BC的匝数为次侧绕组AD的匝数为，记，二次的绕组ad、bc的匝数为，则一。可得电流输入输出关系[4]：把一次侧绕组电流用相电流替换，即为：式中，为从三相端子流进变压器的电流。输出端口电压方程为：图6 YNd11接线牵引变压器两供电臂输出电压波形从电压输出波形中可以得到α供电臂电压波形超前β供电臂电压波形120°，在对称阻性负载下，两臂电流输出波形幅值相同，相位相差120°，满足理论值。 2.2 Vv接线牵引变压器 Vv 接线牵引变压器是由两个单相牵引变压器并联而成，仿真模型如图7所示.在仿真模型中牵引变压器T1和T2的原、次边变比设置为110kV/27.5kV。对，于

【财务管理】经济管理系统建模上机实验报告新

班级：工商09-1 姓名：杜莹 学号：0910200105 时间： 辽宁工程技术大学上机实验报告

1SPSS的研究报告 1.1数据缺失值处理的过程（软件操作截图） 依次点击：transform—replacemisingvalues进入界面 将工作负荷导入NewVariable，并命名“组织支持缺失值”点击OK，即得到结果 1.2数据项目分析的方法和过程 a.量表中的反向题重新计分 依次点击：transform—recodeintosamevariables，进入界面

1.3数据信度检验的方法和过程 1.4数据效度检验的方法和过程 1.5数据相关分析的方法和过程 1.6数据回归分析的方法和过程 1.7对数据回归分析结果的解释 SPSS操作：Analyze→Regression→Linear 以情感耗竭为因变量，角色冲突为自变量进行一元线性回归分析模型总体参数表（ModelSummary）:AjustedRSquare是校正的判定系数，用来说明角色冲突进入回归方程，可解释情感耗竭的22.4%的变异量。判定系数越接近于1，表明回是平方和占总离差平方和的比例越大，用X的变动解释Y值变动的那部分就越多，回归效果就越好。 回归效果的方差分析表（ANOVA）：P＜0.05，则拒绝H0，即回归方程的回归效果显著。 回归系数及显著性检验表（Coefficients）：BETA为标准回归系数，T检验的P＞0.05，则接受原假设，即相应系数与0无显著差异；若P＜0.05，则拒绝原假设，即相应系数显著异于0。相应系数是否显著异于0也可看系数的95%的置信区间，若区间包含0，则相应系数为0。根据上述分析结果，写出回归方程。情感耗竭＝0.432角色冲突+2.081 1.8体会、收获和软件操作经验 经过运用spss对数据进行处理以及过程中遇到的问题，熟练地掌握了spss对于数据的各种操作过程，了解了如何将初始数据经过处理后去除非价值数据，变成于我们有用的有价值的数据，为决策做出一定的辅助作用。 通过对细微数据的处理和关注，能帮助我们培养细心，认真，严谨的科学精神，为以后的科研项目工作打下坚实的基础。

系统建模与仿真考试题

1.信息时代认识世界（科学研究）的三种方法是：理论研究、（_实验研究_）、（__ 仿真___）。 2.根据系统状态随时间变化是连续性还是间断性的，可将系统划分为（_连续系统_）、 （__离散系统__）。 3.系统仿真中的三个基本概念是系统、（__模型_）、仿真。 4.拟对某系统进行研究，首先要对系统作出明确的描述，即确定系统各个要素：实体、 属性、活动、（__状态_）、（_事件___）。 ?阶段性知识测试 5.系统仿真有三个基本的活动，即系统建模、仿真建模和（__仿真实验__），联系这 三个活动的是系统仿真的三要素，即系统、模型和计算机（硬件和软件）。 6.系统仿真的一般步骤是：（1）调研系统，明确问题、（2）（___设立目标，收集数据 __）、（3）建立仿真模型、（4）编制程序、（5）运行模型，计算结果、（6）（_统计分析，进行决策__） ?阶段性知识测试 7.仿真软件发展经历了四个阶段（1）高级程序语言阶段；（2）仿真程序包、初级仿 真语言阶段；（3）商业化仿真语言阶段；(4) （_一体化建模与仿真环境_）阶段。 8.常用的仿真软件有Arena、Automod、MATLAB、Promodel、（__WITNESS______）、 （______FLEXSIM___）。 9.求解简单系统问题的“原始”方法是（___解析解决____），借助（___实验__）可大大 提高该方法的效率和精度。 ?阶段性知识测试 10.排队系统可简化表示为A/B/C/D/E。其中A为到达模式；B为（服务模式）、C为服 务台数量、D为系统容量；E为排队规则。 11.常见的排队规则有：先到先服务、后到后服务、优先级服务、最短处理时间优先服 务、随机服务等。请以连线方式将下列排队规则名称的中英文对照起来。 先进先出FIFO 后进先出LIFO 随机服务SIRO 最短处理时间优先SPT 优先级服务PR ?阶段性知识测试 12.模型中，习惯称实体为成分。成分可分为主动成分和被动成分。请问排队系统中的 随机到达的顾客属于（主动）成分（主动/被动）。 13.事件是改变系统状态的瞬间变化的事情。一般指活动的开始和结束。事件可分为必 然事件（主要）、条件事件（次要）、系统事件。其中（______）一般不出现在将来事件表中（FEL）。 14.活动是具有指定长度的持续时间，其开始时间是确定。排队系统主要活动有 （_______）和服务活动。 ?阶段性知识测试 15.仿真时钟表示仿真时间的变量。Witness仿真系统中仿真钟用系统变量（TIME）表 示。 仿真策略，也称仿真算法。离散事件系统适用的仿真策略有（_事件调度法_）、活动扫描法、进程交互法、三阶段法等。 16.建立输入数据模型需要4个步骤：（1）从现实系统收集数据；（2）（_确定输入数据

MATLAB Simulink系统建模与仿真 实验报告

MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************

实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：

1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：

系统建模与仿真实验报告

实验1 Witness仿真软件认识 一、实验目的 熟悉Witness 的启动；熟悉Witness2006用户界面；熟悉Witness 建模元素；熟悉Witness 建模与仿真过程。 二、实验内容 1、运行witness软件，了解软件界面及组成； 2、以一个简单流水线实例进行操作。小部件（widget）要经过称重、冲洗、加工和检测等操作。执行完每一步操作后小部件通过充当运输工具和缓存器的传送带（conveyer）传送至下一个操作单元。小部件在经过最后一道工序“检测”以后，脱离本模型系统。 三、实验步骤 仿真实例操作: 模型元素说明：widget 为加工的小部件名称；weigh、wash、produce、inspect 为四种加工机器，每种机器只有一台；C1、C2、C3 为三条输送链；ship 是系统提供的特殊区域，表示本仿真系统之外的某个地方； 操作步骤： 1：将所需元素布置在界面：

2：更改各元素名称： 如; 3:编辑各个元素的输入输出规则：

4: 运行一周（5 天*8 小时*60 分钟=2400 分钟），得到统计结果。5:仿真结果及分析： Widget： 各机器工作状态统计表：

分析：第一台机器效率最高位100%，第二台机器效率次之为79%，第三台和第四台机器效率低下，且空闲时间较多，可考虑加快传送带C2、C3的传送速度以及提高第二台机器的工作效率，以此来提高第三台和第四台机器的工作效率。 6：实验小结： 通过本次实验，我对Witness的操作界面及基本操作有了一个初步的掌握，同学会了对于一个简单的流水线生产线进行建模仿真，总体而言，实验非常成功。

系统建模与仿真习题3及答案

系统建模与仿真习题三及答案 1.已知系统 )24(32)(21+++=s s s s s G 、2 103)(2+-=s s s G 求G 1(s)和G 2(s)分别进行串联、并联和反馈连接后的系统模型。 解： clc;clear; num1=[2 3]; den1=[1 4 2 0]; num2=[1 -3]; den2=[10 2]; G1=tf(num1,den1); G2=tf(num2,den2); Gs1=series(G1,G2) Gp1=parallel(G1,G2) Gf=feedback(G1,G2) 结果： Transfer function: 2 s^2 - 3 s - 9 ------------------------------ 10 s^4 + 42 s^3 + 28 s^2 + 4 s Transfer function: s^4 + s^3 + 10 s^2 + 28 s + 6 ------------------------------ 10 s^4 + 42 s^3 + 28 s^2 + 4 s Transfer function: 20 s^2 + 34 s + 6 -------------------------------- 10 s^4 + 42 s^3 + 30 s^2 + s – 9 2.某双闭环直流电动机控制系统如图所示：

利用feedback( )函数求系统的总模型。 解： 模型等价为： 编写程序： clc;clear; s=tf('s'); G1=1/(0.01*s+1); G2=(0.17*s+1)/(0.085*s); G3=G1; G4=(0.15*s+1)/(0.051*s); G5=70/(0.0067*s+1); G6=0.21/(0.15*s+1); G7=(s+2)/s; G8=0.1*G1; G9=0.0044/(0.01*s+1); sys1=feedback(G6*G7,0.212); sys2=feedback(sys1*G4*G5,G8*inv(G7)); sys=G1*feedback(sys2*G2*G3,G9) 结果： Transfer function:

MATLAB实现通信系统仿真实例

补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样

物流系统建模与仿真-考前复习题资料-共12页

物流系统建模与仿真考前复习题 1、名词解释（5*4分） （1）系统：系统是由若干可以相互区别、相互联系而又相互作用的要素所组成，在一定的阶层结构形成中分布，在给定的环境约束下，为达到整体的目的而存在的有机集合体。 （2）物流系统模型：物流系统模型是对物流系统特征要素、有关信息和变化规律的一种抽象表达，描述了系统各要素之间的相互关系、系统与环境之间的相互作用，以反映系统的某些本质。 （3）系统仿真：应用数学模型、相应的实用模型的装置、计算机系统、部分实物的仿真系统，对某一给定系统进行数学模拟、半实物模拟、实物模拟，以便分析、设计、研究这种给定系统；或者利用这种仿真训练给定系统的专业人员。 （4）离散事件系统：指系统状态在某些随机时间点上发生离散变化的系统。离散事件动态系统，本质上属于人造系统 （4）实体：实体是描述系统的三个基本要素(实体、属性、活动)之一。在离散事件系统中的实体可分为两大类：临时实体及永久实体。在系统中只存在一段时间的实体叫临时实体。这类实体由系统外部到达系统，通过系统，最终离开系统。临时实体按一定规律不断地到达(产生)，在永久实体作用下通过系统，最后离开系统，整个系统呈现出动态过程。 （5）事件：事件就是引起系统状态发生变化的行为。从某种意义上说，这类系统是由

事件来驱动的。在一个系统中，往往有许多类事件，而事件的发生一般与某一类实体相联系，某一类事件的发生还可能会引起别的事件发生，或者是另一类事件发生的条件等，为了实现对系统中的事件进行管理，仿真模型中必须建立事件表，表中记录每一发生了的或将要发生的事件类型和发生时问，以及与该事件相联的实体的有关属性等。 （6）仿真时钟：仿真钟用于表示仿真时间的变化。离散事件动态系统的状态是在离散时间点上发生变化的，并且由于引起状态变化的事件发生时间的随机性，仿真钟的推进步长是随机的。如果两个相邻发生的事件之间系统状态不发生任何变化，则仿真钟可以跨过这些“不活动”周期。从一个事件发生时刻推进到下一事件发生时刻，仿真钟的推进呈跳跃性，推进速度具有随机性。 （7）事件调度法：仿真模型中的时间控制部件用于控制仿真钟的推进。在事件调度法中，事件表按事件发生时间先后顺序安排事件。时间控制部件始终从事件表中选择具有最早发生时问的事件记录，然后将仿真钟修改到该事件发生时刻。对每一类事件，仿真模型有相应的事件子程序。每一个事件记录包含该事件的若干个属性，其中事件类型是必不可少的，要根据事件类型调用相应的事件子程序。在事件子程序中，处理该事件发生时系统状态的变化，进行用户所需要的统计计算；如果是条件事件，则应首先进行条件测试，以确定该事件是否确能发生。如果条件不满足，则推迟或取消该事件。该事件子程序处理完后返回时问控制部件。 （8）进程交互法：一个进程包含若干个有序事件及有序活动。进程交互法采用进程描述系统，它将模型中的主动成分所发生的事件及活动按时间顺序进行组合，从而形成进程表，一个成分一旦进入进程，它将完成该进程的全部活动。 （9）连接：通过对象之间的连接定义仿真模型的流程，模型中对象之间是通过端口来

MATLAB的建模和仿真

课程设计说明书 题目：基于Matlab的IIR滤波器设计与仿真班级：2012 级电气五班 姓名：王璐 学号：201295014178 指导教师：张小娟 日期：2015年 1 月12日

课程设计任务书

基于MATLAB的IIR滤波器设计与仿真 前言 数字信号处理（digital signal processing，DSP）是从20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是把信号用数字或符号表示的序列，通过计算机或通用（专用）信号处理设备，用数字的数值计算方法处理（例如滤波、变换、压缩、增强、估计、识别等），以达到提取有用信息便于应用处理的目的。数字信号处理系统有精度高、灵活性高、可靠性高、容易大规模集成、时分复用、可获得高性能指标、二维与多维处理等特点。正是由于这些突出的特点，使得它在通信、语音、雷达、地震测报、声呐、遥感、生物医学、电视、仪器中得到愈来愈广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，Digital Filter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（Infinite Impulse Response）滤波器和有限冲激响应FIR（Finite Impulse Response）滤波器。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来结算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的有点，使MATLAB成为一个强大的数学软件，在新的版本中也加入了对C，FORTRAN，C++，JA V A的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 1 数字滤波器概述 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：Y(eωj)=X(eωj)H(eωj) 其中Y(eωj)、X(eωj)分别是数字滤波器的输出序列和输入序列的频域特性（或称为

系统建模仿真实验一

生态平衡建模实验 姓名：冯雪 系别：自动化系 学号：SA14157014

生态平衡建模实验 1、目标 通过此实验了解系统动力学的仿真方法，学会用所学的建模理论来对实际问题进行建模，并对实际问题进行分析。对此生态平衡系统进行仿真实验，改变不同的控制参数，分析实验的结果，得出系统保持生态平衡的条件，为决策者决策提供理论分析基础。 2、原理 通过此次实验了解系统动力学的仿真方法，学会用所学的建模理论来对实际问题进行建模，并对实际问题进行分析。对此生态平衡系统进行仿真实验，改变不同的控制参数，分析实验的结果，得出系统保持生态平衡的条件，为决策者决策提供理论分析基础。系统动力学(System Dynamics)是美国麻省理工学院J.W 福雷斯特(JayW Forrester)教授创立的一门新兴学科。它按照自身独特的方法论建立系统的动态模型，并借助于计算机进行仿真，以处理行为随时间变化的系统的问题。系统动力学首先强调系统性的观点，以及联系、发展、运动的观点，是研究复杂系统，诸如:社会、经济、环境、人口、生态平衡、产业发展等的有效工具。系统动力学的研究对象主要是社会经济系统。社会经济系统的范围十分广泛，凡是涉及到人类的社会活动和经济活动的系统都属于社会系统。诸如本文要研究的人口系统、资源系统、环境系统、经济系统、科技系统、能源系统，都属于社会经济系统。系统动力学的基础是通过实验方法认识系统的行为，为管理决策者提供决策的依据。系统动力学仿真的基本步骤是： （1）明确建模目的 一般来说，系统动力学对社会系统进行仿真实验的主要目的是认识和预测系统的结构和设计最佳参数，为制定合理的政策提供依据。这一步的工作包括观察系统、专家咨询、收集数据资料等，在涉及具体对象系统时，应根据其要求，仿真目的有所侧重。 （2）确定系统边界 系统动力学是将研究对象视为一个系统来处理的。系统是一个相对的概念，相对于所研究问题的实质和建模的目的而言。一旦所要研究的问题的实质和建模的目的已经确定，系统也就确定了，其边界应该是清晰的和唯一的。确定了系统边界之后，才能确定系统的内生变量和外生变量。内生变量是由系统内部反馈结构决定的变量，外生变量是由影响环境因素确定的变量。系统动力学认为系统的行为是基于系统内部的种种因素而产生的，并假定系统的外部因素不给系统的行为以本质的影响，也不受系统内部因素的控制。 （3）因果关系分析 通过因果关系分析，要明确系统内部各要素之间的因果关系，并用表示因果关系的反馈回路来描述。所谓反馈是指：系统中某要素的增加，使受它影响的系统其他要素也发生变化（增加或者减少）。反馈环分为正反馈和负反馈，而正反馈环使系统表现为增长的行为，负反馈使系统表现为收敛的行为。系统动力学认为反馈环是构造系统的第一层次，其多少是系统复杂程度的标志。观察实际系统获得的信息首先用于这一层次。任意两个系统要素从因果关系来看必然是正因果

生产系统建模与仿真试卷(A卷)

上海海洋大学试卷 姓名：学号：专业班名： 一．简述题（共40分） 1.什么是事件？在单通道排队系统中，哪两个典型事件影响系统的状态？这两个典型事件分别发生时，可能会改变系统哪些状态？（5分） 事件是指引起系统状态发生变化的行为或者事情 在单通道派对系统中的典型事件是：顾客到达和服务结束 顾客到达发生，系统可能会由闲开始变为忙，可能引起队长发生变化 服务结束，系统的状态可能有忙变为闲，可能引起队长发生变化 2.分析FMS（柔性制造系统）中的实体、状态、事件和活动。要求每一项写出2个。（8分） 实体：机床、工件 状态：空闲、加工 事件：工件到达、加工结束 活动：工件到达与工件加工开始这之间的一段事件是一个活动

3．在排队模型中，假定用链表来存放排队等待服务的顾客。链表中只有“到达时间”这样的单属性，当前CLOCK ＝10，已用空间表和可用空间表的情形见下图1，并且任何时候队列中的顾客数不会超过4位。若已知排队系统中依次发生的事件如下表1。 请根据表1中列出的事件画出CLOCK ＝15，CLOCK ＝20，CLOCK ＝25时的已用空间表和可用空间表的情形（注意：画出的图形中必须标上行号）。（8分）

4．库存系统仿真中有哪4种类型的事件？当这4种事件同时发生时，系统如何处理4种事件？（4分） 1 货物到达 2 顾客需求 3 仿真结束 4 月初清库 5．请问输入数据分析的基本步骤有哪些，并简述各个步骤的基本内容？（6分） 输入数据收集 分布的识别 参数估计 拟合度检验 6．在稳态仿真中，哪两种方法能够提高仿真结果的精度？（4分） 重复运行次数和增加运行长度

Matlab通信系统建模与仿真例题源代码-第三章

% ch3example1A.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=buttord(f_p,f_s,R_p,R_s, 's'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=butter(n, Wn, 's'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example1B.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=ellipord(f_p,f_s,R_p,R_s,'s'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=ellip(n,R_p,R_s,Wn,'s'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example2A.m f_N=8000; % 采样率 f_p=2100; f_s=2500; R_p=3; R_s=25; % 设计要求指标 Ws=f_s/(f_N/2); Wp=f_p/(f_N/2); % 计算归一化频率 [n, Wn]=buttord(Wp,Ws,R_p,R_s); % 计算阶数和截止频率 [b,a]=butter(n, Wn); % 计算H(z) figure(1); freqz(b,a, 1000, 8000) % 作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率)

系统建模与仿真习题答案forstudents

第一章习题 1-1什么是仿真？它所遵循的基本原则是什么？ 答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识，统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。 它所遵循的基本原则是相似原理。 1-2在系统分析与设计中仿真法与解析法有何区别？各有什么特点？ 答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，计算。它是一种纯物理意义上的实验分析方法，在对系统的认识过程中具有普遍意义。由于受到理论的不完善性以及对事物认识的不全面性等因素的影响，其应用往往有很大局限性。 仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。 1-3数字仿真包括那几个要素？其关系如何？ 答: 通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与计算机。由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模问题；将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验。 1-4为什么说模拟仿真较数字仿真精度低？其优点如何？。 答：由于受到电路元件精度的制约和容易受到外界的干扰，模拟仿真较数字仿真精度低 但模拟仿真具有如下优点： （1）描述连续的物理系统的动态过程比较自然和逼真。 （2）仿真速度极快，失真小，结果可信度高。 （3）能快速求解微分方程。模拟计算机运行时各运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。 （4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进

《系统建模与仿真》复习题样例 江苏大学

《系统建模与仿真》复习题样例：考试内容主要但不 限于如下内容 一、单项选择题(每题.5分,共32题) 1、下列哪个图标表示输送链Conveyor元素（C）。----序号17 A、 B、 C、 D、 2、某条生产线生产产品A，生产速率为1件/3分钟，生产的产品将送入仓库Buf 存储起来，假设生产线产出的第一件A在仿真时刻3，则运行至仿真时刻60，统计进入Buf的零件A的数量可以使用下面的函数(B)。----序号507 A、NPARTS(Buf) B、NPARTS(A) C、NPARTS2(Buf,A,1) D、APARTS(Buf) 3、一次能处理多个部件，即n个部件输入n个部件输出的是：（B ）。----序号218 A、单处理机 B、批处理机 C、装配机 D、生产机 4、在模型中有一属性元素process_time，表示不同的零件在某一机器上所需要的加工时间，那么，在机器详细设计中，对机器的加工时间cycle time栏中应输入（）----序号144 A、process_time B、process_time（） C、match D、cycle time 5、对缓冲器（buffer）中几个缓冲区用矩形框框起来的可视化设计，其所需要使用的可视化属性是（B ）。----序号134 A、name B、rectangle C、patch D、part queue 6、零件（part）到达系统的时间间隔规律在零件详细设计对话框的（ C）中进行设置。----序号148 A、type B、first arrival C、inter arrival D、to 7、有3个零件nut一次性进入系统缓冲区buf1中等待机器加工，机器加工该零部件的时间为3分钟，则计算Bmaxtime（buf1）的结果是（）。----序号261 A、3 B、6 C、8 D、9 8、可以用于机器（machine）输入（from）规则的是（ A）。----序号80 A、pull B、push C、send D、take 9、下列哪个图标表示时间序列曲线Timeseries元素（A ）。----序号21 A、 B、 C、 D、

MATLAB软件应用及仿真设计

电气专业方向综合课程设计任务书

目录 1.单相半波可控整流系统 (1) 1.1晶闸管的仿真 (1) 1.1.1晶闸管模型 (1) 1.1.2晶闸管参数及其设置 (1) 1.2单相半波可控整流电路的仿真 (2) 1.2.1电路图及工作原理 (2) 1.2.2建立仿真模型 (2) 1.2.3模型参数简介与设置 (3) 1.2.4 仿真结果 (5) 2.晶闸管三相桥式整流系统的仿真 (9) 2.1电路图及工作原理 (9) 2.2建立仿真模型 (9) 2.3模型参数简介与设置 (10) 2.4仿真结果 (12) 3.降压斩波电路（Buck变换器） (14) 3.1可关断晶闸管（GTO）的仿真 (14) 3.1.1可关断晶闸管模型 (14) 3.1.2可关断晶闸管参数及其设置 (14) 3.2 Buck变换器的仿真 (15) 3.2.1电路图及工作原理 (15) 3.2.2建立仿真模型 (15) 3.2.3模型参数简介与设置 (16) 3.4 仿真结果 (17) 4.升压斩波电路（Boost变换器） (18) 4.1绝缘栅双极型晶体管（IGBT）的仿真 (18) 4.1.1绝缘栅双极型晶体管模型 (18) 4.1.2绝缘栅双极型晶体管参数及其设置 (18)

4.2 Boost变换器的仿真 (19) 4.2.1电路图及工作原理 (19) 4.2.2建立仿真模型 (19) 4.2.3模型参数简介与设置 (19) 4.3仿真结果 (20) 5.相位控制的晶闸管单相交流调压器带系统的仿真 (21) 5.1电路图及工作原理 (21) 5.2建立仿真模型 (21) 5.3模型参数设置 (21) 5.4 仿真结果 (22) 6.课程设计总结 (24) 参考文献 (25)