计算方法——拉格朗日插值多项式

【问题描述】考虑[0.0,1.2]内的函数y=f(x)=cos(x)。利用多个（2,3,4等）节点构造拉格朗日插值多项式。

【输入形式】在屏幕上依次输入在区间[0.0,1.2]内的一个值x*，构造插值多项式后求其P(x*)值，和多个节点的x坐标。

【输出形式】输出插值多项式系数矩阵，拉格朗日系数多项式矩阵和P(x*)值（保留小数点后6位有效数字）。

【样例1输入】

0.3

0 0.6 1.2

【样例1输出】

[-0.40043538 -0.05084608 1. ]

[[ 1.38888889 -2.5 1. ]

[-2.77777778 3.33333333 -0. ]

[ 1.38888889 -0.83333333 0. ]]

P2(0.3)=0.948707

【评分标准】根据输入得到的输出准确

import numpy as np
import math
xx=input()
xx=float(xx)
Enter = input()
Enter = Enter.split(" ")
n= len(Enter)
X=np.zeros((1,n))
Y=np.zeros((1,n))
A=np.zeros((n,n))
X[0]=Enter
fxx=[]
fxxx=np.poly1d([0,0])for i in range(0,n):Y[0][i]=math.cos(X[0][i])for i in range(0,n):fxx.append(np.poly1d([0,1]))for i in range(0,n):v=np.poly1d([0,1])for j in range(0,n):if j!=i :fxx[i]=np.polymul(fxx[i],np.poly1d([1,-X[0][j]]))/(X[0][i]-X[0][j])v=np.polymul(v,np.poly1d([1,-X[0][j]]))/(X[0][i]-X[0][j])A[i]=v.cfor i in range(0,n):fxxx=fxx[i]*Y[0][i]+fxxxprint(fxxx.c)
print(A)
print('P{}({})={:.6f}'.format(n-1,xx,fxxx(xx)))'''''
0.3
0 0.6 1.2
'''''

计算方法——拉格朗日插值多项式相关推荐

拉格朗日多项式图像matlab,matlab实现拉格朗日函数，拉格朗日插值多项式
%拉格朗日插值多项式利用矩阵求解 x=1:0.2:3;%已知数据点x坐标向量:x y=sin(x);%已知数据点x坐标向量:y x1=1.1:0.2:3.1;%插值点的x坐标:x1 L=zeros( ...
求拉格朗日多项式matlab,拉格朗日插值多项式积分求圆周率近似Matlab实现
Lagrange 插值多项式积分求圆周率近似摘要: 公式1:y1=4/(1+x^2) 公式2:y2=4*sqrt(1-x^2) 分别对公式1.公式2求其拉格朗日插值多项式,再对其求0-1上的定积分来 ...
拉格朗日插值多项式的理解与C++实现
这篇记录一下我对拉格朗日插值的理解以及代码的实现. 线性插值: 如图给出了两个点(x0,y0),(x1,y1),假定x1>x0, 现在要估计某个给定x的函数值,其中x∈(x0,x1). 线性插 ...
数值计算方法-拉格朗日插值法
要求: 编写拉格朗日插值函数,输入为(x,f,order1) 其中x代表节点数据,为n*2矩阵 f代表是否分段 oreder1:若分段,oreder1代表函数阶次,只能为1或者2 若不分段,则函数阶次 ...
matlab做三次拉格朗日插值多项式_从零开始的matlab学习笔记——（24）曲线拟合...
matlab应用--求极限,求导,求积分,解方程,函数绘图,三维图像,拟合函数....更多内容尽在个人专栏:matlab学习才发现这个数据插值的返回值不是个函数表达式,是一系列点集,不能使用poly ...
matlab拉格朗日插值多项式代码
个人用,大家也可以直接复制,直接调用即可 function yy = lag(x,y,xx)m=length(x);n=length(y);if m~=n,error('')ends=0;for i= ...
拉格朗日Lagrange插值多项式
分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow 也欢迎大家转载本篇文章.分享知识,造福人民,实现我们中华民族伟大复兴! 本文转载 ...
【计算方法】拉格朗日插值法
概念设f(x)在N+1个点(x0,y0)-(xn,yn)处的值已知,其中值xk在区间[a,b]上,xk互不相同,满足a≤x0<x1<-<xn≤b,yk=f(xk).求任一插值点x对 ...
【例题+习题】【数值计算方法复习】【湘潭大学】（二）
第二章:函数基本逼近(一)--插值逼近目录第二章:函数基本逼近(一)--插值逼近写在前面的话知识点(重点) 拉格朗日插值公式 Neville 插值公式牛顿插值公式拉格朗日插值多项式和牛顿插 ...

计算方法——拉格朗日插值多项式

计算方法——拉格朗日插值多项式

计算方法——拉格朗日插值多项式相关推荐

最新文章

热门文章