b-spline python_SciPyTutorial-一元B样条插值
12. Scipy Tutorial-一元样条插值
前边使用过线性插值linear、最邻近插值nearest以及三次插值cubic,这些插值凸点点多,不够光滑,本章研究样条插值spline,实现尽量让插值函数的曲线显得更光滑(收敛性)一些。样条插值是使用一种名为样条的特殊分段多项式进行插值的形式。由于样条插值可以使用低阶多项式样条实现较小的插值误差,这样就避免了使用高阶多项式所出现的龙格现象,所以样条插值得到了流行。
scipy.interpolate包里提供了两个函数splev和splrep共同完成(B-样条:贝兹曲线(又称贝塞尔曲线))插值,和之前一元插值一步就能完成不同,样条插值需要两步完成,第一步先用splrep计算出b样条曲线的参数tck,第二步在第一步的基础上用splev计算出各取样点的插值结果。
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import splev, splrep
def f(x):
return x ** 2 + 10 * np.sin(x) + 1
x = np.linspace(0, 10, 20)
y = f(x)
plt.plot(x, y, '*-')
plt.show()
x2 = np.linspace(0, 10, 200)
spl = splrep(x, y, s = 0)
y2 = splev(x2, spl)
plt.plot(x, y, 'o', x2, y2)
plt.show()
b-spline python_SciPyTutorial-一元B样条插值相关推荐
- matlab spline边界约束,三阶样条插值(一阶导数边界条件) matlab程序
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function S=Threch1(X,Y,dy0,dyn,xi) % X 为已知数据的横坐标 %Y 为已知数据的纵坐标 %xi 插值点处的横坐标 %S ...
- matlab指令vpa(j10),MATLAB指令集(最新整理)
<MATLAB指令集(最新整理)>由会员分享,可在线阅读,更多相关<MATLAB指令集(最新整理)(6页珍藏版)>请在人人文库网上搜索. 1.主题词意义主题词意义主题词意义主题 ...
- matlab练习_MATLAB教程-台大郭彦甫-第十四节,含练习答案
14-回归与内插 一.Polynomial curve fitting(多项式曲线拟合) (一)Simple Linear Regression(简单线性回归) 1.A bunch of data p ...
- 秦九韶算法matlab实验报告,数值分析上机实验报告.doc
实验报告一 题目: (绪论) 非线性方程求解及误差估计 摘要:非线性方程的解析解通常很难给出,因此线性方程的数值解法就尤为重要.本实验采用两种常见的求解方法二分法.Newton法和改进的Newton法 ...
- matlab已知xyz拟合标定map,matlab练习题
<matlab练习题>由会员分享,可在线阅读,更多相关<matlab练习题(16页珍藏版)>请在人人文库网上搜索. 1.1. 在命令窗口输入语句:a=1,2,3;4,5,6;7 ...
- 亮温模型--相关概念
亮温模型–相关概念 本博文是亮温模型(Brightness Temperature Model)学习的第1篇论文,主要记录亮温模型学习的相关内容. 目录 亮温模型--相关概念 1.亮温模型相关 ...
- MATLAB数据分析(插值运算和曲线拟合)
MATLAB数据分析(插值运算和曲线拟合) MATLAB数据分析 在数学研究和工程应用中,我们常常需要对某些环节进行反复的测验或测试,以获取更多的数据支持最优决策,但是由于环境条件的限制,我们往往只能 ...
- 【美赛学习记录】模型
美赛学习记录-2022年2月7日 代码! 线性回归 数据拟合 插值 最优化求极值 插值 ARIMA 复杂网络实验 模型验证 K-Fold Cross-validation k折交叉验证 [基础模型] ...
- Matlab数据分析与计算,进程线程面试题总结
③ mean(A,dim):当dim为1时,该函数等同于mean(A):当dim为2时,返回一个列向量,其第i个元素是A的第i行的算术平均值. 求中值的函数是median,其调用方法和mean函数相同 ...
- 关于MATLAB的学习记录(纯入门用)
其实大二时有修过MATLAB这门课,但最终还是流于应付考试的表面. 今年终于决定试试美赛,于是乎又捡起来这门课,学习内容是按照bilibili上面中南大学的视频讲解,再加上自己编写代码和先前学习的印象 ...
最新文章
- 人工神经外网络中为什么ReLu要好过于Tanh和Sigmoid Function?
- Spring Boot 2.6.1 发布:为 Spring Cloud 2021铺路!
- LeetCode 17. Letter Combinations of a Phone Number
- 简单说说你对Java内存模型的原子性的了解吧?
- CyclicBarrier(回环栅栏)
- NS3官方参考资料manual和tutorial
- 用了Python,老板再也不用担心我写不了CUDA了!
- Jquery取得iframe中元素的几种方法(转载)
- 详细解读八大无线网络安全技术利弊
- php imagick gif,PHP基于php_imagick_st-Q8.dll实现JPG合成GIF图片的方法
- 9大最佳工程施工项目管理系统
- 领跑衫获奖感言 课程总结
- matlab text函数
- python查看list的shape_列表list、数组np.array等的len,size,shape操作
- 凑热闹,列个非主流书单:(1)分析与解决问题、通用管理(2)技术管理/CTO...
- 酷比魔方iwork1x i30双系统版,重装单Ubuntu 20.04系统记录(重力传感器、触摸屏、启动菜单、声卡问题、优化充电慢)
- pandas使用to_feather函数将dataframe保存为feather文件(需要依赖pyarrow包)提升大文件读取效率、pandas使用read_feather函数读取feather文件
- 杰里之BQB认证出现的一些问题分析篇
- Nutch安装第四天,进入正题,Nutch2.4的配置和编译
- UltraEdit编辑器中关键字自动补齐功能设置