![](/img/trans.png)
[英]Decomposition of a Toeplitz matrix into sum of two circulant matrices (matlab/octave)
[英]Sum matrix error octave
我從形狀為N·2·M的矩陣開始,將其稱為x。 然后,我有另一個形狀為M的theta矩陣。以下代碼給我以下八度誤差:不符合參數(op1為2x2x2,op2為2x2x2)。
x(n+1,:,:) = x(n+1,:,:) + [cos(theta(:)) sin(theta(:))]
怎么了
PS。 這是完整的代碼:
M=30;
N=32;
for n=1:1:N
z = [1.0 0.0]
x(1,1,1:M) = z(1);
x(1,2,1:M) = z(2);
x(n+1,:,:) = x(n,:,:) + randn([1,2,M])
theta(:) = randn([M,1]);
x(n+1,:,:) = x(n+1,:,:) + [cos(theta(:)) sin(theta(:))]
end
size([cos(theta(:)) sin(theta(:))])
為30x2,而size(x(n+1,:,:))
為1x2x30
他們不匹配。
我不知道這在數學上是否是您所需要的(即,如果要加在一起的索引是應該的),但是您可以通過以下方法解決此問題:
x(n+1,:,:) = x(n+1,:,:) + reshape([cos(theta(:)) sin(theta(:))].',1,2,M);
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.