Info

This question is closed. Reopen it to edit or answer.

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

1 view (last 30 days)
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 Comments
John D'Errico
John D'Errico on 26 Nov 2021
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.

Answers (0)

Community Treasure Hunt

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

Start Hunting!