声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1510|回复: 0

[绘图技巧] 转子跌落仿真程序

[复制链接]
发表于 2013-9-13 22:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
QQ图片20130913222420.jpg 求解这个微分方程。这里Fn,Ft都含有xy的函数
我需要画出横坐标x(-0.2,0.2),纵坐标y(-0.2,0.2)的图。可以我画出之后图是空白的。希望大侠能帮小弟看一下程序问题
M函数为 Guiji_fun.m
function f=Guiji_fun1(t,y)
         g=9.8;
         k=3.378e8;
         b=10/9;
         c=1000;
         m=1.1;
         Jr=1.9*10^(-4);
         s=1.25*10^(-4);
         r=5.875*10^(-3);
         u=1;
         F(1)=k*(sqrt(y(1)^2+y(2)^2)-s)^b+c*((y(2)*y(1)+y(3)*y(4))/(sqrt((y(1))^2)+(y(3))^2));
         F(2)=u*F(1);%F(t)
f=[y(2);-(F(1)*cos(atan(y(2)/y(1)))+F(2)*sin(atan(y(2)/y(1))))/m;y(4);-g+(F(2)*cos(atan(y(2)/y(1)))-F(1)*sin(atan(y(2)/y(1))))/m;y(6);-(F(2)*r)/Jr];


窗口执行
>> y0=[0,0,0,0,0,6000];
>> [t,y]=ode45(@Guiji_fun1,[-0.2,0.2],y0);
>> plot(y(:,1),y(:,3))
Warning: Imaginary parts of complex X and/or Y arguments ignored



回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 11:20 , Processed in 0.068170 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表