基于MatlabSimulink建模与仿真关系

基于MatlabSimulink建模与仿真关系
基于MatlabSimulink建模与仿真关系

基于Matlab/Simulink建模与仿真关系

摘要:用MATLAB/Simulink来设计数字逻辑电路给设计者开辟了一种全新的设计方法的可能性,使用Simulink模块进行以代替传统的原理图进行数字逻辑电路的分析可以使初学者能更形象直观地看到电路运行时的波形变化,可以大大提高学习效率。故利用MATLAB仿真软件在日常的数字电路的研究学习中有着越来越重要的作用。

用MATLAB/Simulink来设计数字逻辑电路给设计者开辟了一种全新的设计方法的可能性。使用Simulink模块进行设计、调试、和仿真,以代替传统的原理图进行数字逻辑电路的分析。使初学者能更形象直观地看到电路运行时的波形变化,可以大大提高学习效率。

仿真结果表明,利用Matlab进行数字电路的设计、调试,结果直观、省时省力。它不仅能用来仿真本论文中提到的电路,而且能广泛地应用于其它数字电路的仿真,是数字电路设计、调试的有效工具。关键词:数字逻辑电路;仿真;MATLAB、Simulink;

引言:TLAB的名称源自Matrix Laboratory,它是一种科学计算软件,专门以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作。Simulink 是基于MATLAB的框图设计环境,可以用来对各种动态系统进行建模、分析和仿真,它的建模范围广泛,可以针对任何能够用数学来描述的系统进行建模。

MATLAB程序设计

MATLAB有两种工作方式:一种是交互式的命令行工作方式;另一种是M文件的程序工作方式。在前一种工作方式下,MATLAB被当做一种高级数学演算纸和图形表现器来使用,MATLAB提供了一套完整的而易于使用的编程语言,为用户提供了二次开发的工具,下面主要介绍MATLAB控制语句和程序设计的基本方法。

Matlab是当今最优秀的科技应用软件之一,它以强大的科学计算与可视化功能、简单易用、开放式可扩展环境,特别是所附带的30多种面向不同领域的工具箱支持,使得它在许多科学领域中成为计算机辅助设计和分析、算法研究和应用开发的基本工具和首选平台。Matlab具有其他高级语言难以比拟的一些优点,编写简单,编程效率高,易学易懂,因此Matlab语言也被通俗地称为演算纸式的科学算法语言。在控制、通信、信号处理及科学计算等领域中,Matlab 都被广泛地应用,已经被认可为能够有效提高工作效率、改善设计手段的工具软件。

Simulink仿真

Simulink是一种图形化仿真软件, 是基于MATLAB语言环境下实现动态系统建模、仿真的一个集成环境, 其中存储了大量系统模型。调用模块库中的各个模块, 进行有关信号连接, 就可完成所需模型的仿真。Simulink的组成很有层次性, 它的下面有很多子库, 如连续子库、非连续子库、离散子库、数学函数子库、信号源子库、信号输出子库及附加模块库等, 各个子库下面含有该类型的众多模块,

复制这些模块到模型窗口中就可实现各种仿真。可按照其有关性质在相应的子库中寻找所需模块。

Simulink是MATLAB软件的扩展,它是实现动态系统建模和仿真的一个软件包,它与MATLAB语言的主要区别在于,其与用户交互接口是基于Windows的模型化图形输入,其结果是使得用户可以把更多的精力投入到系统模型的构建,而非语言的编程上。

使用Simulink进行仿真一般分为两步:首先需要在仿真模型编辑窗口中搭建好自己的仿真模型,设置好具体模型参数和仿真参数;然后,就可以开始仿真,Simulink将根据用户搭建的模型,模拟系统在设定环境下的具体行为。

组合逻辑电路仿真

组合逻辑电路由最基本的“与门”电路、“或门电路”和“非门”电路组成,其输出值仅依赖于其输入变量的当前值,

与输入变量的过去值无关—即不具记忆和存储功能。常见的组合逻辑电路有数据选择器、编码器、译码器、数据分配器、数据比较器等典型的中规模组合逻辑器件。组合逻辑电路的结构具有以下特点:输出输入之间没有反馈延迟通路;电路中不含有记忆功能的元件。

