龙格库塔原理详解及解微分方程组的实现

资源包含龙格库塔法的实现原理，利用从原理上matlab实现四阶龙格库塔法解微分方程组，微分方程组为CPG震荡图

%% 函数定义

%此处函数为CPG震荡模型
function [yy]=equation（tvalue）
    alpha=5;%x轴方向收敛速度
    beta=10;%y轴方向收敛速度
    a=150; %椭圆长轴
    b=100; %椭圆短轴
    omiga=pi;% 振荡器频率，越小收敛的越快
    
    x=value（1）;
    y=value（2）;
    
    yy（11）=alpha*（1-x*x/a/a-y*y/b/b）*x+omiga*a/b*y;
    yy（12）=beta*（1-x*x/a/a-y*y/b/b）*y-omiga*b/a*x;
end



%%调用举例
% [T1F1]=runge_kutta（@equation[0;0.2]0.01010）;
% plot（F1（:1）F1（:2）‘.‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-12-31 20:01  龙格库塔\
     文件         475  2018-12-31 19:50  龙格库塔\equation.m
     文件         528  2018-12-31 19:48  龙格库塔\runge_kutta.m
     文件      406528  2018-12-31 19:59  龙格库塔\第六章-常微分方程数值解（龙格-库塔法）.ppt