24小时热门版块排行榜    

CyRhmU.jpeg
查看: 319  |  回复: 3

驴肉火烧

金虫 (小有名气)

[求助] 编程问题

各位大侠能帮我看看我的matlab程序吗,运行不出来
%建立方程,这是一个分段函数,t1是分界点,v和s是关于x的函数
function dy=f(x,y)
global I1 I2 m2 m3 F s g R k1 k2 v s t1
if x     dy=2*sqrt(1/2*m2*v^2+1/2*m3*v^2+1/2*k2*s^2+(m2+m3)*g*s-F*s-1/2*k1*(R-R*cos(y))/(100+50));
else
    dy=2*sqrt(1/2*m2*v^2+1/2*k2*s^2+1/2*k3*s^2+m2*g*s-1/2*k1*(R-R*cos(y))/(I1+I2))
end
%编程
global I1 I2 m2 m3 F s g R k1 k2 v s t1
I1=100
I2=50
m2=1.5
m3=0.7
F=100
s=0.008
g=9.8
R=0.1
v=0.00003*x^4-0.0002*x^3+0.008*x^2
s=0.000006*x^5 -0.00005*x^4 + 0.003*x^3
t1=10
[x,y]=ode45('f',[0 15],0);
plot(x,y)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

驴肉火烧

金虫 (小有名气)

我编程那段少发了k1=600000
k2=45000
2楼2011-10-07 21:28:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖


驴肉火烧(金币+1): 谢谢你,是我的方程写的有问题,少了对k3的定义 2011-10-08 09:48:18
sunyang1988(金币+1): 谢谢提示 2011-10-11 17:28:35
再好好看看你的积分函数吧,是不是有误
3楼2011-10-08 09:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

驴肉火烧

金虫 (小有名气)

引用回帖:
3楼: Originally posted by change0618 at 2011-10-08 09:26:00:
再好好看看你的积分函数吧,是不是有误

好的,我再看看
4楼2011-10-08 09:31:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 驴肉火烧 的主题更新
信息提示
请填处理意见