How do I create a for loop for the this example?

Question

Relly Syam 2021-4-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/813150-how-do-i-create-a-for-loop-for-the-this-example

回答： Walter Roberson 2021-4-26

 % 
clear 
clc 
syms E0 E1 E2 E3 c0 c1 c2 c3 t r0 r1 r2 r3 K0 K1 K2 K3 format long
c0=0 c1=1/3; c2=2/3; c3=1; r0=0; r1=1; r2=2; r3=3;
E0=[euler(r0,sym(c0)) euler(r0,c1) euler(r0,c2) euler(r0,c3)] 
E1=[euler(r1,c0) euler(r1,c1) euler(r1,c2) euler(r1,c3)] 
E2=[euler(r2,c0) euler(r2,c1) euler(r2,c2) euler(r2,c3)] 
E3=[euler(r3,c0) euler(r3,c1) euler(r3,c2) euler(r3,c3)] 
E=[E0' E1' E2' E3']
K0=[int(euler(r0,t), 0,c0) int(euler(r0,t), 0,c1) int(euler(r0,t), 0,c2) int(euler(r0,t), 0,c3)] 
K1=[int(euler(r1,t), 0,c0) int(euler(r1,t), 0,c1) int(euler(r1,t), 0,c2) int(euler(r1,t), 0,c3)] 
K2=[int(euler(r2,t), 0,c0) int(euler(r2,t), 0,c1) int(euler(r2,t), 0,c2) int(euler(r2,t), 0,c3)] 
K3=[int(euler(r3,t), 0,c0) int(euler(r3,t), 0,c1) int(euler(r3,t), 0,c2) int(euler(r3,t), 0,c3)] 
K=[K0' K1' K2' K3']
Ek=E-K 
F=[6*c0-3*c0^2;6*c1-3*c1^2;6*c2-3*c2^2;6*c3-3*c3^2] 
InvInv_Ek=inv(Ek) 
C=inv(Ek)*F 
Ua=E*C
Ue=[6*c0 6*c1 6*c2 6*c3]'

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Walter Roberson 2021-4-26

https://www.mathworks.com/help/symbolic/sym.euler.html

The function accepts vectors for both arguments.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2021-4-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/813150-how-do-i-create-a-for-loop-for-the-this-example#answer_684760

在 MATLAB Online 中打开

format long g
syms E0 E1 E2 E3 c0 c1 c2 c3 t r0 r1 r2 r3 K0 K1 K2 K3
c0=0; c1=1/3; c2=2/3; c3=1; r0=0; r1=1; r2=2; r3=3;
EvalAt = [c0, c1, c2, c3];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3], EA).', EvalAt, 'uniform', 0);
K = horzcat(ktemp{:}).'
K = 4×4
                         1                      -0.5                         0                      0.25
                         1        -0.166666666666667        -0.222222222222222          0.12037037037037
                         1         0.166666666666667        -0.222222222222222         -0.12037037037037
                         1                       0.5                         0                     -0.25