数学建模——插值算法

news2024/9/23 13:27:52

概念:数模比赛中,常常需要根据有已知的函数点进行数、模型处理和分析,而有时候现有的数据是极少的,不足以支撑分析的进行,这时就需要使用一些数学的方法,“模拟产生“一些新的但又比较靠谱的值来满足需求,这就是插值的作用。

一维插值问题:

通过已有的点和对应的值,构造函数将其要插入的横坐标但如得到插入点(纵坐标)

插值法的定义:

插值法的分类:

插值法原理:

拉格朗日插值法与牛顿插值法:

(这两种方法只做了解,为的是引出龙格现象和后续我们mtalab中真正常用到两种方法,所以下面只放一些相关的ppt)

 

 

 

龙格现象:

(直白点就是插入的次数越多不已经数据就越精确,反而会产生龙格现象:即放插入值多是图像更加偏离)

 

(所以我们在不熟悉曲线运动的趋势时少用多次插值)

如何解决龙格现象

分段插值即:在要插入点x(横坐标)左右各找一点,将其对应的值连接成一个线段,x在线段上对应的值y即是我们要插入的点。这样故好的那束次数过多,导致数据不精确

如何使得插值法更加精确呢?

A.埃尔米特(Hermite)插值

在matlab中用pchip函数来表示分段三次埃尔米特插值法,pchip(x,y,new_x)中x和y可以是一个n列一行的矩阵,new_x的长度是插入后的长度。

Plot是绘图函数,plot(x,y,‘线段或点的形状’),

Eg

三次样条插值:

 

Matlab中三次样条插值和分段三次埃尔比特插值法相只是函数名字不同,传参和返回值都是一样的

 分段三次埃尔米特插值法和三次样条插值法绘图后的比较(我么用的是用这两个插值函数购构建sinx函数图像):

不难看出三次样条插值法更光滑

插值法解决生活实例的例题Mtalab代码的实现与讲解:

本题是给出我们1 3 5 7 9 11 13 15周的数据让我们利用插值法将1-15周的数据补全,且用matlab构建每个指标的图像。(只需要补全一号池的数据即可)

其中所用到的函数的讲解:

figure(num)函数:在同一个脚本文件里面,需要给每个图像编号,否则只会显示最后一个图像。(如果我们每个图像在单独的一个面板且展示的时候需要都展示出来,因此我们需要使用到figure函数。)

eg:

legend函数:标注图像中的线段/点对应的插值方法和显示的位置等。

Legend(‘名字’,‘location(坐标一般写位置是需要先写这个)’,‘方向(eg:SouthEast东南)’)

Eg:

 代码实现及其注释讲解:

clear;clc
load data.mat;
x=Z(1,:);
[n,m]=size(Z)%求出x的行(数据个数)和列(周数)
data_name={'周数','轮虫','溶氧','COD','水温','PH值','盐度','透明度','总碱度','氯离子','透明度','生物量'};
disp(['共有' num2str(n-1) '个指标进行插值']);
disp('正在进行一号池三次样条插值,请等待:');
%设置一个矩阵P,保存插入好的数据
%1.开辟一个矩阵(大小为n-1行和15列)
P=zeros(n-1,15);
for i=2:n
new_x=1:15;%插入后的总列数
y=Z(i,:);
%将每次插入好的数据保存下来
pi=spline(x,y,new_x);
%绘图
%a.设置图形的大小
%subpolt(x,y,i)的作用:
%我们一共要绘制11个图,为了方面观看需要将这11个图放在同一个面板上观察,而subpolt的作用就是绘制这样一个面板的函数
%前两个参数x,y表示这个面版最多能方4*3个图,一共四行,每行最多方3个,i则代表了位置,从左->右,上->下依次排序1-x*y
subplot(4,3,i-1);
%b.传入数据构建图形
plot(x,y,'ro',new_x,pi,'-');
%c.设置横坐标的长度(1-15),纵坐标不进行改变
%axis函数设置图形的横纵坐标范围
axis([0 15 -inf inf])
%d.设置图形的名称(每个图形标记好相对应的指标名称)
%ylable为图像设置名称
ylabel(data_name{i})
%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
%将每次保存的插入好的输入依次保存到一个大的矩阵当中
P(i-1,:)=pi;
end
%为每组图中的线段和点标记好所用的插入方法:
legend('原始数据','三次样条插值数据','Location','SouthEast')
%周数插入到P矩阵的上方
P=[1:15;P];


