Need help in coding this summation, please.

1 次查看（过去 30 天）
Luciano Campos
Luciano Campos 2016-3-11
评论： Ced 2016-3-12
Hi,
Can anybody help me with coding this summation into matlab, please:
where J is (3x6), A is (6x6) and Phi is (3x3x9).
I've tried with this code:
% for ii=1:hor
for mm=1:ii
G(:,:,ii)=sum(kron(J*A'^(ii-mm),Phi(:,:,mm)),3);
end
end %
but it doesn't seem to work, as I get incorrect results. I think the mistake is in the sum function, but I don't know other way of coding the summation (cumsum function doesn't work either).
I'd really appreciate some help here, as I'm completely stuck.
best,
  3 个评论
显示 1更早的评论
Luciano Campos
Luciano Campos 2016-3-11
Hi Ced, Thanks for your comment. I tried what you suggest, but it still won't work. The original formula is correct. If you develop the original summation you get:
I think the problem is that I don't know how to sum up consecutive terms as neither sum nor cumsum function seem to work.
Thanks for your help,
best,
Ced
Ced 2016-3-11
编辑：Ced 2016-3-11
I just found it strange to start G indexing a 1, but Phi indexing at 0.
Yeah, I messed up, sorry. Your original equation looks correct to me, but yes, you sum up too much. With sum, for each mm, you sum up over all of Phi again, which is not what you want.
Assuming you initialize G with zero:
for ii=1:hor
for mm=1:ii
G(:,:,ii) = G(:,:,ii) + kron(J*A'^(ii-mm),Phi(:,:,mm));
end
end

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2016-3-12
Is the x in a circle indicating cross-product or is it indicating exclusive or?
  1 个评论
Ced
Ced 2016-3-12
I believe it stands for the tensor product, hence the kronecker operator.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by