matlab 椭圆参数传递,[平面几何][Matlab] 平面椭圆参数与一般式之间的转换

最新推荐文章于 2023-12-25 08:57:58 发布

weixin_39671621

最新推荐文章于 2023-12-25 08:57:58 发布

阅读量356

点赞数

文章标签： matlab 椭圆参数传递

该博客介绍了如何将椭圆的参数形式转换为一般形式，并给出了MATLAB代码进行验证。首先，通过旋转和平移矩阵，推导出椭圆一般式的系数表达式。然后，从一般式反推出椭圆的长半轴、短半轴和倾角。最后，提供了MATLAB代码来计算这些参数，并绘制椭圆图形。

摘要由CSDN通过智能技术生成

椭圆的一般式为：\[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;