作业帮 > 数学 > 作业

matlab用四阶龙格库塔法解微分方程组;

来源:学生作业帮 编辑:大师作文网作业帮 分类:数学作业 时间:2024/11/10 07:50:05
matlab用四阶龙格库塔法解微分方程组;
dy(1)=-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));
dy(2)=(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);
以上述计算结果再算一下表达式;
x=3000*sin(70*2*pi/360)-y(1)*sin(y(2));
z=3000*cos(70*2*pi/360)-12*t-y(1)*cos(y(2));
(每一组x,z对应一组y1,y2)并画出x,z的图像.
matlab用四阶龙格库塔法解微分方程组;
function dy=test(t,y)
dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);];
[t,y]=ode45('test',[0.01,1],[1,1])
plot(t,y(:,1),t,y(:,2));
x=3000*sin(70*2*pi/360)-y(1).*sin(y(2));
z=3000*cos(70*2*pi/360)-12*t-y(1).*cos(y(2));
plot(t,x,t,z)
最后的这图像不好.你可再看看.方法是没问题的.
你可到MATLAB中文论坛去请教关于MATLAB的问题,那儿高手多的很.