Hi , i want to build a matrix

[Jan]

Okay. What is your question? What is C, A, B and ny? What have you tried so far?

[Husam]

https://www.sciencedirect.com/science/article/abs/pii/S0019057809000111

[Husam]

I'm trying to simulate this article

[Jan]

Fine. You do not expect that the members of the forum read this article to find out, what C, A, B and ny is, hopefully. This is your turn. So please explain the details. Post the current code and ask a specific question.

[Husam]

tanx

Answer 1

[David Hill]

C=rand(1,3);
B=rand(3,1);
A=rand(3,3);
n=size(A,1);
a=eye(n);
y=4;
H=diag(repelem(C*B,n*y));
for k=-1:-1:-n*y+1
  a=a+A^-k;
  H=H+diag(repelem(C*a*B,n*y+k),k);
end
H

H = 12×12

0.9203 0 0 0 0 0 0 0 0 0 0 0 1.7417 0.9203 0 0 0 0 0 0 0 0 0 0 2.7749 1.7417 0.9203 0 0 0 0 0 0 0 0 0 3.9719 2.7749 1.7417 0.9203 0 0 0 0 0 0 0 0 5.3924 3.9719 2.7749 1.7417 0.9203 0 0 0 0 0 0 0 7.0677 5.3924 3.9719 2.7749 1.7417 0.9203 0 0 0 0 0 0 9.0469 7.0677 5.3924 3.9719 2.7749 1.7417 0.9203 0 0 0 0 0 11.3838 9.0469 7.0677 5.3924 3.9719 2.7749 1.7417 0.9203 0 0 0 0 14.1437 11.3838 9.0469 7.0677 5.3924 3.9719 2.7749 1.7417 0.9203 0 0 0 17.4028 14.1437 11.3838 9.0469 7.0677 5.3924 3.9719 2.7749 1.7417 0.9203 0 0

[Husam]

thank u very much