%%?函數定義

%此處函數為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