各位大大好
小弟做代数矩阵时 想要将
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)