利用线性函数做插值

每一段的线性函数：

#Program 0.6 Linear Interploation

import numpy as np

import matplotlib.pyplot as plt

#分段线性插值闭包

def get_line(xn, yn):

def line(x):

index = -1

#找出x所在的区间

for i in range(1, len(xn)):

if x <= xn[i]:

index = i-1

break

else:

i += 1

if index == -1:

return -100

#插值

result = (x-xn[index+1])*yn[index]/float((xn[index]-xn[index+1])) + (x-xn[index])*yn[index+1]/float((xn[index+1]-xn[index]))

return result

return line

xn = [i for i in range(-50,50,10)]

yn = [i**2 for i in xn]

#分段线性插值函数

lin = get_line(xn, yn)

x = [i for i in range(-50, 40)]

y = [lin(i) for i in x]

#画图

plt.plot(xn, yn, 'ro')

plt.plot(x, y, 'b-')

plt.show()

python分段算利润、税收

''' 题目:企业发放的奖金根据利润提成. 利润(I)低于或等于10万元时,奖金可提10%: 利润高于10万元,低于20万元时,低于10万元的部分按10%提成,高于10万元的部分,可提成7.5%: 2 ...

Python 分段利润提成

题目:企业发放的奖金根据利润提成.利润(I)低于或等于10万元时,奖金可提10%:利润高于10万元,低于20万元时,低于10万元的部分按10%提成,高于10万元的部分,可提成7.5%:20万到40万之 ...

Matlab插值函数

x=0:2*pi; y=sin(x); xx=0:0.5:2*pi; %interp1对sin函数进行分段线性插值,调用interp1的时候,默认的是分段线性插值 y1=interp1(x,y,xx) ...

建模算法(八)&mdash&semi;&mdash&semi;插值

插值:求过已知有限个数据点的近似函数 拟合:已知有限个数据点,求近似函数,不要求过已知数据点,只要求在某种意义下在这些点的误差最小 (一)插值方法 一.拉格朗日多项式插值 1.插值多项式 就是做出一个 ...

【数学建模】day04-插值与拟合

关于插值原理,这篇文章里总结过. 插值,是在有限个数据点的情况下,模拟出更多的点来适应实际问题的需要. 拟合,是在已知数据点基础上,以已知点处最小误差为标准,模拟出近似函数. 二者有似,实则不同,ma ...

MATLAB数学实验总结

L1 MATLAB 基础知识 P6 表1-3 数据显示格式 format rat format long P20 表2-5 常用的矩阵函数 zeros(m,n) %零阵 eye(n) %单位阵 one ...

Matlab随笔之插值与拟合(上)

原文:Matlab随笔之插值与拟合(上) 1.拉格朗日插值 新建如下函数: function y=lagrange(x0,y0,x) %拉格朗日插值函数 %n 个节点数据以数组 x0, y0 输入(注 ...

Matlab随笔之插值与拟合(下)

原文:Matlab随笔之插值与拟合(下) 1.二维插值之插值节点为网格节点 已知m x n个节点:(xi,yj,zij)(i=1…m,j=1…n),且xi,yi递增.求(x,y)处的插值z. Matl ...

Python Cookbook&lpar;第3版&rpar;中文版：15&period;21 诊断分段错误

15.21 诊断分段错误¶ 问题¶ 解释器因为某个分段错误.总线错误.访问越界或其他致命错误而突然间奔溃. 你想获得Python堆栈信息,从而找出在发生错误的时候你的程序运行点. 解决方案¶ faul ...

随机推荐

SPL--spl&lowbar;autoload&lowbar;register

spl_autoload_register() : 调用未定义类时,系统会按顺序调用注册到spl_autoload_register()函数的所有函数,而不是调用__autoload函数. 解决问题: ...

黄聪：Discuz&excl;的SEO优化策略二：如何去掉页脚多余的信息

论坛搭建好,首先是把多余的东西都砍掉. 页脚的信息在我看来,都是很多余的信息,如下图: 要怎么消灭掉它们呢? 1.进入 全局 -- 站点信息 2.站点名称改为你的论坛名称,它会出现在内页的标题最末位. ...

擦亮自己的眼睛去看SQLServer之简单Insert&lpar;转&rpar;

摘要:本来是打算先写SQLServer历史的,不过感觉写那部分内容比较难还需要多查些资料.于是调整了下顺序写下简单的Insert语句. 不过感觉写那部分内容比较难还需要多查些资料.于是调整了下顺序写下 ...

tomcat session思考

最近在做统计服务器在线人数时,按照当前服务器session数量来实现,具体方法:实现HttpSessionListener接口来,每次sessionCreate的时候增加1,sessionDestro ...

codevs 1036 商务旅行 (倍增LCA)

/* 在我还不知道LCA之前 暴力跑的SPFA 70分 三个点TLE */ #include #include #include

mysql习题

如图表创建数据库. create table class( cid int auto_increment primary key, caption ) )engine=innodb default c ...

PHP算法学习(5) 位运算

svn地址:svn://gitee.com/zxadmin/live_z 2019年2月14日11:38:46 <?php /* * 位运算学习笔记 * * 1,php所有的数都是有符号的,无法 ...

面试：C&plus;&plus;不可继承类

面试中可能遇到让设计一个无法被继承的类.最简单的实现是将该类的构造函数设置为私有的,然后通过静态成员函数调用私有构造函数实例化对象,这样的类确实不可继承,但是使用起来非常不方便,必须使用静态成员实例化 ...

bootstrap学习参考网站

----https://www.     evget    .com   /article /

&num; 20155337《网络对抗》Web基础

20155337Exp8 Web基础 实践目标 1. 实践内容 (1).Web前端HTML 能正常安装.启停Apache.理解HTML,理解表单,理解GET与POST方法,编写 ...