一、正定矩阵与椭圆的关系

二维空间中椭圆最基本的形式为
x2a2+y2b2=1\frac{x^{2}}{a^{2}} +\frac{y^{2}}{b^{2}} =1a2x2​+b2y2​=1
上面的这个方程写成矩阵的形式为
[xy]T[1a2001b2][xy]=xTAx=1\begin{bmatrix} x\\ y \end{bmatrix}^{T} \begin{bmatrix} \frac{1}{a^{2}} & 0\\ 0& \frac{1}{b^{2}} \end{bmatrix} \begin{bmatrix} x\\ y \end{bmatrix} =x^{T}Ax=1[xy​]T[a21​0​0b21​​][xy​]=xTAx=1
AAA的特征值为
λ1=1a2,λ2=1b2\lambda _{1}=\frac{1}{a^{2}},\space \lambda _{2}=\frac{1}{b^{2}} λ1​=a21​, λ2​=b21​
AAA的归一化特征向量为
μ1=[120],μ2=[012]\mu _{1}=\begin{bmatrix} \frac{1}{\sqrt{2} }\\ 0 \end{bmatrix},\space \mu _{2}=\begin{bmatrix} 0\\ \frac{1}{\sqrt{2} } \end{bmatrix}μ1​=[2​1​0​], μ2​=[02​1​​]
椭圆的长短轴分别沿着矩阵AAA的两个特征向量的方向,而两个与之对应的特征值分别是半长轴和半短轴的长度的平方的倒数,即满足
a=1λ1,b=1λ2a=\frac{1}{\sqrt{\lambda_{1} }},\space b=\frac{1}{\sqrt{\lambda_{2} }}a=λ1​​1​, b=λ2​​1​
更一般的,对于如下方程所示椭圆
ax2+2bxy+cy2=1ax^{2}+2bxy+cy^{2}=1ax2+2bxy+cy2=1
化为矩阵形式是
[xy]T[abbc][xy]=xTAx=1\begin{bmatrix} x\\ y \end{bmatrix}^{T} \begin{bmatrix} a & b\\ b & c \end{bmatrix} \begin{bmatrix} x\\ y \end{bmatrix} =x^{T}Ax=1[xy​]T[ab​bc​][xy​]=xTAx=1

二、椭圆的一般方程

椭圆的一般方程为
[(x−x0)cos⁡θ+(y−y0)sin⁡θ]2a2+[(x−x0)sin⁡θ−(y−y0)cos⁡θ]2b2=1\frac{\left [\left (x-x_{0} \right ) \cos \theta +\left (y-y_{0} \right )\sin \theta \right ]^{2}}{a^{2}} +\frac{\left [\left (x-x_{0} \right )\sin \theta -\left (y-y_{0} \right )\cos \theta \right ]^{2}}{b^{2}} =1 a2[(x−x0​)cosθ+(y−y0​)sinθ]2​+b2[(x−x0​)sinθ−(y−y0​)cosθ]2​=1
其中 (x0,y0)(x_{0},y_{0})(x0​,y0​) 为椭圆中心点,θ为椭圆倾角,椭圆倾角如下图紫色线与黄色线所成锐角

将该椭圆变为矩阵形式为
[x−x0y−y0]T[cos⁡2θa2+sin⁡2θb2cos⁡θsin⁡θa2+cos⁡θsin⁡θb2cos⁡θsin⁡θa2+cos⁡θsin⁡θb2sin⁡2θa2+cos⁡2θb2][x−x0y−y0]=1\begin{bmatrix} x-x_{0}\\ y-y_{0} \end{bmatrix}^{T} \begin{bmatrix} \frac{\cos ^{2}\theta}{a^{2}}+\frac{\sin ^{2}\theta}{b^{2}} & \frac{\cos\theta\sin\theta}{a^{2}}+\frac{\cos\theta\sin\theta}{b^{2}}\\ \frac{\cos\theta\sin\theta}{a^{2}}+\frac{\cos\theta\sin\theta}{b^{2}} & \frac{\sin ^{2}\theta}{a^{2}}+\frac{\cos ^{2}\theta}{b^{2}} \end{bmatrix} \begin{bmatrix} x-x_{0}\\ y-y_{0} \end{bmatrix}=1 [x−x0​y−y0​​]T[a2cos2θ​+b2sin2θ​a2cosθsinθ​+b2cosθsinθ​​a2cosθsinθ​+b2cosθsinθ​a2sin2θ​+b2cos2θ​​][x−x0​y−y0​​]=1

三、绘制一般椭圆的matlab程序

ellipse_l = 5;  %椭圆半长轴
ellipse_s = 2;  %椭圆半短轴
ellipse_theta = (10) / 180 * pi;   %椭圆倾角
ellipse_center = [5, 4.5];  %椭圆中心点坐标%椭圆展开后方程的各个参数
a = ellipse_l^2*sin(ellipse_theta)^2+ellipse_s^2*cos(ellipse_theta)^2;
b = 2*(ellipse_s^2-ellipse_l^2)*sin(ellipse_theta)*cos(ellipse_theta);
c = ellipse_l^2*cos(ellipse_theta)^2+ellipse_s^2*sin(ellipse_theta)^2;
d = -ellipse_l^2*ellipse_s^2;%绘制椭圆图像
figure(1)
fimplicit(@(x,y)a.*(x-ellipse_center(1)).^2+b.*(x-ellipse_center(1)).*(y-ellipse_center(2))+c.*(y-ellipse_center(2)).^2+d,[0 10 0 10]);
hold on
scatter(ellipse_center(1),ellipse_center(2))
hold on