最后的插值数据完整后的结果和利用matlab构建的图像:

 

 博主主要跟着清风数学建模的课程学习,其中里面的一些图片都来源于上课视频的截图。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/749247.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

为什么很多人说不建议学JAVA,说很难找到工作?

前言 在回答这个问题之前,我曾经也看到了类似的其他语言相关的问题: 那么现在问题来了,什么语言都不建议学,那什么是可以学的呢?看多了所谓的不建议学习的文章和回答固然可以了解到该语言的些许缺点和当前环境下的一些…

NTLM协议

0x01 LM Hash加密算法 LM Hash本质是一个DES加密,从Windows Vista和Windows Server 2008开始就禁用了LM Hash。LM Hash明文密码被限定在14位以内,密码超过14位则会停止使用LM Hash。下面介绍一下LM Hash的加密流程,以Pss1234为例&#xff1a…

Element-Plus搭建CMS页面结构 引入第三方图标库iconfont(详细)

Element-Plus组件库使用 element plus组件库是由饿了么前端团队专门针对vue框架开发的组件库,专门用于电脑端网页的。因为里面集成了很多组件,所以使用他可以非常快速的帮我们实现网站的开发。 安装: npm install element-plus --save 引入…

制作Visual Studio离线安装包

vs2015之后官网就不提供离线安装包了,使用离线安装包就需要自己手动制作一个; 以vs2019为例: 先去官网下载在线安装器 官网下载地址:Visual Studio 较旧的下载 - 2019、2017、2015 和以前的版本 (microsoft.com) 展开2019的标签…

表格tr根据内容撑开高度并且每个td继承tr的高度(height:inherit,align-items: stretch)

项目设计图 实际上数据又长又短,例如类型有可能出现2行甚至3行的字数,所以整个行要适配最高的高度,然后其他子元素继承行的高度。 ps:原生的table或者uview的u-table都不怎么好使,所以用view来做 在做的过程中出现高度…

澎湃认证:浪潮信息携手可利邦推出隐私计算一体机解决方案

北京可利邦与浪潮信息联合发布隐私计算一体机解决方案,以浪潮信息inMerge超融合产品为基础设施平台,预装可利邦隐私计算软件平台,通过SRDC一体机进行整机柜交付,助力客户快速上线隐私计算应用系统。该方案已通过浪潮信息澎湃技术认…

Node连接Mongodb数据库

1.初始化 npm init 2.安装mongoose npm i mongoose 3.导入mongoose const mongooserequire("mongoose") 4.连接mongodb服务 mongoose.connect("mongodb://127.0.0.1:27017/user") 说明:mongodb是协议,user是数据库,如果没有会自动创…

Oracle通过函数调用dblink同步表数据方案(全量/增量)

创建对应的包,以方便触发调用 /*包声明*/ CREATE OR REPLACE PACKAGE yjb.pkg_scene_job AS /*创建同步任务*/FUNCTION F_SYNC_DRUG_STOCK RETURN NUMBER;/*同步*/PROCEDURE PRC_SYNC_DRUG_STOCK(RUNJOB VARCHAR2) ; END pkg_scene_job; /*包体*/ CREATE OR REPL…

Dubbo最核心功能——服务暴露的配置、使用及原理

系列文章目录 【收藏向】从用法到源码,一篇文章让你精通Dubbo的SPI机制 面试Dubbo ,却问我和Springcloud有什么区别? 超简单,手把手教你搭建Dubbo工程(内附源码) 文章目录 系列文章目录前言一、服务暴露的…

