Info

此问题已关闭。 请重新打开它进行编辑或回答。

I need assistance in writing a matlab code for this iterative example

1 次查看(过去 30 天)
TRYING TO SOLVE FOR u and finding maximum iterations needed to solve
Q1 = [ 25, 10, 15, 20; 10, 12, 18, 23; 12, 15, 17,19 ];
F1 = [ 6;7;12];
F2 = [ 48;64;72];
iter = 0;
N = 100; % maximum number of iterations
uold = u0;
QA = Q1;
F1 = RF1;
F2 = RFM; % max load capacity
u0 = u1;
dF = 0.2;
[row,col] = size(u);
tol = 0.0001;
syms 'u'
R1 = F2 - F1;
u2 = inv(Q1)*R1;
u3 = u0 + u2;
F3 = F1 + dF;
R2 = F2 - F3; % unbalanced force must be equal to 0 in order to converge
Q2 = F3\u3;
u4 = inv(Q2)*R2;
u5 = u3 + u4;
F4 = F3 + dF;
R3 = F2 - F4;
Q3 = F4\u5;
u6 = inv(Q3)*R3;
u7 = u5 + u6;
F5 = F4 + dF;
R4 = F2 - F5;
Q4 = F5\u7;
u8 = inv(Q4)*R4;
u9 = u7 + u8;
F6 = F5 + dF;
R5 = F3 - F6;
  2 个评论
John D'Errico
John D'Errico 2021-11-26
You do understand that the inverse of a matrix is undefined when the matrix is not square?
Q1 = [ 25, 10, 15, 20; 10, 12, 18, 23; 12, 15, 17,19 ]
Q1 = 3×4
25 10 15 20 10 12 18 23 12 15 17 19
size(Q1)
ans = 1×2
3 4
Therefore no number of iterations will suffice, sicne your code will error out almost immediately.

回答(0 个)

Community Treasure Hunt

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

Start Hunting!

Translated by