时序逻辑电路也由上述基本逻辑门电路组成,但存在反馈回路—它的输出值不仅依赖于输入变量的当前值,也依赖于输入变量的过去值。常见的时序电路有触发器、锁存器等中规模时序器件。由于只分高、低电平,抗干扰力强,精度和保密性佳。广泛应用于计算机、数字控制、通信、自动化和仪表等方面。最基本的有与电路、或电路和

非电路。

收获、体会与建议:

利用MATLAB实现时序逻辑电路的设计与仿真,从仿真结果可以看出它们均可以达到技术指标要求,而且方法简单、快捷,大大减轻了工作量。我对本次设计的结果还比较满意,经过一周的设计、分析,终于圆满的完成此次课程设计。虽然在这个过程中也遇到一些困难,但通过与同学的交流和自己查资料,我对MATLAB有了更深的体会。通过本次设计我感觉到自己对MATLAB的Simulink仿真不是很熟悉,在用到所需元件时都要进行查阅,在今后要进一步增加自己的MATLAB 学习,提高自己MATLAB汇编技能。

引用:【薛定宇. 基于MATLAB/Simulink的系统仿真技术与应用[M]. 清华大学出版社, 2002.】

【陈桂明. 应用MATLAB建模与仿真[M]. 科学出版社, 2001】.

【李谨, 邓卫华. AMESim与MATLAB/Simulink联合仿真技术及应用[J]. 指挥控制与仿真, 2004, 26(5):61-64.】

现代信号处理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 σσ-=--=-+

系统建模与仿真

一、基本概念 1、数字正弦载波调制 在通信中不少信道不能直接传送基带信号,必须用基带信号对载波波形的某些参量进行控制,使得载波的这些参量随基带信号的变化而变化,即所谓数字正弦载波调制。 2、数字正弦载波调制的分类。 在二进制时, 数字正弦载波调制可以分为振幅键控(ASK)、移频键控(FSK)和移相键控(PSK)三种基本信号形式。如黑板所示。 2、高斯白噪声信道 二、实验原理 1、实验系统组成 2、实验系统结构框图

图 1 2FSK信号在高斯白噪声信道中传输模拟框图 各个模块介绍p12 3、仿真程序 x=0:15;% x表示信噪比 y=x;% y表示信号的误比特率,它的长度与x相同FrequencySeparation=24000;% BFSK调制的频率间隔等于24KHz BitRate=10000;% 信源产生信号的bit率等于10kbit/s SimulationTime=10;% 仿真时间设置为10秒SamplesPerSymbol=2;% BFSK调制信号每个符号的抽样数等于2 for i=1:length(x)% 循环执行仿真程序 SNR=x(i);% 信道的信噪比依次取中的元素 sim('project_1');% 运行仿真程序得到的误比特率保存在工作区变量BitErrorRate中 y(i)=mean(BitErrorRate); end hold off% 准备一个空白的图 semilogy(x,y);%绘制的关系曲线图,纵坐标采用对数坐标 三、实验结论

图 4 2FSK信号误比特率与信噪比的关系曲线图 系统建模与仿真(二) ——BFSK在多径瑞利衰落信道中的传输性能 一、基本概念 多径瑞利衰落信道 二、实验原理 1、实验系统组成

Matlab Simulink 仿真步骤

