Info

Questa domanda è chiusa. Riaprila per modificarla o per rispondere.

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

1 visualizzazione (ultimi 30 giorni)
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 Commenti
John D'Errico
John D'Errico il 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.

Risposte (0)

Community Treasure Hunt

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

Start Hunting!

Translated by