我试图为 CG 运行测试用例,我需要生成:
- 对称正定矩阵
- 大小 > 10,000
- 全密
仅使用矩阵索引,必要时使用 1 个向量(如)
条件数小于 1000。
我努力了:
A=rand(N,N)
使用然后生成随机矩阵A'A
使其对称。PD。[这增加了条件数]使用如图所示的矢量方法,但我似乎无法获得
(x,i,j)
可以确保 Sym 和 PD 的功能。
经过大量的实验,我想出了:
a(it,jt) = (vec(it)+vec(jt))/((it-1)^2+(jt-1)^2);
如果
a(it,it) = x(it)
如果
但这是 PD 直到大约 500x500。
- XLATM。[所有的分级和缩放,太难理解了。特别是因为我无法理解底层的线性代数]
有人可以给我一个满足上述要求的 x(向量)和 i,j(索引)函数吗?