Lagrange插值法实验:求拉格朗日插值多项式和对应x的近似值matlab实现(内附代码)

本文主要是介绍Lagrange插值法实验:求拉格朗日插值多项式和对应x的近似值matlab实现(内附代码),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、实验要求

已知函数表:
在这里插入图片描述
求出Lagrange 插值多项式,并计算x=1.2处的y的近似值。

二、MATLAB代码

求解多项式:

X = input('请输入横坐标向量X:\nX=');  % 获取用户输入的横坐标向量
Y = input('请输入纵坐标向量Y:\nY=');  % 获取用户输入的纵坐标向量
m = length(X);  % 获取节点个数
L = ones(m, m);  % 初始化基函数系数矩阵
for k = 1 : m  % 循环计算每个基函数的系数V = 1;  % 初始化基函数的临时变量for i = 1 : mif k ~= i  % 排除当前节点时的计算V = conv(V, poly(X(i))) / (X(k) - X(i));  % 在当前节点处计算基函数的系数endendL(k, :) = V;  % 将基函数的系数存储到矩阵中
end
l = sym(zeros(m, 1));  % 初始化基函数数组
for k = 1 : m  % 遍历基函数系数矩阵l(k) = poly2sym(L(k, :));  % 将基函数系数转化为符号表达式
end
fprintf('基函数为:\n');
for k = 1 : m  % 输出每个基函数的表达式fprintf('q%d(x)=%s\n', k, char(l(k)));
end
P = sym(0);  % 初始化拉格朗日多项式
for k = 1 : m  % 计算拉格朗日多项式的表达式P = P + Y(k) * l(k);
end
fprintf('拉格朗日多项式为:\nP(x)=%s\n', char(P));  % 输出拉格朗日多项式的表达式

求解近似值:

% 定义拉格朗日插值多项式的系数
coeff = [1/12, -29/6, 1/12, 62/3];% 计算 x=1.2 处的 y 的近似值
x = 1.2;
y_approx = polyval(coeff, x);% 显示结果
disp(['当 x = ' num2str(x) ' 时,y 的近似值为 ' num2str(y_approx)]);

三、实验结果

请输入横坐标向量X:
X=[1, 2, 4, 5]
请输入纵坐标向量Y:
Y=[16,12,8,9]
基函数为:
q1(x)=(11x^2)/12 - (19x)/6 - x^3/12 + 10/3
q2(x)=(29x)/6 - (5x^2)/3 + x^3/6 - 10/3
q3(x)=(4x^2)/3 - (17x)/6 - x^3/6 + 5/3
q4(x)=(7x)/6 - (7x^2)/12 + x^3/12 - 2/3
拉格朗日多项式为:
P(x)=x^2/12 - (29*x)/6 + x^3/12 + 62/3
在这里插入图片描述
当 x = 1.2 时,y 的近似值为 13.9507

参考:lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题)

这篇关于Lagrange插值法实验:求拉格朗日插值多项式和对应x的近似值matlab实现(内附代码)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1608

相关文章

Linux下删除乱码文件和目录的实现方式

《Linux下删除乱码文件和目录的实现方式》:本文主要介绍Linux下删除乱码文件和目录的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux下删除乱码文件和目录方法1方法2总结Linux下删除乱码文件和目录方法1使用ls -i命令找到文件或目录

SpringBoot+EasyExcel实现自定义复杂样式导入导出

《SpringBoot+EasyExcel实现自定义复杂样式导入导出》这篇文章主要为大家详细介绍了SpringBoot如何结果EasyExcel实现自定义复杂样式导入导出功能,文中的示例代码讲解详细,... 目录安装处理自定义导出复杂场景1、列不固定,动态列2、动态下拉3、自定义锁定行/列,添加密码4、合并

mybatis执行insert返回id实现详解

《mybatis执行insert返回id实现详解》MyBatis插入操作默认返回受影响行数,需通过useGeneratedKeys+keyProperty或selectKey获取主键ID,确保主键为自... 目录 两种方式获取自增 ID:1. ​​useGeneratedKeys+keyProperty(推

Spring Boot集成Druid实现数据源管理与监控的详细步骤

《SpringBoot集成Druid实现数据源管理与监控的详细步骤》本文介绍如何在SpringBoot项目中集成Druid数据库连接池,包括环境搭建、Maven依赖配置、SpringBoot配置文件... 目录1. 引言1.1 环境准备1.2 Druid介绍2. 配置Druid连接池3. 查看Druid监控

Linux在线解压jar包的实现方式

《Linux在线解压jar包的实现方式》:本文主要介绍Linux在线解压jar包的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux在线解压jar包解压 jar包的步骤总结Linux在线解压jar包在 Centos 中解压 jar 包可以使用 u

c++ 类成员变量默认初始值的实现

《c++类成员变量默认初始值的实现》本文主要介绍了c++类成员变量默认初始值,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录C++类成员变量初始化c++类的变量的初始化在C++中,如果使用类成员变量时未给定其初始值,那么它将被

Qt使用QSqlDatabase连接MySQL实现增删改查功能

《Qt使用QSqlDatabase连接MySQL实现增删改查功能》这篇文章主要为大家详细介绍了Qt如何使用QSqlDatabase连接MySQL实现增删改查功能,文中的示例代码讲解详细,感兴趣的小伙伴... 目录一、创建数据表二、连接mysql数据库三、封装成一个完整的轻量级 ORM 风格类3.1 表结构

基于Python实现一个图片拆分工具

《基于Python实现一个图片拆分工具》这篇文章主要为大家详细介绍了如何基于Python实现一个图片拆分工具,可以根据需要的行数和列数进行拆分,感兴趣的小伙伴可以跟随小编一起学习一下... 简单介绍先自己选择输入的图片,默认是输出到项目文件夹中,可以自己选择其他的文件夹,选择需要拆分的行数和列数,可以通过

Python中将嵌套列表扁平化的多种实现方法

《Python中将嵌套列表扁平化的多种实现方法》在Python编程中,我们常常会遇到需要将嵌套列表(即列表中包含列表)转换为一个一维的扁平列表的需求,本文将给大家介绍了多种实现这一目标的方法,需要的朋... 目录python中将嵌套列表扁平化的方法技术背景实现步骤1. 使用嵌套列表推导式2. 使用itert

Python使用pip工具实现包自动更新的多种方法

《Python使用pip工具实现包自动更新的多种方法》本文深入探讨了使用Python的pip工具实现包自动更新的各种方法和技术,我们将从基础概念开始,逐步介绍手动更新方法、自动化脚本编写、结合CI/C... 目录1. 背景介绍1.1 目的和范围1.2 预期读者1.3 文档结构概述1.4 术语表1.4.1 核