汉克尔矩阵,随机的条目

75(30天)

显示旧的评论

哈贾尔Alshaikh 2023年4月8日在一37

0
链接

这个问题直接联系

https://es.mathworks.com/matlabcentral/answers/1943584-hankel-matrix-random-entries

评论道: 哈贾尔Alshaikh2023年4月9日23:25

我有一个汉克尔矩阵X:

                          第一行= 1:n;
                         
                          col1 = (n + 1): (2 * n - 1);
                         
                          %创建满秩汉克尔矩阵的大小n ^ 2
                         
                          H_full =汉克尔(第一行,col1);
                         
                          %计算低秩近似
                         
                          [U, V] =圣言(H_full);
                         
                          U = U (:, 1: r);
                         
                          S = (1: r, 1: r);
                         
                          V = V (:, 1: r);
                         
                          X = U * * V ';

我想选择随机反对角的,让他们在一组叫做ω

例如我有5×5矩阵,我想选择两个随机反对角线,让里面的索引条目属于这些一组叫做ω。所以在5 * 5矩阵有9反对角线,我想选择第二和第四antidiagonalsω= {(1、2),(2,1),(1,4),(2、3),(2),(4,1)}

11日评论
显示隐藏 10年长的评论

Dyuman Joshi 2023年4月8日于5:05

编辑:Dyuman Joshi 2023年4月8日在38

代码中有未定义的变量。

                                n = 5;
                               
                                第一行= 1:n;
                               
                                col1 = (n + 1): (2 * n - 1);
                               
                                %创建满秩汉克尔矩阵的大小n ^ 2
                               
                                H_full =汉克尔(第一行,col1);

H_full不是大小nxn,随着col1 n - 1的元素。如果你想H_full nxn大小,改变

col1 = (n + 1): (2 * n);

指数的反对角

                                z = rot90(重塑(1:n ^ 2, n, n));
                               
                                % Antidiagonals
                               
                                印第安纳州= 2 [4];
                               
                                % ctr = min(印第安纳州,2 * n-ind);
                               
                                ω=细胞(1,元素个数(印第安纳州));
                               
                                为k = 1:元素个数(印第安纳州)
                               
                                [r、c] = ind2sub ([n n],诊断接头(z,印第安纳州(k) - n));
                               
                                kω{1,}= [r c];
                               
                                结束
                               
                                ω{1}
                               
                                   ans =
                                   2×2
                                  
                                    1 2 2 1

                                ω{2}
                               
                                   ans =
                                   4×2
                                  
                                    1 4 2 3 3 2 4 1

哈贾尔Alshaikh 2023年4月9日在0:14

非常感谢你muchy你工作很好但是我想生成任何随机抽样的方式,而不是印第安纳=[2 - 4],然而,我写了如下运行并不阻止我不没有的原因

                                m = n;
                               
                                mn = m * n;
                               
                                n = 1000;
                               
                                d = 0.4;
                               
                                r = 10;
                               
                                好吧= 0;
                               
                                z = rot90(重塑(1:锰、n, n));
                               
                                而~好吧
                               
                                % Antidiagonals
                               
                                %印第安纳= [1:2 * n - 1]
                               
                                印第安纳州= randperm (n)(2 * 1),圆(((一维)* (2 * n) 1) / 10);
                               
                                % ctr = min(印第安纳州,2 * n-ind);
                               
                                ω=细胞(1,元素个数(印第安纳州));
                               
                                为k = 1:元素个数(印第安纳州)
                               
                                [s、c] = ind2sub ([n n],诊断接头(z,印第安纳州(k) - n));
                               
                                kω{1,}= c [s];
                               
                                结束
                               
                                ω= sub2ind (n [m],年代,c)
                               
                                %包括指数
                               
                                (I, J) = ind2sub (n [m],ω)
                               
                                φ=稀疏(I, J 1 m, n)
                               
                                sumr =全(sum(φ,1))
                               
                                sumc =全(sum(φ,2))
                               
                                如果任何(sumr < r) | |任何(sumc < r)
                               
                                disp (“不”)
                               
                                继续
                               
                                其他的
                               
                                好吧= 1;
                               
                                结束
                               
                                结束

Dyuman Joshi 2023年4月9日在7:08

在定义有一个paranthesis人士印第安纳州的从而得到2的输出

                                n = 1 e3; d = 0.4;
                               
                                % |
                               
                                % v
                               
                                印第安纳州= randperm (n)(2 * 1),圆(((一维)* (2 * n) 1) / 10)
                               
                                   印第安纳州=
                                   1×1999
                                  
                                    1138 523 1858 604 1644 639 170 1933 801 1147 1753 737 754 393 1710 1800 709 449 544 1295 657 444 249 888 1681 1642 1542 1404 1626 1602

                                  ans = 120

                                %纠正
                               
                                印第安纳州= randperm ((2 * n) 1,圆(((一维)* (2 * n) 1) / 10))
                               
                                   印第安纳州=
                                   1×120
                                  
                                    1308 1849 1094 764 1949 1723 445 439 1404 341 753 331 358 846 1377 226 79 1929 155 1786 481 1348 98 419 276 755 1092 1053 664

我不明白这条线

ω= sub2ind (n [m],年代,c)

你为什么盖写ω后循环?

哈贾尔Alshaikh 2023年4月9日23:25

我不是与matlab多好,我不知道你是否明白我想要的正确与否

登录置评。

在回答这个问题。

答案(0)

在回答这个问题。

类别

MATLAB 语言基础知识矩阵和数组创建和连接矩阵

找到更多的在创建和连接矩阵在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的

汉克尔矩阵,随机的条目

11日评论
显示隐藏 10年长的评论

答案(0)

另请参阅

类别

标签

社区寻宝

汉克尔矩阵,随机的条目

11日评论 显示隐藏 10年长的评论

答案(0)

另请参阅

类别

标签

社区寻宝

11日评论
显示隐藏 10年长的评论