声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1239|回复: 4

[编程技巧] 恳请各位高手们救助啊。关于应用ode45计算时,出现了如下结果,

[复制链接]
发表于 2011-10-31 18:59 | 显示全部楼层 |阅读模式

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

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

x
编程一:
clc,clear
Pr=25000;
Pa=101000;
D=0.028;
L0=0.137;
m=1.4;
V0=9.5852e-005;
x0=[0;-1.1514];
t_final=0.018;
[t,y]=ode45('oppo',[0,t_final],x0,[]);
figure(1)
plot(t,y(:,1))
figure(2)
plot(t,y(:,2))
编程二:
function dx=oppo(t,x)
Pr=25000;
Pa=101000;
D=0.028;
L0=0.137;
m=1.4;
V0=1/4*pi*D^2*L0+1/6*pi*D^3;
dx=[x(2)
    (-(Pr+Pa)*(V0/((acos((D-x(1))/D)*D*(D-x(1))+(1/4*pi-1/2*acos((D-x(1))/D))*D^2-(1/2*D-1/2*x(1))*(2*D*x(1)-x(1).^2)^(1/2))*(L0+1/2*pi*x(1))+1/6*pi*(D-x(1)).^3))^m*(1/(1-(D-x(1)).^2/D^2)^(1/2)*(L0+1/2*pi*x(1))+1/2*acos((D-x(1))/D)*D*pi-m/((acos((D-x(1))/D)*D*(D-x(1))+(1/4*pi-1/2*acos((D-x(1))/D))*D^2-(1/2*D-1/2*x(1))*(2*D*x(1)-x(1).^2)^(1/2))*(L0+1/2*pi*x(1))+1/6*pi*(D-x(1)).^3)*((1/(1-(D-x(1)).^2/D^2)^(1/2)*(D-x(1))-acos((D-x(1))/D)*D-1/2*D/(1-(D-x(1)).^2/D^2)^(1/2)+1/2*(2*D*x(1)-(x(1)).^2)^(1/2)-1/2*(1/2*D-1/2*x(1))/(2*D*x(1)-x(1).^2)^(1/2)*(2*D-2*x(1)))*(L0+1/2*pi*x(1))+1/2*(acos((D-x(1))/D)*D*(D-x(1))+(1/4*pi-1/2*acos((D-x(1))/D))*D^2-(1/2*D-1/2*x(1))*(2*D*x(1)-x(1).^2)^(1/2))*pi-1/2*pi*(D-x(1)).^2)*acos((D-x(1))/D)*D*(L0+1/2*pi*x(1)))+Pa*(1/(1-(D-x(1)).^2/D^2)^(1/2)*(L0+1/2*pi*x(1))+1/2*acos((D-x(1))/D)*D*pi))*x(1)/m];
运行后出现了:Warning: Imaginary parts of complex X and/or Y arguments ignored,原函数里有个开根号的,但是手动算过几个数,并不是复数啊,是不是程序本身出问题了呢?高手指点啊……

回复
分享到:

使用道具 举报

 楼主| 发表于 2011-10-31 19:14 | 显示全部楼层
发表于 2011-11-1 00:34 | 显示全部楼层
个人时间/水平有限, 仅建议
若LZ确定没复数, 设些断点或增加些判断式找出那裡出现复数, 再看看是否合理
发表于 2011-12-22 10:13 | 显示全部楼层
我的问题跟楼主类似   不知你的问题解决没有?
发表于 2011-12-24 23:38 | 显示全部楼层
有空试了下LZ的程式, 但根本报错不同!?
Error: File: oppo.m Line: 6 Column: 4
The input character is not valid in MATLAB statements or expressions.

本想试著改, 但LZ的式子那麼长又很乱, 所以放弃了... sorry
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-9 00:37 , Processed in 0.064252 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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