How to get rid of the error: Error using horzcat. Dimensions of matrices being concatenated are not consistent.

Question

Mallouli Marwa 2016-4-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/277760-how-to-get-rid-of-the-error-error-using-horzcat-dimensions-of-matrices-being-concatenated-are-not

评论： Mallouli Marwa 2016-4-7

Hi

When i solve a second order differential equation with constant coefficients, i obtained a solution. But when i change these coefficients by the matrix form i have obtained this error:

 Error using horzcat. Dimensions of matrices being concatenated are not consistent.

And my new code is a function :

function xdot = equadif(t,x)
% Function file for mass with spring and damping.
% Position is first variable, velocity is second variable.
F=[12;23:0];
M=[1,0,0;1,3,0;0,0,2];
cp = 1.38*10^-9;
K = [5,-6,1;-6,2,20;4,2,5];
TETA = [1,1,0;2,0,0;1,0,0];
C =[2,10,0;5,4,0;1,5,0];
R = 10*10^3;
A = [zeros(3,3), eye(3,3), zeros(3,1); -K*M^-1, -C*M^-1, TETA*M^-1; zeros(3,3), TETA/cp, -1/(R*cp)];
B = [0; M^-1 ; 0];
xdot = A*x+B*F;
end
[t,x]=ode45(@equadif,[0,0.01,0.03],[0,1,0]);
plot(t,x(:,1),'b',t,x(:,2),'r',t,x(:,3),'--')

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Muhammad Usman Saleem 2016-4-7

please format your code? difficult to understand

请先登录，再进行评论。

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

Answer 1

Roger Stafford 2016-4-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/277760-how-to-get-rid-of-the-error-error-using-horzcat-dimensions-of-matrices-being-concatenated-are-not#answer_216912

编辑：Roger Stafford 2016-4-7

在 MATLAB Online 中打开

Both the expressions for A and B are invalid. In A the three quantities

zeros(3,3), eye(3,3), zeros(3,1)

have a total of 7 columns, whereas

-K*M^-1, -C*M^-1, TETA*M^-1

have a total of 9 columns. These are incompatible. Also

zeros(3,3), TETA/cp, -1/(R*cp)

have only 7 columns. Moreover zeros(3,3) and TETA/cp have 3 rows whereas -1/(R*cp) has only one row.

In B, the three quantities have 1 column, 3 columns, and 1 column, respectively. These don't match.

Try laying out these various quantities with little dots and you will see that they don't combine in a rectangular pattern.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Mallouli Marwa 2016-4-7

I have changed the matrix A in the right form: TETA = [1;2;1]; A = [zeros(3,3), eye(3,3), zeros(3,1); -M^-1*K, -M^-1*C, -M^-1*TETA; (-1/R)*TETA', zeros(1,3), -1/(R*cp)]; But ther is still an error in B: B = [zeros(3,1); M^-1 ; zeros(1,3)]; How can i change it ??

请先登录，再进行评论。