这个其实就是以t为未知量，把两点坐标代入直线方程求解。
可以直接用solve求解，但得到的结果可能会包含t0自身（有很小的误差），需要将其过滤。
参考代码：
% 常数定义
a=2;b=1;t0=11/18*pi;
x0=a*cos(t0);y0=b*sin(t0);
kx=3;ky=-2;
% 解方程
syms t
t1=solve(ky*(a*cos(t)-x0)-kx*(b*sin(t)-y0));
% 对解进行过滤
t1=double(t1);
t1((a*cos(t1)-x0).^2+(b*sin(t1)-y0).^2
% 绘图
clf
ezplot(a*cos(t),b*sin(t))
hold on;
plot(x0,y0,'ro')
plot(a*cos(t1),b*sin(t1),'mp');
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024