数学建模-最优包衣厚度终点判别法(主成分分析)

💞💞 前言
hello hello~ ,这里是viperrrrrrr~💖💖 ,欢迎大家点赞🥳🥳关注💥💥收藏🌹🌹🌹
💥个人主页:viperrrrrrr的博客
💥 欢迎学习数学建模算法、大数据、前端等知识,让我们一起向目标进发!


基于近红外光谱的肠溶片最优包衣厚度终点判别法

       包衣是将片剂的外表面均匀地包裹上一层衣膜的过程,旨在控制药物在胃肠道中的释放部位和速度,遮盖苦味或不良气味,防潮、避光,改善外观等。然而,包衣膜太薄或太厚都不利于药效,并且包衣终点的判断方法目前存在一定的难度。近红外光谱技术(NIRS)是一种高效、无需试剂、无污染的分析方法,通过近红外光谱仪、化学计量学软件和应用模型,能快速、简便地实现多组分检测。为实现包衣终点的准确判断,对数据进行分析并完成以下问题:

问题一:对药品在不同包衣时间段包衣片剂的近红外光谱进行特征峰提取,选择具有有效信息的波长片段,即波长选择。

问题二:分析药品包衣厚度分类规律,建立合适的模型对药品包衣不同厚度进行划分,给出方法及结果,并进行灵敏度分析。

问题三:对于不同的包衣厚度,通过建立模型分析包衣之间的关联性,判别出最优的包衣厚度。

我们本次主要解决问题一

问题一

包衣是将片剂均匀地包裹衣膜,用于控制药物释放、遮盖苦味等。然而,包衣膜太薄或太厚都不利于药效,并且包衣终点的判断方法目前存在一定的难度。而近红外光谱技术是一种无污染、快速、多组分检测的分析方法,它适用于包衣终点的确定。本文基于一批现有的红外光谱相关数据,建立主成分分析聚类分析Bayes判别等模型和梯度下降算法,实现了片剂包衣最佳终点的判断。建立了主成分分析模型进行特征峰选取,利用主成分分析中广泛使用的降维技术,通过线性变换将高维数据集转化为低维数据集,同时保留数据集中的主要信息。由于不同的特征会有不同的量纲,这可能会影响到模型的性能,因此在应用PCA之前,通常需要对数据进行数据预处理。列出协方差矩阵反映数据集中各特征之间相关性的矩阵。通过求解协方差矩阵的特征向量和特征值,可以得到数据集的主成分。这些主成分是原始特征空间中的线性变换,它们是新的、相互独立的、能够捕捉到数据变动的最大方差的向量。

首先,通过滑动平均滤波法对附件数据进行平滑处理,即滤波。基于附件数据,由图(1)可见原始波长噪音数据较多,会影响最终数据结果[1]。使用5倍滑动滤波进行降噪处理后,由图(2)可见,成功过滤了部分噪声数据,使光谱变得更加平滑。

图1 原始数据图

图2 5倍滑动滤波处理后数据图

5.1.2主成分分析模型

根据问题一的分析,我们建立了主成分分析模型进行了对不同包衣时间段包衣片剂的近红外光谱的特征峰选取。首先,将波长定义为X轴,在不同包衣时间下的片剂包衣、素片为Y 轴,在此基础上进行主成分分析,在进行特征根选取时,发现有两种情况可以选择,第一种情况是选取特征根大于1的成分,可以选取出2个主成分,第二种情况是按照公式(1)算取方差贡献率以选取特征根,可以选取出3个主成分。再算取累计方差贡献率以验证特征峰选取的合理性[2]

(1)

通过minitab(见附录)从中提取的2个主成分的特征值,对第一种情况进行分析,见图3,可知其得分向量(具体见附录)和特征根为:

,方差累计贡献率rat_1=0.99785。说明选取的这两个主成分可以解释99.785%的原数据,具有较强的代表性。再对第二种情况进行分析,可以得到特征根和方差贡献率,方差累计贡献率rat_{2}=0.99897

说明选取的这三个主成分可以解释99.897%的原数据,具有更强的代表性。但是对比第一种情况,代表性并提升幅度过小,综合考虑后,最终在问题一中选取第一种情况,即两个主成分为最终解。

图3 主成分分析碎石图

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

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

相关文章

【Leetcode】【240404】1614. Maximum Nesting Depth of the Parentheses

BGM(?):圣堂之门-阿沁《梵谷的左耳》 Description A string is a valid parentheses string (denoted VPS) if it meets one of the following: It is an empty string “”, or a single character not equal to “(” or “)…

【剑指offr--C/C++】JZ7 重建二叉树

