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

一、实验要求

已知函数表:

求出Lagrange 插值多项式,并计算x=1.2处的y的近似值。

二、MATLAB代码

求解多项式:

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 : m
        if k ~= i  % 排除当前节点时的计算
            V = conv(V, poly(X(i))) / (X(k) - X(i));  % 在当前节点处计算基函数的系数
        end
    end
    L(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));  % 输出拉格朗日多项式的表达式

求解近似值:

matlab 复制代码
% 定义拉格朗日插值多项式的系数
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 - (19 x)/6 - x^3/12 + 10/3

q2(x)=(29x)/6 - (5 x^2)/3 + x^3/6 - 10/3

q3(x)=(4x^2)/3 - (17 x)/6 - x^3/6 + 5/3

q4(x)=(7x)/6 - (7 x^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实现(内附代码及例题)

相关推荐
coding随想1 小时前
JavaScript中的BOM:Window对象全解析
开发语言·javascript·ecmascript
念九_ysl1 小时前
Java 使用 OpenHTMLToPDF + Batik 将含 SVG 遮罩的 HTML 转为 PDF 的完整实践
java·开发语言·pdf
yaoxin5211232 小时前
124. Java 泛型 - 有界类型参数
java·开发语言
liulilittle2 小时前
深度剖析:OPENPPP2 libtcpip 实现原理与架构设计
开发语言·网络·c++·tcp/ip·智能路由器·tcp·通信
88号技师2 小时前
2025年6月一区-田忌赛马优化算法Tianji’s horse racing optimization-附Matlab免费代码
开发语言·算法·matlab·优化算法
勤奋的知更鸟2 小时前
Java 编程之模板方法模式
java·开发语言·模板方法模式
上单带刀不带妹3 小时前
手写 Vue 中虚拟 DOM 到真实 DOM 的完整过程
开发语言·前端·javascript·vue.js·前端框架
-凌凌漆-4 小时前
【Qt】QStringLiteral 介绍
开发语言·qt
程序员爱钓鱼4 小时前
Go语言项目工程化 — 常见开发工具与 CI/CD 支持
开发语言·后端·golang·gin
军训猫猫头4 小时前
1.如何对多个控件进行高效的绑定 C#例子 WPF例子
开发语言·算法·c#·.net