MATLAB基础与应用简明教程 张明等编著 北京航空航天大学出版社(2001.01) MATLAB软件环境是美国New Mexico大学的Cleve Moler博士首创的,全名为MATrix LABoratory(矩阵实验室)。它建立在20世纪七八十年代流行的LINPACK(线性代数计算)和ESPACK(特征值计算)软件包的基础上。LINPACK和ESPACK软件包是从Fortran语言开始编写的,后来改写为C语言,改造过程中较为复杂,使用不便。MA TLAB是随着Windows环境的发展而迅速发展起来的。它充分利用了Windows环境下的交互性、多任务功能语言,使得矩阵计算、数值运算变得极为简单。MA TLAB语言是一种更为抽象的高级计算机语言,既有与C语言等同的一面,又更为接近人的抽象思维,便于学习和编程。同时,它具有很好的开放性,用户可以根据自己的需求,利用MA TLAB提供的基本工具,灵活地编制和开发自己的程序,开创新的应用。 本书重点介绍了MA TLAB的矩阵运算、符号运算、图形功能、控制系统分析与设计、SimuLink仿真等方面的内容。 Chap1 MATLAB入门与基本运算 本章介绍MATLAB的基本概念,包括工作空间;目录、路径和文件的管理方式;帮助和例题演示功能等。重点介绍矩阵、数组和函数的运算规则、命令形式,并列举了可能得到的结果。由于MA TLAB的符号工具箱是一个重要分支,其强大的运算功能在科技领域有特殊的帮助作用。 1.1 MATLAB环境与文件管理 1.2 工作空间与变量管理 1.2.1 建立数据 x1=[0.2 1.11 3]; y1=[1 2 3;4 5 6]建立一维数组x1和二维矩阵y1。分号“;”表示不显示定义的数据。 MATLAB还提供了一些简洁方式,能有规律地产生数组: xx=1:10 %xx从1到10,间隔为1 xx=-2:0.5:1 %xx从-2到1,间隔为0.5 linespace命令等距离产生数组,logspace在对数空间中等距离产生数组。对于这一类命令,只要给出数组的两端数据和维数就可以了。 xx=linespace(d1,d2,n) %表示xx从d1到d2等距离取n个点 xx=logspace(d1,d2,n) %表明xx从10d1到10d2等距离取n个点 1.2.2 who和whos命令 who: 查看工作空间中有哪些变量名 whos: 了解这些变量的具体细节 1.2.3 exist命令 查询当前的工作空间内是否存在一个变量,可以调用exist()函数来完成。 调用格式:i=exist(…A?); 式中,A为要查询的变量名。返回的值i表示A存在的形式: i=1 表示当前工作空间内存在一个变量名为A的矩阵; i=2 表示存在一个名为A.m的文件; i=3 表示MATLAB的工作路径下存在一个名为A.mex的文件;

物流仿真大作业.doc

物流系统仿真 期末作业 题目:Manufacturing System Planning and Scheduling 班级:物流工程131 学号:1311393003 1311393008 姓名:黎宇帆张力夫 日期:2015-09-19 成绩:

