1.基底rho_i既不能为0,不然K奇异;也不能为很小,不然u表达式中系数很大,这样优化结果很烂
解决办法:rho_i中非0元:线性变换到 [Min,Max];
rho_i中0元:当i为奇数,Min;当i为偶数,-Min。
2.基底rho_i在组装成K的时候,有相消的部分,消不干净剩下很小的数,也会使得u表达式中系数很大
解决办法:测试了下,nx,ny,nz都为奇数就没有相消,有偶数的时候,只有18,19,19的情形下可以。这个与力F也应该有关,暂时不测了,都取奇数应该可以。
3.C44到底变不变
a.确定U^T:E:V展开式子是啥,从计算角度确定下
解决:U^T:E:V展开式子是爱因斯坦和。
网友评论