声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2865|回复: 0

[mathematica] mathematica的数值计算精度

[复制链接]
发表于 2010-8-3 07:17 | 显示全部楼层 |阅读模式

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

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

x
最近在做一个有限元分析,发现用mathematica编的程序生成的单元质量和刚度矩阵没有对称,具体如下:
单元刚度矩阵eleKK=Transpose[bK1]. bK2. bK1
其中bK1为8乘240阶,bK2为8乘8阶矩阵。
bK2经验证为对称矩阵,从单元刚度矩阵表达式分析,应该为对称,但eleKK==Transpose[eleKK]为False,很奇怪,大家有没有遇到这种情况,急问!!!!!
附件后缀名改为.nb即可运行。

test.txt

31.27 KB, 下载次数: 2

源程序

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-23 11:34 , Processed in 0.064288 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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