如何检查一个复杂rectagular矩阵的列proportionaly吗?

1视图(30天)
我有一个复杂的矩阵G大小45的x 2,在magitude明智比第一第二列是非常小的。
在G创建的代码列所示代码。这里的附加文件g1和t。
我在寻找一个正确的方法来检查,如果列是成比例的 复杂的 rectagular矩阵G ?
清晰的所有;
关闭所有;
负载g1.mat
负载t.mat
G = (g1 g1。* t);

接受的答案

约翰D 'Errico
约翰D 'Errico 2023年2月22日
编辑:约翰D 'Errico 2023年2月22日
检查两个向量(在你的情况下,只是一个矩阵的列)比例是非常容易的,然后计算比例常数。你需要小心,你不能测试向量之间的相关系数,因为它是微不足道的,有两种截然不同的向量并不成比例,然而,分享一个完美的关联。例如向量x, vectot y = 1 + x的相关系数1,但他们并不成比例。
负载g1.mat
负载t.mat
类属性名称大小字节ans 1 x30 60字符x33 cmdout 66字符g1 45 x1 720双复杂t 45 x1 360双
t
t = 45×1
1.0 e 03 * -0.1735 -0.1824 -0.1939 -0.2054 -0.2143 -0.2183 -0.2166 -0.2094 -0.1986 -0.1867
情节(t)
不幸的是,我认为你误解了能做些什么,甚至比例到底意味着什么。g1和g1。* t彼此并不成正比。你multutiplied g1的元素,不恒定向量t, t不是常数。t / t的元素发生巨大的变化,所以,它甚至改变的迹象。所以g1不是正比于g1。* t。
例如,考虑向量
x = [1:5]”;
y = randn (5、1);
(x, y, y x。*)
ans = 5×3
1.0000 0.9764 0.9764 2.0000 0.6719 1.3439 3.0000 1.0125 3.0375 4.0000 1.7801 7.1205 5.0000 1.5934 7.9669
是否有任何可能的原因决定x, x。* y相互“比例”?不。
鉴于这一切,你会做什么呢?你可能决定制定问题,有一个标量恒定值,这样只有向量g1和g1。不知道向量t * t,我们可以近似的第二个向量与产品g1 *年代?
老实说,如果你知道向量t
s =意味着(t);
可能是你想要做什么。但你可以试着估计标量常数。最好的估计是:由反斜杠符:
g2 = g1。* t;
s_ = g1、g2
s_ e-21i e-05 = -8.2677 + 3.7221
正如你所看到的,因为g1和g2是复杂的估值,少量的虚部爬,但这只是浮点垃圾。
在这是一个小问题,更大的值(级)的g1成为那些有最大影响的结果估计。事实上,我们不会找到s_是一样的意思(t)。
意思是(t)
ans = -5.1546 e-05
然后相关的问题是,是真正的两个向量距离比例?这是一个有效的问题。我们可能会比较规范(g1 * s - g2),然后按比例缩小的g1的规范。
格式长g
规范(g1 * s_ (g2) /规范(g1)
ans =
0.000114243043261004
这将给我们一些meaure我们如何做。这里少数是好的。或者,也许我们可以制定它的人们经常使用,在附近1 R ^ 2是一件好事。所以:
1 -范数(g1 * s_ (g2) /规范(g1)
ans =
0.999885756956739
那么你会说模型是aruably体面,s_是一个体面的近似。最后,您可能决定绘制两个向量。
阴谋(1:元素个数(g1) (* s_ g1, g2),“o”)
警告:虚部复杂的X和Y参数被忽略。
或许问题是你需要欣赏比例意味着什么。

答案(1)

威廉•罗斯
威廉•罗斯 2023年2月22日
相关性是一个衡量两个向量的比例。向量是复杂时,相关性是复杂的。复杂的相关性的大小是在0和1之间,是一个合理的比例。
下面的代码显示每一列的意思,它显示了复杂的相关性和复杂的大小相关。
负载g1.mat
负载t.mat
G = (g1 g1。* t);
流(“G(: 1):意味着= %。3 e + * % .3e \ n ',真正的(平均(G(: 1))),图像放大(平均(G (: 1))))
G(: 1):指= 1.023 e-01 + i * 9.341 e-02
流(“G(:, 2):意味着= %。3 e + * % .3e \ n ',真正的(平均(G(:, 2))),图像放大(平均(G (:, 2))))
G(: 2):指= -3.074 e-06 + i * 2.825 e-05
流(= %的相关性。3 f + * % .3f \ n ',
真正的(corr (G (: 1), G(:, 2))),图像放大(corr (G (: 1), G (:, 2))))
相关我* -0.531 = 0.629 +
流(Abs (corr) = % .3f \ n '、abs (corr (G (: 1), G (:, 2))))
Abs (corr) = 0.823
祝你的工作。

类别

找到更多的在随机数生成帮助中心文件交换

社区寻宝

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

开始狩猎!