Polymath软件的应用题材.pptVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Polymath软件的应用 Polymath软件简介 科学与工程中数学问题的数值求解工具 Polymath软件简介 辅助功能 高级计算器 单位换算 常用物理化学常数查询 结果输出到Excel和Matlab软件 Polymath软件界面 Polymath软件使用方法 应用举例 线性方程组求解 给定一组数据(xi,yi),i 1, 2 , …, m ,做拟合直线p x a + bx , 均方误差为 : 整理得到拟合曲线满足的方程: 练习一 下表为实验测得的某一物性和温度之间的关系数据,表中x为温度数据,y为物性数据。请用线性函数拟合温度和物性之间的关系。 解:设拟合直线 ,并计算得下表: 应用举例 非线性方程组求解 练习二 用PLOYMATH求解下列非线性方程组: 农场主的财产有500只兔子(x)和200只狐狸(y)。用POLYMATH画一条最多到500d的期间内作为时间函数的狐狸和兔子浓度的曲线。 问题的提出 化工设计及化工模拟计算中,有大量的物性参数及各种设备参数。实验测量得到的常常是一组离散数据序列(xi ,yi) 图1-1所示为“噪声” 图1-2所示为无法同时满足某特定的函数 问题的提出 在化学化工中,许多模型也要利用数据拟合技术,求出最佳的模型和模型参数。 如在某一反应工程实验中,我们测得了如下表所示的实验数据: 问题的提出 确定在其他条件不变的情况下,转化率y和温度T的具体关系,现拟用两种模型去拟合实验数据,两种模型分别是: 拟合的标准 向量Q与Y之间的误差或距离有以下几种定义方法: (1)用各点误差绝对值的和表示 (2)用各点误差按绝对值的最大值表示 (3)用各点误差的平方和表示 均方误差 拟合的标准 由于计算均方误差的最小值的原则容易实现而被广泛采用。按均方误差达到极小构造拟合曲线的方法称为最小二乘法。 拟合的标准 实例 实验测得二甲醇(DME)的饱和蒸汽压和温度的关系如下表 : 拟合的标准 实例 如果以直线拟合p a+bt,即拟合函数是一条直线。通过计算均方误差Q a , b 最小值而确定直线方程(见图1-3) 拟合的标准 实例 如果采用二次拟合,通过计算下述均方误差: 拟合的标准 实例 比较图1-3和图1-4以及各自的相关系数和方差可知: 对于DME饱和蒸汽压和温度之间的关系,在实验温度范围内用二次拟合曲线优于线性拟合。 二次拟合曲线具有局限性,由图1-4观察可知,当温度低于-30℃时,饱和压力有升高的趋势,但在拟合的温度范围内,二次拟合的平均绝对偏差又小于一次拟合,故对物性数据进行拟合时,不仅要看在拟合条件下的拟合效果,还必须根据物性的具体性质,判断在拟合条件之外的物性变化趋势,以便使拟合公式在已做实验点数据之外应用。 * * 主要功能 线性方程组(LEQ) 非线性方程组(NLE) 微分方程组(DEQ) 数据拟合(REG) 输出形式:报告、图形和表格 线性方程组 非线性方程组 微分方程组 数据拟合 输入需要求解的方程 输入其中参数的计算方法及参数值 输入自变量的初始值、终值和计算步长或者函数值的初始猜值 在工具栏中选定求解方法、输出格式(图形或表格) 运行程序,获得计算结果 输入窗口 运行结果 方程数 运行按钮 (1) Q a , b 的极小值需满足: 线性方程组求解 或 2 称式 2 为拟合曲线的法方程。 线性方程组求解 可用消元法或克莱姆方法解出方程: 线性方程组求解 34 23 17 24 21 21 14 22 14 y 130 100 72 70 64 62 60 55 42 x 13 17 16 14 12 13 13 12 12 11 10 11 y 40 36 31 30 29 25 23 22 21 16 15 13 x 121 100 121 144 144 … 1156 61640 143 150 176 252 264 … 4420 18913 11 10 11 12 12 … 34 344 13 15 16 21 22 … 130 956 1 2 3 4 5 … 21 Σ x2 xy y x 编号 将数据代入法方程组(1-12)中,得到: 解方程得:a 8.2084 , b 0.1795 。 拟合直线为: 线性方程组求解 线性方程组求解 输入方程 输入初值 运行结果 微分方程组求解 微分方程组求解 练习二 微分方程组求解 原程序 运行报告 图形输出 表格输出 图1-1 含有噪声的数据 图1-2 无法同时满足某特定函数的数据序列 数据拟合与分析 1-2 1-3 1-4 1-5 1-6 序号 温度 ℃ 蒸气压 MPa 1 -23.7 0.101 2 -10 0.174

文档评论(0)

知识宝库 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档