声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2139|回复: 1

[Maple] maple计算问题

[复制链接]
发表于 2008-10-11 11:09 | 显示全部楼层 |阅读模式

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

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

x
> restart:
> omega22:=-k[2]*(k[1]^2*J^2*m+k[1]*J*m*sqrt(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[1]*k[2])-2*S^2*k[1]*sqrt(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[1]*k[2])-2*k[1]*J*k[2]*m^2+4*S^2*k[1]*k[2]*m-k[2]*m^2*sqrt(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[1]*k[2])+k[2]^2*m^3)/(-m*J+S^2)/(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[1]*k[2])^(1/2)/(-k[2]*m+J*k[1]+sqrt(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[1]*k[2])):
> omega[22]:=1/2/(J*m-S^2)*(k[2]*m+J*k[1]-sqrt(k[2]^2*m^2-2*J*k[1]*k[2]*m+J^2*k[1]^2+4*S^2*k[2]*k[1])):
> simplify(omega[22]-omega22);
                                  0
> l:=244;p:=1.29;theta[0]:=Pi/4;h:=10/1000;c[1]:=-18;c[2]:=3.8;c[3]:=6;d[1]:=-1.02;d[2]:=0.99;d[3]:=-0.21;M[L]:=1.7;a:=0.4;dd:=0.028:M:=h*dd*0.9:dh:=dd+h*sin(theta[0]);U:=y;m:=3*l*(M[L]+M)*(2/Pi)^2;J:=(2/Pi)^2*l*(M[L]+M)*(a)^2;mu[1]:=2*m*0.04*sqrt(k[1]/m):mu[2]:=(2*J*0.06*sqrt(k[2]/J)):a[1] := (-1/2*p*U*dh*c[1]-p*U*dh*c[2]*theta[0]-3/2*p*U*dh*c[3]*theta[0]^2)*l;b[1] := (-1/2*p*U*dh^2*d[1]-p*U*dh^2*d[2]*theta[0]-3/2*p*U*dh^2*d[3]*theta[0]^2)*l; k[1]:=3*4*26000/l:k[2]:=864:S:=3*l*(2/Pi)^2*M*dd/2;a[3] := -1/2*p/U*dh*c[3]*l;b[3] := -1/2*p/U*dh^2*d[3]*l;
> simplify(omega[22]-omega22);
                             35.10262596
这是我的程序,蓝色部分是显示的计算结果。可以看到直接omega[22]-omega22进行化简,结果是0,这说明omega[22]和omega22是相等的。但是将各个参数的值都给出来以后,omega[22]-omega22=35.10262595,两者竟然不相等了。有哪位高手给解释一下。
回复
分享到:

使用道具 举报

发表于 2008-10-25 11:51 | 显示全部楼层

回复 楼主 zhanghb 的帖子

中间是不是出现了一些舍入的问题?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-23 17:02 , Processed in 0.080248 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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