用Python作一条已知曲线的等距曲线
创始人
2024-04-24 04:44:49
0

参考资料:

该如何作一条已知曲线的等距曲线? - 知乎

等距线_百度百科


目录

1.等距线

2.数学推导 

3.示例

4.代码与结果


 1.等距线

等距线(equidistant line)亦称平行曲线,一种平面曲线,即由一已知曲线所产生的另一曲线。平行曲线如下图。

2.数学推导

2.1 核心思想: 

从原曲线法线方向偏移所需的距离,就可以得到这些平行曲线

2.2数学推导

对于参数曲线r(t),其切线速度就是其第一导数:

v(t)=r'(t)

归一化后得单位切线矢量:

T(t)=(T_{x}(t),T_{y}(t))=\frac{v(t)}{\left | v(t) \right |}

在二维中,垂直于切线的单位法矢量只需旋转切线90度,可以得到垂直于切线的单位法矢量为

N(t)=(-T_{y}(t),T_{x}(t))

那么,设偏移距离为a,r(t)的两条平行曲线为

s(t)=r(t)\pm aN(t)

3.示例

本文以y=4x^{3}+3x^{2}+2x+1为例,将该方程转为参数方程,得s(t)为,

s(t)=\left\{\begin{matrix}t \\ 4t^{3}+3t^{2}+2t+1 \end{matrix}\right.

对s(t)求导得v(t),

v(t)=\left\{\begin{matrix}1 \\ 12t^{2}+6t+2 \end{matrix}\right.

归一化后得单位切线矢量T(t),

T(t)=\left\{\begin{matrix}\frac{1}{\sqrt{1+(12t^{2}+6t+2)^{2}}} \\ \frac{12t^{2}+6t+2}{\sqrt{1+(12t^{2}+6t+2)^{2}}} \end{matrix}\right.

将T(t)旋转90度得到N(t),

N(t)=\left\{\begin{matrix}-1*\frac{12t^{2}+6t+2}{\sqrt{1+(12t^{2}+6t+2)^{2}}} \\ \frac{1}{\sqrt{1+(12t^{2}+6t+2)^{2}}} \end{matrix}\right.

于是平移距离为a的曲线为,

S(t)=r(t)\pm aN(t)

4.代码与结果

 PS:注意横、纵坐标尺度不一致,想看更明显的效果可以试验一下y=x^2

import numpy as np
import matplotlib.pyplot as pltt = np.linspace(-5,5,40)
x1 = t
y1 = 4*t**3 + 3*t**2 + 2*t + 1x2 = t + 3*-1*(12*t**2+6*t+2)/np.sqrt(1+(12*t**2+6*t+2)**2)
y2 = 4*t**3 + 3*t**2 + 2*t + 1 + 3*1/np.sqrt(1+(12*t**2+6*t+2)**2)x3 = t - 3*-1*(12*t**2+6*t+2)/np.sqrt(1+(12*t**2+6*t+2)**2)
y3 = 4*t**3 + 3*t**2 + 2*t + 1 - 3*1/np.sqrt(1+(12*t**2+6*t+2)**2)plt.plot(x1,y1)
plt.plot(x2,y2)
plt.plot(x3,y3)

 

相关内容

热门资讯

监控摄像头接入GB28181平... 流程简介将监控摄像头的视频在网站和APP中直播,要解决的几个问题是:1&...
Windows10添加群晖磁盘... 在使用群晖NAS时,我们需要通过本地映射的方式把NAS映射成本地的一块磁盘使用。 通过...
protocol buffer... 目录 目录 什么是protocol buffer 1.protobuf 1.1安装  1.2使用...
在Word、WPS中插入AxM... 引言 我最近需要写一些文章,在排版时发现AxMath插入的公式竟然会导致行间距异常&#...
【PdgCntEditor】解... 一、问题背景 大部分的图书对应的PDF,目录中的页码并非PDF中直接索引的页码...
Fluent中创建监测点 1 概述某些仿真问题,需要创建监测点,用于获取空间定点的数据࿰...
educoder数据结构与算法...                                                   ...
MySQL下载和安装(Wind... 前言:刚换了一台电脑,里面所有东西都需要重新配置,习惯了所...
MFC文件操作  MFC提供了一个文件操作的基类CFile,这个类提供了一个没有缓存的二进制格式的磁盘...
有效的括号 一、题目 给定一个只包括 '(',')','{','}'...