一般情况下的椭圆方程相关推荐

  1. 正常情况及匝间短路情况下,三相PMSM电机电压方程推导及对电感的理解

    匝间短路下(interturn short)的电压方程 最近在学习故障检测的知识(fault detection),fault detection可以分为三类方法:based on model/sig ...

  2. 利用指针访问opencv Mat类型的矩阵,以及求椭圆方程的函数

    1. 能利用数组初始化Mat类型矩阵吗? 2.为什么要用数组初始化Mat矩阵,以及什么情况下使用,使用数组初始化Mat矩阵的条件是什么? 3.怎么利用指针访问Mat矩阵,是怎么进行的? 答:是为了使用 ...

  3. 硕博生,在导师完全散养的情况下,如何提升自己的科研能力?(转载)

    链接:https://www.zhihu.com/question/385466539 来源:深度学习与计算机视觉 声明:仅做学术分享,侵删 看到一些文章写的不错,我对此进行转载和收藏,在自己学习的同 ...

  4. 【组合数学】递推方程 ( 通解定义 | 无重根下递推方程通解结构定理 )

    文章目录 一.通解定义 二.无重根下递推方程通解结构定理 一.通解定义 递推方程解的形式 : 满足 H(n)−a1H(n−1)−a2H(n−2)−⋯−akH(n−k)=0H(n) - a_1H(n-1 ...

  5. 标准椭圆和任意椭圆方程之间的变换公式推导

    我们在高中数学中就学习过标准的椭圆方程如下:x2a2+y2b2=1\frac{x^2}{a^2}+\frac{y^2}{b^2}=1a2x2​+b2y2​=1这个方程表示一个中心在坐标原点.长轴延xx ...

  6. unity 使用超椭圆方程对图片进行裁切制作圆角矩形

    今天无意间在网上发现了一个公式,叫做超椭圆方程,效果就是圆角矩形: 原文地址: http://blog.csdn.net/tom_221x/article/details/50924154 原文中是使 ...

  7. python解椭圆方程的例题_《椭圆》方程典型例题20例(含标准答案)

    X 1 + X 2 二 X M 1 一 2 1 a 2 2~ , a 1 1 a 2 <椭圆>方程典型例题 20 例 典型例题一 例 1 椭圆的一个顶点为 A 2,0 ,其长轴长是短轴长的 ...

  8. 工业机器人什么情况下会出现奇点_机器人奇点是什么?他们怎么会像黑洞一样?...

    如果你对科学感兴趣,"奇点"很可能会让你想起黑洞.自从美国LIGO实验室证明了引力波的存在之后,黑洞就被媒体广泛报道,为大众所熟知.根据物理学家推论,在黑洞的中心,存在一个&quo ...

  9. 93---Python 极坐标系下绘制椭圆

    阅读这篇之前,推荐阅读-90-Python 直角坐标系下绘制椭圆形(超链接点击跳转). 直角坐标系下椭圆方程的表达式为: x 2 a 2 + y 2 b 2 = r 2

最新文章

  1. UltraTextEditor
  2. RabbitMQ负载均衡(4)——LVS
  3. (2.15)备份与还原--使用作业备份、清理过期备份、清理历史记录、事务日志是否备份过...
  4. Flink的时间语义和Watermark
  5. jq发送动态变量_「系统架构」Nginx调优之变量的使用(3)
  6. msgpack在java和go中序列化的区别
  7. php开启ziparchive类,php压缩解压文件ZipArchive类的方法使用教程
  8. 计算机二级考试vb内容,计算机二级VB考试考什么内容?(一)
  9. 削峰填谷 matlab,风电与电动汽车协同并网调度环境模型
  10. Winfrom+Devexpress 之GridControl ConvertTo TileView 制作精美的卡片式看板
  11. 提莫攻击的中毒持续时间
  12. 远程连接腾讯云MySQL数据库
  13. 正则表达式检测电子邮件_电子邮件中的垃圾邮件检测
  14. 秒变配色高手!怎么都不会错的6条网页设计配色原则
  15. K-means最优K值计算(利用SSE)
  16. python列表过滤的方法
  17. 揭秘换肤技术:实现系统滚动条换肤功能
  18. Go语言if条件判断
  19. 西电保研至东南经验贴
  20. CentOS 6.5 本地yum源配置

热门文章

  1. 全球首家!苹果市值达 3 万亿美元,AR 和自动驾驶是下一重点?
  2. c语言stl用法,C STL快速入门!学习使我快乐
  3. Android7.0以上 安装Ca证书
  4. 怎么升级计算机配置,旧电脑配置升级攻略,看完再决定要不要配置新电脑
  5. notify()会立刻释放锁么?
  6. Google及Facebook第三方登录问题,“将你登入到此应用时出错,请稍后再试”
  7. 工业级静默活体检测开源算法技术解析记录
  8. 计算用户浏览时长(微信浏览器)
  9. 如何开启系统打印机服务器,windows7系统如何开启打印机服务
  10. 7人团队敏捷管理手册