二重のfor文による行列への出力

上段右辺のインデックス付けされているclaやcd、clは56×41のデータがそれぞれ入っており、その他の文字は定数その他の文字は定数値です。c11～c33などはA,B行列の各要素でありc21以降はこの場では省略しております。二重のfor文によってA,B行列を56×41通りだしたいのですがA,B行列もインデックス付けされてしまいうまくいきません。

for j=1:size(cla,2)

for l=1:size(cla,1)

c11(l,j)=-ro.*v.*s(l,j).*cd(l,j)./m;

c12(l,j)=ro.*v.^2.*s(l,j)./(2.*m).*(-cd(l,j)+cl(l,j)./57.3);

c14=-g/57.3;

A(l,j)=[c11(l,j) c12(l,j) 0 c14;

c21(l,j) c22(l,j) c23 0;

c31(l,j) c32(l,j) c33(l,j) 0;

0 0 1 0];

B(l,j)=[0;

d21(l,j);

d31(l,j);

0];

end

end

理想は4×4のA,B行列が56×41通り出せることです。

Answer 1

[Atsushi Ueno]

「配列 A(l,j)=***」を「cell 配列 A{l,j}=***」に変更すると理想の形になります。

下記は1×4のベクトル(中身は適当な値)を 56×41 通り出した例です。

for j = 1:41
    for l = 1:56
        A{l,j} = [l+j l+j l j];
    end
end
A

