1、积分的数值求解函数 积分 : (求面积 )trapz(x , y)给定数据点 x和 y， 计算y=f(x)下的梯形面积积分。精度更高的求积分函数 quad,QUADL;近似求原函数并做图mm2302二维积分 :Dblquad三维积分 : Triplequad 微分 : 具有其内在的敏感性 ,故对于被积分的数据使通过实验获得时 ,避免使用 .应该相对数据做 ( ),然后对 ( )进行微分 . Polyfit Polyderdiff(x)数组元素间的差值 (diff(y)./diff(x)近似求微分 )，结果分析 mm2305,mm2306 结论二维 微分 :梯度 gradient 主要用作图形数
2、据可视化(mm2307)quiver曲率 del2(mm2308)微分的数值求解函数微分方程 (ODE)（ 1）1.IVP格式2.IVP微分方程解法程序 t , y=ode45(fname , tspan , yo)用 4阶 /5阶龙格 -库塔算法解微分方程组 other 6 ODE solvers: ODE23, ODE113, ODE15S, ODE23S, ODE23T, ODE23TB 设置选项 :OPTIONS=ODESET (name1,value1,), ODEGET ；微分方程 (ODE)（ 2） output functions: ODEPLOT,ODEPHAS2,ODEPH
3、AS3, ODEPRINT； evaluating solution: DEVAL； ODE examples:RIGIDODE,BALLODE, ORBITODE边界值问题 (BVP)和偏微分方程问题基本用法 （连续时间系统仿真 1）1、微分方程变形： ydot=ldefile(t,y)m文件的编写 function ydot=vdpol(t,y,mu) %VDPOL van der Pol equation. % Ydot=VDPOL(t,Y) % Ydot(1) = Y(2) % Ydot(2) = mu*(1-Y(1)2)*Y(2)-Y(1) % mu = 2 %ydot = y(2,
4、:); mu*(1-y(1,:)2)*y(2,:)-y(1,:); ydot = y(2); mu*(1-y(1)2)*y(2)-y(1);2、 求解命令的编写 2401.doc基本用法 （连续时间系统仿真 2）例 :微分方程的数值解法是在科学与工程计算中经常遇到的问题。著名的 Lorenz 模型的状态方程表示为：若令 ，且假设初值为 ，e 为一个小常数，洛伦兹 .doc离散系统仿真 Filter ConSIMULINK动态仿真SIMULINK动态仿真是 MATLAB软件包之一，用于可视化的动态系统仿真，适用于连续、离散系统，也适用于线性和非线性系统。1） 启动 SIMULINK， 进入 SI
5、MULINK窗口 ；2）在 SIMULINK窗口下，借助模块库，创建系统 框图模型并调整模块参数。3）两种方式 （ SIMULINK模型窗口、 matlab command 窗口 sim） 下，设置仿真参数后， 启动仿真。4）输出仿真结果：利用 scope模块，利用变量返 回值、利用 to workspace模块 。SIMULINK动态仿真 (2) 模块平衡点 （ trim)和线性化 (linmod、linmod2、 dlinmod等） 子系统创建及封装 S-函数 非线性系统优化设计作业 2:系统时间响应和动态仿真应用 4阶龙格 -库塔法和使用 simulink软件包仿真程序，求如下系统的响应 。
