python线性插值

2024-03-20 02:28
文章标签 python 线性插值

本文主要是介绍python线性插值,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

假设有一个一维数组,但是此数组中只有部分位置上有值,其它位置数据缺失,现在想用线性插值的方法将其填充。

示例代码:

import numpy as np# 假设你有一个长度为171的数组,名为full_data,其中有13个数据点
# 用 None 表示缺失的数据点
full_data = [None] * 171# 假设你已经有了13个数据点的索引和值
known_indices = [11, 21, 31, 41, 51, 61, 71, 81, 91, 101, 121, 141, 161]
known_values = [0.207, 0.217, 0.184, 0.177, 0.162, 0.312,0.253, 0.317, 0.333, 0.336, 0.352, 0.288, 0.295]# 使用 NumPy 的interp函数进行线性插值
# 注意:interp函数要求已知数据点的索引必须是单调递增的
# 这里假设已知索引是单调递增的
full_indices = np.arange(len(full_data))
full_data = np.interp(full_indices, known_indices, known_values)# 打印结果
print(full_data)

结果展示:

[0.207   0.207   0.207   0.207   0.207   0.207   0.207   0.207   0.2070.207   0.207   0.207   0.208   0.209   0.21    0.211   0.212   0.2130.214   0.215   0.216   0.217   0.2137  0.2104  0.2071  0.2038  0.20050.1972  0.1939  0.1906  0.1873  0.184   0.1833  0.1826  0.1819  0.18120.1805  0.1798  0.1791  0.1784  0.1777  0.177   0.1755  0.174   0.17250.171   0.1695  0.168   0.1665  0.165   0.1635  0.162   0.177   0.1920.207   0.222   0.237   0.252   0.267   0.282   0.297   0.312   0.30610.3002  0.2943  0.2884  0.2825  0.2766  0.2707  0.2648  0.2589  0.2530.2594  0.2658  0.2722  0.2786  0.285   0.2914  0.2978  0.3042  0.31060.317   0.3186  0.3202  0.3218  0.3234  0.325   0.3266  0.3282  0.32980.3314  0.333   0.3333  0.3336  0.3339  0.3342  0.3345  0.3348  0.33510.3354  0.3357  0.336   0.3368  0.3376  0.3384  0.3392  0.34    0.34080.3416  0.3424  0.3432  0.344   0.3448  0.3456  0.3464  0.3472  0.3480.3488  0.3496  0.3504  0.3512  0.352   0.3488  0.3456  0.3424  0.33920.336   0.3328  0.3296  0.3264  0.3232  0.32    0.3168  0.3136  0.31040.3072  0.304   0.3008  0.2976  0.2944  0.2912  0.288   0.28835 0.28870.28905 0.2894  0.28975 0.2901  0.29045 0.2908  0.29115 0.2915  0.291850.2922  0.29255 0.2929  0.29325 0.2936  0.29395 0.2943  0.29465 0.2950.295   0.295   0.295   0.295   0.295   0.295   0.295   0.295   0.295  ]

线性插值绘图代码:

import matplotlib.pyplot as plt
import numpy as np# 数据
x = [11, 21, 31, 41, 51, 61, 71, 81, 91, 101, 121, 141, 161]
y = [0.207, 0.217, 0.184, 0.177, 0.162, 0.312,0.253, 0.317, 0.333, 0.336, 0.352, 0.288, 0.295]
x1 = np.arange(0, 171)
y1 = [0.207, 0.207, 0.207, 0.207, 0.207, 0.207, 0.207, 0.207, 0.207,0.207, 0.207, 0.207, 0.208, 0.209, 0.21, 0.211, 0.212, 0.213,0.214, 0.215, 0.216, 0.217, 0.2137, 0.2104, 0.2071, 0.2038, 0.2005,0.1972, 0.1939, 0.1906, 0.1873, 0.184, 0.1833, 0.1826, 0.1819, 0.1812,0.1805, 0.1798, 0.1791, 0.1784, 0.1777, 0.177, 0.1755, 0.174, 0.1725,0.171, 0.1695, 0.168, 0.1665, 0.165, 0.1635, 0.162, 0.177, 0.192,0.207, 0.222, 0.237, 0.252, 0.267, 0.282, 0.297, 0.312, 0.3061,0.3002, 0.2943, 0.2884, 0.2825, 0.2766, 0.2707, 0.2648, 0.2589, 0.253,0.2594, 0.2658, 0.2722, 0.2786, 0.285, 0.2914, 0.2978, 0.3042, 0.3106,0.317, 0.3186, 0.3202, 0.3218, 0.3234, 0.325, 0.3266, 0.3282, 0.3298,0.3314, 0.333, 0.3333, 0.3336, 0.3339, 0.3342, 0.3345, 0.3348, 0.3351,0.3354, 0.3357, 0.336, 0.3368, 0.3376, 0.3384, 0.3392, 0.34, 0.3408,0.3416, 0.3424, 0.3432, 0.344, 0.3448, 0.3456, 0.3464, 0.3472, 0.348,0.3488, 0.3496, 0.3504, 0.3512, 0.352, 0.3488, 0.3456, 0.3424, 0.3392,0.336, 0.3328, 0.3296, 0.3264, 0.3232, 0.32, 0.3168, 0.3136, 0.3104,0.3072, 0.304, 0.3008, 0.2976, 0.2944, 0.2912, 0.288, 0.28835, 0.2887,0.28905, 0.2894, 0.28975, 0.2901, 0.29045, 0.2908, 0.29115, 0.2915, 0.29185,0.2922, 0.29255, 0.2929, 0.29325, 0.2936, 0.29395, 0.2943, 0.29465, 0.295,0.295, 0.295, 0.295, 0.295, 0.295, 0.295, 0.295, 0.295, 0.295]# 绘制折线图
plt.scatter(x, y, color='red', label='pre-interpolation')
plt.plot(x, y, color='pink')
# dashes设置虚线与虚线之间的间隔
plt.plot(x1, y1, alpha=1, linestyle='--',color='blue', dashes=(10, 10), label='post-interpolation')# 添加标题和标签
plt.xlabel('X-axis')
plt.ylabel('Y-axis')plt.legend()# 显示图形
plt.show()

