Question

如何将此代码转换为for循环，以便我可以添加更多力量和时刻？这是很多代码，它基本上计算单位向量并找到x方向，y方向和z方向上的力的总和，并对时刻做同样的事情。

这仅适用于一个特定情况，因此我想将其置于for循环中，因此当应用新力时，它们将被考虑在内。

任何正确方向的帮助或推动都将是一个巨大的帮助！提前致谢。

load ('inputdata.mat')

A= zeros(6,6);
B= zeros (6,1);


% Calculate Force in x,y,z direction, sums them and puts them in the B 
% vector positions 1,2,3 

NF1=((ForceDir(1,2))^2+(ForceDir(1,3))^2+(ForceDir(1,4))^2)^(1/2);

F1i=(ForceDir(1,2)/NF1)* ForceDir(1,1);

F1j=(ForceDir(1,3)/NF1)* ForceDir(1,1);

F1k=(ForceDir(1,4)/NF1)*ForceDir(1,1);

NF2= ((ForceDir(2,2))^2+(ForceDir(2,3))^2+(ForceDir(2,4))^2)^(1/2);

F2i=(ForceDir(2,2)/NF2)* ForceDir(2,1);

F2j=(ForceDir(2,3)/NF2)* ForceDir(2,1);

F2k=(ForceDir(2,4)/NF2)* ForceDir(2,1);

B(1,1)= F1i+F2i; %External Force in x

B(2,1)= F1j+F2j; %External Force in y

B(3,1)= F1k+F2k; %External Force in z

% Calculate Moments in x,y,z direction and puts them in B vector 
% position 3,4,5

unitvectorForce1=[F1i F1j F1k];

MomentByForce1= cross(ForceCoor(1,:),unitvectorForce1);

unitvectorForce2= [F2i F2j F2k];

MomentByForce2= cross(ForceCoor(2,:),unitvectorForce2);

MomentNormal=((DirMoment(1,2))^2+(DirMoment(1,3))^2+(DirMoment(1,4))^2)^(1/2);

M1i= (DirMoment(1,2)/MomentNormal)*DirMoment(1,1);

M1j= (DirMoment(1,3)/MomentNormal)*DirMoment (1,1);

M1k= (DirMoment(1,4)/MomentNormal)*DirMoment(1,1);

unitvectorMoment= [M1i M1j M1k];

B(4,1)= MomentByForce1(1)+MomentByForce2(1)+ unitvectorMoment(1)

B(5,1)= MomentByForce1(2)+MomentByForce2(2)+ unitvectorMoment(2)

B(6,1)= MomentByForce1(3)+MomentByForce2(3)+ unitvectorMoment(3)

Answer 1

如果我理解正确，这应该可以解决问题：

for ii = 1:N
  NF(ii) = ((ForceDir(ii,2))^2+(ForceDir(ii,3))^2+(ForceDir(ii,4))^2)^(1/2);
  Fi(ii) = ForceDir(ii,2) / NF(ii) * ForceDir(ii,1);
  Fj(ii) = ForceDir(ii,3) / NF(ii) * ForceDir(ii,1);
  Fk(ii) = ForceDir(ii,4) / NF(ii) * ForceDir(ii,1);
  %% Etc
end

转换为for循环

1 个答案: