[討論] 在特定條件下運算代數矩陣

看板MATLAB作者 (WANG)時間10年前 (2015/04/07 01:38), 10年前編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
各位大大好 小弟做代數矩陣時 想要將 kx=sym('kx') ky=sym('ky') kz=sym('kz') 這三個代數 賦予kx^2+ky^2+kz^2=1的性質 要怎麼在程式中加入此項條件呢? 有高手可幫忙嗎? 感激不盡 我打的程式碼如下 ============================================================================== kx=sym('kx') ky=sym('ky') kz=sym('kz') sida=sym('sida') cosa=kx/((kx^2+ky^2)^0.5) sina=ky/((kx^2+ky^2)^0.5) cosb=kz/((kx^2 + ky^2 + kz^2)^(0.5)) sinb=(kx^2+ky^2)^(0.5)/((kx^2 + ky^2 + kz^2)^(0.5)) coss=cos(sida) sins=sin(sida) Rza=[ cosa sina 0; -sina cosa 0; 0 0 1] Ryb=[ cosb 0 -sinb; 0 1 0; sinb 0 cosb] Rzs=[coss -sins 0; sins coss 0; 0 0 1] R=inv(Rza)*inv(Ryb)*Rzs*Ryb*Rza simplify(R) -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 125.230.106.226 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1428341911.A.92F.html ※ 編輯: achoodrumer (125.230.106.226), 04/07/2015 01:39:45
文章代碼(AID): #1L8iINal (MATLAB)
文章代碼(AID): #1L8iINal (MATLAB)