拉格朗日三次插值公式_差值函数

拉格朗日三次插值公式_差值函数第一部分:问题分析(1)实验题目:拉格朗日插值算法具体实验要求:要求学生运用拉格朗日插值算法通过给定的平面上的n个数据点,计算拉格朗日多项式Pn(x)的值,并将其作为实际函数f(x)的估计值。用matlab编写拉格朗日插值算法的代码,要求代码实现用户输入了数据点(xi,f(xi))、插值点之后,程序能够输出插值点对应的函数估值。(2)实验目的:让同学们进一步掌握拉格朗日插值算法的原理以及运算过程,并且通过matlab编程培养实际的上机操作能力和代码能力。第二部分:数学原理要估计任一点..

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

第一部分:问题分析

(1)实验题目:拉格朗日插值算法

具体实验要求:要求学生运用拉格朗日插值算法通过给定的平面上的n个数据点,计算拉格朗日多项式Pn(x)的值,并将其作为实际函数f(x)的估计值。用matlab编写拉格朗日插值算法的代码,要求代码实现用户输入了数据点(xi,f(xi))、插值点之后,程序能够输出插值点对应的函数估值。

(2)实验目的:让同学们进一步掌握拉格朗日插值算法的原理以及运算过程,并且通过matlab编程培养实际的上机操作能力和代码能力。

第二部分:数学原理 

拉格朗日三次插值公式_差值函数 要估计任一点ξ,ξ≠xi,i=0,1,2,…,n,则可以用Pn(ξ)的值作为准确值f(ξ)的近似值,此方法叫做“插值法”。

称式(*)为插值条件(准则),含xi(i=0,1,…,n)的最小区间[a,b],其中a=min{x0,x1,…,xn},b=max{x0,x1,…,xn}。

具体实现形式:

拉格朗日三次插值公式_差值函数

 第三部分:程序设计流程

(1)langrange插值函数(被调用者):

拉格朗日三次插值公式_差值函数

(2)执行函数(面向用户:调用者):

拉格朗日三次插值公式_差值函数

 第四部分:代码实现

拉格朗日插值函数实现:

(每次最外层的for循环,打印一次插值点以及其期望值)

function[] = Lagrange(x,f,x0)
%得出数据点的个数
n = length(x) ;
%得出插值点的个数
m = length(x0);
    for i = 1:m
    D = x0(i);
    y = 0.0;
    for k = 1:n
    %置l(x)的初值
    l = 1.0;
        for j = 1:n
            %迭代表示出l(xj)
            %记住,第k个是取不到的!!!!
            if j~=k
            l = l*(D-x(j))/(x(k)-x(j));
            end
        end
        %Pn(x)
        y = y + l*f(k);
    end
        xx = num2str(D,'%.4f');
        y =num2str(y,'%.4f');
        disp('f(x)的近似值点坐标为:');
        disp(['(',xx,',',y,')']);
    end
end

Jetbrains全家桶1年46,售后保障稳定

 用户调用部分: 

(1)情形一:

(实现结果:给定插值函数的区间,给定原函数f(x),给定插值点–>实现拉格朗日估值的计算)

注意:

1.拉格朗日坐标点的x坐标向量,通过区间加步长实现:所以用户只需修改区间端点,程序便可以根据输入的区间等分次数n,自行合成

2.通过syms x将原函数f(x)表示出来(由用户自己输入),程序便可以根据输入的f(x)自动计算出对应数据节点的函数值 

%根据插值要求,自行算出所需要输出的坐标值
n = input("插值等分的次数n=");
X = -5:(10/n):5;
m = length(X);%获取节点的个数
Y = zeros(1,m);%创造全0向量为插值点函数值做铺垫
syms x;
%在此输入测试的函数f(x)
f(x) = (1+x^2)^(-1);
for i = 1:m
    Y(i) = f(X(i));
end
x0 = input('请输入插值点x0数组:');
Lagrange(X,Y,x0);

(2)情形二:

(实现结果:针对给定的无规律离散数据点,且未给定对应的原函数f(x))

%输入数据点的坐标
x = input("输入数据点x坐标的向量:");
y = input("输入数据点y坐标的向量:");

%输入所要估值的插值点的x坐标向量
x0 = input('请输入插值点x0数组:');

Lagrange(X,Y,x0);

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/230923.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 再次发布一个jq的下拉框式日期选择插件 dateSelector ~

    再次发布一个jq的下拉框式日期选择插件 dateSelector ~

  • python彩色文字输出

    python彩色文字输出

  • java h2 数据库_H2数据库介绍「建议收藏」

    java h2 数据库_H2数据库介绍「建议收藏」一、H2数据库简介1、H2数据库是一个开源的关系型数据库。H2是一个嵌入式数据库引擎,采用java语言编写,不受平台的限制,同时支持网络版和嵌入式版本,有比较好的兼容性,支持相当标准的sql标准,支持集群2、提供JDBC、ODBC访问接口,提供了非常友好的基于web的数据库管理界面二、在Java中操作H2数据库1、以嵌入式(本地)连接方式连接H2数据库这种连接方式默认情况下只允许有一个客户端连接到…

    2022年10月12日
  • android四种启动模式_android:windowSoftInputMode

    android四种启动模式_android:windowSoftInputMode(1)添加头文件:#include(2)在特定驱动结构体中添加early_suspend结构:#ifdefCONFIG_HAS_EARLYSUSPENDstructearly_suspendearly_suspend;#endif(3)在驱动probe函数中注册相关early_suspend结构体:#ifdefCONFIG_HAS_EARLYSUSPEND

  • pkl文件是什么_python pkl文件

    pkl文件是什么_python pkl文件这里只介绍关于字典类型和列表类型的数据的存储,其他格式的请君自行探索。1、字典类型importpickle#首先导入这个库,没有安装的话,自行百度,很简单dict_data={“name”:[“张三”,”李四”]}withopen(“dict_data.pkl”,’wb’)asfo:#将数据写入pkl文件pickle.dump(dict_data,fo)with…

  • discuz php接口文档,Discuz二次开发手册.doc[通俗易懂]

    discuz php接口文档,Discuz二次开发手册.doc[通俗易懂]Discuz二次开发手册Discuz文件说明,有助于discuz爱好者,进行自己的开发,在这里提供方便admincp.php——后台系统设置主程序文件,一般只处理菜单的显示的访问权限,不处理管理控制。ajax.php——论坛模板的ajax判断及数据返回都在这里进行attachment.php——附件文件,仅仅处理附件下载的功能。announcement.php——论坛公告的显示,一般很少改con…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号