matlab曲线 投影,如何将空间曲线投影到三个坐标平面上？

clc;clear;

syms u1 t1 v1 bt2 alf th1 th2

%syms th1

dt=pi/180;

v1=3*dt;

i=16/15;

bet1=21.5*dt;

u1=70;

for alf=2.5*dt;

th1=0:pi/60:2*pi;

st=-th1;

t1=atan(-i*sin(alf)*sin(st)/(i*(sin(alf)*cos(st)*sin(bet1)+cos(alf)*cos(bet1))-cos(bet1)));

x=-u1.*sin(v1).*(cos(t1).*((4190228212462809*sin((16*th1)/15).*sin(alf))/4503599627370496+(1650574788106375.*cos((16*th1)/15).*sin(th1))/4503599627370496-(1650574788106375*sin((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496)-sin(t1).*(cos((16*th1)/15).*cos(th1)+sin((16*th1)/15).*cos(alf).*sin(th1)))-u1.*cos(v1).*((1650574788106375.*sin((16*th1)/15).*sin(alf))/4503599627370496-(4190228212462809*cos((16*th1)/15).*sin(th1))/4503599627370496+(4190228212462809*sin((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496);

y=u1.*sin(v1).*(cos(t1).*((1650574788106375*sin((16*th1)/15).*sin(th1))/4503599627370496-(4190228212462809*cos((16*th1)/15).*sin(alf))/4503599627370496+(1650574788106375*cos((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496)-sin(t1).*(sin((16*th1)/15).*cos(th1)-cos((16*th1)/15).*cos(alf).*sin(th1)))-u1.*cos(v1).*((1650574788106375*cos((16*th1)/15).*sin(alf))/4503599627370496+(4190228212462809*sin((16*th1)/15).*sin(th1))/4503599627370496&#