RK45(dydt, trange, yinit, tol): 类似MATLAB的ode45，RK45适用于八度音程，因为...
5星
- 浏览量: 0
- 大小:None
- 文件类型：None
简介：
RK45函数模仿了MATLAB中的ode45算法，在八度系统中求解微分方程。它通过Runge-Kutta法提供变步长积分，确保结果达到指定的容差(tol)，从初始条件(yinit)出发，沿给定时间范围(trange)演进(dydt定义系统动态)。
函数像 ode45 一样求解常微分方程。
dydt = 一阶微分方程组
trange = 自变量的范围
yinit = 初始条件
tol =误差绝对公差，默认值=1e-6
例子：
[t,y] = RK45(@(t,y)[y(2); (1-y(1)^2)*y(2)-y(1)],[0 20],[2 ; 0]);
subplot (1,2,2);
plot(t,y,o)
title(van der Pol 方程 (\mu = 1) with RK45);
xlabel(时间 t);
ylabel(解决方案 y);
legend(y_1，y_2);
[t2，y2] = ode45(@(t,y)[y(2); (1-y(1)^2)*y(2)-y(1)], [0 20], [2 ; 0]);
subplot (1,2,1);
plot(t2,y2,-);
全部评论 (0)
还没有任何评论哟~
