Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2008-07-08
年龄: 44
帖子: 2
声望力: 0 ![]() |
![]()
程序如下:
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 请教,这是为什么? 数值积分怎么会出现复数? |
![]() |
![]() |
主题工具 | |
显示模式 | |
|
|
![]() |
||||
主题 | 主题作者 | 版面 | 回复 | 最后发表 |
[求助]急求钻石滤波器设计 | 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 |