A = 56x41 cell array
  Columns 1 through 10

    {[   2 2 1 1]}    {[   3 3 1 2]}    {[   4 4 1 3]}    {[   5 5 1 4]}    {[   6 6 1 5]}    {[   7 7 1 6]}    {[   8 8 1 7]}    {[   9 9 1 8]}    {[ 10 10 1 9]}    {[ 11 11 1 10]}
    {[   3 3 2 1]}    {[   4 4 2 2]}    {[   5 5 2 3]}    {[   6 6 2 4]}    {[   7 7 2 5]}    {[   8 8 2 6]}    {[   9 9 2 7]}    {[ 10 10 2 8]}    {[ 11 11 2 9]}    {[ 12 12 2 10]}
    {[   4 4 3 1]}    {[   5 5 3 2]}    {[   6 6 3 3]}    {[   7 7 3 4]}    {[   8 8 3 5]}    {[   9 9 3 6]}    {[ 10 10 3 7]}    {[ 11 11 3 8]}    {[ 12 12 3 9]}    {[ 13 13 3 10]}
    {[   5 5 4 1]}    {[   6 6 4 2]}    {[   7 7 4 3]}    {[   8 8 4 4]}    {[   9 9 4 5]}    {[ 10 10 4 6]}    {[ 11 11 4 7]}    {[ 12 12 4 8]}    {[ 13 13 4 9]}    {[ 14 14 4 10]}
    {[   6 6 5 1]}    {[   7 7 5 2]}    {[   8 8 5 3]}    {[   9 9 5 4]}    {[ 10 10 5 5]}    {[ 11 11 5 6]}    {[ 12 12 5 7]}    {[ 13 13 5 8]}    {[ 14 14 5 9]}    {[ 15 15 5 10]}
    {[   7 7 6 1]}    {[   8 8 6 2]}    {[   9 9 6 3]}    {[ 10 10 6 4]}    {[ 11 11 6 5]}    {[ 12 12 6 6]}    {[ 13 13 6 7]}    {[ 14 14 6 8]}    {[ 15 15 6 9]}    {[ 16 16 6 10]}
    {[   8 8 7 1]}    {[   9 9 7 2]}    {[ 10 10 7 3]}    {[ 11 11 7 4]}    {[ 12 12 7 5]}    {[ 13 13 7 6]}    {[ 14 14 7 7]}    {[ 15 15 7 8]}    {[ 16 16 7 9]}    {[ 17 17 7 10]}
    {[   9 9 8 1]}    {[ 10 10 8 2]}    {[ 11 11 8 3]}    {[ 12 12 8 4]}    {[ 13 13 8 5]}    {[ 14 14 8 6]}    {[ 15 15 8 7]}    {[ 16 16 8 8]}    {[ 17 17 8 9]}    {[ 18 18 8 10]}
    {[ 10 10 9 1]}    {[ 11 11 9 2]}    {[ 12 12 9 3]}    {[ 13 13 9 4]}    {[ 14 14 9 5]}    {[ 15 15 9 6]}    {[ 16 16 9 7]}    {[ 17 17 9 8]}    {[ 18 18 9 9]}    {[ 19 19 9 10]}
    {[11 11 10 1]}    {[12 12 10 2]}    {[13 13 10 3]}    {[14 14 10 4]}    {[15 15 10 5]}    {[16 16 10 6]}    {[17 17 10 7]}    {[18 18 10 8]}    {[19 19 10 9]}    {[20 20 10 10]}
    {[12 12 11 1]}    {[13 13 11 2]}    {[14 14 11 3]}    {[15 15 11 4]}    {[16 16 11 5]}    {[17 17 11 6]}    {[18 18 11 7]}    {[19 19 11 8]}    {[20 20 11 9]}    {[21 21 11 10]}
    {[13 13 12 1]}    {[14 14 12 2]}    {[15 15 12 3]}    {[16 16 12 4]}    {[17 17 12 5]}    {[18 18 12 6]}    {[19 19 12 7]}    {[20 20 12 8]}    {[21 21 12 9]}    {[22 22 12 10]}
    {[14 14 13 1]}    {[15 15 13 2]}    {[16 16 13 3]}    {[17 17 13 4]}    {[18 18 13 5]}    {[19 19 13 6]}    {[20 20 13 7]}    {[21 21 13 8]}    {[22 22 13 9]}    {[23 23 13 10]}
    {[15 15 14 1]}    {[16 16 14 2]}    {[17 17 14 3]}    {[18 18 14 4]}    {[19 19 14 5]}    {[20 20 14 6]}    {[21 21 14 7]}    {[22 22 14 8]}    {[23 23 14 9]}    {[24 24 14 10]}
    {[16 16 15 1]}    {[17 17 15 2]}    {[18 18 15 3]}    {[19 19 15 4]}    {[20 20 15 5]}    {[21 21 15 6]}    {[22 22 15 7]}    {[23 23 15 8]}    {[24 24 15 9]}    {[25 25 15 10]}
    {[17 17 16 1]}    {[18 18 16 2]}    {[19 19 16 3]}    {[20 20 16 4]}    {[21 21 16 5]}    {[22 22 16 6]}    {[23 23 16 7]}    {[24 24 16 8]}    {[25 25 16 9]}    {[26 26 16 10]}

  Columns 11 through 19

    {[ 12 12 1 11]}    {[ 13 13 1 12]}    {[ 14 14 1 13]}    {[ 15 15 1 14]}    {[ 16 16 1 15]}    {[ 17 17 1 16]}    {[ 18 18 1 17]}    {[ 19 19 1 18]}    {[ 20 20 1 19]}
    {[ 13 13 2 11]}    {[ 14 14 2 12]}    {[ 15 15 2 13]}    {[ 16 16 2 14]}    {[ 17 17 2 15]}    {[ 18 18 2 16]}    {[ 19 19 2 17]}    {[ 20 20 2 18]}    {[ 21 21 2 19]}
    {[ 14 14 3 11]}    {[ 15 15 3 12]}    {[ 16 16 3 13]}    {[ 17 17 3 14]}    {[ 18 18 3 15]}    {[ 19 19 3 16]}    {[ 20 20 3 17]}    {[ 21 21 3 18]}    {[ 22 22 3 19]}
    {[ 15 15 4 11]}    {[ 16 16 4 12]}    {[ 17 17 4 13]}    {[ 18 18 4 14]}    {[ 19 19 4 15]}    {[ 20 20 4 16]}    {[ 21 21 4 17]}    {[ 22 22 4 18]}    {[ 23 23 4 19]}
    {[ 16 16 5 11]}    {[ 17 17 5 12]}    {[ 18 18 5 13]}    {[ 19 19 5 14]}    {[ 20 20 5 15]}    {[ 21 21 5 16]}    {[ 22 22 5 17]}    {[ 23 23 5 18]}    {[ 24 24 5 19]}
    {[ 17 17 6 11]}    {[ 18 18 6 12]}    {[ 19 19 6 13]}    {[ 20 20 6 14]}    {[ 21 21 6 15]}    {[ 22 22 6 16]}    {[ 23 23 6 17]}    {[ 24 24 6 18]}    {[ 25 25 6 19]}
    {[ 18 18 7 11]}    {[ 19 19 7 12]}    {[ 20 20 7 13]}    {[ 21 21 7 14]}    {[ 22 22 7 15]}    {[ 23 23 7 16]}    {[ 24 24 7 17]}    {[ 25 25 7 18]}    {[ 26 26 7 19]}
    {[ 19 19 8 11]}    {[ 20 20 8 12]}    {[ 21 21 8 13]}    {[ 22 22 8 14]}    {[ 23 23 8 15]}    {[ 24 24 8 16]}    {[ 25 25 8 17]}    {[ 26 26 8 18]}    {[ 27 27 8 19]}
    {[ 20 20 9 11]}    {[ 21 21 9 12]}    {[ 22 22 9 13]}    {[ 23 23 9 14]}    {[ 24 24 9 15]}    {[ 25 25 9 16]}    {[ 26 26 9 17]}    {[ 27 27 9 18]}    {[ 28 28 9 19]}
    {[21 21 10 11]}    {[22 22 10 12]}    {[23 23 10 13]}    {[24 24 10 14]}    {[25 25 10 15]}    {[26 26 10 16]}    {[27 27 10 17]}    {[28 28 10 18]}    {[29 29 10 19]}
    {[22 22 11 11]}    {[23 23 11 12]}    {[24 24 11 13]}    {[25 25 11 14]}    {[26 26 11 15]}    {[27 27 11 16]}    {[28 28 11 17]}    {[29 29 11 18]}    {[30 30 11 19]}
    {[23 23 12 11]}    {[24 24 12 12]}    {[25 25 12 13]}    {[26 26 12 14]}    {[27 27 12 15]}    {[28 28 12 16]}    {[29 29 12 17]}    {[30 30 12 18]}    {[31 31 12 19]}
    {[24 24 13 11]}    {[25 25 13 12]}    {[26 26 13 13]}    {[27 27 13 14]}    {[28 28 13 15]}    {[29 29 13 16]}    {[30 30 13 17]}    {[31 31 13 18]}    {[32 32 13 19]}
    {[25 25 14 11]}    {[26 26 14 12]}    {[27 27 14 13]}    {[28 28 14 14]}    {[29 29 14 15]}    {[30 30 14 16]}    {[31 31 14 17]}    {[32 32 14 18]}    {[33 33 14 19]}
    {[26 26 15 11]}    {[27 27 15 12]}    {[28 28 15 13]}    {[29 29 15 14]}    {[30 30 15 15]}    {[31 31 15 16]}    {[32 32 15 17]}    {[33 33 15 18]}    {[34 34 15 19]}
    {[27 27 16 11]}    {[28 28 16 12]}    {[29 29 16 13]}    {[30 30 16 14]}    {[31 31 16 15]}    {[32 32 16 16]}    {[33 33 16 17]}    {[34 34 16 18]}    {[35 35 16 19]}

  Columns 20 through 28

    {[ 21 21 1 20]}    {[ 22 22 1 21]}    {[ 23 23 1 22]}    {[ 24 24 1 23]}    {[ 25 25 1 24]}    {[ 26 26 1 25]}    {[ 27 27 1 26]}    {[ 28 28 1 27]}    {[ 29 29 1 28]}
    {[ 22 22 2 20]}    {[ 23 23 2 21]}    {[ 24 24 2 22]}    {[ 25 25 2 23]}    {[ 26 26 2 24]}    {[ 27 27 2 25]}    {[ 28 28 2 26]}    {[ 29 29 2 27]}    {[ 30 30 2 28]}
    {[ 23 23 3 20]}    {[ 24 24 3 21]}    {[ 25 25 3 22]}    {[ 26 26 3 23]}    {[ 27 27 3 24]}    {[ 28 28 3 25]}    {[ 29 29 3 26]}    {[ 30 30 3 27]}    {[ 31 31 3 28]}
    {[ 24 24 4 20]}    {[ 25 25 4 21]}    {[ 26 26 4 22]}    {[ 27 27 4 23]}    {[ 28 28 4 24]}    {[ 29 29 4 25]}    {[ 30 30 4 26]}    {[ 31 31 4 27]}    {[ 32 32 4 28]}
    {[ 25 25 5 20]}    {[ 26 26 5 21]}    {[ 27 27 5 22]}    {[ 28 28 5 23]}    {[ 29 29 5 24]}    {[ 30 30 5 25]}    {[ 31 31 5 26]}    {[ 32 32 5 27]}    {[ 33 33 5 28]}
    {[ 26 26 6 20]}    {[ 27 27 6 21]}    {[ 28 28 6 22]}    {[ 29 29 6 23]}    {[ 30 30 6 24]}    {[ 31 31 6 25]}    {[ 32 32 6 26]}    {[ 33 33 6 27]}    {[ 34 34 6 28]}
    {[ 27 27 7 20]}    {[ 28 28 7 21]}    {[ 29 29 7 22]}    {[ 30 30 7 23]}    {[ 31 31 7 24]}    {[ 32 32 7 25]}    {[ 33 33 7 26]}    {[ 34 34 7 27]}    {[ 35 35 7 28]}
    {[ 28 28 8 20]}    {[ 29 29 8 21]}    {[ 30 30 8 22]}    {[ 31 31 8 23]}    {[ 32 32 8 24]}    {[ 33 33 8 25]}    {[ 34 34 8 26]}    {[ 35 35 8 27]}    {[ 36 36 8 28]}
    {[ 29 29 9 20]}    {[ 30 30 9 21]}    {[ 31 31 9 22]}    {[ 32 32 9 23]}    {[ 33 33 9 24]}    {[ 34 34 9 25]}    {[ 35 35 9 26]}    {[ 36 36 9 27]}    {[ 37 37 9 28]}
    {[30 30 10 20]}    {[31 31 10 21]}    {[32 32 10 22]}    {[33 33 10 23]}    {[34 34 10 24]}    {[35 35 10 25]}    {[36 36 10 26]}    {[37 37 10 27]}    {[38 38 10 28]}
    {[31 31 11 20]}    {[32 32 11 21]}    {[33 33 11 22]}    {[34 34 11 23]}    {[35 35 11 24]}    {[36 36 11 25]}    {[37 37 11 26]}    {[38 38 11 27]}    {[39 39 11 28]}
    {[32 32 12 20]}    {[33 33 12 21]}    {[34 34 12 22]}    {[35 35 12 23]}    {[36 36 12 24]}    {[37 37 12 25]}    {[38 38 12 26]}    {[39 39 12 27]}    {[40 40 12 28]}
    {[33 33 13 20]}    {[34 34 13 21]}    {[35 35 13 22]}    {[36 36 13 23]}    {[37 37 13 24]}    {[38 38 13 25]}    {[39 39 13 26]}    {[40 40 13 27]}    {[41 41 13 28]}
    {[34 34 14 20]}    {[35 35 14 21]}    {[36 36 14 22]}    {[37 37 14 23]}    {[38 38 14 24]}    {[39 39 14 25]}    {[40 40 14 26]}    {[41 41 14 27]}    {[42 42 14 28]}
    {[35 35 15 20]}    {[36 36 15 21]}    {[37 37 15 22]}    {[38 38 15 23]}    {[39 39 15 24]}    {[40 40 15 25]}    {[41 41 15 26]}    {[42 42 15 27]}    {[43 43 15 28]}
    {[36 36 16 20]}    {[37 37 16 21]}    {[38 38 16 22]}    {[39 39 16 23]}    {[40 40 16 24]}    {[41 41 16 25]}    {[42 42 16 26]}    {[43 43 16 27]}    {[44 44 16 28]}

  Columns 29 through 37

    {[ 30 30 1 29]}    {[ 31 31 1 30]}    {[ 32 32 1 31]}    {[ 33 33 1 32]}    {[ 34 34 1 33]}    {[ 35 35 1 34]}    {[ 36 36 1 35]}    {[ 37 37 1 36]}    {[ 38 38 1 37]}
    {[ 31 31 2 29]}    {[ 32 32 2 30]}    {[ 33 33 2 31]}    {[ 34 34 2 32]}    {[ 35 35 2 33]}    {[ 36 36 2 34]}    {[ 37 37 2 35]}    {[ 38 38 2 36]}    {[ 39 39 2 37]}
    {[ 32 32 3 29]}    {[ 33 33 3 30]}    {[ 34 34 3 31]}    {[ 35 35 3 32]}    {[ 36 36 3 33]}    {[ 37 37 3 34]}    {[ 38 38 3 35]}    {[ 39 39 3 36]}    {[ 40 40 3 37]}
    {[ 33 33 4 29]}    {[ 34 34 4 30]}    {[ 35 35 4 31]}    {[ 36 36 4 32]}    {[ 37 37 4 33]}    {[ 38 38 4 34]}    {[ 39 39 4 35]}    {[ 40 40 4 36]}    {[ 41 41 4 37]}
    {[ 34 34 5 29]}    {[ 35 35 5 30]}    {[ 36 36 5 31]}    {[ 37 37 5 32]}    {[ 38 38 5 33]}    {[ 39 39 5 34]}    {[ 40 40 5 35]}    {[ 41 41 5 36]}    {[ 42 42 5 37]}
    {[ 35 35 6 29]}    {[ 36 36 6 30]}    {[ 37 37 6 31]}    {[ 38 38 6 32]}    {[ 39 39 6 33]}    {[ 40 40 6 34]}    {[ 41 41 6 35]}    {[ 42 42 6 36]}    {[ 43 43 6 37]}
    {[ 36 36 7 29]}    {[ 37 37 7 30]}    {[ 38 38 7 31]}    {[ 39 39 7 32]}    {[ 40 40 7 33]}    {[ 41 41 7 34]}    {[ 42 42 7 35]}    {[ 43 43 7 36]}    {[ 44 44 7 37]}
    {[ 37 37 8 29]}    {[ 38 38 8 30]}    {[ 39 39 8 31]}    {[ 40 40 8 32]}    {[ 41 41 8 33]}    {[ 42 42 8 34]}    {[ 43 43 8 35]}    {[ 44 44 8 36]}    {[ 45 45 8 37]}
    {[ 38 38 9 29]}    {[ 39 39 9 30]}    {[ 40 40 9 31]}    {[ 41 41 9 32]}    {[ 42 42 9 33]}    {[ 43 43 9 34]}    {[ 44 44 9 35]}    {[ 45 45 9 36]}    {[ 46 46 9 37]}
    {[39 39 10 29]}    {[40 40 10 30]}    {[41 41 10 31]}    {[42 42 10 32]}    {[43 43 10 33]}    {[44 44 10 34]}    {[45 45 10 35]}    {[46 46 10 36]}    {[47 47 10 37]}
    {[40 40 11 29]}    {[41 41 11 30]}    {[42 42 11 31]}    {[43 43 11 32]}    {[44 44 11 33]}    {[45 45 11 34]}    {[46 46 11 35]}    {[47 47 11 36]}    {[48 48 11 37]}
    {[41 41 12 29]}    {[42 42 12 30]}    {[43 43 12 31]}    {[44 44 12 32]}    {[45 45 12 33]}    {[46 46 12 34]}    {[47 47 12 35]}    {[48 48 12 36]}    {[49 49 12 37]}
    {[42 42 13 29]}    {[43 43 13 30]}    {[44 44 13 31]}    {[45 45 13 32]}    {[46 46 13 33]}    {[47 47 13 34]}    {[48 48 13 35]}    {[49 49 13 36]}    {[50 50 13 37]}
    {[43 43 14 29]}    {[44 44 14 30]}    {[45 45 14 31]}    {[46 46 14 32]}    {[47 47 14 33]}    {[48 48 14 34]}    {[49 49 14 35]}    {[50 50 14 36]}    {[51 51 14 37]}
    {[44 44 15 29]}    {[45 45 15 30]}    {[46 46 15 31]}    {[47 47 15 32]}    {[48 48 15 33]}    {[49 49 15 34]}    {[50 50 15 35]}    {[51 51 15 36]}    {[52 52 15 37]}
    {[45 45 16 29]}    {[46 46 16 30]}    {[47 47 16 31]}    {[48 48 16 32]}    {[49 49 16 33]}    {[50 50 16 34]}    {[51 51 16 35]}    {[52 52 16 36]}    {[53 53 16 37]}

  Columns 38 through 41

    {[ 39 39 1 38]}    {[ 40 40 1 39]}    {[ 41 41 1 40]}    {[ 42 42 1 41]}
    {[ 40 40 2 38]}    {[ 41 41 2 39]}    {[ 42 42 2 40]}    {[ 43 43 2 41]}
    {[ 41 41 3 38]}    {[ 42 42 3 39]}    {[ 43 43 3 40]}    {[ 44 44 3 41]}
    {[ 42 42 4 38]}    {[ 43 43 4 39]}    {[ 44 44 4 40]}    {[ 45 45 4 41]}
    {[ 43 43 5 38]}    {[ 44 44 5 39]}    {[ 45 45 5 40]}    {[ 46 46 5 41]}
    {[ 44 44 6 38]}    {[ 45 45 6 39]}    {[ 46 46 6 40]}    {[ 47 47 6 41]}
    {[ 45 45 7 38]}    {[ 46 46 7 39]}    {[ 47 47 7 40]}    {[ 48 48 7 41]}
    {[ 46 46 8 38]}    {[ 47 47 8 39]}    {[ 48 48 8 40]}    {[ 49 49 8 41]}
    {[ 47 47 9 38]}    {[ 48 48 9 39]}    {[ 49 49 9 40]}    {[ 50 50 9 41]}
    {[48 48 10 38]}    {[49 49 10 39]}    {[50 50 10 40]}    {[51 51 10 41]}
    {[49 49 11 38]}    {[50 50 11 39]}    {[51 51 11 40]}    {[52 52 11 41]}
    {[50 50 12 38]}    {[51 51 12 39]}    {[52 52 12 40]}    {[53 53 12 41]}
    {[51 51 13 38]}    {[52 52 13 39]}    {[53 53 13 40]}    {[54 54 13 41]}
    {[52 52 14 38]}    {[53 53 14 39]}    {[54 54 14 40]}    {[55 55 14 41]}
    {[53 53 15 38]}    {[54 54 15 39]}    {[55 55 15 40]}    {[56 56 15 41]}
    {[54 54 16 38]}    {[55 55 16 39]}    {[56 56 16 40]}    {[57 57 16 41]}

[雄太]

回答ありがとうございます。同様にやってみたところA,B行列を56×41通り出すことはできました。しかし、このA,B行列はシステムツールボックスの制御系設計に使うつもりでしたが、以下のようにセルでは使えないとのエラーが出ました。

"lqr" コマンドは、クラス "cell" のモデルに対して使用できません。

どうにか対処することはできないでしょうか、もしこちらにも知見をお持ちでしたら伺いたいです。

[Atsushi Ueno]

lqr 関数が cell 配列を受け付けない、即ち一度に複数の方程式を計算出来ないので、行列を作った時と同様に for 文を回す必要があります。for 文と同様の繰り返し処理を行う cellfun 関数を使う手もあります。

[雄太]

ありがとうございます。cellfun関数でうまくいけました。