rewrite A=[31] matrix to [33][31(var)]=A

Question

shahin hashemi 2017-12-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/372016-rewrite-a-3-1-matrix-to-3-3-3-1-var-a

评论： Star Strider 2017-12-9

采纳的回答： Star Strider

在 MATLAB Online 中打开

dear all

i have 3*1 matrix like matrix E below

syms qd td ed q teta
E =[ ed*cos(q)*sin(teta) - cos(teta)*(qd*cos(q)*sin(teta) + td*cos(teta)*sin(q)) - td*sin(q)*sin(teta)^2 ; cos(teta)*(td*cos(teta)*sin(q) - qd*sin(q)*sin(teta)) + ed*sin(q)*sin(teta) + td*cos(q)*sin(teta)^2  ;  ed*cos(teta) + cos(q)*sin(teta)*(qd*cos(q)*sin(teta) + td*cos(teta)*sin(q)) - sin(q)*sin(teta)*(td*cos(teta)*sin(q) - qd*sin(q)*sin(teta))]

that i want to rewrite it in [3*3]*[qd;td;ed]=[E] this form

is there any Command that can help me to calculate this [3*3] matrix by matlab

answer should be like this:

[                 -cos(q)*cos(teta)*sin(teta),                        - cos(teta)^2*sin(q) - sin(q)*sin(teta)^2, cos(q)*sin(teta)]
[                 -cos(teta)*sin(q)*sin(teta),                          sin(q)*cos(teta)^2 + cos(q)*sin(teta)^2, sin(q)*sin(teta)]
[ cos(q)^2*sin(teta)^2 + sin(q)^2*sin(teta)^2, cos(q)*cos(teta)*sin(q)*sin(teta) - cos(teta)*sin(q)^2*sin(teta),        cos(teta)]

i really apreciated if you could help me

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2017-12-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/372016-rewrite-a-3-1-matrix-to-3-3-3-1-var-a#answer_295448

在 MATLAB Online 中打开

Doing it by matrix division (the mldivide function) did not work because the rows are not consistent.

The only solution I can come up with is:

syms qd td ed q teta
E =[ ed*cos(q)*sin(teta) - cos(teta)*(qd*cos(q)*sin(teta) + td*cos(teta)*sin(q)) - td*sin(q)*sin(teta)^2 ; cos(teta)*(td*cos(teta)*sin(q) - qd*sin(q)*sin(teta)) + ed*sin(q)*sin(teta) + td*cos(q)*sin(teta)^2  ;  ed*cos(teta) + cos(q)*sin(teta)*(qd*cos(q)*sin(teta) + td*cos(teta)*sin(q)) - sin(q)*sin(teta)*(td*cos(teta)*sin(q) - qd*sin(q)*sin(teta))];
v = [qd;td;ed];
for k1 = 1:size(E,1)
    C{k1} = coeffs(E(k1), v);
end
Out = [C{1}; C{2}; C{3}]