各位大大好
小弟做代數矩陣時 想要將
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)