为啥矩阵乘法显示矩阵乘法维度不正确？是什么原因？

Uqzalq

为啥矩阵乘法显示矩阵乘法维度不正确？是什么原因？代码如下：
/ \$ b: z) \2 a5 u( h/ Q
% A/ u  D7 t6 H% `function sys = mdlDerivatives( t,x,u)
global uc ur
- T% L" k; V/ ?% }* b0 Oyd1 = u(1);
& L) k- |# n$ A* D' n( |2 g0 x. ]yd2 = u(2);
5 a) t$ D! Z! T- idyd1 = cos(t);
. \) h" i' ^+ m4 I) m, p. P( z- Gdyd2 = cos(t);
ddyd1 =- sin(t) ;
. M$ K: q0 f2 |5 g8 ^ddyd2 =- sin(t);
x1= u(3);x2= u(4);
% E8 v. S% W7 \6 Y/ Vx3 = u(5);x4 = u(6);
5 Y: [6 I- D5 W% ye1 = yd1 - x1;
e2 = yd2 - x3;
3 P9 i' f# \9 g; g' Kde1 = dyd1 - x2;
% \" s0 l# f! j" I" _1 ude2 = dyd2 - x4;
nmn1 = 20 ; nmn2 = 20 ;
beta11 = nmn1 ;beta21 = nmn2;
K0 = 15 *eye(2);
( p% A3 k- d# z6 T3 |" I  h9 rv1 = ddyd1 + beta11 * de1;
v2 = ddyd2 + beta21 * de2 ;
# i6 L3 ~" ?1 |: vv=[v1 v2]';
s1=de1+nmn1*e1;
s2=de2+nmn2*e2;
& E7 p9 p  A0 T# J7 us=[s1 s2]';
% Contro1 1aw(26)
6 k2 G7 k: s" LepO = 0.1;
Ip= eye(2);
  I- U, p2 @2 ?: ^rou0 = 0.2;
$ s* ~5 @* N' {# P0 E& G3 w! }6 ~! v/ Tu1(1) = exp( -1/2* ((x(1) +1.25)/0.6)^2);
u1(2) = exp( -1/2 * (x(1)/0.6)^2);
' |) z; C, t$ y  yu1(3) = exp( - 1/2 * ((x(1)-1.25)/0.6)^2);
" j5 H0 W, ]3 a( yu2(1) =exp( - 1/2* ((x(2) + 1.25)/0.6)^2) ;
u2(2) = exp( - 1/2 * (x(2)/0.6)*2) ;
u2(3) =exp( - 1/2* ((x(2) - 1.25)/0.6)*2) ;
u3(1) =exp( -1/2* ((x(3) + 1.25)/0.6)^2) ;
  u7 p& G0 f! w2 J9 Z% v+ Q  |  i, ^u3(2) = exp( - 1/2 *(x(3)/0.6)^2) ;
6 R% Q' J. w4 U- o- Du3(3) =exp( -1/2* ((x(3)-1.25)/0.6)^2);
& {, M) r8 ]  `' }6 Y) Mu4(1) =exp( - 1/2* ((x(4) + 1.25)/0.6)^2) ;
u4(2) = exp( - 1/2 * (x(4)/0.6)^2);
( T0 B/ Y6 Y: ?( H% }' g* wu4(3) = exp( -1/2* ((x(4)-1.25)/0.6)^2) ;
FS1 =0;
FS2 = zeros(3^4,1);
for L1= 1:3
    for L2 = 1:3
) U% U, J: _/ g: z' A        for L3 = 1:3
            for L4 = 1:3
                FS2(3^3 *(L1-1) + 3^2 *(L2 - 1) + 3*(L3 - 1) + L4) = u1(L1) * u2(L2) * u3(L3) *u4(L4);
                FS1= FS1+ u1(L1) * u2(L2)* u3(L3) * u4(L4);
/ S2 D" O/ E9 x, W8 Q            end
        end
    end
end
FS=FS2/FS1;
w_f1 = FS;w_f2 = FS;
w_g11 = FS;w_g12 = FS;
w_g21 = FS;w_g22 = FS;
thtaf1 = x(1:3^4);
thtaf2 = x(3^4 + 1:3^4 *2);
thtag11 = x(3^4 * 2 + 1:3^4 * 3 );
thtag12 = x(3^4 * 3 + 1:3^4 * 4);
thtag21 = x(3^4*4 + 1:3^4 * 5);
thtag22 = x(3^4*5 + 1:3^4* 6);
- ~# |6 \6 b' e5 ~- o0 w# ?- Sf1 = w_f1'*thtaf1 ;
8 o( C# F, t  b( rf2 = w_f2'*thtaf2;
g11 = w_g11'*thtag11;
2 _1 T% U4 e' P% a$ ]0 y1 n3 A/ Ng12 = w_g12'*thtag12;
g21 = w_g21 *thtag21 ;
% [( ]+ E) G- l8 V( H  pg22 = w_g22 *thtag22 ;
' W! ^3 f/ l3 ]0 Z9 w' O, Dep_f = x(3^4 *6 +1:3^4 *6 +2) ;
7 W) w* R- x6 \* A! b9 hep_g11 = x(6*3^4 + 3);
ep_g12 = x(6 * 3*4 + 4);
ep_g21 = x(6 * 3*4 + 5) ;
, q' V' I- v9 C- M/ w5 |7 y6 Yep_g22 = x(6 * 3*4 + 6);
ep_g = [ep_g11 ep_g12 ; ep_g21 ep_g22];
ep_u = x( 3^4 * 6 + 7);
" D8 k9 \5 t1 D% E( U( W& q  o* _  fdelta = x(3^4* 6 + 8);
& t) v" z- k2 u2 u$ V' EF= [f1 f2]';
( P# \5 I: i, W2 ~G= [g11 g12;g21 g22];
7 h# V, H7 ]; H) j' Q0 o- |k1 = 20;
uc = G'*inv(ep0 *Ip + G*G')*(-F + v+K0*s+ k1*G*s);
uO = epO*inv(epO * Ip + G*G') * (-F + v + K0*s + k1*G*s);
) k$ @" E$ I, R3 ]ur = s * abs(s') * (ep_f + ep_g * abs(uc) + ep_u' * abs(u0))/(norm(s)^2 + de1ta);
xite_f1= 0.5;
xite_f2 = 0.5;
$ Z0 K& D1 t" m0 D$ }; Oxite_g11 = 0.5;
xite_g12 = 0.5;
xite_g21 = 0.5;
6 @3 p7 [; Y" ^4 N9 Qxite_g22 = 0.5;
xite0 = 0.001;
%F
for i = 1:1:3^4
& S$ n+ O6 Z/ i8 N8 u; b    sys(i) =- xite_f1 * w_f1(i) * s(1);
! V8 m& [; s, Q8 E$ S3 aend
* x# t. P+ p. z4 W! Wfor i= 3^4 + 1:1:2*3^4
# \! G' v. R7 v- l( M8 k& z    sys(i)=-xite_f2 * w_f2(i-3^4) * s(2);
5 e2 O: X6 x+ ?. w% F4 iend
%G
for i=2*3^4+1:1:3*3^4
/ O+ S4 H6 v! g7 L    sys( i) =- xite_g11*w^g11(i-2*3^4) * s(1)*(uc(1) - k1*s(1));
" W6 v+ I1 W% s' H) @+ [& Jend
+ {% }; F! O. `+ \for i = 3*3^4 + 1:1:4*3^4
    sys(i) =- xite_g12 * w_g12(i-3*3^4) * s(1) * (uc(2)-k1*s(2));
end
for i = 4*3^4 + 1:1:5*3^4
    sys(i) =-xite_g21*w_g21(i - 4*3^4) * s(2) * (uc(1) - k1*s(1));
end
for i=5*3^4 + 1:1:6*3^4
' l; Z, \( Y/ Z# Y* W    sys(i) =- xite_g22 * w_g22(i - 5 *3^4) * s(2) * (uc(2) - k1*s(2));
/ c& Q2 j1 D; ?6 Q; H0 L. nend
, s7 c  }: z9 `7 G: _gama0 = 0.001;

sys(6 * 3^4 + 1) = gama0*(abs(s(1))) ; % ep_f1
sys(6 * 3^4 + 2) = gama0* (abs(s(2)));% ep_f2
+ t( U! ~' I2 P  r: @%ep_g
sys(6 *3^4 + 3)        = gama0 * (abs(s(1)))        * abs(uc(1) );        % ep_g1 1
sys(6 * 3^4 + 4)= gama0 * (abs(s(1)))        * abs(uc(2)) ;        % ep_g12
$ R, ]  u, D, E4 isys(6 * 3^4 + 5)= gama0 * (abs(s(2)))        * abs(uc(1)) ;        % ep_g21
sys(6 * 3^4 + 6)= gama0 * (abs(s(2)))        * abs(uc(2)) ;        % ep_g22
% ep_u
sys(6 * 3^4 + 7) = gama0*(abs(s')) *abs(u0);
) u# i' A1 `7 B$ D- x% de1ta
0 s- _/ e& R, f6 N& \$ zsys(6 * 3^4 + 8) =- xite0* (abs(s') * (ep_f + ep_g * abs(uc) + ep_u' *abs(uO)))/(norm(s)^2 + de1ta) ;
function sys = mdlOutputs(t,x,u)
8 w3 P# m* R" u! |global uc ur
' o0 d7 C' M! W/ _" D7 oif t == 0
uc = zeros(2,1);
: c* V# ^. p* t( d6 y/ fur = zeros(2,1);
end
ut = uc + ur;
, H, \' U. A  wsys(1) = ut(1);
sys(2)=ut(2);

CCxiaom

帮顶一下

piday123

代码太长，别人也不是你课题组的成员，猜不出你每一步都是想实现具体什么步骤的。这种情况下，一般都是建议你自己一步一步测试到底哪里有问题。
$ u) b0 z6 q/ Q5 ^比如先把function的那行加个百分号变成注释，然后自己随便生成个测试用的数组如 t = 0 : 1000; x = 1 : 1000; u = 2 : 1000; 之类的开始逐行运行，看哪步报错。
比如在这个测试数组下，这几行
2 G$ l9 e" d* ^# c

• f1 = w_f1'*thtaf1 ;
• f2 = w_f2'*thtaf2;
• g11 = w_g11'*thtag11;
• g12 = w_g12'*thtag12;
  2 Z: Y3 @" f, E  J4 R6 i

* y- x* e! y( ~' K
维度就不对，而且一查就发现加'转置导致的。
后面运行还会出现 ep0 这种未定义的量被调用来进行计算的情况，我就不一一列举了，你自己完全可以逐行排查的。
还是那句话，这里的坛友不是你课题组的成员，顶多用测试数据看看哪里运行可能会出错，具体怎么改，只有设置全部算法步骤的你本人最清楚。

6 W( r" y& ?$ y2 M8 G
  `7 A% V& w! s& t* p7 `6 S

Uqzalq

piday123 发表于 2020-7-8 19:20
! [$ R" z9 ^3 ~; p8 x代码太长，别人也不是你课题组的成员，猜不出你每一步都是想实现具体什么步骤的。这种情况下，一般都是建议 ...

2 D: n3 D- K; k% M, w代码是有点长，谢谢啦，我一步一步测试找问题

xiaogegepcb

看看问题咋解决的