% right-hand side of arenstorf problem function yDot = arenstorf(t,y) global mu muHat % unpack y
u1 = y(1);
u2 = y(2);
u1Dot = y(3);
u2Dot = y(4); % define variables of convenience
r1 = ((u1+mu)^2+u2^2)^(3/2);
r2 = ((u1-muHat)^2+u2^2)^(3/2); u1DDot = u1 + 2*u2Dot - muHat*(u1+mu)/r1 - mu*(u1-muHat)/r2;
u2DDot = u2 - 2*u1Dot - muHat*u2/r1 - mu*u2/r2; % pack up derivatives
yDot = [u1Dot u2Dot u1DDot u2DDot]';
相关文章
- Matlab b AudioToolBox实时运算
- matlab的帮助文档切换成中文(求助贴)
- 【MATLAB】最简单最可靠的方法!!再次续更!离线下载runtime后,matlab无法识别runtime的问题(matlaba)
- 解决MATLAB对VS高版本出现“错误使用 mex未找到支持的编译器或 SDK。”的问题
- 【matlab】 matlab实现计算两张图片的相似度
- Matlab配置C++/VS等编译环境(mex -setup 找不到编译器问题解决)
- Matlab 机器人工具箱 机械臂工作空间仿真
- 《DSP using MATLAB》示例Example5.9
- 利用matlab/simulink绘制bode图,以及用origin绘制反正切函数
- Matlab 帮助文档不能打开解决办法