电力系统数字仿真作业
电力系统及其自动化研硕10-13 韩暘
1 . Matlab潮流计算作业
本程序的功能是用牛顿——拉夫逊法进行潮流计算。其中：
B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳；5、支路的变比；6、支路首端处于K 侧为1，1侧为0
B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值；4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量；6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；
3为PV节点；
程序如下：
n=input('请输入节点数:n=');
nl=input('请输入支路数:nl=');
isb=input('请输入平衡母线节电号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入由支路参数形成的矩阵:B1=');%变压器侧为1，否则为0
B2=input('请输入各节点参数形成的矩阵:B2=');
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);
% % %---------------------------------------------------
for i=1:nl %支路数
if B1(i,6)==0 %左节点处于1侧
p=B1(i,1);q=B1(i,2);
else %左节点处于K侧
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元
Y(q,p)=Y(p,q); %非对角元
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧
end
%求导纳矩阵
disp('导纳矩阵Y=');
disp(Y)
e(1)=1.00;e(2)=1.00;e(3)=1.10;e(4)=1.05;
f(1)=0;f(2)=0;f(3)=0;f(4)=0;
G=real(Y);
B=imag(Y);
%设S=P+Qj；
S(1)=-0.30-0.18j;
S(2)=-0.55-0.13j;
S(3)=0.50;
P=real(S);
Q=imag(S);
V(3)=1.10;
V(4)=1.05+0j;
k=0;
precision=1;
N1=3;
while precision > 0.00001
for m=1:N1
for n=1:4
Pt(n)=(e(m)*(G(m,n)*e(n)-B(m,n)*f(n))+f(m)*(G(m,n)*f(n)+B(m,n)*e(n))); end
dP(m)=P(m)-sum(Pt);
end
for m=1:3
for n=1:4
Qt(n)=(f(m)*(G(m,n)*e(n)-B(m,n)*f(n))-e(m)*(G(m,n)*f(n)+B(m,n)*e(n))); end
dQ(m)=Q(m)-sum(Qt);
end
Vt=e(3)*e(3)+f(3)*f(3);
dV=(V(3)^2)-Vt;
for m=1:3
for n=1:4
Bi(n)=G(m,n)*f(n)+B(m,n)*e(n);
Ai(n)=G(m,n)*e(n)-B(m,n)*f(n);
end
end
for m=1:3
for n=1:3
if m==n
H(m,n)=sum(Ai)-G(m,n)*e(m)-B(m,n)*f(m);
N(m,n)=-sum(Bi)+B(m,n)*e(m)-G(m,n)*f(m);
J(m,n)=sum(Bi)+B(m,n)*e(m)-G(m,n)*f(m);
L(m,n)=-sum(Ai)+G(m,n)*e(m)+B(m,n)*f(m);
J(3,3)=-2*e(3);
L(3,3)=-2*f(3);
else
H(m,n)=-G(m,n)*e(m)-B(m,n)*f(m);
N(m,n)=B(m,n)*e(m)-G(m,n)*f(m);
J(m,n)=B(m,n)*e(m)-G(m,n)*f(m);
L(m,n)=G(m,n)*e(m)+B(m,n)*f(m);
J(3,1)=0;
L(3,1)=0;
J(3,2)=0;
L(3,2)=0;
end
end
end
%生成雅克比矩阵。
%对角线上的元素
for m=1:3
Ja(2*m-1,2*m-1)=H(m,m);
Ja(2*m-1,2*m)=N(m,m);
Ja(2*m,2*m-1)=J(m,m);
Ja(2*m,2*m)=L(m,m);
end
%非对角线上的元素
for m=1:3
for n=1:3
if m==n
else
Ja(2*m-1,2*n-1)=H(m,n);
Ja(2*m-1,2*n)=N(m,n);
Ja(2*m,2*n-1)=J(m,n);
Ja(2*m,2*n)=L(m,n);
end
end
end
PQ=[dP(1);dQ(1);dP(2);dQ(2);dP(3);dV];
dU=-inv(Ja)*PQ;
precision=max(abs(dU));
for n=1:3
e(n)=e(n)+dU(2*n-1);
f(n)=f(n)+dU(2*n);
end
k=k+1;
for n=1:4
Un(n)=e(n)+f(n)*j;
end
disp('第')
k-1
disp('次迭代')
disp('各节点电压')
Un
disp('迭代过程中节点不平衡量的变化情况') PQ
end
disp('计算平衡节点的功率')
for m=1:4
In(m)=Y(m,4)*Un(m);
end
S=Un(4)*sum(conj(In))
输入数据为：
B1=[1 2 0.1+0.4i 0.03056i 1 0;
2 4 0.08+0.4i 0.02826i 1 0;
3 1 0.3i 0 0.909 1 0;
4 1 0.12+0.5i 0.0384i 1 0]
B2=[1 2 0.1+0.4i 0.03056i 1 0;
2 4 0.08+0.4i 0.02826i 1 0;
3 1 0.3i 0 1. 1 1;
4 1 0.12+0.5i 0.0384i 1 0]
仿真结果为：
请输入节点数:n=4
请输入支路数:nl=>> 4
请输入平衡母线节电号:isb=4
请输入误差精度:pr=0.00001
请输入由支路参数形成的矩阵:B1=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 0.9091 0;4 1 0.12+0.5i 0.0384i 1 0]
请输入各节点参数形成的矩阵:B2=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 1.1 1;4 1 0.12+0.5i 0.0384i 1 0]
导纳矩阵Y=
1.0421 - 8.2428i -0.5882 +
2.3529i 0 +
3.6666i -0.4539 + 1.8911i
-0.5882 + 2.3529i 1.0690 - 4.7274i 0 -0.4808 + 2.4038i
0 + 3.6666i 0 0 - 3.3333i 0
-0.4539 + 1.8911i -0.4808 + 2.4038i 0 0.9346 - 4.2616i
第
ans =
次迭代
各节点电压
Un =
0.9938 - 0.0118i 0.9771 - 0.1127i 1.1000 + 0.1234i 1.0500
迭代过程中节点不平衡量的变化情况
PQ =
-0.2773
-0.0509
-0.5260
0.0196
0.5000
第
ans =
1
次迭代
各节点电压
Un =
0.9848 - 0.0106i 0.9590 - 0.1088i 1.0929 + 0.1249i 1.0500 迭代过程中节点不平衡量的变化情况
PQ =
0.0003
-0.0046
0.0020
-0.0650
0.0028
-0.0152
第
ans =
2
次迭代
各节点电压
Un =
0.9844 - 0.0085i 0.9585 - 0.1083i 1.0924 + 0.1289i 1.0500 迭代过程中节点不平衡量的变化情况
PQ =
0.0011
-0.0010
-0.0026
-0.0001
0.0065
-0.0001
第
ans =
3
次迭代
各节点电压
Un =
0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1289i 1.0500 迭代过程中节点不平衡量的变化情况
PQ =
-0.0004
0.0015
-0.0002
0.0004
0.0004
-0.0000
第
ans =
4
次迭代
各节点电压
Un =
0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1290i 1.0500
迭代过程中节点不平衡量的变化情况
PQ =
1.0e-004 *
-0.2215
-0.0211
0.3567
-0.3126
0.0816
-0.0000
第
ans =
5
次迭代
各节点电压
Un =
0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1290i 1.0500
迭代过程中节点不平衡量的变化情况
PQ =
1.0e-005 *
0.0014
0.0954
-0.2909
0.8859
0.0524
-0.0000
计算平衡节点的功率
S =
0.3679 + 0.2647i
2 . Simulink仿真作业
供电系统如图所示，各元件参数如下：
线路L：长50km，x=0.4Ω/km；
变压器T：
假定供电点电压为106.5kV，保持恒定，空载运行时变压器低压母线发生三相短路。
搭建仿真电路图如下：
仿真结果如下：
图示为短路电流的波形图。
现代工程控制理论 实验报告 实验名称：控制系统数字仿真技术 实验时间： 2015/5/3 目录 一、实验目的 (2) 二、实验容 (3)
三、实验原理 (3) 四、实验方案 (6) 1、分别离散法； (6) 2、整体离散法； (7) 3、欧拉法 (9) 4、梯形法 (9) 5、龙格——库塔法 (10) 五、实验结论 (11) 小结： (14) 一、实验目的 1、探究多阶系统状态空间方程的求解； 2、探究多种控制系统数字仿真方法并对之进行精度比较；
二、 实验容 1、 对上面的系统进行仿真，运用分别离散法进行分析； 2、 对上面的系统进行仿真，运用整体离散法进行分析； 3、 对上面的系统进行仿真，运用欧拉法进行分析； 4、 对上面的系统进行仿真，运用梯形法进行分析； 5、 对上面的系统进行仿真，运用龙泽——库塔法进行分 析； 6、 对上面的几种方法进行总计比较，对他们的控制精度分 别进行分析比较； 三、 实验原理 1、 控制系统状态空间方程整体离散法的求解； 控制系统的传递函数一般为 x Ax Bu Y Cx Du ? =+=+ 有两种控制框图简化形式如下： KI 控制器可以用框图表示如下：
惯性环节表示如下： 高阶系统(s)(1)n K G T = +的框图如下 对于上面的框图可以简写传递函数 x Ax Bu Y Cx Du ? =+=+ 根据各环节间的关系可以列写出式子中出现的系数A 、B 、C 和D ，下面进行整体离散法求传递函数的推导
00 ()0 ...*()...()(t)(0)...*(t)(0)(t)(0)()(0)At At At At At t t At t t A AT t AT A At t t At At A At A t x Ax Bu e e x e Ax e Bu d e x dt Bue dt dt e x Bue dt e x x Bue d e x x e e Bue d x x e Bue d t KT x kT x e τ ττ τττττ ? -? -----------=+=+=?=?=+=+?=+==????? ?①①得②③ ③得令()0 (1)(1)[(1)]0 (1)[(1)]0 ...(1)[(1)](0)...*(1)()(1)T (1)()()() ,kT A kT A kT k T A k T A k T AT k T AT A k T kT T T AT At AT At AT Bue d t K T x k T x e Bue d e x k e x k Bue d k t x k e x k e Budt e x k e Bdt u k e ττττττ τ?-+?++-++-+=++=+-+-=+-=+=+=+?Φ=? ? ? ??④ 令⑤ ⑤④得令令0 (1)()(1) T At m m e Bdt x k x k x k Φ=+=Φ?+Φ?+?得 这样，如果知道系数，就可以知道高阶系统的传递函数和状态空间方程。 2、 在控制系统的每一个环节都加一个采样开关，构成分别 离散法求解系统的状态空间方程； 采样开关其实是一个零阶保持器
工业电力系统动态建模和仿真分析 (Industrial power system dynamic modeling and simulation analysis) 一、概述 工业电力系统： 大型电力系统复杂性：本身有发电机、电动机 中型工业电力系统：即使无发电机，也包括大量中压电动机 意义、内容： 1、确定通过动态建模与仿真分析验证： 1、机组的暂态稳定（极限切除时间） 2、特定的大容量电动机的电压稳定 3、校验电流电压型保护的定植 4、确定低频减载与孤网运行 二、介绍原件与组成： （一）、同步电机实用模型： 1、意义：对于dq0坐标下同步电机方程，如果单独考虑与定子d绕组、q绕组相独立的零轴绕组，则在计及d,q,f,D,Q5个绕组的电磁过渡过程（以绕组磁链或电流为状态量）以及转子机械过渡过程（以ω及δ为状态量）时，电机为七阶模型。对于一个含有上百台发电机的多机电力系统，若再加上其励磁系统、调速器和原动机的动态方程，则将会出现“维数灾”给分析计算带来极大的困难。因此在实际工程问题中，常对同步电机的数学模型作不同程度的简化，以便在不同的场合下使用。 2、对派克方程中的转子变量 若,则 可用定子侧等效量取代原来的转子量，得到用这些实用等效量表示的同步电机实用方程。原派克方程中的定子量，保留易测量及计算的和及和，而消去和两个变量。 3、三阶实用模型 其简单而又能计算励磁系统动态，因而广泛的应用于精度要求不十分高，但仍需计及励磁系统动态的电力系统动态分析中，较适用于凸极机。 模型导出基于： （1）、忽略定子d绕组、q绕组的暂态，即定子电压方程中取P=P=0 （2）、在定子电压方程中，设ω≈（p.u.）在速度变化不大的过渡过程中，其引起的误差很小。 （3）、忽略D绕组、Q绕组，其作用可在转子运动方程补入阻尼项来近似考虑。 及以下三个定子侧等效实用变量： 为消除转子励磁绕组的变量 、 定子励磁电动势 电机(q轴)空载电动势 电机瞬变电动势 （二）、励磁系统数学模型： 描述同步发电机励磁系统(包括励磁调节器)物理过程的数学方程。是电力系统机电暂态过程数学模型的重要组成部分，主要应用于电力系统稳定计算。
上海市高等教育自学考试 工业自动化专业（独立本科段）（B080603）控制系统数字仿真 (02296) 自学考试大纲 上海交通大学自学考试办公室编上海市高等教育自学考试委员会组编 2013年
I、课程的性质及其设置的目的和要求 （一）本课程的性质与设置的目的 “控制系统数字仿真”是利用数字计算进行各种控制系统分析、设计、研究的有力工具，是控制系统工程技术人员必须掌握的一门技术。 本课程是工业自动化专业的专业课程，也是一门理论和实际紧密结合的课程。 通过本课程的学习，学生能掌握系统仿真的基本概念、基本原理及方法；掌握基本的仿真算法及能用高级编程语言在微机上编程实现，学会使用常用的仿真软件。为学习后继课程、从事工程技术工作、科学研究以及开拓性技术工作打下坚实的基础。 （二）本课程的基本要求 1．要求掌握系统、模型、仿真的基本概念，这是学好仿真这门课程的概念基础。 2．掌握常用的连续系统数学仿真算法及能用某种高级编程语言上机实现。 3．初步掌握利用微机来分析、设计、研究控制系统的方法与仿真技术。 （三）本课程与相关课程的联系 先修课程：自动控制原理、现代控制理论基础、高级编程语言。
II、课程内容与考核目标 第1章概论 （一）学习目的和要求 通过本章学习，了解系统的概念，系统的分类方法及特点，仿真的应用目的。了解模型的基本概念，熟悉模型的分类方法及特点。掌握仿真的基本概念，仿真的分类方法及特点。熟悉仿真的一般步骤，仿真技术的应用，熟悉计算机仿真的三要素及基本活动。 （二）课程内容 第一节系统、模型与仿真 1.系统 2.模型 3.仿真 4.仿真科学与技术的发展沿革 第二节系统仿真的一般知识 1.相似理论 2.基于相似理论的系统仿真 3.系统仿真的类型 4.系统仿真的一般步骤 第三节仿真科学与技术的应用 1.仿真在系统设计中的应用 2.仿真在系统分析中的应用 3.仿真在教育与训练中的应用 4.仿真在产品开发及制造过程中的应用 第四节当前仿真科学与技术研究的热点 1.网络化仿真技术 2.复杂系统/开放复杂巨系统的建模与仿真
控制系统数字仿真题库 填空题 1.定义一个系统时，首先要确定系统的；边界确定了系统的范围，边界以外对系统的作用称为系统的，系统对边界以外环境的作用称为系统的。 1.定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以外环境的作用称为系统的输出。 2．系统的三大要素为：、和。 2．系统的三大要素为：实体、属性和活动。 3．人们描述系统的常见术语为：、、和 3．人们描述系统的常见术语为：实体、属性、事件和活动。 4．人们经常把系统分成四类，分别为：、、和 4．人们经常把系统分成四类，它们分别为：连续系统、离散系统、采样数据系统和离散-连续系统。 5、根据系统的属性可以将系统分成两大类：和。 5、根据系统的属性可以将系统分成两大类：工程系统和非工程系统。 6．根据描述方法不同，离散系统可以分为： 和。 6．根据描述方法不同，离散系统可以分为：离散时间系统和离散事件系统。 7. 系统是指相互联系又相互作用的的有机组合。 7. 系统是指相互联系又相互作用的实体的有机组合。 8．根据模型的表达形式，模型可以分为和数学模型二大类，期中数学模型根据数学表达形式的不同可分为二种，分别为：和。8．根据模型的表达形式，模型可以分为物理模型和数学模型二大类，期中数学模型根据数学表达形式的不同可分为二种，分别为：静态模型和动态模型。 9．连续时间集中参数模型的常见形式为有三种，分别为：、和。 9．连续时间集中参数模型的常见形式为有三种，分别为：微分方程、状态方程和传递函数。 10、采用一定比例按照真实系统的样子制作的模型称为，用数学表达式来描述系 统内在规律的模型称为。 10、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统 内在规律的模型称为数学模型。 11．静态模型的数学表达形式一般是方程和逻辑关系表达式等，而动态模型的数学表达形式一般是方程和方程。 11．静态模型的数学表达形式一般是代数方程和逻辑关系表达式等，而动态模型的数
某某大学 《电力系统建模及仿真课程设计》总结报告 题目：基于MATLAB的电力系统短路故障仿真于分析 姓名 学号 院系 班级 指导教师
摘要：本次课程设计是结合《电力系统分析》的理论教学进行的一个实践课程。 电力系统短路故障，故障电流中必定有零序分量存在，零序分量可以用来判断故障的类型，故障的地点等，零序分量作为电力系统继电保护的一个重要分析量。运用Matlab电力系统仿真程序SimPowerSystems工具箱构建设计要求所给的电力系统模型，并在此基础上对电力系统多中故障进行仿真，仿真波形与理论分析结果相符，说明用Matlab对电力系统故障分析的有效性。实际中无法对故障进行实验，所以进行仿真实验可达到效果。 关键词：电力系统；仿真；短路故障；Matlab；SimPowerSystems Abstract: The course design is a combination of power system analysis of the theoretical teaching, practical courses. Power system short-circuit fault, the fault current must be zero sequence component exists, and zero-sequence component can be used to determine the fault type, fault location, the zero-sequence component as a critical analysis of power system protection. SimPowerSystems Toolbox building design requirements to the power system model using Matlab power system simulation program, and on this basis, the power system fault simulation, the simulation waveforms with the theoretical analysis results match, indicating that the power system fault analysis using Matlab effectiveness. Practice can not fault the experiment, the simulation can achieve the desired effect. Keywords: power system; simulation; failure; Matlab; SimPowerSystems - 1 - 目录 一、引言 ............................................ - 3 -
仿真计算 1、在PSCAD中建立典型的同步发电机模型，对同步发电机出口三相短路进行仿真研究。要求： （1）运行“同步发电机短路”模型，截取定子三相短路电流波形，并对波形进行分析，验证与理论分析中包含的各种分量是否一致； 图一同步发电机短路模型
图二、定子三相短路电流 定子三相短路电流中含有直流分量和交流分量，其中周期分量会衰减。三相短路电流直流分量大小不等，但衰减规律相同，均按指数规律衰减，衰减时间常数为Ta，由定子回路电阻和等值电感决定，大约在0.2s。交流分量也按指数规律衰减，它包括两个衰减时间常数，分为次暂态过程、暂态过程和稳态过程。 （2）修改电抗参数Xd（Xd’，X’’d），增加或者减小，截取定子三相电流，并与第一步结果对比分析； 图一是Xd`=0.314 p.u，Xd``=0.280 p.u情况下的定子电流波形；图二是Xd`=0.514 p.u, Xd``=0.280 p.u情况下的定子电流波形。显然，随着Xd`的增大定子的电流在减少。
图三、定子三相短路电流 （3）修改时间常数Td（Td’，T’’d），增加或者减小，截取定子三相电流，并与第一步结果对比分析。 参数Td’=6.55s ，Td”=0.039s时定子电流如图一所示；当参数变为Td’=3.55s ，Td”=0.039s是定子电流如图三所示，显然
图四、定子三相短路电流 2、利用暂态仿真软件对下面的简单电网进行建模，对模型中各元件参数进行详细说明，并进行短路计算。将故障点的电流电压波形及线路M端的电流电压波形、相量图粘贴到课程报告上。 要求：
（1）短路类型为①三相故障；②A相接地；③BC两相故障。 （2）两端系统电势夹角取15o δ=。 （3）故障点设置为线路MN中点（25km处）。 （4）仿真结果包括M、N两侧和短路点处的三相电压、电流的瞬时值波形和短路发生后时刻的三相电压、电流相量图。 三、课程学习心得 通过本课程的学习，你有哪些体会和心得，请写出来。可以从以下几个方面考虑，但不局限于这些方面：通过课程你学到了哪些知识；学会了哪些方法；对电力系统的认识；对课程的建议等。 课程的开始复习了一下简单的电力系统稳态分析部分，然后就进行了课程的重点就是电力系统的暂态分析，其中包括PARK变换、标么值下的磁链方程和电压方程、同步发电机各种电势的表达式、发电机阻抗的概述、（次）暂态电抗和（次）暂态电势、发电机三相短路电流、对称分量法、叠加定理、电力系统简单故障分析。学习了几种电力系统分析中的方法，例如分析同步发电机短路时PARK变换将静止三相坐标系的量转化为旋转坐标系dq0的量，还有分析不对称故障时对称分量法转化到相对简单的对称故障分析中。
控制系统数字仿真 题型举例与总复习 一、填空题 A类基本概念题型 1、系统是指相互联系又相互作用的实体的有机组合。 2、定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以为环境的作用称为系统的输出。 3、系统的三大要素为：实体、属性和活动。 4、根据系统的属性可以将系统分成两大类：工程系统和非工程系统。 5、相似原理用于仿真时，对仿真建模方法的三个基本要求是稳定性、准确性和快速性。 6、根据模型种类不同，系统仿真可分为三种：物理仿真、数字仿真和半实物仿真。 7、按照系统模型特征分类，仿真可分为连续系统仿真及离散事件系统仿真两大类。 8、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统内在规律的模型称为数学模型。 9、计算机仿真是指将模型在计算机上进行试验的过程。 10、系统仿真的三个基本活动是系统建模、仿真建模和仿真试验，计算机仿真的三个要素为：系统、模型与计算机。 11、如果某数值计算方法的计算结果对初值误差和计算误差不敏感，则称该计算方法是稳定的。 12、数值积分法步长的选择应遵循的原则为计算稳定性及计算精度。 13、采样数值积分方法时有两种计算误差，分别为截断误差和舍入误差。 14、三阶隐式啊达姆氏算法的截断误差为O(?4)，二阶龙格-库塔法的局部截断误差为O(?3)，四阶龙格-库塔法的局部截断误差为O(?5)。 15、在判定数值积分方法的稳定域时，使用的测试方程为y?=μy。 16、龙格-库塔法的基本思想是用几个点上函数值的线性组合来避免计算函数的高阶导数，提高数值计算的精度。 17、连续系统仿真中常见的一对矛盾为计算速度和计算精度。 18、离散相似法在采样周期的选择上应当满足采样定理。 19、保持器是一种将离散时间信号恢复成连续信号的装置，零阶保持器能较好地再现阶跃信号，一阶保持器能较好地再现斜坡信号。 20、实际信号重构器不可能无失真地重构信号，具体表现为信号重构器会对被重构的信号产生相位的滞后和幅度的衰减。 21、一般将采样控制系统的仿真归类为连续系统仿真。 22、在控制理论中，由系统传递函数来建立系统状态方程的问题被称为“实现问题”。 23、常用的非线性环节包括:饱和非线性、失灵非线性、迟滞回环非线性。
开题报告 电气工程及其自动化 基于MATLAB/simulink的船舶电力系统建模与故障仿真 一、综述本课题国内外研究动态，说明选题的依据和意义 1、本课题国内外研究动态 船舶电力系统是一个独立的、小型的完整电力系统，主要由电源设备、配电系统和负载组成。船舶电站是船上重要的辅助动力装置，供给辅助机械及全船所需电力。它是船舶电力系统的重要组成部分，是产生连续供应全船电能的设备。船舶电站是由原动机、发电机和附属设备（组合成发电机组）及配电板组成的。最近几年，船舶电站采用电子技术、计算机控制技术，实现船舶电站自动化和船舶电站的全自动控制，实现无人值班机舱。船舶自动化技术正朝着微机监控、全面电气、综合自动化方向发展。船舶电站运行的可靠性、经济性及其自动化程度对保证船舶的安全运营具有极其重要的意义。 国外的某些造船业发达的国家在二十世纪中叶就着手船舶电力系统领域的探索，在船舶电力系统稳态、暂态过程等方面进行了细致的研究。近些年来，挪威挪控公司困.R.co咖l)、英国船商公司(TRANSS)、德国西门子公司(SIEMENS)、-日本三菱公司(MITSUBISHD等大公司开始进行船舶电力系统的建模与分析方面的研究工作。国内针对船舶电力系统的研究起步相对较晚，虽然取得了一定成果，但在理论先进性、系统完整性等方面还存在一定差距，这也在一定程度上导致了目前国产船电设备与世界主要造船国家船电设备存在一定差距、装船率偏低等一系列问题。 目前,国内外最常用的建模软件有四中：分别是：matlab、lingo、Mathematica 和SAS。MATLAB用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。Matlab开发效率高，自带很多数学计算函数，对矩阵支持好。Lingo可以用于求解非线性规划，也可以用于一些线性和非线性方程组的求解等，功能十分强大，是求解优化模型的最佳选择。Mathematica是一款科学计算软件，很好地结合了数值和符号计算引擎、图形系统、编程语言、文本系统、和与其他应用程序的高级连接。SAS 是一个模块化、集成
电 力 系 统 仿 真 作 业 论 文 电控学院 电气0903 刘娟 0906060301
离散可编程三相电压源PLL和可变频率正序电压和功率测量 the Discrete 3-Phase Programmable Voltage Source PLL and Variable-Frequency Positive-Sequence Voltage and Power Measurements 线路图： 线路结构： 一个25KV，100MVA的短路等效电路网络给一个5MW，5Mvar的负载供电。电源的内部电压通过离散的三相可编程电压源装置来提供。三相电压电流测量装置用来检测三个负载电压和电流。 离散的三相PLL装置用来测量频率，也产生一个基于频率变化的系统电压信号。PLL用来驱动两个测量装置，并把变化的频率考虑在内。其中一个用来计算正序负载电压的标幺值，另外一个用来计算负载的有功和无功功率。这两个装置和PLL必须初始化，以保证初始处在稳态。 PLL和两个测量装置分别在Extras/Discrete in the Control Block 和 Extras/Discrete Measurements中可以找到。 整个系统（包括网络，PLL和测量装置）以50us的采集时间来离散。仿真时间4.0秒，仿真参数ode45(Dormand-Prince)。
基本原理： PLL的概念 PLL其实就是锁相环路，简称为锁相环。许多电子设备要正常工作，通常需要外部的输入信号与内部的振荡信号同步，利用锁相环路就可以实现这个目的。锁相环路是一种反馈控制电路，简称锁相环（PLL）。目前锁相环主要有模拟锁相环，数字锁相环以及有记忆能力（微机控制的）锁相环。 PLL的特点 锁相环的特点是：利用外部输入的参考信号控制环路内部振荡信号的频率和相位。因锁相环可以实现输出信号频率对输入信号频率的自动跟踪，所以锁相环通常用于闭环跟踪电路。锁相环在工作的过程中，当输出信号的频率与输入信号的频率相等时，输出电压与输入电压保持固定的相位差值，即输出电压与输入电压的相位被锁住，这就是锁相环名称的由来。 PLL的组成 锁相环通常由鉴相器（PD）、环路滤波器（LF）和压控振荡器（VCO）三部分组成，锁相环组成的原理框图如下图所示。 锁相环中的鉴相器又称为相位比较器，它的作用是检测输入信号和输出信号的相位差，并将检测出的相位差信号转换成uD（t）电压信号输出，该信号经低通滤波器滤波后形成压控振荡器的控制电压uC（t），对振荡器输出信号的频率实施控制。
《控制系统数字仿真》课程 大作业 姓名： 学号： 班级： 日期： 同组人员：
目录 一、引言 (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上机 题目：基于MATLAB的电力系统自动重合闸 所在学院： 学科专业： 学号： 学生姓名： 指导教师： 二零一五年四
摘要 分析了单相自动重合闸的工作特性，并利用MATLAB软件搭建了220kv电力系统的自动重合闸的仿真模型，模拟系统发生单相接地、三相相间短路故障，断路器跳闸后自动重合闸的工作过程。 关键词：电力系统自动重合闸MATLAB 短路故障
目录 1 引言 (1) 2 模型中主要模块的选择和参数 (2) 2.1同步发电机模块 (2) 2.2 变压器模块 (2) 2.3 输电线路模块 (3) 2.3.1 150km线路模块 (3) 2.3.2 100km线路模块 (4) 2.1 电源模块 (5) 2.3 负载模块 (6) 2.3.1 三相串联RLC负载Load1 (6) 2.3.2 三相串联RLC负载Load4 (7) 2.4 断路器模块 (8) 2.5 测量模块 (9) 2.6 显示模块 (9) 2.7 其他模块 (9) 2.8 仿真参数设置 (10) 3 仿真结果及波形分析 (10) 3.1 线路单相重合闸 (10) 3.2 线路三相重合闸 (12) 总结 (13) 参考文献 (14)
基于Matlab的电力系统自动重合闸 1 引言 随着技术的发展，电力系统的规模越来越复杂。从实际条件与安全角度考虑，不太可能进行电力系统科研实验，因而电力系统数字仿真成为了电力系统研究、规划和设计的重要手段。电力系统仿真软件如BPA，EMTP，PSCAD/ EMTDC ，NETOMAC，PSASP，MATLAB等，正向着多功能，具有更高的可移植性方向发展。其中在MATLAB 中，电力系统模型可以在Simulink环境下直接搭建，Simulink电力系统元件库中有多种多样的电气模块，电力系统大多数元件都包含。其中，可以直接调用。电力系统大部分故障是瞬时性故障，因此采用自动重合闸后，电力系统发生瞬时性故障时供电的连续性、系统的稳定性得到很大的提高。此外，自动重合闸有效纠正由于断路器或继电保护误动作引起的误跳闸。 本文以MATLAB为工具，对简单系统的线路单相重合闸和线路三相重合闸进行分析与研究。 1.1 仿真模型的设计和实现 电力系统正常运行时可以认为是三相对称的，即电压、电流对称，且具有正弦波形。下图为理想情况下220kv电力系统的模型。 图 1 220kv电力系统模型
控制系统数字仿真题库 一、填空题 1. 定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以为环境的作用称为系统的输出。 2．系统的三大要素为：实体、属性和活动。 3．人们描述系统的常见术语为：实体、属性、事件和活动。 4．人们经常把系统分成四类，它们分别为：连续系统、离散系统、采样数据系统和离散-连续系统。 5、根据系统的属性可以将系统分成两大类：工程系统和非工程系统。 6．根据描述方法不同，离散系统可以分为：离散时间系统和离散事件系统。 7. 系统是指相互联系又相互作用的实体的有机组合。 8．根据模型的表达形式，模型可以分为物理模型和数学模型二大类，其中数学模型根据数学表达形式的不同可分为二种，分别为：静态模型和动态模型。 9、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统内在规律 的模型称为数学模型。 10．静态模型的数学表达形式一般是代数方程和逻辑关系表达式等，而动态模型的数学表达形式一般是微分方程和差分方程。 11．系统模型根据描述变量的函数关系可以分类为线性模型和非线性模型。 12 仿真模型的校核是指检验数字仿真模型和数学模型是否一致。 13．仿真模型的验证是指检验数字仿真模型和实际系统是否一致。 14．计算机仿真的三个要素为：系统、模型与计算机。 15．系统仿真的三个基本活动是系统建模、仿真建模和仿真试验。 16．系统仿真根据模型种类的不同可分为：物理仿真、数学仿真和数学-物理混合仿真。 17．根据仿真应用目的的不同，人们经常把计算机仿真应用分为四类，分别为： 系统分析、系统设计、理论验证和人员训练。 18．计算机仿真是指将模型在计算机上进行实验的过程。 19. 仿真依据的基本原则是:相似原理。 20. 连续系统仿真中常见的一对矛盾为计算速度和计算精度。 21．保持器是一种将离散时间信号恢复成连续信号的装置。 22．零阶保持器能较好地再现阶跃信号。 23. 一阶保持器能较好地再现斜坡信号。 24. 二阶龙格-库塔法的局部截断误差为O()。 25．三阶隐式阿达姆斯算法的截断误差为:O()。
昆明学院 《电力系统建模及仿真》课程设计 任务书 适用于：电气工程及其自动化专业 （电气工程方向） 自动控制与机械工程学院电子电气教研室 2015年6月
一、课程设计的目的 该课程设计是在完成《电力系统分析》的理论教学之后安排的一个实践教学环节。其目的在于巩固和加深对电力系统潮流和短路电流计算基本原理的理解，学习和掌握应用计算机进行电力系统设计和计算的方法，培养学生独立分析和解决问题的能力。 二、课程设计的基本要求 掌握电力系统等值模型和参数计算，以及潮流和短路计算的基本原理，学会应用计算机计算系统潮流分布和短路电流的方法。 三、课程设计选题原则 该课程设计是根据电力系统分析课程内容，结合实际工程和科研的电力系统网络进行系统的潮流和短路电流计算。 四、课程设计的任务及要求 1、基本要求 （1）用Matlab中Simulink组件的SimPowerSystems工具箱构建设计要求所给的电力系统模型，在所给电力系统中K处选取不同故障类型（三相短路、单相接地短路、两相短路、两相接地短路进行仿真，比较仿真结果，给出自己的结论。（电力系统接线图见附录1，选做一题） （2）基于Matlab/Simulink，搭建附录2所示电力网络模型，并进行潮流计算。 2、课程设计论文编写要求 纸张A4、要求书写整齐，字数不少于2000字。 （1）封面包括：《电力系统建模与仿真课程设计》总结报告、专业、班级、学号、姓名、指导教师（具体格式附后） （2）论文包括目录、摘要、正文、参考文献、心得体会等。 要求：画出完整电路图、参数标注清楚；按照具体项目要求，完成仿真内容并记录仿真结果，给出自己的结论。 五、时间分配 1、查阅资料、熟悉Matlab中Simulink组件的SimPowerSystems工具箱（1天）； 2、基于Matlab/Simulink的电力系统短路故障的仿真与分析（3天）；
1-1什么是仿真？它所遵循的基本原则是什么？ 答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识，统汁数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。 它所遵循的基本原则是相似原理。 1-2在系统分析与设计中仿真法与解析法有何区別？各有什么特点？ 答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，il?算。它是一种纯物理意义上的实验分析方法，在对系统的认识过程中具有普遍意义。由于受到理论的不完善性以及对事物认识的不全而性等因素的影响，其应用往往有很大局限性。 仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。 1-3数字仿真包括那几个要素？其关系如何？ 答：通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与让算机。由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模问题：将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验。 1-4为什么说模拟仿真较数字仿真精度低？其优点如何？o 答：由于受到电路元件精度的制约和容易受到外界的下?扰，模拟仿真较数字仿真精度低 但模拟仿真具有如下优点： （1）描述连续的物理系统的动态过程比较自然和逼真。 （2）仿真速度极快，失真小，结果可信度髙。 （3）能快速求解微分方程。模拟汁算机运行时0运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。 （4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进行非实时仿真。 （5）易于和实物相连。 1-5什么是CAD技术？控制系统CAD可解决那些问题？ 答：CAD技术，即计算机辅助设计（Computer Aided Design），是将计算机高速而精确的计算能力, 大容量存储和数据的能力与设讣者的综合分析，逻辑判断以及创造性思维结合起来，用以快速设计进程，缩短设计周期，提髙设计质量的技术。 控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制理论。此外，自适应控制，自校正控制以及最优控制等现代控制测略都可利用CAD技术实现有效的分析与设计。 1-6什么是虚拟现实技术？它与仿真技术的关系如何？ 答：虚拟现实技术是一种综合了计算机图形技术，多媒体技术，传感器技术，显示技术以及仿真技术等多种学科而发展起来的高新技术。 1-7什么是离散系统？什么是离散事件系统？如何用数学的方法描述它们？ 答：本书所讲的“离散系统”指的是离散时间系统，即系统中状态变量的变化仅发生在一组离散时刻上的系统*它一般采用差分方程.离散状态方程和脉冲传递函数来描述。 离散事件系统是系统中状态变量的改变是由离散时刻上所发生的事件所驱动的系统。这种系统的输入输出是随机发生的，一般采用概率模型来描述。 1-8如图1-16所示某卫星姿态控制仿真实验系统，试说明： （1）若按模型分类，该系统属于那一类仿真系统？ （2）图中“混合汁算机”部分在系统中起什么作用？ （3）与数字仿真相比该系统有什么优缺点？ 答：（1）按模型分类，该系统属于物理仿真系统“ （2）混合计算机集中了模拟仿真和数字仿真的优点，它既可以与实物连接进行实时仿真，计算一些复杂函数，又可以对控制系统进行反复迭代讣算。其数字部分用来模拟系统中的控制器，而模拟部分用于模拟控制对象。（4）与数字仿真相比，物理仿真总是有实物介入，效果逼真，精度高，具有实时性与在线性的特点, 但其构成复杂，造价较髙，耗时过长，通用性不强。
控制系统数字仿.. 交卷时间：2016-04-01 21:13:58 一、单选题 1. (2分) 列出工作内存中的变量名称以及细节，只需在命令窗口输入________。 ? A. what ? B. who ? C. echo on ? D. whose 得分：0知识点：控制系统数字仿真作业题 答案D解析 2. (2分) 在Simulink中，运行系统仿真的工具栏图标为 ? A. ? B. ? C. ? D. 得分：0知识点：控制系统数字仿真作业题 答案C解析 3. (2分) 设A＝[0 2 3 4;1 3 5 0]，B＝[1 0 5 3;1 5 0 5]则A>=B的结果为________。
? A. ? B. ? C. ? D. 得分：0知识点：控制系统数字仿真作业题 答案B解析 4. (2分) 若B=[3 2 7 4 9 6 1 8 0 5]，则B([end-3:end])为________。 ? A. 3 7 1 ? B. 3 2 7 4 9 9 4 7 2 3 ? C. 3 4 ? D. 1 8 0 5 得分：0知识点：控制系统数字仿真作业题 答案D解析 5. (2分) 执行以下指令之后E,F的值分别为________。 A=[1 2 3; 4 5 6]; B=[3 4 5; 7 8 9]; C=3; E = A+B; F = B+C ? A. E=[4 5 6;7 8 9] F=[6 7 8;10 11 12] ? B. E=[6 7 8;10 11 12] F=[4 5 6;7 8 9] ? C. E=[4 5 6;7 8 9 F=[6 4 5;10 8 9] ? D. E=[4 5 6;7 8 9] F=[3 4 8;7 8 12]
控制系统数字仿真实验报告 班级：机械1304 姓名：俞文龙 学号： 0801130801
实验一数字仿真方法验证1 一、实验目的 1.掌握基于数值积分法的系统仿真、了解各仿真参数的影响； 2.掌握基于离散相似法的系统仿真、了解各仿真参数的影响； 3.熟悉MATLAB语言及应用环境。 二、实验环境 网络计算机系统（新校区机电大楼D520），MATLAB语言环境 三实验内容 （一）试将示例1的问题改为调用ode45函数求解，并比较结果。 实验程序如下； function dy = vdp(t,y) dy=[y-2*t/y]; end [t,y]=ode45('vdp',[0 1],1); plot(t,y); xlabel('t'); ylabel('y');
（二）试用四阶RK 法编程求解下列微分方程初值问题。仿真时间2s ，取步长h=0.1。 ?????=-=1 )0(2y t y dt dy 实验程序如下: clear t0=0; y0=1; h=0.1; n=2/h; y(1)=1; t(1)=0; for i=0:n-1 k1=y0-t0^2; k2=(y0+h*k1/2)-(t0+h/2)^2; k3=(y0+h*k2/2)-(t0+h/2)^2;
k4=(y0+h*k3)-(t0+h)^2; y1=y0+h*(k1+2*k2+2*k3+k4)/6; t1=t0+h; y0=y1; t0=t1; y(i+2)=y1; t(i+2)=t1; end y1 t1 figure(1) plot(t,y,'r'); xlabel('t'); ylabel('y'); （三）试求示例3分别在周期为5s的方波信号和脉冲信号下的响应，仿真时间20s，采样周期Ts=0.1。
如图所示为一无穷大功率供电的三相对称系统，短路发生前系统处于稳定运行状态。假设a 相电流为)sin(i |0|0 ?αω-+=t (1-1) 式中，222|0|m )'()'(L L R R U I m +++=ω,) '()'(arctan 0R R L L ++=ω? 假设t=0s 时刻，f 点发生三相短路故障。此时电路被分成俩个独立回路。由无限大电源供电的三相电路，其阻抗由原来的)'()'(L L j R R +++ω突然减小为L j R ω+。由于短路后的电路仍然是三相对称的，依据对称关系可以得到a 、b 、c 相短路全电流的表达式 式中，22m )(L R U I m ω+=为短路电流的稳态分量的幅值。 短路电流最大可能瞬时值称为短路电流的冲击值，以m i 表示。冲击电流主要用于检验电气设备和载流导体在短路电流下的受力是否超过容许值，即所谓的动稳定度。由此可得冲击电流的计算式为 式中，im K 称为冲击系数，即冲击电流值对于短路电流周期性分量幅值的倍数；αT 为时间常数。 短路电流的最大有效值m I 是以最大瞬时值发生的时刻（即发生短路经历约半个周期）为中心的短路电流有效值。在发生最大冲击电流的情况下，有 短路电流的最大有效值主要用于检验开关电器等设备切断短路电流的能力。 无穷大功率电源供电系统仿真模型构建 假设无穷大功率电源供电系统如图所示，在0.02s 时刻变压器低压母线发生三相短路故障，仿真其短路电流周期分量幅值和冲击电流的大小。线路参数为;km 17.0,km 4.0,5011Ω=Ω==r x km L 变压器额定容量A MW S N ?=20，电压U s %=10.5，短路损耗KW P s 135=?,空载损耗KW P 220=?，空载电流I 0%=0.8，变比11110=T K ，高低压绕组均为Y 形联结；并设供点电压为110KV 。其对应的Simulink 仿真模型如图所示 无穷大功率电源供电系统的Simulink 仿真图 在Simulink 仿真图中各模块名称及提取路径见表。 表 仿真电路中各模块名称及提取路径
辽宁工程技术大学 电力系统继电保护综合训练一 设计题目方向性电流保护的建模与仿真 指导教师刘健辰 院（系、部）电气与控制工程学院 专业班级电网13—1班 学号1305080116 姓名苏小平 日期2017/01/05
智能电网系综合训练标准评分模板
电力系统继电保护析综合训练一 任务书 本次综合训练目的在于通过对双侧电源电力系统的方向电流保护的建模与仿真，巩固和运用所学到的方向电流保护理论知识，掌握Matlab 仿真软件的使用方法，培养学生分析问题和解决问题的能力。 双侧电源电力系统结构图如下： 系统基本参数如下，线路长度和短路点位置见后面的班级数据表。 电源：o 11510kV M E =∠ ,o 1050kV N E =∠ ,o ,,0.22673.13s M s N Z Z ==∠Ω 线路：LGJ-240/40型架空线，单位正序阻抗o 10.45173.13/km z =∠Ω。 设计要求： 利用Matlab/Simulink 建立仿真模型，完成仿真计算，分析仿真结果。 设计说明书内容： 1、 任务书 2、 电网相间短路的方向电流保护原理 3、 利用Matlab/Simulink 建立仿真模型 4、 设置故障，完成仿真计算 5、 分析仿真结果 6、 重新设置两侧电源参数，分析对方向电流保护的影响
说明：1）1~7组每组3人；第8组4人。 2）将自己姓名填入表中
目录 一、综合训练目的 (1) 二、电网相间短路的方向电流保护原理 (1) 三、中性点不接地系统故障特征 (1) 四、仿真 (1) （1）、仿真模型 (2) （2）、设置故障、观察故障特征 (2) 五、结果分析 (3) 六、重新设置两侧电源参数，分析对方向电流保护的影响 (3) 参考文献 (4)
吉大16春学期《控制系统数字仿真》在线作业一答案 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII
吉大16春学期《控制系统数字仿真》在线作业一 一、单选题（共 15 道试题，共 30 分。） 1. 数值积分法中，其计算精度p=2的算法是（）。 A. 欧拉法 C. 四阶—龙格库塔法 D. 以上都不是 满分：2 分 2. i=2; a=2i;b=2*i;c=2*sqrt(-1);程序执行后；a, b, c的值分别是（）。 A. a=4,b=4,c=2.0000i B. a=4,b=2.0000i, c=2.0000i C. a=2.0000i, b=4,c=2.0000i D. a=2.0000i,b=2.0000i,c=2.0000i 满分：2 分 3. 下列哪条指令是求矩阵的行列式的值（）。 A. inv B. diag C. det D. eig 满分：2 分 4. CAD软件中我们一般都用（）单位来做图以达到最佳的效果。 A. 米 B. 厘米 C. 毫米 D. 分米 满分：2 分 5. 在CAD命令输入方式中以下不可采用的方式有（）。
A. 点取命令图标 B. 在菜单栏点取命令 C. 用键盘直接输入 D. 利用数字键输入 满分：2 分 6. 角度x=[30 45 60]，计算其正弦函数的运算为（）。 A. SIN（deg2rad(x)） B. SIN(x) C. sin(x) D. sin(deg2rad(x)) 满分：2 分 7. 绘制系统奈氏曲线的命令是（）。 A. step B. pzmap C. nyquist D. sgrid 满分：2 分 8. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为（）。 A. a'*b B. a.*b C. a*b D. a-b 满分：2 分 9. 下列哪个变量的定义是不合法的（）。 A. abcd-3 B. xyz_3 C. abcdef D. x3yz