仿真作业2013-GPSS

仿真作业一、（排队系统仿真）

1．某小邮电所的顾客到达时间间隔为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/Simulink 仿真实验 实验一 过程控制系统建模 ............................................................................................................. 1 实验二 PID 控制 ............................................................................................................................. 2 实验三 串级控制 ............................................................................................................................. 6 实验四 比值控制 ........................................................................................................................... 13 实验五 解耦控制系统 . (19) 实验一 过程控制系统建模 指导内容：（略） 作业题目一： 常见的工业过程动态特性的类型有哪几种？通常的模型都有哪些？在Simulink 中建立相应模型，并求单位阶跃响应曲线。 作业题目二： 某二阶系统的模型为2 () 22 2n G s s s n n ?ζ??= ++，二阶系统的性能主要取决于ζ，n ?两个参数。试利用Simulink 仿真两个参数的变化对二阶系统输出响应的影响，加深对二阶 系统的理解，分别进行下列仿真： （1）2n ?=不变时，ζ分别为0.1, 0.8, 1.0, 2.0时的单位阶跃响应曲线； （2）0.8ζ=不变时，n ?分别为2, 5, 8, 10时的单位阶跃响应曲线。

信号与系统的MATLAB仿真

信号与系统的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)'); % 加注波形标题 运行结果：

信号与系统

信号与系统 单项选择题 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)

《机械系统动力学仿真分析软件》

| 论坛社区 《机械系统动力学仿真分析软件》(MSC.ADAMS.2005.R2)R2 资源分类： 软件/行业软件 发布者： Coolload 发布时间： 2005-12-18 20:22 最新更新时间： 2005-12-19 07:04 浏览次数： 14548 实用链接： 收藏此页 eMule资源 下面是用户共享的文件列表，安装eMule后，您可以点击这些文件名进行下载 [机械系统动力学仿真分析软件].[$u]MSC.ADAMS.2005.R2.rar201.2MB [机械系统动力学仿真分析软 295.4MB 件].MSC_ADAMS_V2005_ISO-LND-CD1.iso [机械系统动力学仿真分析软185.0MB

件].MSC_ADAMS_V2005_ISO-LND-CD2.bin [机械系统动力学仿真分析软 6.5KB 件].Msc.Adams.v2005.Iso-Lnd-Cd1-Crack.rar 全选480.4MB eMule主页下载eMule使用指南如何发布 中文名称：机械系统动力学仿真分析 软件 英文名称：MSC.ADAMS.2005.R2 版本：R2 发行时间：2005年12月15日 制作发行：美国MSC公司 地区：美国 语言：英语 简介： [通过安全测试] 杀毒软件：Symantec AntiVirus 版本： 9.0.0.338 病毒库：2005-12-16 共享时间：10:00 AM - 24:00 PM(除 非线路故障或者机器故障) 共享服务器：Razorback 2.0 [通过安装测试]Windows2000 SP4 软件版权归原作者及原软件公司所 有，如果你喜欢，请购买正版软件

信号与系统习题答案(教学参考)

《信号与系统》复习题 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．求下列微分方程描述的系统冲激响应和阶跃响应。

运动控制系统仿真作业

运动控制系统仿真作业 利用Matlab解运动控制系统习题 习题2-5在转速、电流双闭环调速系统中，两个调节器均采用PI调节器。当系统带额定负载运行时，转速反馈线突然断线，系统重新进入稳态后，电流调节器的输入偏差电压是否为零？为什么？ 解：（一）结合电流、转速调节器的设计建立转速、电流双闭环调速系统模型。设有某晶闸管供电的双闭环直流调速系统，整流装置采用三相桥式电路，基本数据如下：直流电动机：220V，136A，1460r/min，e C=0.132V2min/r，允 许过载倍数λ=1.5； 晶闸管装置放大系数s K=40； 电枢回路总电阻R=0.5Ω； 时间常数l T=0.03s，m T=0.18s； 电流反馈系数β=0.05V/A(≈10V/1.5N I)； 转速反馈系数α=0.007V2min/r(≈10V/N n)。 设计要求：设计电流调节器，要求电流超调量5%iσ=。设计转速调节器，要求转速无静差，空载起动到额定转速时的转速超调量10%iσ=，并检验转速超调量的要求能否得到满足。 1.设计电流调节器 1）确定时间常数 ①整流装置滞后时间常数s T。三相桥式电路的平均失控时间s

T=0.0017s。②电流滤波时间常数oi T。取oi T=0.002s。 ③电流环小时间常数之和￡i T。按小时间常数近似处理，取￡i s oi T T T=+=0.0037s。 2）选择电流调节器结构 根据设计要求10%iσ=，并保证稳态电流无差，可按典型I型系统设计电流调节器。电流环控制对象是双惯性的，因此可用PI型电流调节器，其传递函数为 (1)()i i ACR i K s W s s ττ+=检查对电源电压的抗扰性能： ￡i l T T=0.030.0037s s=8.11，由表1可知，各项指标都是可以接受的。 电流调节器超前时间常数：i l Tτ==0.03s。 电流环开环增益：要求10%iσ=时，根据表2可知，￡i I K T =0.5，因此 1￡i0.50.5135.10.0037I K s T s -===于是，ACR的比例系数为 ￡i135.10.030.5 1.013400.05 I i i K R K Tτ??===?4）校验近似条件

