开学已经是第二周了,我的《微分几何》也上课两周了,进度比较慢,现在才讲到平面曲线的曲率。在平面曲线$\boldsymbol{t}(t)=(x(t),y(t))$某点上可以找出单位切向量。

$$\boldsymbol{t}=\left(\frac{dx}{ds},\frac{dy}{ds}\right)$$

其中$ds^2 =dx^2+dy^2$,将这个向量逆时针旋转90度之后,就可以定义相应的单位法向量$\boldsymbol{n}$,即$\boldsymbol{t}\cdot\boldsymbol{n}=0$。

常规写法

让我们用弧长$s$作为参数来描述曲线方程,$\boldsymbol{t}(s)=(x(s),y(s))$,函数上的一点表示对$s$求导。那么我们来考虑$\dot{\boldsymbol{t}}$,由于$\boldsymbol{t}^2=1$,对s求导得到

$$\boldsymbol{t}\cdot\dot{\boldsymbol{t}}=0$$

也就是说$\dot{\boldsymbol{t}}$与$\boldsymbol{t}$垂直,由于只是在平面上,所以$\dot{\boldsymbol{t}}$与$\boldsymbol{n}$平行。即

$$\dot{\boldsymbol{t}}=\kappa \boldsymbol{n}$$

类似地,有$\dot{\boldsymbol{n}}$与$\boldsymbol{t}$平行。并且对$\boldsymbol{t}\cdot\boldsymbol{n}=0$求导得到

$$\dot{\boldsymbol{t}}\cdot\boldsymbol{n}+\boldsymbol{t}\cdot\dot{\boldsymbol{n}}=0$$

将$\dot{\boldsymbol{t}}=\kappa \boldsymbol{n}$代入上式得到

$$\dot{\boldsymbol{n}}=-\kappa \boldsymbol{t}$$

$\kappa$被称为曲线在该点的曲率。

复数表示

以上是教科书的标准写法,但事实上,研究平面曲线的最方便的工具还是复数。将$\boldsymbol{r}(s)$用一个带参数的复数表示$z(s)$,那么上面的两式可以写成更简洁的一个式子

$$\ddot{z}(s)=i\kappa (s) \dot{z}(s) $$

这样写的好处还在于,任意给出曲率函数$\kappa (s) $,我们就可以求出对应的曲线

$$z(s)=\int e^{i\int \kappa (s)ds}ds $$

这是简洁而有效的。

另外,不妨设$dz=ds e^{i\phi}$,那么

$$\dot{z}=e^{i\phi}$$

自然地

$$\ddot{z}=e^{i\phi}\left(i\dot{\phi}\right)$$

所以曲率可以表示为

$$\kappa=\dot{\phi}$$

各种坐标

利用它可以很方便地推导出各种坐标系下的曲率表达式,如曲线为一般的参数方程$(x(t),y(t))$时,用函数加一撇表示对t求导,有$ds=\sqrt{x'(t)^2+y'(t)^2}dt,\phi=\arctan\left(\frac{y'(t)}{x'(t)}\right)$,那么