制造系统规划与调度 翻译 2.1引言 现代生产调度工具是非常强大的,提供了广阔的范围内调整工具的行为的真实过程要求的选项和参数。 然而,更多的选项的存在,它就在实践中找到的工具的最佳配置更加困难。 即专家们经常无法预测的多种可能性的影响。 测试甚至一小部分在现实中可能的配置,对实际生产过程的影响可能需要几个月的时间,可能会严重降低整体性能。 因此,这样的试验在实践中是不可行的。 优化的生产调度仿真模型比使用真正的过程更安全,更便宜,更快,更容易测试。为了在一个中等规模的制造公司充分使用先进的调度工具的优势,找到它的一个最佳的规则和参数的优化配置。 模块化仿真模型的整个业务的制造系统和生产过程中阳极氧化阶段是建立以测试不同的调度配置的影响。调度工具的配置测试和优化进行了离线使用的仿真模型。实际生产过程不受干扰,可以非常快速、低成本的找到最优配置。 2.2问题描述 位于英国的一个中型制造商,生产一系列的不同的小压铝零件和一系列大批量的其他面向消费者的产品。典型的应用包括香水的喷雾组件和哮喘患者的分配器。这是一个高度竞争的行业,成功取决于是否能实现高效率和低成本制造。所以生产调度是非常重要的。 在过去,该公司安装的软件工具可以支持生产过程中的各个区域调度。全面提高公司绩效,增加产量和减少产品的交货时间,他们计划建立自动电抗器的供应链规划服务器–总调度系统协调当地所有的业务和生产区。为了提供最好的解决方案,调度工具供应商,预优国际(https://www.360docs.net/doc/893627067.html,)决定使用模拟求解调度工具的优化配置。 问题是建立一个仿真工具,它将接受的到来客户订单和生产订单排序以满足这些需求。一个重要的地方是模型的生产过程本身,以确保它的主要阶段的最佳时刻加载。阳极氧化阶段是整个生产过程中特别重要的,因此,它必须是非常详细的模拟,以测试到整体订单的交货时间可以通过阳极氧化过程阶段优化减少到什么程度。 在这种情况下的研究主要目标是以下几个: (1)为了确定公司模型间的相关业务和生产过程和确定订单和交货时间, (2)在规划部门分析和优化业务流程,为了处理传入的需求和规划生产订单。 (3)测试的整体生产时间,提高灵敏度,特别是确定是否引入特定排序规则的生产订单将减少在阳极氧化处理阶段总的处理时间。

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机电系统仿真大作业

一曲柄滑块机构运动学仿真 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初始值设置

《生产系统建模与仿真》教学大纲

《生产系统建模与仿真》教学大纲 (理论课程) 开课系(部):工程学院课程编号:010396 课程类型:专业课总学时:48 学分:3 适用专业:工业工程开课学期:2014-2015学年第一学期 先修课程:概率论与数理统计、C语言程序设计、系统工程导论 一、课程简述 《生产系统建模与仿真》是面向工程实际的应用型课程,是工业工程系的主导课程之一。学生通过本课程的学习能够初步运用仿真技术来发现生产系统中的关键问题,并通过改进措施的实现,提高生产能力和生产效率。 本课程具有较强的理论性,同时具有较强的实践性和应用性,能够有效增强学生的系统仿真理论基础,提高学生对系统仿真、分析工作的适应性,培养其开发创新能力。 本课程的教学目标是培养学生的设计能力、创新能力和工程意识。课程以制造型生产企业为核心,通过理论教学和实践环节相结合,阐述了离散事件系统建模与仿真技术在生产企业分析中的基本原理和方法。其容涉及计算机仿真技术在生产系统分析中的作用和原理、仿真软件的介绍,重点介绍排队系统、库存系统、加工系统以及输入、输出数据分析。本课程的目的是要求学生通过学习、课堂教育和上机训练,能了解如何运用计算机仿真技术模拟生产系统的布置和调度管理;并熟悉和掌握计算机仿真软件的基本操作和能够实现的功能;使学生了解计算机仿真的基本步骤。 二、课程要求 (一)教学方法 1、启发式课堂讨论 针对关键知识点、典型题和难题,通过教师提问,鼓励学生回答问题或请到讲台前做题,并请其他学生评判或提出不同的答案或不同的解决方法。目的是加强学生自主学习的能力和判断能力,培养主动思考的习惯,启发学生的探索精神。 2、重视在教学中加强知识演进的逻辑规律的讲解 提高学生的逻辑思维能力,培养学生分析问题、解决问题的能力。 3、加强计算机辅助设计、分析 将Flexsim仿真软件引入教学中。应用计算机辅助设计、分析,能方便的改变系统

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”,参数设置如下图:

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; %%根据抽样时间间隔进行抽样

系统建模与仿真习题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:

simulink-matlab仿真教程

simulink matlab 仿真环境教程 Simulink 是面向框图的仿真软件。 演示一个Simulink 的简单程序 【例1.1】创建一个正弦信号的仿真模型。 步骤如下: (1) 在MATLAB 的命令窗口运行simulink 命令,或单击工具栏中的图标,就可以打开Simulink 模块库浏览器 (Simulink Library Browser) 窗口,如图1.1所示。 (2) 单击工具栏上的图标或选择菜单“File ”——“New ”——“Model ”,新建一个名为“untitled ”的空白 模型窗口。 (3) 在上图的右侧子模块窗口中,单击“Source ”子模块库前的“+”(或双击Source),或者直接在左侧模块和工具箱栏单击Simulink 下的Source 子模块库,便可看到各种输入源模块。 (4) 用鼠标单击所需要的输入信号源模块“Sine Wave ”(正弦信号),将其拖放到的空白模型窗口“untitled ”,则“Sine Wave ”模块就被添加到untitled 窗口;也可以用鼠标选中“Sine Wave ”模块,单击鼠标右键,在快捷菜单中选择“add to 'untitled'”命令,就可以将“Sine Wave ”模块添加到untitled 窗口,如图1.2 所示。 图7.1 Simulink 界面

(5) 用同样的方法打开接收模块库“Sinks”,选择其中的“Scope ”模块(示波器)拖放到“untitled”窗口中。 (6) 在“untitled”窗口中,用鼠标指向“Sine Wave”右侧的输出端,当光标变为十字符时,按住鼠标拖向“Scope”模块的输入端,松开鼠标按键,就完成了两个模块间的信号线连接,一个简单模型已经建成。如图1.3所示。 (7) 开始仿真,单击“untitled”模型窗口中“开始仿真”图标,或者选择菜单“Simulink”——“Start”,则仿真开始。双击“Scope”模块出现示波器显示屏,可以看到黄色的正弦波形。如图1.4所示。 (8) 保存模型,单击工具栏的图标,将该模型保存为“Ex0701.mdl”文件。 1.2 Simulink的文件操作和模型窗口 1.2.1 Simulink的文件操作 1. 新建文件 新建仿真模型文件有几种操作: ?在MATLAB的命令窗口选择菜单“File”“New”“Model”。 图7.2 Simulink界面 图7.3 Simulink模型窗口 图7.4 示波器窗口

曾华艳组离散事件系统仿真大作业

新疆财经大学实验报告 课程名称:物流管理综合实验 实验项目名称:系统建模与仿真 学号: 2013104059 姓名:曾华艳 班级:物流管理11-1 指导教师:林秋平 2014年 6月 2日

新疆财经大学实验报告

《铁路局联通营业厅排队仿真分析实验报告》 一、实验目的 (一)通过对铁路局联通营业厅运作的观察,建立计算机仿真全过程,对营业厅运作进行数据采集、建模和仿真分析,为联通营业厅提出改进和优化方案的建议。 (二)通过这次实验活动,全面了解计算机仿真技术在物流领域、生产制造领域等离散事件系统中的应用,理解仿真技术如何辅助管理人员进行决策。 (三)通过分组合作的形式,提供一种系统仿真工作中常见的团队协作方式的实践体验,培养协调工作、共同完成任务的能力。 二、系统描述 人们进入联通营业厅,首先要通过取票系统拿到自己的号,先在等待区等待叫号系统报自己的号。一共有2个服务台,2个服务台同时工作,哪个服务台叫到几号,拿这个号码的人就去哪个服务台,叫号系统按顺序叫号,2个服务台叫号不会发生重复现象。我们组决定针对铁路局联通营业厅叫号排队办理业务的过程进行研究,因此我们采集了仿真模型相关数据。记录了每位顾客到达时间、等待时间和离开时间。将收集的数据整理,录入excel中,并计算出了顾客的到达时间间隔和被服务时间,再利用flexsim建立仿真模型进行仿真分析与优化。 三、小组分工 (一)本组成员 1.组长:曾华艳 2.组员:晁芙蓉、陈磊、阿尔孜姑丽、宗泽宁、张振恒 (二)小组分工 1.调查收集数据和模型优化:全体成员 2.数据录入:晁芙蓉、张振恒、阿尔孜姑丽 3.数据处理:宗泽宁、阿尔孜姑丽 4.仿真模型建立与分析:陈磊、曾华艳 5.实验报告:曾华艳、晁芙蓉、宗泽宁 6.PPT 制作:张振恒、陈磊

大作业题目

控制系统仿真大作业 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要素,后合成解

第一章 系统建模与仿真概述

第一章系统建模与仿真概述 系统:系统是由两个以上相互区别或相互作用的单元有机的结合在起来,完成某一功能的综合体。 系统的特征:1.系统的整体性 2.系统的层次性 3.系统的相关系 4.系统的目的性 5.系统对环境的适应性系统: 模型:模型是对系统的特征要素,有关信息和变化规律的一种抽象表述、它反映 了系统某些本质属性,描述了系统各要素间的相互关系,系统与环境之间的相互 作用。 模型的意义:1.客观实体系统很难做试验,或者根本不能做实验。 2.对象问题虽然可以做试验,但是利用模型更便于理解。 3.模型易于操作,利用模型的参数变化来了解现实问题的本质和规 律更加经济方便。 系统模型的种类:抽象模型和形象模型 抽象模型:数学模型图形模型计算机模型概念模型 形象模型:模拟模型实体模型 建立模型的步骤: 1.弄清问题,掌握实际情况 2.搜集资料 3.确定因素之间的关系 4.构造建模 5.求解模型 6.检验模型的正确性 系统建模预防针的一般方法和步骤(P17) 仿真的发展趋势:建模方法面对对象仿真分布交互仿真人工智能与 计算机仿真虚拟现实仿真 Internet网上仿真 第二章商贸物流系统建模与仿真 商贸流通在社会经济中的地位与作用:1,商贸流通是连接生产和消费的纽带; 2,商贸流通对生产具有反作用; 3,商贸流通是国民经济现代化的支柱。 商贸活动的内容: 1,商流,对象物所有权转移的活动称为商流。 2,物流,是指事物从供给方向需求方的转移。

3,资金流,主要是指资金流的转移过程,包括付款,转账等过程,是 整个商贸活动的目的。 4,信息流,指商品信息的提供,商品促销信息,技术支持,售后服务 等内容,也包括诸如询单价,报单价,付款通知单,转账通知单等商业贸易单证以及交易 方的支付能力和支付信誉。 预测:所谓预测就是人们对某一不确定的或未知事件的表述。 预测的作用:从变化的事物中找出使事物发生变化的变化的固有规律,寻找和研究各种变化现象的背景及其演变的逻辑关系,从而去揭示事物未来的面貌。 判断预测方法:一,部门负责人评判预测法;二,销售人员估计法;三,德尔菲法;四, 历时类比法。 德尔菲法:依靠技术专家小组背靠背景来判断,来代替面对面的会议,是不同专家将分歧的幅度和理由都能够表达出来,经过客观分析以求达到客观规律的一致意见。 时间序列预测技术:一,移动平均预测法(计算题p30例2); 二,指数平均预测法。 DRP:是分销需求计划的简称,它是MRP原理和技术在流通领域中的应用。该技术主要解决分销物资的应用和调度问题,其基本目标是合理进行分销物资和资源配置,以达到既有效 地满足市场需求优势的配置费用最省的目的。 *DRP的基本概念 1.库存:指仓库或物流中心实际存在的物资数量。 2.安全库存:为便于生产经营活动正常进行,防止因需求货供应的波动 引起缺货或停工待料,经常在仓库各项目保持一定数量的计划库存量, 成为安全库存。 3.期初和期末库存:指在论述的时间段开始和结束时本单位的实际库存。 4.进货提前期:指从发出订货到所定货物运回并入库所需要的时间长度。 5.送货提前期:指从接收订单到货物送到用户手中并接收入库的时间长度。 6.在途物资:指供应商已经接受订单备货,但尚未来到本单位入库的物资。 7.订货批量:指一次订货所订的物资数量。 8.时间周期:就是根据实际需要划分的时间段信息,如一日,周,月划分。 9.计划期:是指DRP尽心运算的整个时间段,可能是一个月,一个季度 或一年,他可划分为几个计划周期。 10.物流中心:从事物流活动的具有完善的信息网络的场所或组织。 BOD简介:B OD是MRP中物料清单BOM的概念和结构在分销领域的运用,它同BOM在产品结构树中连接各零件和成品一样,在供应方和各个需求方之间架起了一座沟通的桥梁。 DRP在分销网络中的运作原理(p43DRP原理图)

电力电子电路建模与分析大作业要点

西安理工大学 研究生课程论文/研究报告 课程名称:电力电子系统建模与分析 任课教师: 完成日期:2016 年7 月 5 日 专业:电力电子与电力传动 学号: 姓名: 同组成员: 成绩:

题目要求 某用户需要一直流电源,要求:直流输出24V/200W,输出电压波动及纹波均<1%。用户有220V交流电网(±10%波动变化)可供使用: (1) 设计电源主电路及其参数; (2) 建立电路数学模型,获得开关变换器传函模型; (3) 设计控制器参数,给出控制补偿器前和补偿后开环传递函数波特图,分 析系统的动态和稳态性能; (4) 根据设计的控制补偿器参数进行电路仿真,实现电源要求; (5) 讨论建模中忽略或近似因素对数学模型的影响,得出适应性结论(量化 性结论:如具体开关频率、具体允许扰动幅值及频率等)。 主要工作 本次设计主要负责电源主电路及其参数的的设计,以及建立电路数学模型并获得开关变换器传函模型这两部分内容,具体如下: (1) 本次设计电源主电路及其参数,采用从后向前的逆向设计思想。首先根据系统输出要求,设计了后级DC/DC型Buck电路的参数。接着设计了前级不控整流电路以及工频变压器的参数。考虑到主电路启动运行时的安全性,在主电路中加入了软启动电路; (2) 本次DC/DC变换器的建模并没有采用传统的状态空间平均方法,而是采用更为简单、直观的平均开关建模方法,建立了Buck变换器小信号交流模型。最后,推到出了开关变换器的传递函数模型,并给出了Buck电路闭环控制框图。

1 设计主电路及其参数 1.1主电路设计 根据题目要求,系统为单相交流220V/50Hz 输入,直流24V/200W 输出。对于小功率单相交流输入的场合,由于二极管不控整流电路简单,可靠性高,产生的高次谐波较少,广泛应用于不间断电源(UPS)、开关电源等场合。所以初步确定本系统主电路拓扑为:前级AC-DC 电路为电源经变压器降压后的二极管不控整流,后级DC-DC 电路为Buck 斩波电路,其中Buck 电路工作在电感电流连续模式(CCM ),前后级之间通过直流母线和直流电容连接在一起。系统主电路结构如图1-1所示。 AC 220V/50Hz L C 1 C 2R D S 图1-1 系统主电路结构图 1.2主电路参数设计 本次设计电源主电路参数,采用从后向前的逆向设计思想。先对后级DC/DC 型Buck 电路的参数进行设计,接着对前级不控整流电路以及工频变压器的参数进行设计。下面分别对后级的Buck 电路和前级经变压器降压后的不控整流电路各参数进行分析设计。 1.2.1 输出电阻计算 根据系统电路参数:220,50;24;200i o U V Hz U V P W ===,可计算: 输出电流: /200/248.33O O I P U W V A ==≈ (1-1) 负载等值电阻: /24/8.33 2.88O O R U I V A ==≈Ω (1-2)

微纳系统仿真大作业

1. 用有限差分法和有线元方法把以下问题变成数值方程,并说明两种方法的异同: 2(,)0x y ??= 边界条件: (,0)(,1)0;(0,)(1,)1; x x y y ????==== 解: (1) 有限差分法 2(,)0x y ??= 即为20xx yy u u u ?=+=,其中(),u x y ?= 将定义域等分,步长均为h ,则 ()()()()()()222,,,,4,u x h y u x h y u x y h u x y h u x y u O h h ++-+++--?=+ 略去高阶无穷小,由20xx yy u u u ?=+=得 ()()()()(),,,,4,0u x h y u x h y u x y h u x y h u x y ++-+++--= 定义域离散,离散点为(),i j x y ,则上式可化为 ()()()()(),,,,4,0i j i j i j i j i j u x h y u x h y u x y h u x y h u x y ++-+++--= 定义域为01,01x y ≤≤≤≤,因为本题只是为了说明原理,故将其简单等分为33?单元,节点编号从()()0,02,2到。 则 1,1,,1,1,40i j i j i j i j i j u u u u u +-+-+++-= 其位移矢量为

()()()()()()()()()0,00,10,21,01,11,22,02,12,24110141101410141101411011410141011410140u u u u u u u u u ????-??????????-????????????-??????-????????????=-??????-????????????-????????-??????????-?????????? 由边界条件 (,0)(,1)0;(0,)(1,)1; x x y y ????==== 可知()()()()0,00,22,02,2,,,u u u u 有冲突,可以区位均值即()()()()0,00,22,02,20.5u u u u ====,而()()1,11,20u u ==,()()0,12,10u u ==,带入化简即可求得()1,1u (2) 有限元法 使用有限元法的计算流程为: 求解区域离散化; 构造插值函数形成分段光滑的坐标函数系; 用 Ritz 方法求解微分方程 对2(,)0x y ??=构造函数 ()221,2x y dxdy x x ???????????∏=+???? ? ????????????? ?? 首先将整个区域离散为三角形的子区域如下图, 三角形微小子区域中的值由三 角形节点值的插值结果表示,即 (),i i j j k k x y N N N ????=++

2015数学建模选修大作业

中华女子学院 成绩2014 — 2015学年第二学期期末考试 (论文类) 论文题目数学建模算法之蒙特卡罗算法 课程代码1077080001 课程名称数学建模 学号130801019

姓名陈可心 院系计算机系 专业计算机科学与技术 考试时间2015年5月27日 一、数学建模十大算法 1、蒙特卡罗算法 该算法又称随机性模拟算法,是通过计算机仿真来解决问题的算法,同时可以通过模拟可以来检验自己模型的正确性,是比赛时必用的方法。接下来本文将着重介绍这一算法。 2、数据拟合、参数估计、插值等数据处理算法 比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用Matlab作为工具。 3、线性规划、整数规划、多元规划、二次规划等规划类问题 建模竞赛大多数问题属于最优化问题,很多时候这些问题可以用数学规划算法来描述,通常使用Lindo、Lingo软件实现。这个也是我们数学建模选修课时主要介绍的问题,所以对这方面比较熟悉,也了解了Lindo、Lingo软件的基本用法。 4、图论算法 这类算法可以分为很多种,包括最短路、网络流、二分图等算法,涉及到图论的问题可以用这些方法解决,上学期数据结构课程以及离散数学课程中都有介绍。它提供了对很多问题都很有效的一种简单而系统的建模方式。

5、动态规划、回溯搜索、分治算法、分支定界等计算机算法 这些算法是算法设计中比较常用的方法,很多场合可以用到竞赛中 6、最优化理论的三大非经典算法:模拟退火法、神经网络、遗传算法 这些问题是用来解决一些较困难的最优化问题的算法,对于有些问题非常有帮助,但是算法的实现比较困难,需慎重使用。 7、网格算法和穷举法 网格算法和穷举法都是暴力搜索最优点的算法,在很多竞赛题中有应用,当重点讨论模型本身而轻视算法的时候,可以使用这种暴力方案,最好使用一些高级语言作为编程工具。 8、一些连续离散化方法 很多问题都是实际来的,数据可以是连续的,而计算机只认的是离散的数据,因此将其离散化后进行差分代替微分、求和代替积分等思想是非常重要的。 9、数值分析算法 如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比如方程组求解、矩阵运算、函数积分等算法就需要额外编写库函数进行调用。10、图象处理算法 赛题中有一类问题与图形有关,即使与图形无关,论文中也应该要不乏图片的,这些图形如何展示以及如何处理就是需要解决的问题,通常使用Matlab进行处理。 二、蒙特卡罗方法 2.1算法简介 蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,1946年,美国拉斯阿莫斯国家实验室的三位科学家John von Neumann,Stan Ulam 和 Nick

MATLAB电路仿真实例

题14.14 图(a)所示电路,已知 V )2cos(15S t u =二端口网络阻抗参数矩阵 Ω?? ????=46j 6j 10Z 求ab 端戴维南等效电路并计算电压o u 。 u -+o u 图题14.14 (一)手动求解: 将网络N 用T 型电路等效,如图(b)所示 S U +-o U 等效阻抗 Ω=-+-?+ -=4.6j615j6j6)15(6j 6j 4i Z 开路电压 V 2j302 15j6j6105j6OC =?∠?+-+=U V 1482 18.3j46.42j3j4j4Z j4OC o ?∠=+?=?+=U U i

所以 )1482cos(18.3o ?+=t u V (二)Matlab 仿真: ⒈分析:本次仿真需输入各阻抗Zl 、Z1、Z2、Z3、Z4以及激励源Us 的参数值,仿真结果需输出开路电压Uoc 、等效阻抗Zi 以及电感两端电压U0的幅值和相位信息,并绘制Uoc ,U0的值随时间变化的波形曲线。其中各元件与原图的对应关系如下图所示: ⒉编辑M 文件的源程序如下: clear %清空自定义变量 z1=4-6j;z2=6j;z3=10-6j;z4=5;us=15*exp(j*0);zl=4j;%输入各元件参数 zi=z1+(z2*(z3+z4)/(z2+z3+z4));%等效阻抗zi 的计算表达式uoc=us*z2/(z2+z3+z4);%开路电压uoc 的计算表达式u0=zl/(zi+zl)*uoc;%电感两端电压uo 的计算表达式disp('The magnitude of zi is'); %在屏幕上显示“The magnitude of zi is ”disp(abs(zi)) %显示等效阻抗zi 的模disp('The phase of zi is'); %在屏幕上显示“The phase of zi is ”disp(angle(zi)*180/pi)%显示等效阻抗zi 的辐角 disp('The magnitude of uoc is'); %在屏幕上显示“The magnitude of uoc is ” disp(abs(uoc))%显示开路电压uoc 的模

相关文档
最新文档