txrtxs 发表于 2010-8-3 07:17

mathematica的数值计算精度

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