四阶五级matlab,微分方程数值解法matlab(四阶龙格

微分方程的数值解法四阶龙格—库塔法 (The Fourth-Order Runge－Kutta Method) 常微分方程(Ordinary differential equations, ODE) 初值问题---给出初始值边值问题---给出边界条件一.解ODE的基本机理: 3. 根据式(2.2)编写计算导数的M函数文件-ODE文件例题1：著名的Van der Pol方程四阶Runge-Kutta公式三. Runge-Kutta 法解Van der Pol 方程的Matlab 程序结构主程序：RK_vanderpol.m 子程序：RK_sub.m(函数文件) 解法2：采用Runge_Kutta法编程计算主程序：RK_vanderpol.m t0=0; tN=20; y0=[0.25; 0]; h=0.001; t = t0 : h : tN; N = length (t); j = 1; for i = 1 : N t1 = t0 + h; K1 = RK_sub(t0, y0); K2 = RK_sub(t0 + h/2, y0 + h*K1/2); K3 = RK_sub(t0 + h/2, y0 + h*K2/2); K4 = RK_sub(t0 + h, y0 + h*K3); y1 = y0 + (h/6)*(K1 + 2*K2 + 2*K3 + K4); yy1(j)=y1(1); yy2(j)=y1(2); t0=t1; y0=y1; j=j+1; end subplot (121), plot (t, yy1, t, yy2); grid subplot (122), plot (yy2, yy1); grid 四. Matlab对应命令：ode23,ode45 五. 动力学系统的求解 1. 动力学方程 3. Matlab 程序(主程序：ZCX) t0;Y0;h;N;P0,w; %输入初始值、步长、迭代次数、初始激励力； for i = 1 : N t1 = t0 + h P=[P0*sin(w*t0);0.0;0.0] %输入t0时刻的外部激励力 K1 = ZCX_sub (t0, Y0, P ) P= %输入 (t0+h/2) 时刻的外部激励力 K2 = ZCX_sub (t0 + h/2, Y0 + hK1/2, P ) K3 = ZCX_sub (t0 + h/2, Y0 + hK2/2, P ) P= %输入 (t0+h) 时刻的外部激励力 K4 = ZCX_sub (t0 + h, y0 + hK3, P) Y1 = y0 + (h/6) (K1 + 2K2 + 2K3 + K4) t1, Y1 (输出 t1, y1) next i 输出数据或图形 Matlab 程序(子程序：ZCX_sub.m) 例题2：三自由度质量弹簧系统 * * 与初值常微分方程解算有关的指令 ode23 ode45 ode113 ode23t ode15s ode23s ode23tb 2. 把高阶方程转换成一阶微分方程组 1. 列出微分方程初始条件令 (2.1) (2.2) (2.3) 例：著名的Van der Pol方程令降为一阶初始条件把t,Y作为输入宗量，把作为输出宗量 %M function file name: dYdt.m function Yd = f (t, Y) Yd = f (t,Y) 的展开式例Van der Pol方程 %M function file name: dYdt.m function Yd = f (t, Y) Yd=zeros(size(Y)); 4. 使编写好的ODE函数文件和初值供微分方程解算指令(solver)调用 Solver解算指令的使用格式 [t, Y]=solver (‘ODE函数文件名’, t0, tN, Y0, tol); ode45 输出宗量形式说明： t0：初始时刻；tN：终点时刻Y0：初值； tol：计算精度 % 主程序

四阶五级matlab,微分方程数值解法matlab(四阶龙格—库塔法).ppt相关推荐

matlab非线性数值解法,Matlab非线性方程数值解法（2）
实验目的用Matlab实现非线性方程的Aitken加速法.牛顿法.简化牛顿法和牛顿下山法. 实验要求给出Aitken加速法.牛顿法.简化牛顿法和牛顿下山法算法. 用Matlab实现Aitken加速 ...
matlab非线性数值解法,matlab 数值分析非线性方程与方程组的数值解法
feval function x=fun(a,b) x=a+b; 调用 1.feval(@fun,a,b); 2.feval('fun',a,b); 作为参数时 function x=main(fun ...
[常微分方程的数值解法系列五] 龙格-库塔(RK4)法
龙格-库塔法简介基本思想具体方法一阶二阶求解参数特殊二阶三阶高阶步长选择例子在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其 ...
二阶龙格库塔公式推导_[常微分方程的数值解法系列五] 龙格-库塔(RK4)法
在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主要是应用于一些简单和特殊的微分方程求解中,对于 ...
有确定项微分方程的matlab程序,微分方程的数值解法matlab四阶龙格—库塔法课件...
<微分方程的数值解法matlab四阶龙格-库塔法课件>由会员分享,可在线阅读,更多相关<微分方程的数值解法matlab四阶龙格-库塔法课件(36页珍藏版)>请在人人文库网上搜索 ...
常微分方程数值解matlab欧拉,matlab 常微分方程数值解法源程序代码
matlab 常微分方程数值解法源程序代码所属分类:其他开发工具:matlab 文件大小:16KB 下载次数:41 上传日期:2019-02-13 11:03:29 上传者:XWLYF 说明 ...
一节双曲型方程基于MATLAB的求解,双曲方程基于matlab的数值解法
<双曲方程基于matlab的数值解法>由会员分享,可在线阅读,更多相关<双曲方程基于matlab的数值解法(9页珍藏版)>请在人人文库网上搜索. 1.双曲型方程基于MATLAB ...
四阶龙格库塔法的基本思想_数值常微分方程-欧拉法与龙格-库塔法
大三时候在跳蚤市场闲逛,从一位数学院的学长那里买了一些闲书,最近翻出来刚好有李荣华.刘播老师的<微分方程数值解法>和王仁宏老师的<数值逼近>,结合周善贵老师的<计算物理& ...
微分方程数值解法(PID仿真用一阶被控对象库PLC算法实现)
微分方程除极特殊情况外,大部分不可能求出它的精确解,只能用各种近似方法得到满足一定精度的近似解,之前大家学过的教程知道微分方程有级数解法和Picard逐步逼近法,这些方法可以给出解的近似表达式,称为近 ...
微分方程数值解法（实际应用）
微分方程数值解法(实际应用) 假设有这样一种情况,关于人行走的,从A到B地,行走速度越来越慢,速度是关于时间1/(t+A)的函数,现在求当s和时间t的关系: 从微分的角度看,v=s'=1/(t+A) ...

四阶五级matlab,微分方程数值解法matlab(四阶龙格—库塔法).ppt

四阶五级matlab,微分方程数值解法matlab(四阶龙格—库塔法).ppt相关推荐

最新文章

热门文章