python离散点切向量的计算与可视化
问题
给定离散点如下
t = np.arange(0,2,0.1)
y = t**2-2*t+1
如图可见是一个参数方程为f(t)=t2−2t+1f(t) = t^2-2t+1f(t)=t2−2t+1的二次函数,计算曲线在这些离散点上的切向量。
计算方法
方法一
对曲线上某点Pi(xi,yi)P_{i}(x_{i},y_{i})Pi(xi,yi),我们将Pi+1−Pi−12\frac{P_{i+1}-P_{i-1}}{2}2Pi+1−Pi−1认为是曲线在点PiP_{i}Pi上的切向量。
u = np.gradient(t)
v = np.gradient(y)
fig = plt.figure()
plt.scatter(t,y)
plt.quiver(t,y,u,v)
plt.show()
方法二
采用三次样条曲线对离散点进行拟合,再求其一阶导,作为切向量。
from scipy.interpolate import splprep, splev
tck,myu = splprep([t,y],s=0)
u,v = splev(myu,tck,der=1)
fig = plt.figure()
plt.scatter(t, y)
plt.quiver(t,y,u,v)
plt.show()
python离散点切向量的计算与可视化相关推荐
- python单自由度振动计算-数据可视化2
参考 单自由度振动理论推导参考: 结果 输出位移由3部分组成: 代码 import matplotlib.pyplot as plt import numpy as np from moviepy.v ...
- Python语言程序设计基础科学计算与可视化小练习
提示:这里可以添加系列文章的所有文章的目录,目录需要自己手动添加 例如:第一章 Python 机器学习入门之pandas的使用 提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档 一些有 ...
- python三维矩阵可视化_科学计算三维可视化---Mlab基础(基于Numpy数组的绘图函数)...
Mlab了解 Mlab是Mayavi提供的面向脚本的api,他可以实现快速的三维可视化,Mayavi可以通过Mlab的绘图函数对Numpy数组建立可视化. 过程为: .建立数据源 .使用Filter( ...
- python计算商品总价_GitHub - ideaOzy/data_analysis: 基于Python的南京二手房数据采集及可视化分析...
基于Python的南京二手房数据采集及可视化分析 1 内容简介 首先通过爬虫采集链家网上所有南京二手房的房源数据,并对采集到的数据进行清洗:然后,对清洗后的数据进行可视化分析,探索隐藏在大量数据背后的 ...
- Python计算并可视化商品批发时商家收益和顾客节省情况
问题描述: 假设某商品进价49元,售价75元,现在商场搞活动,顾客每多买一个就给优惠1%.现在要求计算并可视化顾客购买num个该物品时商家收益和顾客节省情况的折线图.并标记商家收益最大的批发数量和收益 ...
- python 科学计算三维可视化笔记(第三周 高级进阶)
python 科学计算三维可视化笔记 第三周 高级进阶 python 科学计算三维可视化笔记 第三周 高级进阶 一.Mayavi 入门 1. Mayavi 库的基本元素 2. 快速绘制实例 3. Ma ...
- 【数据可视化应用】Python反距离权重(IDW)插值计算及可视化绘制
本文我们将介绍IDW(反距离加权法(Inverse Distance Weighted)) 插值的Python计算方法及插值结果的可视化绘制过程.主要涉及的知识点如下: IDW简介 自定义Python ...
- python 克里金空间插值_Python克里金(Kriging)插值计算及可视化绘制
前面两篇推文我们分别介绍了使用Python和R进行IDW(反距离加权法) 插值的计算及结果的可视化过程,详细内容可见如下: 本期推文,我们将介绍如何使用Python进行克里金(Kriging)插值计算 ...
- python 节气计算_Python可视化-二十四节气与生日间隔天数统计
作为日常可视化作品之一,今天分析一下自97年开始,60年内每次生日与春分.雨水两个节气的相隔天数,看看哪一年会有重合的运气. 立春,十二节令之首,旧称正月节:雨水,十二中气立始.立春后,继之雨水,东风 ...
- AI算法工程师 | 03人工智能基础-Python科学计算和可视化(三)Pandas
文章目录 Python 之 数据处理分析模块 Pandas 一.Pandas 开发环境搭建 二.Pandas 数据类型 1. Series 对象创建 2. DataFrame 对象创建 3. 获取 S ...
最新文章
- PLM在企业中的实际价值与意义
- ASP.NET MVC 控制器激活(二)
- 每日一皮:据说程序员都喜欢这些东西...
- java中sql去除游标_java.sql.SQLException:-ORA-01000:已超过最大打开游标
- XML内部DTD约束 Day24
- windows cmd编辑文本
- Linux部署之批量自动安装系统之DHCP篇
- 【PHP】百度贴吧自动/云签到脚本(PC 端)
- svn执行reflash/cleanup报错wc.db解决办法
- COM 组件创建实例失败,原因是出现以下错误: 8001010a解决办法
- Android JSON解析库Gson和Fast-json的使用对比和图书列表小案例
- facade 门面模式和mediator 调停者
- 上万元游戏拼多多7块搞定 PICO防不住
- MiL.k x Bithumb x Yanolja宣布建立三方商务及市场营销合作关系
- linux风扇转速,ubuntu系统调节GPU风扇转速
- 串行、并行、并发,别再傻傻分不清了!
- Java开发人员必须掌握的两个Linux魔法工具(四)
- 第三天:ptyhon基础知识
- python:拼多多订单接口api
- 于众目睽睽之下隐藏图像:深度隐写术