Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2008-07-08, 11:01   #1
pizizhou
初级会员
 
注册日期: 2008-07-08
年龄: 44
帖子: 2
声望力: 0
pizizhou 正向着好的方向发展
默认 [求助]不定积分算出来的结果为复数

程序如下:
syms t a b aa bb cc H u;
a=4200;b=2400;H=10200;
k=0.4829;
cc=0.449874264; %倾角gama
u=28.9;%偏心矩
rou=sqrt((a*sin(t))^2+(b*cos(t))^2);
L1=0.199533834*cos(t)+0.02077658*(cos(t))^3+0.0052831*(cos(t))^5;
L2=L1+0.001801055*(cos(t))^7+0.000708066*(cos(t))^9+0.000299304*(cos(t))^11;
Lt=a*0.013970421*t*180/pi-a*sin(t)*L2; %利用近似公式计算弧长
kt=H/4-k*Lt;

%m1aa0=-kt*a*sin(t)*cos(cc)/rou+a*cos(t)*sin(cc);
%mqaa60=(0.5*rou+u*(a^2+b^2)/(2*a*b)+u^2/rou)*sin(2*t);

m1aa0=-kt*a*sin(t)*cos(cc)/rou;
mqaa60=u^2/rou*sin(t)*cos(t);

m1aq7=int(m1aa0*mqaa60,t,0,pi/2);
vpa(m1aq7,12)


此时积分结果为复数:-122.360887054-.437256737137e-9*i
奇怪的是,假如把其中的m1aa0,mqaa60换成%后面的表达式,则计算结果为一实数:381974.875637

请教,这是为什么?
数值积分怎么会出现复数?
pizizhou 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码


相似的主题
主题 主题作者 版面 回复 最后发表
[求助]急求钻石滤波器设计 ailin MATLAB论坛 0 2008-09-05 23:53
【求助】怎么把Scope的接口增加2个或者更多啊 jerry1 MATLAB论坛 3 2008-04-09 23:03
[求助]请高手帮忙看一下这个关于微分方程组的程序,错在哪里~跪谢 柠檬小巫 MATLAB论坛 0 2008-04-07 20:58
【资料】Matlab教程之keyman系列四:Matlab处理优化问题 keyman MATLAB论坛 6 2007-11-29 11:17


所有时间均为北京时间。现在的时间是 12:55


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.