用matlab求直线和椭圆的交点坐标!>> syms a b t0 t kx ky>> %已知椭圆上的一点(a*cos(t0),b*sin(t0)),过该点的矢量的方向为(kx,ky),求该向量与椭圆的另一个交点坐标(a*cos(t),b*sin(t)),即求出旋转
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/26 09:18:06
用matlab求直线和椭圆的交点坐标!>> syms a b t0 t kx ky>> %已知椭圆上的一点(a*cos(t0),b*sin(t0)),过该点的矢量的方向为(kx,ky),求该向量与椭圆的另一个交点坐标(a*cos(t),b*sin(t)),即求出旋转
用matlab求直线和椭圆的交点坐标!
>> syms a b t0 t kx ky
>> %已知椭圆上的一点(a*cos(t0),b*sin(t0)),过该点的矢量的方向为(kx,ky),求该向量与椭圆的另一个交点坐标(a*cos(t),b*sin(t)),即求出旋转角t;
>> %显然t与t0,a,b,kx,ky有关(t的定义见图示)
用matlab求直线和椭圆的交点坐标!>> syms a b t0 t kx ky>> %已知椭圆上的一点(a*cos(t0),b*sin(t0)),过该点的矢量的方向为(kx,ky),求该向量与椭圆的另一个交点坐标(a*cos(t),b*sin(t)),即求出旋转
这个其实就是以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<eps)=[]
% 绘图
clf
ezplot(a*cos(t),b*sin(t))
hold on;
plot(x0,y0,'ro')
plot(a*cos(t1),b*sin(t1),'mp');