信号与系统的MATLAB仿真

成绩
课程设计说明书(计算书、论文)
题 目 信号与系统的 MATLAB 仿真
课 程 名 称 院 （系）
信号与系统 电子通信工程学院
专 业 班 级 学 生 姓 名 学 号
设 计 地 点 指 导 教 师
设计起止时间：
年
月 日
至
年
月
日

1．
课程设计应达到的目的
(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
?求零、极点 ?并绘图 ?冲激响应

（4）求解差分方程：y(n)-y(n-1)-2y(n-2)=f(n) ?f(n)=( 1 )nε (n)
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代表二阶

车辆系统动力学仿真大作业(带程序)

Assignment Vehicle system dynamics simulation 学院：机电学院 专业：机械工程及自动化 姓名： 指导教师：

The model we are going to analys: The FBD of the suspension system is shown as follow:

According to the New's second Law, we can get the equation: 2 )()(221211mg z z c z z k z m --+-=???? 221212)()(z k mg z z c z z k z m w +-----=? ??? 0)()()()(222111222111=-++--+-++--+? ? ? ? ? ? ? ?w w w w z L z k z L z k z L z c z L z c z m χχχχ 0)()()()(2222111122221111=-++----++---? ? ? ? ? ? ? ?w w w w z L z L k z L z L k z L z L c z L z L c J χχχχχ d w w w w Q z L z k z L z c z m ,111111111)()(-=------? ? ? ? ?χχ d w w w w Q z L z k z L z c z m ,222222222)()(-=-+--+-? ????χχ When there is no excitation we can get the equation: 2)()(221211mg z z c z z k z m --+-=???? 2 21212)()(z k mg z z c z z k z m w +-----=? ??? Then we substitude the data into the equation, we write a procedure to simulate the system: Date: ???? ?? ??? ??==?==?===MN/m 0.10k m 25.1s/m kN 0.20MN/m 0.1m kg 3020kg 2100kg 3250w 2l c k I m m by w b

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

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

目录 一、引言 (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的信号与系统仿真及应用

本科毕业（论文） 题 目 （中、英文 ） 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

控制系统计算机仿真作业

兰州理工大学 《控制系统计算机仿真》 上机报告Ⅰ 院系：电气工程与信息工程学院 班级：14级自动化3班 姓名：孙悦 学号：1405220323 时间：2017年6月15日 电气工程与信息工程学院

《控制系统计算机仿真》上机实验任务书Ⅰ（2017） 一、上机实验内容及要求 1.matlab软件 要求利用课余时间熟悉掌握matlab软件的基本数值运算、基本符号运算、基本程序设计方法及常用的图形命令操作。 2.各章节仿真实验内容及要求 具体实验内容及要求请详见上机实验报告。 二、上机实验时间安排及相关事宜 1.依据课程教学大纲要求，上机实验学时共16学时，学生须在每次上机之前 做好相应的准备工作，以确保在有限的机时内完成仿真实验要求的内容； 2.实验完成后按规定完成相关的仿真实验报告； 3.仿真实验报告请按有关样本制作并A4打印，侧面装订，作为成绩评定的一 部分。 自动化系《控制系统计算机仿真》课程组 2017年3月

一、Matlab 基础操作 1-1用MATLAB 语言求下列系统的状态方程、传递函数、零极点增益和部分分式形式的模型参数，并分别写出其相应的数学模型表达式： （1）2450351024 247)(234 23+++++++=s s s s s s s s G 程序如下: num=[7,24,24] den=[10,35,50,24] [A,B,C,D]=tf2ss(num,den) 系统的状态方程： A = -3.5000 -5.0000 -2.4000 1.0000 0 0 0 1.0000 0 B = 1 0 0 C = 0.7000 2.4000 2.4000 D = 零极点增益形式： [Z,P,K]=tf2zp(num,den) Z = -1.7143 + 0.6999i -1.7143 - 0.6999i

实验一信号与系统仿真实验

实验一 信号与系统仿真实验 希望同学们根据实验任务要求事先做好预习，上机实验完成后应写出相应的实验报告(要求附程序与仿真结果)。 一、 实验目的 了解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后，

《控制系统仿真》课程论文模板

基于全转速范围内的直接转矩控制调速系统的研究摘要：建立了一种包含弱磁控制的直接转矩调速控制系统仿真模型，实现了感应电动机全转速范围内的速度控制，既可以实现同步转速以下的恒转矩调速，又可以实现弱磁范围内的恒功率调节。仿真结果显示系统的调速性能良好，在换向调速过程中，无抖动和超调现象，很好地实现了电动机调速的四象限特性，弱磁控制过程平稳。对研究全转速范围内直接转矩调速控制具有较好的参考价值。 关键词：直接转矩控制；调速；弱磁；仿真 Research of Direct Torque Control Speed Adjustment System within Full-speed Range Abstract: A simulation system of direct torque control speed adjustment system is introduced, it include field-weakening control. The simulation system realized induction-motor speed control within full-speed range, constant torque adjustment less than foundation speed and constant power adjustment greater than foundation speed. The results of simulation show that system performance is excellent. The ripple and overshoot of speed are reduced; the four-quadrant characteristic of induction-motor varying speed is achieved successfully; the field-weakening control is very smooth. The simulation model has better reference meaning for researching direct torque control system within full-speed range. Keyword: direct torque control, adjustment speed, flux weakening, simulation 1．引言 直接转矩控制系统具有控制结构简单、动态响应快等特点，它在很大程度上解决了矢量控制中计算复杂、特性易受电动机参数变化的影响、实际性能难以达到理论分析结果的一些重要技术问题[1]。直接转矩控制采用双滞环控制策略，电机在运行中转矩及定子磁链脉动较大，影响了电机运行的稳定性，所以研究的焦点大都集中在减小转矩及定子磁链脉动上[2-3]。为了提高定子磁链的估计精度，改善直接转矩控制系统的动静态性能等，提出了多种行之有效的减小脉动方法[4-6]。在直接转矩的调速控制方面，有同步转速以下恒转矩调速的研究[7]，也有基于弱磁范围内的速度控制等[8]，但基于全速度范围内的控制仿真研究还不多见[9，10]。本文建立了一种直接转矩控制的感应电动机变频调速系统仿真模型，可以实现电动机全转速范围内的速度控制。所谓全转速范围，是指电动机调速系统既可以实现同步转速以下的恒转矩调速，又可以实现弱磁范围内的恒功率调节，而且能实现转速的四象限运行特性。 2．直接转矩控制调速系统的建模 直接转矩控制（Direct Torque Control，DTC），是基于动态模型的双闭环控制系统。外环采用转速闭环，控制转速的大小；内环采用磁链和转矩滞环比较控制，产生转矩的快速动态响应。 2.1 磁链与转矩控制 要进行磁链和转矩滞环的闭环控制，必须取得异步电动机的定子磁链和电磁转矩，磁链观测和转矩观测通常采用的数学表达 式为： dt R t i t u t s s s s ) )( )( ( )(- =? ψ（1） ) ( 2 3 β α α β ψ ψ s s s s e i i P T- =（2） 式中：s(t)为定子磁链，u s(t)为定子相电压，i s(t)为定子相电流，R s为定子电阻，

信号与系统作业答案

) 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

控制系统仿真大作业

控制系统仿真 实验报告 专业班级：自动F0903 姓名：罗新勇 学号： 200948280311 指导教师：张杰

实验一、熟悉MATLAB 环境及矩阵、数组的数学计 算 一、 实验目的 1、熟悉启动和退出Matlab 的方法； 2、熟悉Matlab 命令窗口的组成； 3、掌握建立矩阵的方法； 二、 实验内容： 1、帮助命令 使用help 命令，查找 sqrt （开方）函数的使用方法； 2、先求下列表达式的值，然后显示Matlab 工作空间的使用情况并保存全部变量。 .3,9.2,8.2,...,8.2,9.2,0.3,2 3.0ln )3.0sin(2 )3(545.0212),1log(21)2(185sin 2)1(3.03.032 220 1---=+++-= ?? ? ? ??-+=++=+=-a a a e e z i x x x z e z a a 其中 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

?? ???=<≤+-<≤-<≤=5.2:5.0:0,3 2,1221, 110,)4(22 2 4t t t t t t t t z 其中 提示：用逻辑表达式求分段函数值。 （1）z1=2*sin(85/180*pi)/(1+(exp(1))^2) z1 = 0.2375 （2）x=[2 1+2i;-0.45 5] x = 2.0000 1.0000 + 2.0000i -0.4500 5.0000 z2=0.5*log(x+sqrt(1+x^2)) z2 = 0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i （3）a=-3.0:0.1:3.0 a = Columns 1 through 5 -3.0000 -2.9000 -2.8000 -2.7000 -2.6000 Columns 6 through 10 -2.5000 -2.4000 -2.3000 -2.2000 -2.1000 Columns 11 through 15 -2.0000 -1.9000 -1.8000 -1.7000 -1.6000 Columns 16 through 20 -1.5000 -1.4000 -1.3000 -1.2000 -1.1000 Columns 21 through 25 -1.0000 -0.9000 -0.8000 -0.7000 -0.6000 Columns 26 through 30 -0.5000 -0.4000 -0.3000 -0.2000 -0.1000 Columns 31 through 35 0 0.1000 0.2000 0.3000 0.4000 Columns 36 through 40 0.5000 0.6000 0.7000 0.8000 0.9000 Columns 41 through 45 1.0000 1.1000 1.2000 1.3000 1.4000 Columns 46 through 50 1.5000 1.6000 1.7000 1.8000 1.9000 Columns 51 through 55 2.0000 2.1000 2.2000 2.3000 2.4000 Columns 56 through 60 2.5000 2.6000 2.7000 2.8000 2.9000 Column 61