仿真作业一、（排队系统仿真）
某小邮电所的顾客到达时间间隔为18+6分钟的均匀分布，邮局职员对每个顾客的服务时间服从16+4分钟的均匀分布，系统为只有一个服务员的等待制排队系统，试对100个顾客做仿真，并求：
（1）系统的平均等待队长；
（2）顾客在系统中的平均等待时间；
（3）系统服务员忙的概率
要求：利用GPSS语言在计算机上编制仿真程序实现仿真，同时对仿真结果进行分析。
2．有一个理发店只有一个理发师，顾客到达间隔服从λ=0.2人/分的负指数分布，服务员对顾客的理发时间服从μ=0.4人/分，假设理发店等待的顾客座位只有三个，若到达的顾客发现所有这些座位已被占领，则不等待便离去。试仿真100个顾客，并求：
（1）系统服务员忙的概率；
（2）100个顾客中由于座位已被占而没有接受服务离开系统的顾客数。
要求：利用GPSS语言在计算机上编制仿真程序实现仿真，同时对仿真结果进行分析。
注：参数为λ=0.2人/分的负指数分布可用系统自带的负指数分布随机数产生，函数为：GENERATE (EXPONENTIAL( 1, 0, 1/λ ))
3．某汽车加油站只有一台加油器，加油的汽车到达间隔时间为100+10秒的均匀分布。到达加油站后先排队等候再加油，加油服务时间为80+8秒的均匀分布。加油后司机要去缴费窗口排队等待缴费，缴费时间为50+30秒的均匀分布。试对100个顾客做仿真，并求：
（1）加油站与缴费窗口的利用率；
（2）加油站与缴费处顾客的平均等待时间。
要求：利用GPSS语言在计算机上编制仿真程序实现仿真，同时对仿真结果进行分析。
4．在第3题中，若加油站有两台效率相同的加油器，而加油的汽车排一个队(省去缴费环节)。设加油的汽车到达间隔时间为40+10秒的均匀分布。试对100个顾客做仿真，并求：各个加油器的利用率和汽车的平均等待时间。
要求：利用GPSS语言在计算机上编制仿真程序实现仿真，同时对仿真结果进行分析。
5.（选作）某诊所有2个位医生，普通病人以每隔7±5分钟的速率到达，每次看病需10±6分钟，急诊病人以每隔60±40分钟速率到达，每次看病需25±15分钟，急诊具有优先级，可以优先服务，请仿真100小时，估算急诊病人和普通病人的平均排队等待时间。要求：利用GPSS语言在计算机上编制仿真程序实现仿真，同时对仿真结果进行分析。
nGDOU-B-11-112广东海洋大学学生实验报告书（学生用表) 课程名称课程号学院(系)信息学院 专业班级 学生姓名学号 实验地点04002 实验日期 实验一连时间信号的MATLAB表示 和连续时间LTI系统的时域分析 一、实验目的 1．掌握MATLAB产生常用连续时间信号的编程方法，并熟悉常用连续时间信号的波形和特性； 2．运用MATLAB符号求解连续系统的零输入响应和零状态响应； 3．运用MATLAB数值求解连续系统的零状态响应； 4．运用MATLAB求解连续系统的冲激响应和阶跃响应； 5．运用MATLAB卷积积分法求解系统的零状态响应。 二、实验原理 1. 连续信号MATLAB实现原理 从严格意义上讲，MATLAB数值计算的方法并不能处理连续时间信号。然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能够被MATLAB处理，并且能较好地近似表示连续信号。
MATLAB提供了大量生成基本信号的函数。比如常用的指数信号、正余弦信号等都是MATLAB的内部函数。为了表示连续时间信号，需定义某一时间或自变量的范围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图。 三、实验内容 1．实例分析与验证 根据以上典型信号的MATLAB函数，分析与验证下列典型信号MATLAB程序，并实现各信号波形图的显示，连续信号的图形显示使用连续二维图函数plot()。 (1) 正弦信号：用MATLAB命令产生正弦信号2sin(2/4) ππ+，并会出时间0≤t≤3的波形图。 程序如下： K=2;w=2*pi;phi=pi/4; t=0:0.01:3; ft=K*sin(w*t+phi); plot(t,ft),grid on; axis([0,3,-2.2,2.2]) title('正弦信号')
信号与系统的MATLAB 仿真 一、信号生成与运算的实现 1.1 实现)3(sin )()(π±== =t t t t S t f a )(sin )sin()sin(sin )()(t c t t t t t t t S t f a '=' '== ==πππ π ππ m11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果： 1.2 实现)10() sin()(sin )(±== =t t t t c t f ππ m12.m t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数 plot(t,f); % 绘制sinc(t)的波形 运行结果： 1.3 信号相加：t t t f ππ20cos 18cos )(+= m13.m syms t; % 定义符号变量t f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：
1.4 信号的调制：t t t f ππ50cos )4sin 22()(+= m14.m syms t; % 定义符号变量t f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.5 信号相乘：)20cos()(sin )(t t c t f π?= m15.m t=-5:0.01:5; % 定义时间范围向量 f=sinc(t).*cos(20*pi*t); % 计算函数f(t)=sinc(t)*cos(20*pi*t) plot(t,f); % 绘制f(t)的波形 title('sinc(t)*cos(20*pi*t)'); % 加注波形标题 运行结果：
Flexsim应用案例示例 示例一港口集装箱物流系统仿真 （根据：肖锋，基于Flexsim集装箱码头仿真平台关键技术研究，武汉：武汉理工大学硕士学位论文，2006改编） 1、港口集装箱物流系统概述与仿真目的 1.1港口集装箱物流系统概述 1.2港口集装箱物流系统仿真的目的 2、港口集装箱物流系统的作业流程 2.1港口集装箱物流系统描述 2.2港口集装箱物流系统作业流程 2.3港口集装箱物流系统离散模型分析 3、港口集装箱物流系统仿真模型 3.1港口集装箱物流系统布局模型设计 3.2港口集装箱物流系统设备建模 3.3港口集装箱物流系统仿真 4、仿真运行及数据分析 4.1仿真运行及数据处理 4.2仿真数据的结果分析 小结与讨论 示例二物流配送中心仿真 （根据：XXX改编） 1、物流配送中心概述与仿真目的 1.1物流配送中心简介 1.2仿真目的 2、配送中心的作业流程描述 2.1配送中心的功能 2.2配送中心的系统流程
3、配送中心的仿真模型 3.1配送中心的仿真布局模型设计 3.2配送中心的设备建模 3.3配送中心的仿真 4、仿真运行及数据分析 4.1仿真运行及数据处理 4.2仿真数据结果分析 4.3系统优化 小结与讨论 “我也来编书”示例 示例一第X章排队系统建模与仿真学习要点 1、排队系统概述 2、排队系统问题描述 3、排队系统建模 4、排队系统仿真 5、模型运行与结果分析 小结 思考题与习题（3-5题） 参考文献 1、李文锋，袁兵，张煜.2010.物流系统建模与仿真（第6章） 北京：科学出版社 2、王红卫，谢勇，王小平，祁超.2009.物流系统仿真（第6章） 北京：清华大学出版社 3、马向国，刘同娟.2012.现代物流系统建模、仿真及应用案例（第5章）
信号与系统 单项选择题 1、 ( ) 1. D. x(t) 2. -x(t) 3. x(0) 4. -x(0) 2、设是带限信号， rad/s，则对进行均匀采样的最大间隔为( ) 1. 0.2s 2. 0.5s 3. 0.1s 4. 0.3s 3、下列信号中属于数字信号的是（）。 1. 2. 3. 4. 4、设系统输入输出关系为y(t)=x(t)cos(t) ，则系统为（）。 1.因果稳定
2.非因果稳定 3.因果不稳定 4.非因果不稳定 5、关于无失真传输的充要条件，下列叙述中正确的是（）。 1.系统的幅频特性为常数 2.系统的相频特性与频率成正比 3. 4. 6、 1. 0 2. 1 3.无穷大 4.不存在 7、 1. 2. 1 3. 4.无法确定 8、关于数字频率，下列表达中错误的是（） 1.数字频率的高频为π附近
2.数字频率的低频为0和2π附近 3.数字频率为模拟频率对采样频率归一化的频率 4.数字频率的单位为Hz 9、 1. 2. 3. 4. 10、关于三个变换之间的关系，下列叙述错误的是（）。 1.若原信号收敛，虚轴上的拉氏变换就是傅里叶变换 2. s域的左半平面映射到z域的单位圆内部
3.从s域到z域的映射是单值映射 4. s域的右半平面映射到z域的单位圆外部 11、关于信号的分解，下列叙述正确的是（） 1.傅里叶级数是一致性意义下的正交分解 2.任意普通信号可分解为冲激函数的叠加，可用卷积形式来描述 3.信号能分解为实分量和虚部分量，故可对信号进行滤波 4.由于信号的可分解性，故在时域中可用冲激响应来表征系统12、 1. 2 2. 4 3. -2 4. -4 13、 1. 2. 3. 4. 14、关于稳定性的描述，下列叙述中错误的是（）。
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
知识点2 1. 结合具体制造系统或服务系统，分析离散事件动态系统的基本特征。 2. 什么叫“状态空间爆炸”？产生状态空间爆炸的原因是什么？它给系统性能分析带来哪些 挑战？ 3. 常用的离散事件系统建模方法有哪些，它们是如何分类的？ 4. 什么是马尔可夫特性？它在离散事件系统建模与分析中有什么作用？ 5. 根据功能不同，仿真模型（程序）可以分为哪三个层次？分析三个层次之间的关系。 6. 分析事件调度法、活动循环法、进程交互法和消息驱动法等仿真调度方法的特点，在分 析每种调度方法基本原理的基础上，阐述几种仿真调度方法之间的区别与联系，并绘制每种仿真调度方法的流程图。 7. 结合具体的离散事件系统，如银行、理发店、餐厅、超市、医院、作业车间等，采用事 件调度法、活动循环法或进程交互法分析建立此类系统的仿真模型，试分析仿真模型中的建模元素以及仿真调度流程。 8. 从系统描述、建模要点、仿真时钟推进机制等层面，比较事件调度法、活动循环法和进 程交互法的异同之处。 9. 什么叫仿真时钟，它在系统仿真中有什么作用？什么叫仿真时钟推进机制？常用的仿真 时钟推进机制有哪些？它们的主要特点是什么，分别适合于怎样的系统？ 10.结合具体的离散事件系统，分析若采用固定步长时间推进机制、下次事件时间推进机制 或混合时间推进机制时，分别具有哪些优点和缺点，以图形或文字等形式分析时钟推进流程。 11.什么叫仿真效率？什么叫仿真精度？分析影响仿真效率和仿真精度的因素？ 12.从仿真效率和仿真精度的角度，分析和比较三种仿真时钟推进机制的特点，并分析三种 仿真时钟推进机制分别适合于什么样的系统？ 13. 什么是蒲丰投针试验？绘制蒲丰投针试验原理图，通过推导蒲丰投针试验中针与任一直 线相交的概率，分析采用随机投针试验方法来确定圆周率π的原理。 14. 按照蒲丰投针试验的条件和要求，完成投针试验，在统计投针次数、针与直线的相交次 数的基础上，求解π的估计值，并以报表或图形等形式表达试验结果。具体要求如下： ①自行确定针的长度、直线之间的距离。 ②投针10次、20次、30次、40次、50次、…、100次、…、200次、…，分别计算针 与直线相交的概率、π的估计值。 ③以一随机变量描述上述试验结果，并通过编程或采用商品化软件，以图形、报表等形 式表示投针试验结果，分析其中的规律，并给出结论。 ④写出试验报告。 ⑤在熟悉投针试验原理的基础上，编制投针试验仿真程序，动态运行投针试验的过程。15．什么是蒙特卡洛仿真？它有什么特点，蒙特卡洛仿真应用的基本步骤是什么？ 16．采用C或C++等语言，分别编写产生均匀分布、正态分布、指数分布以及威布尔分布的伪随机数序列，通过改变每种分布中参数的数值，分析不同参数数值对随机数值的影响；通过对所产生的伪随机数分布区间的统计、分析和绘图，检验伪随机数的特性及其数值特征。 17. 对于制造系统而言，库存有哪些作用和功能？ 18. 在制造企业中，库存大致可以分成四种类型。简要论述四种库存的名称和功能。 19. 什么是安全库存、订货提前期？确定安全库存和订货提前期时分别需要考虑哪些因素？ 20. 什么叫“订货点法”？要确定订货点，需要哪些条件？订货点法适合于怎样的库存系统？
电力系统分析课程报告姓名 ******* 学院自动化与电气工程学院 专业控制科学与工程 班级 ******* 指导老师 ******* 二〇一六年五月十三
一、同步发电机三相短路仿真 1、仿真模型的建立 选取三相同步发电机模型，以三相视图表示。励磁电压和原动机输入转矩Ef 与Tm均为定常值，且发电机空载。当运行至时，发电机发生三相短路故障。同步发电机三相短路实验仿真模型如图1所示。 图1 同步发电机三相短路实验仿真模型 2、发电机参数对仿真结果的影响及分析 衰减时间常数Ta对于直流分量的影响 三相短路电流的直流分量大小不等，但衰减规律相同，均按指数规律衰减，衰减时间常数为Ta，由定子回路的电阻和等值电感决定（大约）。pscad同步发电机模型衰减时间常数Ta对应位置如图3所示（当前Ta=）。 图3 同步发电机模型参数Ta对应位置
1）Ta=时，直流分量的衰减过程（以励磁电流作为分析）如图4所示。 图4 Ta=发生短路If波形 2）Ta=时，直流分量的衰减过程（以励磁电流作为分析）如图5所示。 图5 Ta=发生短路If波形 短路时刻的不同对短路电流的影响 由于短路电流的直流分量起始值越大，短路电流瞬时值就越大，而直流分量的起始值于短路时刻的电流相位有关，即直流分量是由于短路后电流不能突变而产生的。 Pscad模型中对短路时刻的设置如图6所示 图6 Pscad对于短路时刻的设置 1）当在t=时发生三相短路，三相短路电流波形如图7所示。 图7 t=时三相短路电流波形 2）当在t=时发生三相短路，三相短路电流波形如图8所示。 图8 t=6时三相短路电流波形 Xd、Xd`、Xd``对短路电流的影响 1) Xd的影响 Pscad中对于Xd的设置如图9所示： 图9 Pscad对于D轴同步电抗Xd的设置 下面验证不同Xd时A相短路电流的稳定值。 i.Xd=（标幺制，下同）时，仿真波形如图10所示 图10 Xd=时A相短路电流波形 ii.Xd=10时，仿真波形如图11所示 图11 Xd=时A相短路电流波形 2）Xd`的影响 在Pscad中暂态电抗Xd`的设置如图13所示： 图13 Pscad对于暂态电抗Xd的设置 下面验证不同Xd`时A相短路电流的暂态过程。 i.Xd`=时A相短路电流的波形如图14所示：
、系统、模型和仿真三者之间具有怎样的相互关系 答：系统是研究的对象，模型是系统的抽象，仿真通过对模型的实验以达到研究系统的目的。 、通过因特网查阅有关蒲丰投针实验的文献资料，理解蒙特卡罗方法的基本思想及其应用的一般步骤。 答：蒲丰投针实验内容是这样的：在平面上画有一组间距为a的平行线，将一根长度为L（L （1）实体流图 (2)活动循环图 、以第二章中图2-5所示的并行加工中心系统为对象，建立Petri 网模型。 3214所示Petri 网模型的运行过程，并将分析结果同例3-5相比较。 、任取一整数作为种子值，采用第三题中得到的随机数发生器生成随机数序列的前200项数据，并对其统计性能进行检验。 解：由第3题可得到一个随机数发生器： a=5 b=9 c=3 m=512 取种子值，生成的随机数序列前200项数据如下： n n 5000032458 4 t t P t P P P P t P (2)t3发 生后 t t P t P P P P t P (3)t2发 生后 (4)t1不能 发生 t t P t P P P P t P (5)t4发 生后 《信号与系统》复习题 1． 已知f(t)如图所示，求f(-3t-2)。 2． 已知f(t)，为求f(t0-at)，应按下列哪种运算求得正确结果？（t0和a 都为正值） 3．已知f(5-2t)的波形如图，试画出f(t)的波形。 解题思路：f(5-2t)?? ???→?=倍 展宽乘22/1a f(5-2×2t)= f(5-t) ??→?反转f(5+t)??→?5 右移f(5+t-5)= f(t) 4．计算下列函数值。 （1） dt t t u t t )2(0 0--?+∞ ∞-） （δ （2） dt t t u t t )2(0 --?+∞ ∞ -） （δ （3） dt t t e t ?+∞ ∞ --++） （2)(δ 5．已知离散系统框图，写出差分方程。 解：2个延迟单元为二阶系统，设左边延迟单元输入为x(k) 左○ ∑:x(k)=f(k)-a 0*x(k-2)- a 1*x(k-1)→ x(k)+ a 1*x(k-1)+ a 0*x(k-2)=f(k) (1) 右○ ∑: y(k)= b 2*x(k)- b 0*x(k-2) (2) 为消去x(k)，将y(k)按（1）式移位。 a 1*y(k-1)= b 2* a 1*x(k-1)+ b 0* a 1*x(k-3) (3) a 0*y(k-2)= b 2* a 0*x(k-2)-b 0* a 0*x(k-4) (4) (2)、（3）、（4）三式相加：y(k)+ a 1*y(k-1)+ a 0*y(k-2)= b 2*[x(k)+ a 1*x(k-1)+a 0*x(k-2)]- b 0*[x(k-2)+a 1*x(k-3)+a 0*x(k-4)] ∴ y(k)+ a 1*y(k-1)+ a 0*y(k-2)= b 2*f(k)- b 0*f(k-2)═>差分方程 6．绘出下列系统的仿真框图。 )()()()()(10012 2t e dt d b t e b t r a t r dt d a t r dt d +=++ 7．判断下列系统是否为线性系统。 （2） 8．求下列微分方程描述的系统冲激响应和阶跃响应。 成绩 1． （4）求解差分方程：y(n)-y(n-1)-2y(n-2)=f(n) ?f(n)=( 1 )nε (n) 本科毕业（论文） 题 目 （中、英文 ） in The Signal System 分类 号 学号 密级 公开 学校代码 1107044431 TN911.6 基于MATLAB 的信号系统仿真及应用 The Application of MATLAB in The Signal System 工科 作者姓名 指导教师 学科门类 专业名称 电气工程及其自动化 提交论文日期 成绩评定 二零一五年五月 摘要 当前的科学信息技术正在日新月异的高速发展,而通过应用数字信号处理的方法，已成为一个非常重要的技术手段被广泛应用在通信、音频和图像、遥感，视频等领域。为了更好地了解信号与系统的基本理论和掌握其方法,从而更好地理解和掌握数字信号处理的理论知识,因此在实验过程中我们就需要通过MATLAB 计算机辅助设计平台。 本论文主要探究MATALB在信号与系统中的连续信号和离散信号中的应用，主要从连续和离散两方面入手，进一步掌握信号系统中的相关知识。同时引进计算机软件—MATLAB,对信号系统二阶系统的时域和频域分析，通过它在计算机上对程序进行仿真，阐述信号与系统理论应用与实际相联系。以此激发学习兴趣，变被动接受为主动探知，从而提升学习效果，培养主动思维，学以致用的思维习惯，也可以让人们进一步了解MATLAB软件 关键词：采样定理；MATLAB；信号与系统；抽样定理 Abstract Current, the rapid development of science and information technology are changing and through the application of digital signal processing method, has become a very important technology is widely used in communication, audio and video, remote sensing, video, etc. In order to better understand the basic theory of signal and system, and grasp the method, to better understand and master the theoretical knowledge of digital signal processing, so we need in the process of experiment by MATLAB computer aided design platform. This thesis mainly explores MATALB in signal and system, the application of discrete and continuous signals, mainly from the two aspects of the continuous and discrete, further to master relevant knowledge of signal system. Introduction of computer software - MATAB at the same time, the signal system of second order system time domain and frequency domain analysis, through its d on program on computer simulation, signal and system theory associated with the actual application. To stimulate interest in learning, change passive accept to active detection, so as to improve learning effect, active thinking, to practice habits of thinking, also can let people learn more about MATLAB software. Key words:Sampling theorem; MATLAB; Signals and systems; The sampling theorem 物流系统仿真 期末作业 题目：Manufacturing System Planning and Scheduling 班级：物流工程131 学号：1311393003 1311393008 姓名：黎宇帆张力夫 日期：2015-09-19 成绩： 制造系统规划与调度 翻译 2.1引言 现代生产调度工具是非常强大的，提供了广阔的范围内调整工具的行为的真实过程要求的选项和参数。 然而，更多的选项的存在，它就在实践中找到的工具的最佳配置更加困难。 即专家们经常无法预测的多种可能性的影响。 测试甚至一小部分在现实中可能的配置，对实际生产过程的影响可能需要几个月的时间，可能会严重降低整体性能。 因此，这样的试验在实践中是不可行的。 优化的生产调度仿真模型比使用真正的过程更安全，更便宜,更快，更容易测试。为了在一个中等规模的制造公司充分使用先进的调度工具的优势，找到它的一个最佳的规则和参数的优化配置。 模块化仿真模型的整个业务的制造系统和生产过程中阳极氧化阶段是建立以测试不同的调度配置的影响。调度工具的配置测试和优化进行了离线使用的仿真模型。实际生产过程不受干扰，可以非常快速、低成本的找到最优配置。 2.2问题描述 位于英国的一个中型制造商，生产一系列的不同的小压铝零件和一系列大批量的其他面向消费者的产品。典型的应用包括香水的喷雾组件和哮喘患者的分配器。这是一个高度竞争的行业，成功取决于是否能实现高效率和低成本制造。所以生产调度是非常重要的。 在过去，该公司安装的软件工具可以支持生产过程中的各个区域调度。全面提高公司绩效，增加产量和减少产品的交货时间，他们计划建立自动电抗器的供应链规划服务器–总调度系统协调当地所有的业务和生产区。为了提供最好的解决方案，调度工具供应商，预优国际（https://www.doczj.com/doc/821068077.html,）决定使用模拟求解调度工具的优化配置。 问题是建立一个仿真工具，它将接受的到来客户订单和生产订单排序以满足这些需求。一个重要的地方是模型的生产过程本身，以确保它的主要阶段的最佳时刻加载。阳极氧化阶段是整个生产过程中特别重要的，因此，它必须是非常详细的模拟，以测试到整体订单的交货时间可以通过阳极氧化过程阶段优化减少到什么程度。 在这种情况下的研究主要目标是以下几个： （1）为了确定公司模型间的相关业务和生产过程和确定订单和交货时间， （2）在规划部门分析和优化业务流程，为了处理传入的需求和规划生产订单。 （3）测试的整体生产时间，提高灵敏度，特别是确定是否引入特定排序规则的生产订单将减少在阳极氧化处理阶段总的处理时间。 一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环 3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置 M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置 实验一 信号与系统仿真实验 希望同学们根据实验任务要求事先做好预习，上机实验完成后应写出相应的实验报告(要求附程序与仿真结果)。 一、 实验目的 了解MA TLAB 的基本使用方法和编程技术，以及Simulink 平台的建模与动态仿真方法 ，进一步加深对课程所学内容的理解。 二、 实验项目 1．信号的分解与合成，观察Gibbs 现象。 2．信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 3．信号的频谱分析，观察信号的频谱波形。 4．系统函数的形式转换。 5．用Simulink 平台对系统进行建模和动态仿真。 三、 实验仪器 计算机一人一台；安装Matlab/Simulink 数值仿真软件平台。 四、 实验内容 1、以周期为T ，脉冲宽度为12T 的周期性矩形脉冲为例研究Gibbs 现象。 提示：已知周期方波信号的傅里叶级数系数a k 的表达式如下： π ωπωωk T k a T T a T e a t x k m m k t jk k )sin(22)(101000==== ∑-= 试画出x (t )的波形图（分别取m 等于1，3，7，19，79，T =4T 1），观察Gibbs 现象，通过对不同m 取值的合成波形观察，体会有限项合成信号与原信号的不同，同时，理解函数能量大部分集中在傅里叶级数系数a k 的第一对零点之内的道理 2、求卷积并画图 （1）已知：)2()1()(1---=t u t u t x ，)3()2()(2---=t u t u t x 求：)()()(21t x t x t y *==？并画出其波形。 （2）已知某离散系统的输入和冲激响应分别为：]5,3,2,1,5,3,4,1[][=n x ，]2,4,0,4,2,4[][=n h 。求系统的零状态响应，并绘制出系统的响应图。 提示：求卷积可用),(21x x conv ；画图可用subplot 、plot 和stem 。 3、求频谱并画图 实验指导书 实验1 MATLAB基本操作 一、实验目的 1．熟悉MATLAB实验环境，练习MATLAB命令、m文件、Simulink的基本操作。 2．利用MATLAB编写程序进行矩阵运算、图形绘制、数据处理等。 3．利用Simulink建立系统的数学模型并仿真求解。 二、实验原理 MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。MATLAB有3种窗口，即：命令窗口（The Command Window）、m-文件编辑窗口（The Edit Window）和图形窗口（The Figure Window）。Simulink 是MATLAB的一个部件，它为MATLAB用户提供了一种有效的对反馈控制系统进行建模、仿真和分析的方式。而Simulink另外又有Simulink模型编辑窗口。 1．命令窗口（The Command Window） 当MATLAB启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“>>”后面输入交互的命令，这些命令就立即被执行。 在MATLAB中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m”为后缀，所以称为m文件。 2．m文件编辑窗口（The Edit Window） 我们可以用m文件编辑窗口来产生新的m文件，或者编辑已经存在的m文件。在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m文件编辑窗口；选择菜单“File/Open”就可以打开一个已经存在的m文件，并且可以在这个窗口中编辑这个m文件。 3．图形窗口（The Figure Window） 图形窗口用来显示MATLAB程序产生的图形。图形可以是2维的、3维的数据图形，也可以是照片等。 4．Simulink的启动方式： 启动运行Simulink有两种方式：（1）在Command window中，键入simulink，回车。（2）单击工具栏上Simulink图标。启动Simulink后， 系统仿真II 大作业 自06A-2 赵众源 06101010215 1、考虑如下的微分方程： （1）试在Matlab 环境下采用Euler 法编制其仿真程序，给出在[0，10]秒区间上y 的变化曲线。 （2）试在Matlab —Simulink 环境下搭建仿真模型，给出在[0，10]秒区间上y 的变化曲线。 （3）考虑带有输入的情况： 在Matlab —Simulink 环境下进行仿真，给出在[0，10]区间上y 的变化曲线。 （1） 在M 文件中创建Euler 函数： function [t,y]=euler(odefun,tspan,y0,h) t=tspan(1):h:tspan(2);y(1)=y0; for i=1:length(t)-1 y(i+1)=y(i)+h*feval(odefun,t(i),y(i)); end t=t';y=y'; 在命令窗口得到y 的变化曲线： odefun=inline('2*y*sin(t)-abs(y^3)','t','y'); >> [t,y]=euler(odefun,[0,10],1,0.01); plot(t,y) 3 ()(,)2()sin |()|,(0)1dy t f t y y t t y t y dt ==-=3 ()(,)2()sin |()|(),(0)1,()2cos 2dy t f t y y t t y t u t y u t t dt ==-+== 0.2 0.4 0.6 0.8 1 1.2 1.4 （2）试在Matlab —Simulink 环境下搭建仿真模型，给出在[0，10]秒区间上y 的变化曲线。 （3）考虑带有输入的情况： 在Matlab —Simulink 环境下进行仿真，给出在[0，10]区间上y 的变化曲线。 3 ()(,)2()sin |()|(),(0)1,()2cos 2dy t f t y y t t y t u t y u t t dt ==-+== 仿真计算 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的量，还有分析不对称故障时对称分量法转化到相对简单的对称故障分析中。 ) t t )]( 2 3 --t δt 解 t t t 1.27 （a ）)2()2()(t x t x t y -+-= ① 因为)2()2()0(x x y +-= ，在0=t 的输出与前后时刻的输入都有关，所以系统是记 忆的。 ② 已知)2()2()(111t x t x t y -+-=，)2()2()(222t x t x t y -+-=。当 )()(012t t x t x -=时， )2()2()(01012t t x t t x t y --+--=，而)2()2()(010101t t x t t x t t y +-+--=-， 所以：)()(012t t y t y -≠。因而系统是时变的。 ③已知)2()2()(111t x t x t y -+-=，)2()2()(222t x t x t y -+-=， )2()2()(333t x t x t y -+-=， 当)()()(213t x t x t x +=时，)]2()2([)]2()2([)(21213t x t x t x t x t y -+-+-+-= 所以)()()(213t y t y t y +=，因而系统是可加的。 当)()(12t ax t x =时，)()2()2()(1112t ay t ax t ax t y =-+-=，因而系统是齐次的。 综合系统的可加性与齐次性，所以系统是线性的。 ④因为)2()2()0(x x y +-= ，在0=t 的输出与2=t 的输入也有关，所以系统是非因果 的。 ⑤若+∞ <≤B t x ) (，即输入有界，则： +∞ <≤-+-≤-+-=B t x t x t x t x t y 2)2()2()2()2()(，即输出有界。 所以系统是稳定的。 （d ） ?? ?≥-+<=0),2()(0 ,0)(t t x t x t t y ①)2()0()0(-+=x x y ，即)0(y 与2,0-==t t 的输入有关， ∴系统是记忆系统。 ②令) ()(01t t x t x -= ，则 ?? ?≥--+-<=???≥-+<=0),2()(0,00),2()(0 ,0)(00111t t t x t t x t t t x t x t t y 而 ?? ?≥--+-<=???≥---+-<-=-0000 00000),2()(,00),2()(0,0)(t t t t x t t x t t t t t t x t t x t t t t y ) ()(01t t y t y -≠∴，系统是时变的。 ③令)()(1t ax t x =，则) (0),2()(0 ,0)(1t ay t t ax t ax t t y =???≥-+<=，所以，系统是齐次的。 已知 ???≥-+<=0),2()(0,0)(111t t x t x t t y ，?? ?≥-+<=0),2()(0 ,0)(222t t x t x t t y信号与系统习题答案(教学参考)
信号与系统的MATLAB仿真
课程设计说明书(计算书、论文)
题 目 信号与系统的 MATLAB 仿真
课 程 名 称 院 （系）
信号与系统 电子通信工程学院
专 业 班 级 学 生 姓 名 学 号
设 计 地 点 指 导 教 师
设计起止时间：
年
月 日
至
年
月
日
课程设计应达到的目的
(1)熟悉 Matlab 软件的运行环境 (2)掌握采用 Matlab 软件程序实现信号与系统分析的方法 (3)掌握正确的编程过程和仿真分析 (4)总结对比软件仿真与硬件实验的区别及特点 2．课程设计题目及要求 《信号与系统》课程设计选题主要是要体现本课程的主要教学 内容中的重点部分，同时要求选题能过反映出信号仿真的代表性， 系统分析的应用性， 灵活性， 并且能与原本理论教学中繁琐的数学 计算相比较， 体现出软件计算的方便快捷性， 本课程设计主要包括 四个小设计部分，分别是： (1)信号的产生与简单运算：产生一个方波周期为 4π ，t[0 50]。
(2)?求解微分方程：y"(t)+3y'(t)+2y(t)=2e-2 ε(t)求 yzs； ?求卷积：e-2 ε (t)*e-3 ε (t)
t t
t
(3)求 H (s) ?
2s 2 ? 1 s 3 ? 4s 2 ? 6s ? 9
?求零、极点 ?并绘图 ?冲激响应
3
?f(n)=δ (n)
3.课程设计思路 利用信号与系统中的 matlab 常用命令集求解微分方程，并利用结 果和绘图命令绘图。
4.课程设计原理 设计原理 (1)设计一个简单程序能实现方波信号的生成。 利用Matlab软件的信号处理工具箱(Signal Processing Toolbox)中的专用函数产生 信号并绘出波形。
（2） ?对于求方程的零状态响应，即是求解常微分方程。Matlab 解常微分方程式的语法是 dsolve('equation','condition')，其中equation代表常微分方程式即 y'=g(x,y)， 且须以Dy代表一 微分项y''，condition则为初始条件。 ?利用MATLAB中conv命令求解卷积。 阶微分项y' D2y代表二阶基于MATLAB的信号与系统仿真及应用
物流仿真大作业.doc
matlab机电系统仿真大作业
实验一信号与系统仿真实验
信号与系统仿真实验指导书
系统仿真II大作业
PSCAD的电力系统仿真大作业3
信号与系统作业答案
相关主题
文本预览