牛顿插值java_java实现牛顿插值法

import java.util.Scanner;

public class Newton_interpolation {

/*拷贝向量*/

private static void copy_vector(double from[],double to[]){

int k=from.length;

int k2=to.length;

if(k!=k2){

System.out.println("the two vector's length is not equal!");

System.exit(0);

}

for(int i=0;i

to[i]=from[i];

}

/*牛顿插值法*/

private static double[] Newton_inter_method(double[] X,double[] Y,double X0[]){

int m=X.length;

int n=X0.length;

double[] Y0=new double[n];

double[] cp_Y=new double[m];

for(int i1=0;i1

double t=0;

int j=0;

copy_vector(Y, cp_Y);

int kk=j;

/*求各级均差*/

while(kk

kk=kk+1;

for(int i2=kk;i2

cp_Y[i2]=(cp_Y[i2]-cp_Y[kk-1])/(X[i2]-X[kk-1]);

}

/*求插值结果*/

double temp=cp_Y[0];

for(int i=1;i<=m-1;i++){

double u=1;

int jj=0;

while(jj

u*=(X0[i1]-X[jj]);

jj++;

}

temp+=cp_Y[i]*u;

}

Y0[i1]=temp;

}

return Y0;

}

public static void main(String[] args) {

/*输入插值点横纵坐标*/

System.out.println("Input number of interpolation point:");

Scanner scan=new Scanner(System.in);

int m=scan.nextInt();

System.out.println("Input number of test point:");

int n=scan.nextInt();

double X[]=new double[m];

double Y[]=new double[m];

double X0[]=new double[n];

System.out.println("Input the elements of X:");//已知插值点

for(int i=0;i

X[i]=scan.nextDouble();

}

System.out.println("Input the elements of Y:");//已知插值点的函数值

for(int i=0;i

Y[i]=scan.nextDouble();

}

System.out.println("Input the elements of X0:");//需要求的插值点的横坐标标值

for(int i=0;i

X0[i]=scan.nextDouble();

}

牛顿插值java_java实现牛顿插值法相关推荐

matlab差商chaS,牛顿插值MATLAB算法.doc
MATLAB程序设计期中作业 --编程实现牛顿插值成员:刘川(P091712797) 签名＿＿＿＿＿汤意(P091712817) 签名＿＿＿＿＿王功贺(P091712799) 签名＿＿＿＿＿ ...
数值分析之拉格朗日插值、牛顿插值、分段线性插值实现
1.拉格朗日插值法考虑全局信息的比较经典的插值方法,编程简单,计算量大. #coding=utf-8 from matplotlib import pyplot as pltdef Lg(data, ...
1月16日：拉格朗日中值定理，罗尔定理，柯西中值，拉格朗日插值，牛顿插值，重心插值，拉格朗日乘子法的证明
拉格朗日中值定理 https://www.bilibili.com/video/BV117411E7kx?from=search&seid=17921778669593975548 拉格朗日中 ...
函数插值的python实现——拉格朗日、牛顿插值
函数插值的python实现--拉格朗日.牛顿插值 1. 拉格朗日(Larange)插值实现代码 2. 牛顿(Larange)插值 2.1 牛顿插值多项式的基本形式 2.2 牛顿均差插值多项式 (1) ...
牛顿插值与正交多项式曲线拟合的Matlab实现
一.相关知识提要 1.Newton插值法 2. 正交多项式的方法作曲线拟合二.代码实现 1.Newton插值法 function [] = Newton_Chazhi(x ...
Hermite插值是牛顿插值的极限情形
Hermite插值可以看作牛顿插值的极限状况.为什么可以这么说呢?我们来看一个实例: 构造一个三次多项式 $p_3$ 使得 $p_3(0)=0$,$p_3(1)=1,p_3'(0)=1,p_3'(1) ...
matlab全域基函数,多项式函数插值：全域多项式插值（一）单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...
全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...
插值问题（拉格朗日插值、牛顿插值）
agui_lagrange.m: function f=agui_lagrange(x0,y0,x) % x0为节点向量,y0为节点上的函数值,x为插值点,f为返回插值 n=length(x0);m= ...
拉格朗日插值和牛顿插值的龙格现象
文章目录一.实验目的二.实验设备信息三.实验内容 (一)拉格朗日插值多项式 (二)牛顿插值多项式四.实验步骤 (一)拉格朗日插值函数实现 (二)牛顿插值函数实现 (三)观察拉格朗日插值和牛顿插 ...
重读微积分（六）：差商与牛顿插值
一理解极限二微分学本系列所有代码皆用R语言完成. 4 差商与牛顿插值如果采取间隔不等的采样,差商会变得稍显复杂,对于x0,x1,-,xnx_0,x_1,\ldots,x_nx0,x1,- ...

牛顿插值java_java实现牛顿插值法

牛顿插值java_java实现牛顿插值法相关推荐

最新文章

热门文章