| 查看: 1314 | 回复: 0 | ||
[求助]
如何在matlab中计算对称矩阵的高精度本征函数？
|
|
matlab函数eig的矩阵算法是QZ算法，我用此算法计算一个密集矩阵（阶数仅为32X32）时得到的本征矢量中各个值的量级相差较多，从零点几到E-20多都有分布（如图所示）。 我知道由于浮点数会导致计算误差，但是大概由于eig函数进行了多次浮点数运算，最终导致无法把得到的本征矢量的真实值与误差值分开。 其实只有一次计算倒也无所谓，关键在于我做的程序还要将这些矢量不断迭代，继续求解，最终导致运算结果不准确。 所以，我想问的是： 1. matlab是否有程序包，可以用比较精确地算法，如jacobi-davidson算法，来得到较为精确地本征矢量，至少可以与浮点数运算的误差明显分开 2. 是否可以调整eig函数内部的计算精度解决这个问题？（单纯用eps只能用于单次计算，但是对于eig函数这种调用不好用） 3. 是否用c或者fortran会好一些？  1.png |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有34人回复
-
心脉受损
已经有8人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
博士申请都是内定的吗？
已经有8人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
关于MATLAB中的fft函数
已经有14人回复- 刚刚接触matlab，需要编写三阶矩阵，但其中有方程，不知道该怎么处理，求大神指导
已经有4人回复
- 病态矩阵的精确解求解
已经有17人回复
- matlab 将 数值或者矩阵copy到粘贴板上用什么命令
已经有4人回复
- 在matlab中如何提高矩阵求逆的运算精度？
已经有5人回复
- 我用matlab求矩阵的特性向量矩阵，怎么出错了？
已经有13人回复
- 求助matlab中含有三角矩阵的矩阵方程求解
已经有3人回复
- 用matlab计算两个向量之间的夹角结果出现复数
已经有4人回复
- （急急急）求助matlab大型矩阵求本征值 eigs命令使用后，提示错误
已经有7人回复
- fortran中如何调用matlab中求解特征值函数？
已经有9人回复
- 稀疏矩阵的处理
已经有12人回复
- 用matlab里的nlinfit如何拟合分段函数
已经有8人回复
- 怎么找出lapack函数及函数的使用方法
已经有8人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- Matlab中如何自己定义一个新函数，然后计算时可以调用？
已经有12人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- 菜鸟求助~matlab怎么样让一个1xn的向量中的每个元素是个矩阵？
已经有11人回复
- 【求助】关于matlab的高维数组
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人