idea乱码解决方式大汇总

目录 idea版本: 解决方法: 一、基本方法 1. File -> Settings -> Editor 2. 二、Maven乱码解决方法 三、运行时乱码解决方法 四、因为以前乱设置导致的乱码 idea版本: 解决方法: 一、基本方法 1. File -> Setti…

ETHERCAT主站网关转DEVICENET连接支持ethercat总线的PLC

大家好,今天要和大家分享一款自主研发的通讯网关——远创智控YC-ECTM-DNT。这款产品可是解决了不同协议设备数据交换的麻烦问题,让我们一起来看看它的神奇之处吧! 这款通讯网关有什么特别的呢?首先,它可以连接DEVICENE…

vue2 element-ui el-cascader地址省市区分开单独写

使用 npm 或 yarn 安装 element-china-area-data 包&#xff1a; npm install element-china-area-data 在你的代码中导入 element-china-area-data import { regionData } from element-china-area-data let that; 完整代码 <template><div><el-form ref&quo…

Zookeeper概述及部署

Zookeeper概述及部署 一、Zookeeper 定义二、Zookeeper 特点三、Zookeeper 数据结构四、Zookeeper 应用场景五、Zookeeper选举机制● 第一次启动选举机制● 非第一次启动选举机制 六、部署 Zookeeper 集群1.安装前准备2.安装 Zookeeper 一、Zookeeper 定义 Zookeeper是一个开源…

leetcode:LCP 06. 拿硬币(python3解法)

难度&#xff1a;简单 桌上有 n 堆力扣币&#xff0c;每堆的数量保存在数组 coins 中。我们每次可以选择任意一堆&#xff0c;拿走其中的一枚或者两枚&#xff0c;求拿完所有力扣币的最少次数。 示例 1&#xff1a; 输入&#xff1a;[4,2,1] 输出&#xff1a;4 解释&#xff1a…

mysql-分页数据重复

背景说明 分页查询不同页出现重复数据&#xff0c;底层实现都是使用limit select * from tt1 limit 100,50 第101条开始的50条数据 select * from tt1 limit 100 前100条记录 select * from tt1 limit 100 offset 50 从51条开始&#xff0c;显示后面的100条 原因分析 以前碰…

力扣 53. 最大子数组和

题目来源&#xff1a;https://leetcode.cn/problems/maximum-subarray/description/ C题解&#xff08;来源代码随想录&#xff09;&#xff1a;难点在于判断什么时候更新count。 是当nums[i]为正数吗&#xff1f;那要nums没有正数怎么办&#xff1b; 是当nums[i]比之前子数组…

网络监控的关键指标

网络监控是 IT 的支柱&#xff0c;对于防止可能困扰您的业务的计划外中断至关重要。这就是为什么投资网络监控解决方案是一个明智而安全的举动。特别是随着远程工作成为常态&#xff0c;监控工具使公司能够从世界任何地方监控网络&#xff0c;并有助于在整个过程中保持最佳性能…

小程序:页面跳转闪屏

自己的笔记&#xff0c;随手记录。扛精走开。 1、问题描述 进入页面&#xff0c;是一个组件&#xff0c;通过路由传参判断是由哪个页面进入&#xff0c;不同的页面拿的已选值不一样&#xff0c;需要回显值&#xff0c;在编辑数据。此时会出现一个问题&#xff0c;A页面中进来…

运输问题案例

案例1 运输问题 某部门有3个生产同类型产品的产地&#xff0c;生产的产品由4个销售点出售&#xff0c;各工厂的生产量、各销售点的销售量以及各工厂到各销售点的单位运价&#xff08;元/吨&#xff09;如表1所示&#xff0c;求最佳调运方案&#xff1f; 表1 运输信息 产地\销地…

general 未设置cookie的Secure标志位

解决方案: 在配置文件中增加相应配置即可完成配置&#xff1a; <system.web><httpCookies httponlyCookies"true" requireSSl"true"/></system.web>添加完成后如下所示: