matlab生成不同的椭圆,[平面几何][Matlab] 平面椭圆参数与一般式之间的转换

最新推荐文章于 2022-08-30 21:27:03 发布

aniKuso

最新推荐文章于 2022-08-30 21:27:03 发布

阅读量440

点赞数

文章标签： matlab生成不同的椭圆

椭圆的一般式为：\[A{x^2} + Bxy + C{y^2} + Dx + Ey + F = 0\]

椭圆的参数为：长半轴 $a$短半轴 $b$ 椭圆中心 $(x_{0},y_{0})$ 倾角为 $\theta$ (定义逆时针为正，长轴与x正方向的夹角)

1.由椭圆参数转化为一般式：

推导过程：

椭圆 $C[3*3]$ ,中心在原点，长轴与x轴重合，经过旋转矩阵 ${R} =f({\theta})$ , 平移矩阵 ${T}$$ =$$g$$(x_{0},y_{0})$,

后得到

$C^{'}=T^{T}*R^{T}*C*R*T$

ps:关于旋转矩阵R和平移矩阵T的定义看上篇博文直角坐标系下点/曲线平移与旋转的矩阵计算

即 $H({\theta},x_{0},y_{0},a,b)=A{x^2} + Bxy + C{y^2} + Dx + Ey + F $

对应相等可以得到：

①$A = \frac{{{{\cos }^2}\theta }}{{{a^2}}} + \frac{{{{\sin }^2}\theta }}{{{b^2}}}$

② $B = 2 \cdot \sin \theta \cdot \cos \theta \cdot (\frac{1}{{{a^2}}} - \frac{1}{{{b^2}}})$

③ $C = \frac{{{{\cos }^2}\theta }}{{{b^2}}} + \frac{{{{\sin }^2}\theta }}{{{a^2}}}$

④ $D = - 2 \cdot [{x_0} \cdot (\frac{{{{\cos }^2}\theta }}{{{a^2}}} + \frac{{{{\sin }^2}\theta }}{{{b^2}}}) + {y_0} \cdot \sin \theta \cdot \cos \theta \cdot (\frac{1}{{{a^2}}} - \frac{1}{{{b^2}}})]$

⑤ $E = - 2 \cdot [{x_0} \cdot \sin \theta \cdot \cos \theta \cdot (\frac{1}{{{a^2}}} - \frac{1}{{{b^2}}}) + {y_0} \cdot (\frac{{{{\cos }^2}\theta }}{{{b^2}}} + \frac{{{{\sin }^2}\theta }}{{{a^2}}})]$

⑥$F = \frac{{{{({x_0} \cdot \cos \theta + {y_0} \cdot \sin \theta )}^2}}}{{{a^2}}} + \frac{{{{({x_0} \cdot \sin \theta - {y_0} \cdot \cos \theta )}^2}}}{{{b^2}}} - 1$

matlab推导过程

+验证

clc

syms a b theta x0 y0

% 公式推导

C = [1/a.^2 0 0;

0 1/b.^2 0;

0 0 -1;];

Rot = [cos(theta) sin(theta) 0;

-sin(theta) cos(theta) 0;

0 0 1;];

T = [1 0 -x0;