Runge-Kutta+Method 源代码

% Duffing 吸引子三维相空间图，这里用四阶 Runge-Kutta 法得到微方程的离散序列 % 时间步长 h = 0.01, 而不是1 % 方程表达式 % dx/dt = y % dy/dt = -delta*y - a*x*(1+x^2)…

% Lorenz 吸引子三维相空间图，这里用四阶 Runge-Kutta 法得到微方程的离散序列 % 时间步长 h = 0.01, 而不是1 % 方程表达式 % dx/dt = sigma*(y-x) % dy/dt = r*x - y - x*z …

function x=body1(x,u,v) % 4th order Runge-Kutta solution of differential equation in 'bodydot' persistent delta if nargin==1, …

clear all; % Calculate responses of a 3-story shear building subject to White Noise Excitation on Floor n=3; %Input data …

function [tout,yout] = rk_4(odefile,t0,h,th,y0) %RK_4 实现了 4 阶 Runge-Kutta 定步长微分方程求解。 % % [tout,yout] = rk_4(odefile,t0,h,th,y0…

function [tout,yout] = rk_4(odefile,t0,h,th,y0) %RK_4 实现了 4 阶 Runge-Kutta 定步长微分方程求解。 % % [tout,yout] = rk_4(odefile,t0,h,th,y0…

function [tout,yout] = rk_4(odefile,t0,h,th,y0) %RK_4 实现了 4 阶 Runge-Kutta 定步长微分方程求解。 % % [tout,yout] = rk_4(odefile,t0,h,th,y0…

/****************************************************************************** * …

%% orbit - Program to compute the orbit of a comet. clear all; help orbit; % Clear memory and print header %% * Set initial …

%RK4 4th order Runge-Kutta integration % % Syntax: % [x,Y] = rk4(f,dt,x,[P1,P2,P3,Y]) % % In: % f - Name of function in form…