$$\frac{d\phi}{ds}=\frac{\frac{y''(t)}{x'(t)}-\frac{y'(t)x''(t)}{[x'(t)]^2}}{1+\left(\frac{y'(t)}{x'(t)}\right)^2}\div \left(\frac{ds}{dt}\right)$$

代入整理易得

$$\kappa=\frac{y''(t) x'(t)-x''(t) y'(t)}{[x'(t)^2+y'(t)^2]^{3/2}}$$

在极坐标下,设$r=f(\theta)$,则$z=f(\theta)e^{i\theta}$,那么

$$dz=\left(\frac{d f}{d \theta}+i f\right)e^{i\theta}d\theta$$

所以

$$ds=\sqrt{f^2+\left(\frac{d f}{d \theta}\right)^2}d\theta$$

而$\phi=\arctan\frac{f}{\left(\frac{d f}{d \theta}\right)}+\theta$,那么

$$\frac{d\phi}{ds}=\left[\frac{1-\left(\frac{d^2 f}{d \theta^2}\right) f/\left(\frac{d f}{d \theta}\right)^2}{1+f^2/\left(\frac{d f}{d \theta}\right)^2}+1\right]\div \left(\frac{d s}{d \theta}\right)$$

代入整理得

$$\kappa=\frac{2\left(\frac{d f}{d \theta}\right)^2+f^2-\left(\frac{d^2 f}{d \theta^2}\right)f}{\left[\left(\frac{d f}{d \theta}\right)^2+f^2\right]^{3/2}}$$

三维空间有没有类似方便的东西呢?我也正在思考^_^

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

打赏

微信打赏

支付宝打赏

因为网站后台对打赏并无记录,因此欢迎在打赏时候备注留言。你还可以点击这里或在下方评论区留言来告知你的建议或需求。

如果您需要引用本文,请参考:

苏剑林. (2014, Mar 04). 《平面曲线的曲率的复数表示 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/2403

det曲线_平面曲线的曲率的复数表示相关推荐

  1. det曲线_如何画说话人识别/确认的DET曲线?

    谢邀了. 我画DET曲线也是用Matlab代码的,先给我用的工具链接: DETware_v2-1-tar.gz 这工具只需要运行Eval_Spkr_Det.m这个文件即可,里面有画DET曲线的部分,同 ...

  2. det曲线_一文说透机器学习的主流评价指标

    人工智能或者机器学习的各种文章里面各种指标,还有中英文的,例如FPR.误检率,更有甚者,各种曲线,如ROC.DET等,眼花缭乱,它们之间到底是啥关系? 如果读了这篇文章,懂了这些指标的含义和它们间的关 ...

  3. MATLAB平面曲线按曲率上色三维曲面颜色栏链接自定义属性(二维曲线第三维用颜色表示,三维曲面第四维用颜色表示)

    平面曲线按曲率上色(曲线第三维用颜色表示) 先上效果图 正弦曲线按曲率上色 自由曲线按误差上色 自由曲面按误差上色 曲线第三维用颜色表示实现方法: x=linspace(-10,10,100); y= ...

  4. 模式识别的评价方法:ROC曲线, DET曲线, FPPW, FPPI

    转载自:http://blog.csdn.net/pb09013037/article/details/48949037 因个人在做模式识别相关的工作,模式识别算法最终的性能评价是关键.但苦于网上很难 ...

  5. 模式识别分类器评价指标之DET曲线

    DET(Detection Error Tradeoff )曲线是对二元分类系统误码率的曲线图,绘制出错误拒绝率FRR(False Reject Rate)与错误接受率(False Accept Ra ...

  6. det曲线_11565 P-R、ROC、DET 曲线及 AP、AUC 指标全解析(上)

    机器学习中,最「简单」的一种任务就是二分类任务了.比如,说话人验证(speaker verification)就是一个二分类任务:判断一句话是否是给定说话人说的.再比如,信息检索(informatio ...

  7. 光滑曲线_微分几何笔记(4) —— 二维三维空间中曲线的曲率以及环绕数

     本篇文章我们从一般化的 空间回到我们生活的 空间,看看低维空间中的曲线有哪些性质,主要计算下在非弧长参数下的曲线,曲率挠率的一般表达式. 最后引入环绕数的概念,讲讲怎么数曲线转了多少圈. 4.1 ...

  8. 光滑曲线_曲线的曲率

    曲率(curvature)就是针对曲线上某个点的切线方向角对弧长的转动率, 是通过微分来定义的,表明曲线偏离直线的程度. 数学上表明曲线在某一点的弯曲程度的数值. 曲率越大,表示曲线的弯曲程度越大.曲 ...

  9. cad抛物线曲线lisp_曲线的转弯半径和曲率 - AutoLISP/Visual LISP 编程技术 - CAD论坛 - 明经CAD社区 - Powered by Discuz!...

    ;;;============================================================= ;;; 一般平面曲线参数方程的曲率离心公式 ;;; 功能: 获取曲线上 ...

最新文章

  1. 【C++】Google Protocol Buffer(protobuf)详解(一)
  2. 扩增子统计绘图7三元图
  3. Pycharm的项目文件名是红色的原因及解决办法
  4. [调试技巧] 如何利用WinDbg找出程序崩溃的位置
  5. windows2003+PHP-5.3.8+IIS6.0的安装配置
  6. C#关于读写INI文件
  7. (1)stm32开发之标准库的介绍
  8. python实现装饰器_python装饰器的实现
  9. 运维真的被云革命掉了吗?
  10. three.js 坐标系、camera位置属性、点、线、面
  11. 为何需要实施数据治理
  12. C 线性表的链式存储实现及插入、删除等操作示例
  13. 直接读取Excel文件数据
  14. 2018年视频云服务市场格局进入整合阶段,阿里云视频云位居市场竞争力领导者的位置... 1
  15. MySQL数据库regdate_第十五章 MySQL 数据库
  16. Python爬虫实战+数据分析+数据可视化(前程无忧招聘信息)
  17. java实现12306查票_java抓取12306信息实现火车余票查询示例
  18. html table最小宽度,table宽度比tbody多1
  19. 百度UEditor编译器中获取HTML内容和纯文本,设置UEditor编辑器的内容
  20. html盒子在始终在左侧,div盒子在一行方法(左中右结构 CSS布局)

热门文章

  1. java jpa 异步编程_异步处理时的JPA
  2. 访问nfs_通过NFS访问编年引擎
  3. 工程模式和抽象工厂模式_功能工厂模式
  4. slf4j绑定器_用于ADFLogger的SLF4J绑定–缺少的部分
  5. 使用Google Cloud Storage托管您的Maven工件
  6. activemq优先级_ActiveMQ消息优先级:工作原理
  7. ssh 看apache_使用Apache KeyedObjectPool的ssh连接池
  8. Eclipse MicroProfile:您需要了解的5件事
  9. java btrace_BTrace for Java应用程序简介
  10. 带有Javaslang的Java 8中的功能数据结构