I have a loop command in Maple, do you know how to write it in Matlab?

1 次查看（过去 30 天）
Mehdi
Mehdi 2022-9-15
回答： Walter Roberson 2022-9-15
Ran in:
I have a loop command in Maple, do you know how to write it in Matlab?
for k to 11 do
X[k+1] := eval(V-G1 . EQ, Equate(V, X[k])) :
end do
t=1;
syms p__1 tau__1 tau__2 tau__3
V = [tau__1; tau__2 ;tau__3; p__1]
V = 
EQ =[(1519010779748915*p__1)/9671406556917033397649408 + (37082872076992779027*tau__1)/4611686018427387904 - (10269089160577095*tau__2)/2361183241434822606848 - (2560198925000835947*tau__3)/37778931862957161709568 + 0.02632670329369235494973966905587*sin(25*t);
(4636910109289096185*tau__2)/576460752303423488 - (10269089160577095*tau__1)/2361183241434822606848 - (3843333523122587*p__1)/2417851639229258349412352 + (65247335564564535*tau__3)/4722366482869645213696 - 0.023998204867523233580749613411132*sin(25*t);
(65247335564564535*tau__2)/4722366482869645213696 - (2560198925000835947*tau__1)/37778931862957161709568 - (8969816766329343*p__1)/1208925819614629174706176 + (580819293998174151*tau__3)/72057594037927936 + 0.023995357692118428094048561849008*sin(25*t);
(6650483120930606653*p__1)/11258999068426240 + (1519010779748915*tau__1)/9671406556917033397649408 - (3843333523122587*tau__2)/2417851639229258349412352 - (8969816766329343*tau__3)/1208925819614629174706176]
EQ = 
G = jacobian(EQ,V);
G1 = inv(G);
X=[zeros(4,1)];
X=sym(X);
X(1:4 , 1)=[-0.0032740549266896589425167620590331*sin(25*t);
0.0029834606253022400734802395833198*sin(25*t);
-0.0029769443257114050228215642498688*sin(25*t);
-0.000000000000028494667434775680457822160620761*sin(25*t)]
X = 

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

采纳的回答

Walter Roberson
Walter Roberson 2022-9-15
syms V
for k = 1 : 11
X(k+1) = subs(V - G1 * EQ, V, X(k));
end

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by