How would I solve the following system of equations with a for loop?

Question

Alexandr 2023-11-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2054574-how-would-i-solve-the-following-system-of-equations-with-a-for-loop

移动：Torsten 2023-11-30

z1 = 0;
-3*z1 + 4*z2 - z3 = 0;
z.m-2 - 4*z_m-1 + 6*z.m - 4*z.m+1 + z.m+2 = (-mu * g - f_m) * (delta_x^4)/(E*I);
m = 3:M-2;
%M is a number that determines the size of the array
z.M-2 - 4*z.M-1 + 3*z.M = 0;
z.M = 0;

These are the five equations I need to solve. I tried using Ax = b method but I can't figure out how to get it to work with the third equation. Below is my code.

A = zeros(5,M); %Array A is every constant before each elongation:z from equations 1-5. Rows correspond to each equation, columns correspond to each z.
for m = 3:M-2
    A(3,(m-2):(m+2)) = [1,-4,6,-4,1]; % modifying array for the values of z from equation 3.
    
end
A(1,1) = 1; % modifying array for the values of z from equation 1.
A(2,1:3) = [-3,4,-1]; % modifying array for the values of z from equation 2.
A(4,(M-2):M) = [1,-4,3]; % modifying array for the values of z from equation 4.
A(5,M) = 1; % modifying array for the values of z from equation 5.
b3 = (-mu * g - f((round((M + 1)/2)))) * ((delta_x^4)/(E * I)); % b is a column vector of the left hand side of equations 1-5
b = [0;0;b3;0;0];
z = A\b; % Mx1 column vector z is found by using right side division.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Torsten 2023-11-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2054574-how-would-i-solve-the-following-system-of-equations-with-a-for-loop#answer_1363424

移动：Torsten 2023-11-30

在 MATLAB Online 中打开

Your matrix is not 5xM, but MxM.

The definition of A must read

A = zeros(M); %Array A is every constant before each elongation:z from equations 1-5. Rows correspond to each equation, columns correspond to each z.
A(1,1) = 1; % modifying array for the values of z from equation 1.
A(2,1:3) = [-3,4,-1]; % modifying array for the values of z from equation 2.
for m = 3:M-2
    A(m,(m-2):(m+2)) = [1,-4,6,-4,1]; % modifying array for the values of z from equation 3.
end
A(M-1,(M-2):M) = [1,-4,3]; % modifying array for the values of z from equation 4.
A(M,M) = 1; % modifying array for the values of z from equation 5.