结果展示:
在这里插入图片描述
从图片可以看出,线性插值即是将两点连成直线,插值的数据即是直线上的数据。

这篇关于python线性插值的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/828001

相关文章

Python版本信息获取方法详解与实战

《Python版本信息获取方法详解与实战》在Python开发中,获取Python版本号是调试、兼容性检查和版本控制的重要基础操作,本文详细介绍了如何使用sys和platform模块获取Python的主... 目录1. python版本号获取基础2. 使用sys模块获取版本信息2.1 sys模块概述2.1.1

一文详解Python如何开发游戏

《一文详解Python如何开发游戏》Python是一种非常流行的编程语言,也可以用来开发游戏模组,:本文主要介绍Python如何开发游戏的相关资料,文中通过代码介绍的非常详细,需要的朋友可以参考下... 目录一、python简介二、Python 开发 2D 游戏的优劣势优势缺点三、Python 开发 3D

Python函数作用域与闭包举例深度解析

《Python函数作用域与闭包举例深度解析》Python函数的作用域规则和闭包是编程中的关键概念,它们决定了变量的访问和生命周期,:本文主要介绍Python函数作用域与闭包的相关资料,文中通过代码... 目录1. 基础作用域访问示例1:访问全局变量示例2:访问外层函数变量2. 闭包基础示例3:简单闭包示例4

Python实现字典转字符串的五种方法

《Python实现字典转字符串的五种方法》本文介绍了在Python中如何将字典数据结构转换为字符串格式的多种方法,首先可以通过内置的str()函数进行简单转换;其次利用ison.dumps()函数能够... 目录1、使用json模块的dumps方法:2、使用str方法:3、使用循环和字符串拼接:4、使用字符

Python版本与package版本兼容性检查方法总结

《Python版本与package版本兼容性检查方法总结》:本文主要介绍Python版本与package版本兼容性检查方法的相关资料,文中提供四种检查方法,分别是pip查询、conda管理、PyP... 目录引言为什么会出现兼容性问题方法一:用 pip 官方命令查询可用版本方法二:conda 管理包环境方法

基于Python开发Windows自动更新控制工具

《基于Python开发Windows自动更新控制工具》在当今数字化时代,操作系统更新已成为计算机维护的重要组成部分,本文介绍一款基于Python和PyQt5的Windows自动更新控制工具,有需要的可... 目录设计原理与技术实现系统架构概述数学建模工具界面完整代码实现技术深度分析多层级控制理论服务层控制注

pycharm跑python项目易出错的问题总结

《pycharm跑python项目易出错的问题总结》:本文主要介绍pycharm跑python项目易出错问题的相关资料,当你在PyCharm中运行Python程序时遇到报错,可以按照以下步骤进行排... 1. 一定不要在pycharm终端里面创建环境安装别人的项目子模块等,有可能出现的问题就是你不报错都安装

Python打包成exe常用的四种方法小结

《Python打包成exe常用的四种方法小结》本文主要介绍了Python打包成exe常用的四种方法,包括PyInstaller、cx_Freeze、Py2exe、Nuitka,文中通过示例代码介绍的非... 目录一.PyInstaller11.安装:2. PyInstaller常用参数下面是pyinstal

Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题

《Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题》在爬虫工程里,“HTTPS”是绕不开的话题,HTTPS为传输加密提供保护,同时也给爬虫带来证书校验、... 目录一、核心问题与优先级检查(先问三件事)二、基础示例:requests 与证书处理三、高并发选型:

Python中isinstance()函数原理解释及详细用法示例

《Python中isinstance()函数原理解释及详细用法示例》isinstance()是Python内置的一个非常有用的函数,用于检查一个对象是否属于指定的类型或类型元组中的某一个类型,它是Py... 目录python中isinstance()函数原理解释及详细用法指南一、isinstance()函数