搜索
1.变步长(Variable—Step)求解器可以选择的变步长求解器有：ode45，ode23，ode113，odel5s，ode23s和discret．缺省情况下，具有状态的系统用的是ode45；没有状态的系统用的是discrete．1)ode45基于显式Runge—Kutta(4，5)公式，Dormand—Prince对．它是—个单步求解器(solver)。也就是说它在计算y(tn)时，仅仅利用前一步的计算结果y(tn-1)．对于大多数问题．在第一次仿真时、可用ode45试一下．2)ode23是基于显式Runge—Kutta(2，3)．Bogackt和Shampine对．对于宽误差容限和存在轻微刚性的系统、它比ode45更有效一些．ode23也是单步求解器．3)odell3是变阶Adams-Bashforth—MoultonPECE求解器．在误差容限比较严时，它比ode45更有效．odell3是一个多步求解器，即为了计算当前的结果y(tn），不仅要知道前一步结果y(tn-1)，还要知道前几步的结果y(tn-2)，y(tn-3)，…；4)odel5s是基于数值微分公式(NDFs)的变阶求解器．它与后向微分公式BDFs(也叫Gear方法)有联系．但比它更有效．ode15s是一个多步求解器，如果认为一个问题是刚性的，或者在用ode45s时仿真失败或不够有效时，可以试试odel5s。odel5s是基于一到五阶的NDF公式的求解器．尽管公式的阶数越高结果越精确，但稳定性会差一些．如果模型是刚性的，并且要求有比较好的稳定性，应将最大的阶数减小到2．选择odel5s求解器时，对话框中会显示这一参数．可以用ode23求解器代替。del5s，ode23是定步长、低阶求解器．5)ode23s是基于一个2阶改进的Rosenbrock公式．因为它是一个单步求解器，所以对于宽误差容限，它比odel5s更有效．对于一些用odel5s不是很有效的刚性问题，可以用它解决．6)ode23t是使用“自由”内插式梯形规则来实现的．如果问题是适度刚性，而且需要没有数字阻尼的结果，可采用该求解器．7)ode23tb是使用TR—BDF2来实现的，即基于隐式Runge—Kutta公式，其第一级是梯形规则步长和第二级是二阶反向微分公式．两级计算使用相同的迭代矩阵．与ode23s相似，对于宽误差容限，它比odtl5s更有效．8)discrete(变步长)是simulink在检测到模型中没有连续状态时所选择的一种求解器．2.定步长(Flxed—Step)求解器可以选择的定步长求解器有：ode5，ode4，ode3，ode2，ode1和discrete．1)ode5是ode45的一个定步长版本，基于Dormand—Prince公式．2）ode4是RK4，基于四阶Runge—Kutta公式．3)ode3是ode23的定步长版本，基于Bogacki-Sbampine公式．4)ode2是Heun方法，也叫作改进Euler公式．5)odel是Euler方法．6)discrete(定步长)是不执行积分的定步长求解器．它适用于没有状态的模型，以及对过零点检测和误差控制不重要的模型．3.诊断页（Diagnostics)可以通过选择SimulationParameters对话框的Diagnostics标签来指明在仿真期间遇到一些事件或者条件时希望执行的动作．对于每一事件类型，可以选择是否需要提示消息，是警告消息还是错误消息．警告消息不会终止仿真，而错误消息则会中止仿真的运行．（1）一致性检查一致性检查是一个调试工具．用它可以验证Simulink的0DE求解器所做的某些假设．它的主要用途是确保s函数遵循Simulink内建模块所遵循的规则．因为一致性检查会导致性能的大幅度下阵(高达40％)，所以一般应将它设为关的状态．使用一致性检查可以验证s函数，并有助于确定导致意外仿真结果的原因．为了执行高效的积分运算，Simulink保存一些时间步的结果，并提供给下一时间步使用．例如，某一时间步结束的导数通常可以放下一时间步开始时再使用．求解器利用这一点可以防止多余的导数运算．一致性检查的另一个目的是保证当模块被以一个给定的t(时间)值调用时．它产生一常量输出．这对于刚性求解器(ode23s和odel5s)非常重要，因为当计算Jacobi行列式时．模块的输出函数可能会被以相同的t值调用多次．如果选择了一致性检查，Simulink置新计算某些值，并将它们与保存在内存中的值进行比较，如果这些值有不相同的，将会产生一致性错误．Simulink比较下列量的计算值：1）输出；2）过零点3）导数；4）状态．（2）关闭过零点检测可以关闭一个仿真的过零点检测．对于一个有过零点的模型，关闭过零点检测会加快仿真的速度，但是可能影响仿真结果的精度．这一选项关闭那些本来就有过零点检测的模块的过零点检测．它不能关闭Hircrossing模块的过零点检测．（3）关闭优化I/O存储选择该选项，将导致Simulink为每个模块约I／()值分配单独的缓存，而不是重新利用援存．这样可以充分增加大模型仿真所需内存的数量．只有需要调试模型时才选择该选项．在下列情况下，应当关闭缓存再利用；1)调试一个C-MEXS-函数；2)使用浮点scope或display模块来察看调试模型中的信号．如果缓存再利用打开，并且试图使用浮点scope或display模块来显示缓存已被再利用的信号，将会打开一个错误对话框．（4）放松逻辑类型检验选择该选项，可使要求逻辑类型输入的模块接受双精度类型输入．这样可保证与Simulink3版本之前的模型的兼容性．4.提高仿真性能和精度仿值性能相精度由多种因素决定，包括模型的设计和仿真参数的选择．求解器使用它们的缺省参数值可以使大多数模型的仿真比较精确有效，然而，对于一些模型如果调整求解器相仿真参数将会产生更好的结果．而且，如果对模型的性能比较熟悉，并且将这些信息提供给求解器，得到的仿真效果将会提高。（1）加快仿真速度仿真速度慢的原因有多种．下面列举其中的一些：1)模型中包含有MATLAB的Fcn模块．当模型包含有MATLAB的Fcn模块时，在仿真的每一时间步都会调用MATLAB的解释器，这将大大地减慢仿真的速度．因此应尽可能地使用内建的Fcn模块或者E1ementaryMath模块．2)模型中包含有M文件形式的S函数．M文件形式的S函数也将导致在每一时间步调用MATLAB的解释器．可以考虑将s函数转换为子系统或者c-MEx文件形式的s函数．3)模型中包含有Memory模块．使用Memory模块使得变阶求解器(odel5s和odell3)在每一时间步将阶数设为1阶．4)最大的步长太小．如果改变了最大步长．可以试试重新使用缺省值(auto)运行仿真。5)对精度要求太高．缺省的相对容差（0.1％)通常已经足够了．对于状态值趋于0的模型．如果绝对容差设得太小．仿真时状态值在零点附近会花去太多的时间步．6)时间尺度可能太长．减小时间间隔；7)问题可能是刚性的．而使用的是非刚性求解器．这时可用odel5s试一下；8)模型使用的采样时间相互之间不成倍数关系．相互之间不成倍数的混合采样时间会导致求解器采用足够小的步长，以保证采样时间符合所有的采样时间要求．9)模型包含有代数循环．在每一时间步都会反复计算代数循环，因此这会大大地降低仿真的性能．10）模型中将RandomNumber模块的输出传给了Intergrator模块．对于连续系统，在Sources库中使用Bond—limitedNoise模块．（2）改进仿真精度要检查仿真的精度．仿真运行一段时间以后，减小相对容差到1e-4或者减小绝对容差，并重新运行它．比较两次仿真的结果．如果它们之间没有很大的差别，可以确信结果收敛．如果经过一段时间后，仿真结果变得不稳定，可能是如下原因：1)系统可能不稳定．2)如果使用的是odel5s，可能需要将最大的阶数限制在2阶(求解器稳定的最大阶数)，或者试试用ode23s求解器．如果仿真结果看起来不是很精确，可能是：1)对于一个拥有趋于零的状态值的模型，如果绝对容差设得太大，仿真在零状态值附近花的步数太少．减小绝对容差的大小或者在Integrator对话框中为每—个状态分别调整绝对容差的设定．2）如果减小绝对容差不能有效地提高精度，减少相对容差的大小，减小步长，增加步数．5.处理复数信号在缺省状态下，simulink信号值是实数．尽管如此，模型可以产生相处理具有复数值的信号．可以采取以下任何一种方法在模型中引入复数值信号．1)由顶层瑞口从MATLAB工作空间向模型中装入复数值信号数据．2)在模型中产生一个常数模块，并设置其值为复数．3)产生实数信号分别对应于复数信号的实部和虚部．然后用实-虚复数转换(Real-ImagtoComplexConversion)模块将各部连成复数信号．可以用接受复数的模块来处理复数信号，大多数simulink模块接受复数信导作为输入。一、设置仿真参数和选择解法器设置仿真参数和选择解法器，选择Simulation菜单下的Parameters命令，就会弹出一个仿真参数对话框，它主要用三个页面来管理仿真的参数。Solver页，它允许用户设置仿真的开始和结束时间，选择解法器，说明解法器参数及选择一些输出选项。WorkspaceI/O页，作用是管理模型从MATLAB工作空间的输入和对它的输出。Diagnostics页，允许用户选择Simulink在仿真中显示的警告信息的等级。1、Solver页此页可以进行的设置有：选择仿真开始和结束的时间；选择解法器，并设定它的参数；选择输出项。仿真时间：注意这里的时间概念与真实的时间并不一样，只是计算机仿真中对时间的一种表示，比如10秒的仿真时间，如果采样步长定为0.1，则需要执行100步，若把步长减小，则采样点数增加，那么实际的执行时间就会增加。一般仿真开始时间设为0，而结束时间视不同的因素而选择。总的说来，执行一次仿真要耗费的时间依赖于很多因素，包括模型的复杂程度、解法器及其步长的选择、计算机时钟的速度等等。仿真步长模式：用户在Type后面的第一个下拉选项框中指定仿真的步长选取方式，可供选择的有Variable-step（变步长）和Fixed-step（固定步长）方式。变步长模式可以在仿真的过程中改变步长，提供误差控制和过零检测。固定步长模式在仿真过程中提供固定的步长，不提供误差控制和过零检测。用户还可以在第二个下拉选项框中选择对应模式下仿真所采用的算法。变步长模式解法器有：ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb和discrete。ode45：缺省值，四/五阶龙格－库塔法，适用于大多数连续或离散系统，但不适用于刚性（stiff）系统。它是单步解法器，也就是，在计算y(tn)时，它仅需要最近处理时刻的结果y(tn-1)。一般来说，面对一个仿真问题最好是首先试试ode45。ode23：二/三阶龙格－库塔法，它在误差限要求不高和求解的问题不太难的情况下，可能会比ode45更有效。也是一个单步解法器。ode113：是一种阶数可变的解法器，它在误差容许要求严格的情况下通常比ode45有效。ode113是一种多步解法器，也就是在计算当前时刻输出时，它需要以前多个时刻的解。ode15s：是一种基于数字微分公式的解法器（NDFs）。也是一种多步解法器。适用于刚性系统，当用户估计要解决的问题是比较困难的，或者不能使用ode45，或者即使使用效果也不好，就可以用ode15s。ode23s：它是一种单步解法器，专门应用于刚性系统，在弱误差允许下的效果好于ode15s。它能解决某些ode15s所不能有效解决的stiff问题。ode23t：是梯形规则的一种自由插值实现。这种解法器适用于求解适度stiff的问题而用户又需要一个无数字振荡的解法器的情况。ode23tb：是TR-BDF2的一种实现，TR-BDF2是具有两个阶段的隐式龙格－库塔公式。discrt