一、题目 二、思路及代码 前序遍历:中、左、右。所以前序遍历的第一个节点是树的根节点,第二个节点是左子树的根节点。。。。 中序遍历:左、中、右。树的根节点在中间某处 我们可以根据二者的特点结合一下:对于前序遍历序列{1,2,4…

11-LINUX--信号

一.信号的基本概念 信号是系统响应某个条件而产生的事件,进程接收到信号会执行相应的操作。 与信号有关的系统调用在“signal.h”头文件中有声明 常见信号的值,及对应的功能说明: 信号的值在系统源码中的定义如下: 1. #define …

渗透测试靶机----sec123

渗透测试靶机----sec123 复现打靶,这里先需要搭建好环境 这里还需要将内网中的网站映射到公网中,完全模拟实战渗透测试使用frp轻松实现 这里就搭建好靶机了,准备开始渗透工作 先使用这个网址,扫描看看这里通过扫描,发现这三个端口对应三个网页: http://107.151.243.222:8…

用于无人机小型化设计的高精度温补晶振

用于无人机小型化设计的高精度温补晶振:TG2016SMN和TG2520SMN。无人机的发展可以说是非常的迅速,在安防,农业,交通,电力,直播等领域经常能看到无人机大显身手。无人机的应用场最是非常的广泛,功能更强&…

微软detours代码借鉴点备注

comeasy 借鉴点1 Loadlibray的时间选择 注入库wrotei.dll,为了获取istream的接口,需要loadlibrary,但是在dllmain中是不建议这样做的。因此,动态库在dllmain的时候直接挂载了comeasy.exe的入口 //获取入口 TrueEntryPoint (i…

数据仓库面试总结

文章目录 1.什么是数据仓库?2.ETL是什么?3.数据仓库和数据库的区别(OLTP和OLAP的区别)4.数据仓库和数据集市的区别5.维度分析5.1 什么是维度?5.2什么是指标? 6.什么是数仓建模?7.事实表7.维度表…

SecureCRT通过私钥连接跳板机,再连接到目标服务器

文章目录 1. 配置第一个session(跳板机)2. 设置本地端口3. 设置全局firewall4. 配置第二个session(目标服务器) 服务器那边给了一个私钥,现在需要通过私钥连接跳板机,再连接到目标服务器上 🍉 …

【C++】stack和queue

个人主页 : zxctscl 如有转载请先通知 文章目录 1. stack的介绍和使用1.1 stack的介绍1.2 stack的使用1.3 stack的模拟实现 2. queue的介绍和使用2.1 queue的介绍2.2 queue的使用2.3 queue的模拟实现 3. 容器适配器3.1 概念3.2 STL标准库中stack和queue的底层结构3.…

如何建立一个网页模版

创建一个网页模板涉及的主要步骤如下: 基本HTML结构模板创建: 创建HTML文件: 首先,你需要创建一个新的HTML文件,并在其中编写基本的HTML结构,例如: <!DOCTYPE html> <html lang="zh"> <head><meta charset="UTF-8"

VMware ESXi 6.7

1.浏览器中输入地址&#xff0c;进入管理界面 2.选择 存储 右击浏览&#xff0c;创建新的目录 3.点击 上载 &#xff0c;选择镜像文件 4.等待上载完成 5.点击 虚拟机-新建虚拟机 6.进入新建虚拟机界面 7.进入Windows安装界面 8.安装VMware Tools

python爬虫获取豆瓣前top250的标题(简单)

今天是简略的一篇&#xff0c;简单小实验 import requests from bs4 import BeautifulSoup# 模拟浏览器的构成&#xff08;请求头&#xff09; headers {"User-Agent": "Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Ch…

使用CSS计数器,在目录名称前加上了序号,让目录看起来更加井然有序

目录&#xff08;Text of Contents缩写为TOC&#xff09;&#xff0c;其实就是一篇文章的概要或简述。这好比&#xff0c;去书店买书&#xff0c;先是被这本书的标题所吸引&#xff0c;而后我们才会&#xff0c;翻开这本书目录&#xff0c;看看这本书主要是在讲些什么&#xff…

Restful Web Service

Restful 1.特点 RESTful是一种架构风格&#xff0c;强调简单、轻量级和对资源的状态less操作。RESTful是通过HTTP协议进行通信的。RESTful的应用程序可以调用运行在不同服务器上的服务或函数。RESTful的接口通常使用JSON&#xff0c;但实际上它们都支持多种数据格式。RESTful…

【vue + springboot】切片+断点上传 + 秒传

需要了解透彻请参考&#xff1a;如何实现大文件上传、断点续传、切片上传_断点上传-CSDN博客 说明 切片&#xff1a;根据指定的大小对文件进行切块上传。 断点上传&#xff1a;每次上传完成一个切片后端保存信息&#xff0c;前端每次上传都判断所上传的文件是否存在&#xf…

纪念品分组

题源 贪心做法 &#xff0c;先排序再分组 题目描述 元旦快到了&#xff0c;校学生会让乐乐负责新年晚会的纪念品发放工作。为使得参加晚会的同学所获得 的纪念品价值相对均衡&#xff0c;他要把购来的纪念品根据价格进行分组&#xff0c;但每组最多只能包括两件纪念品&#xff…

1.9 面试经典150题 除自身以外数组的乘积

除自身以外数组的乘积 给你一个整数数组 nums&#xff0c;返回 数组 answer &#xff0c;其中 answer[i] 等于 nums 中除 nums[i] 之外其余各元素的乘积 。 题目数据 保证 数组 nums之中任意元素的全部前缀元素和后缀的乘积都在 32 位 整数范围内。 请 不要使用除法&#xff0…

vue如何搭建项目?

首先我们下载安装router&#xff0c;less / scss&#xff0c;pinia&#xff0c;axios&#xff0c; element / Ant Design Vue&#xff0c;Echarts 以及如何配置别名&#xff0c;自定义组件 npm npm create vitelatest yarn yarn create vite 然后安装依赖 成功后 在src下…

PyTorch搭建Informer实现长序列时间序列预测

目录 I. 前言II. InformerIII. 代码3.1 输入编码3.1.1 Token Embedding3.1.2 Positional Embedding3.1.3 Temporal Embedding 3.2 Encoder与Decoder IV. 实验 I. 前言 前面已经写了很多关于时间序列预测的文章&#xff1a; 深入理解PyTorch中LSTM的输入和输出&#xff08;从i…

python 笔记

文章目录 pdbpdb开始调试pythonpdb设置断点单步执行进入到函数的内部执行到下一个断点或程序结束调用栈查看命令查看当前函数调用堆栈向上一层函数查看调用堆栈查看源代码 importimport 用法 numpy导入numpy模块numpy常用函数np.argmaxnp.sum range生成连续序列生成不连续序列 …