Определите и постройте параметрическую поверхность

Задайте параметрическую поверхность x(u,v), y(u,v), z(u,v) следующим образом.

syms u v
x = cos(u)*sin(v);
y = sin(u)*sin(v);
z = cos(v)*sin(v);

Постройте график поверхности с помощью fsurf.

fsurf(x,y,z)
axis equal

Создайте матрицы вращения

Создайте матрицы 3 на 3 Rx, Ry, и Rz представление поворотов плоскости на угол t о x-, y-, и z-ось, соответственно.

syms t

Rx = [1 0 0; 0 cos(t) -sin(t); 0 sin(t) cos(t)]

Rx = 
$$\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(t\right)& -\sin \left(t\right)\\ 0& \sin \left(t\right)& \cos \left(t\right)\end{array}\right)$$

Ry = [cos(t) 0 sin(t); 0 1 0; -sin(t) 0 cos(t)]

Ry = 
$$\left(\begin{array}{ccc}\cos \left(t\right)& 0& \sin \left(t\right)\\ 0& 1& 0\\ -\sin \left(t\right)& 0& \cos \left(t\right)\end{array}\right)$$

Rz = [cos(t) -sin(t) 0; sin(t) cos(t) 0; 0 0 1]

Rz = 
$$\left(\begin{array}{ccc}\cos \left(t\right)& -\sin \left(t\right)& 0\\ \sin \left(t\right)& \cos \left(t\right)& 0\\ 0& 0& 1\end{array}\right)$$

Вращайте вокруг каждой оси в трёх размерностях

Сначала поверните поверхность вокруг x-ось на 45 степени против часовой стрелки.

xyzRx = Rx*[x;y;z];
Rx45 = subs(xyzRx, t, pi/4);

fsurf(Rx45(1), Rx45(2), Rx45(3))
title('Rotating by \pi/4 about x, counterclockwise')
axis equal

Вращайте вокруг z-ось на 90 степени по часовой стрелке.

xyzRz = Rz*Rx45;
Rx45Rz90 = subs(xyzRz, t, -pi/2);

fsurf(Rx45Rz90(1), Rx45Rz90(2), Rx45Rz90(3))
title('Rotating by \pi/2 about z, clockwise')
axis equal

Вращайте вокруг y-ось на 45 степени по часовой стрелке.

xyzRy = Ry*Rx45Rz90;
Rx45Rz90Ry45 = subs(xyzRy, t, -pi/4);

fsurf(Rx45Rz90Ry45(1), Rx45Rz90Ry45(2), Rx45Rz90Ry45(3))
title('Rotating by \pi/4 about y, clockwise')
axis equal

Масштабирование и вращение

Масштабируйте поверхность на множитель 3 вдоль z-ось. Можно умножить выражение для z по 3, z = 3*z. Более общий подход состоит в том, чтобы создать масштабирующую матрицу, а затем умножить масштабирующую матрицу на вектор координат.

S = [1 0 0; 0 1 0; 0 0 3];
xyzScaled = S*[x; y; z]

xyzScaled = 
$$\left(\begin{array}{c}\cos \left(u\right)\hspace{0.17em}\sin \left(v\right)\\ \sin \left(u\right)\hspace{0.17em}\sin \left(v\right)\\ 3\hspace{0.17em}\cos \left(v\right)\hspace{0.17em}\sin \left(v\right)\end{array}\right)$$

fsurf(xyzScaled(1), xyzScaled(2), xyzScaled(3))
title('Scaling by 3 along z')
axis equal

Поверните масштабированную поверхность вокруг x-, y-, и z-ось на 45 степени по часовой стрелке, по порядку z, затем y, затем x. Матрица вращения для этого преобразования следующая.

R = Rx*Ry*Rz

R = 
$$\begin{array}{l}\left(\begin{array}{ccc}{\cos \left(t\right)}^2& -\cos \left(t\right)\hspace{0.17em}\sin \left(t\right)& \sin \left(t\right)\\ {\sigma }_1& {\cos \left(t\right)}^2-{\sin \left(t\right)}^3& -\cos \left(t\right)\hspace{0.17em}\sin \left(t\right)\\ {\sin \left(t\right)}^2-{\cos \left(t\right)}^2\hspace{0.17em}\sin \left(t\right)& {\sigma }_1& {\cos \left(t\right)}^2\end{array}\right)\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\sigma }_1=\cos \left(t\right)\hspace{0.17em}{\sin \left(t\right)}^2+\cos \left(t\right)\hspace{0.17em}\sin \left(t\right)\end{array}$$

Используйте матрицу поворота, чтобы найти новые координаты.

xyzScaledRotated = R*xyzScaled;
xyzSR45 = subs(xyzScaledRotated, t, -pi/4);

Постройте график поверхности.

fsurf(xyzSR45(1), xyzSR45(2), xyzSR45(3))
title('Rotating by \pi/4 about x, y, and z, clockwise')
axis equal

Проверяйте свойства матрицы вращения R

Матрицы вращения являются ортогональными матрицами. Таким образом, транспонирование R является также его обратной и определяющим R равен 1.

simplify(R.'*R)

ans = 
$$\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right)$$

simplify(det(R))

ans = $$1$$