[英]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.