计算方法——插值法程序实现(一)

例题

给出的函数关系表,分别利用线性插值及二次插值计算的近似值。

|-----------------------------------|----------|----------|----------|----------|----------|
| | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 |
| | 1.105171 | 1.221403 | 1.349859 | 1.491825 | 1.648721 |

参考代码一:Python代码实现(自编码)

python 复制代码
import math
"""
:parameter用于计算插值多项式的系数
"""
def Parameters(data_x,data_y,size):
    parameters=[]
    i = 0
    while i<size:
        j = 0
        temp = 1
        while j<size:
            if i!=j:
                temp*=data_x[i]-data_x[j];
            j += 1
        parameters.append(data_y[i]/temp);
        i += 1
    return parameters
"""
:Calculate用于计算拉格朗日插值公式计算后的值
"""
def Calculate(data_x,parameters,x):
    return_value = 0;
    i = 0
    while i < len(parameters):
        temp = 1
        j = 0
        while j<len(parameters):
            if i!=j:
                temp*=x-data_x[j];
            j+=1
        return_value+=temp*parameters[i]
        i+=1
    return return_value

x1 = [0.2,0.3];
y1 = [1.221403,1.349859];
p1 = Parameters(x1,y1,len(x1));
target_point = 0.27;
real_value = math.exp(target_point);
fit_value = Calculate(x1,p1,target_point)
remainder = real_value-fit_value
print("===============线性拉格朗日插值==================");
print(f"点{target_point}处的真实值为{real_value}")
print(f"点{target_point}处的拟合值为{fit_value}")
print(f'拉格朗日插值余项为{remainder}');
print("===============二次拉格朗日插值==================");
x2 = [0.2,0.3,0.4];
y2 = [1.221403,1.349859,1.491825];
p2 = Parameters(x2,y2,len(x2))
fit_value_2 = Calculate(x2,p2,target_point)
remainder2 = real_value-fit_value_2
print(f"点{target_point}处的拟合值为{fit_value_2}");
print(f'拉格朗日插值余项为{remainder2}');

Python编码计算结果

参考代码二:MATLAB代码实现(自编码)

Matlab 复制代码
%可运行部分代码
clc,clear
format long
x1=[0.2,0.3];
y1=[1.221403,1.349859];
p1=Parameters(x1,y1,length(x1));
point=0.27;
real = exp(point);
disp('点0.27处的真实值为');
disp(real);
fitted1 = Calculate(x1,p1,point);
Remainder1 = real-fitted1;
disp('点0.27处的线性拉格朗日插值结果为');
disp(fitted1);
disp('点0.27处的线性拉格朗日插值余项为');
disp(Remainder1);
x2 = [0.2,0.3,0.4];
y2 = [1.221403,1.349859,1.491825];
p2 = Parameters(x2,y2,length(x2));
fitted2 = Calculate(x2,p2,point);
Remainder2 = fitted2-real;
disp('点0.27处的二次拉格朗日插值结果为');
disp(fitted2);
disp('点0.27处的二次拉格朗日插值余项为');
disp(Remainder2);
format short

Parameters函数

Parameter用于计算插值多项式的系数

Matlab 复制代码
function parameters = Parameters(data_x,data_y,size)
parameters=[];
i=1;
while i<=size
    j=1;
    temp=1;
    while j<=size
        if i ~= j
            temp=(temp*(data_x(i)-data_x(j)));
        end
        % disp(temp);
        j=j+1;
    end
    parameters=[parameters,data_y(i)/temp];
    % disp(parameters);
    i=i+1;
end

Calculate函数

Calculate用于计算拉格朗日插值公式计算后的值

Matlab 复制代码
function return_value = Calculate(data_x,parametres,x)
return_value=0;
i=1;
while i<=length(parametres)
    temp=1;
    j=1;
    while j<=length(parametres)
        if i~=j
            temp = temp*(x-data_x(j));
        end
        j=j+1;
    end
    return_value = return_value+temp*parametres(i);
    i = i+1;
end

MATLAB编码计算结果

参考书目

[1] 李新栋,许文文,张绪浩,任永强. 基于Python的计算方法[M]. 北京:电子工业出版社,2023.

相关推荐
池央5 分钟前
AI性能极致体验:通过阿里云平台高效调用满血版DeepSeek-R1模型
人工智能·阿里云·云计算
我们的五年6 分钟前
DeepSeek 和 ChatGPT 在特定任务中的表现:逻辑推理与创意生成
人工智能·chatgpt·ai作画·deepseek
Yan-英杰7 分钟前
百度搜索和文心智能体接入DeepSeek满血版——AI搜索的新纪元
图像处理·人工智能·python·深度学习·deepseek
Fuweizn9 分钟前
富唯智能可重构柔性装配产线:以智能协同赋能制造业升级
人工智能·智能机器人·复合机器人
小赵起名困难户1 小时前
蓝桥杯备赛1-2合法日期
算法
shichaog1 小时前
腿足机器人之八- 腿足机器人动力学
算法·机器人
weixin_307779131 小时前
Azure上基于OpenAI GPT-4模型验证行政区域数据的设计方案
数据仓库·python·云计算·aws
玩电脑的辣条哥2 小时前
Python如何播放本地音乐并在web页面播放
开发语言·前端·python
taoqick2 小时前
对PosWiseFFN的改进: MoE、PKM、UltraMem
人工智能·pytorch·深度学习
suibian52352 小时前
AI时代:前端开发的职业发展路径拓宽
前端·人工智能