886阶矩阵是小矩阵,根本不是什么大型问题
从你描述的现象来看矩阵元素的模可能太大了,先把A的LU分解算出来,然后把U的对角元做极分解z=re^{it},然后对log(r)求和,e^{it}部分则可直接相乘,这样即使行列式会溢出也能算出来
首先非常感谢您的回答我等会按照你的方法再弄一下!另外我想问一个问题我在进行3000*3000矩阵的逆矩阵转换时我用inv秩是满秩程序出现inv错误提示memory溢出有没有别的方法求逆谢谢!
你得写清楚你用的是什么软件,具体的错误信息是什么双精度3000阶的矩阵只需要72M的内存,如果要把逆矩阵显式求出来也只需要两倍的存贮,远小于现在常见的计算机内存