Polymath软件应用课件

Polymath软件的应用2013年Polymath软件简介科学与工程中数学问题的数值求解工具主要功能线性方程组（LEQ）非线性方程组（NLE）微分方程组（ODE）数据拟合（REG）输出形式：报告、图形和表格Polymath软件简介辅助功能高级计算器单位换算常用物理化学常数查询结果输出到Excel和Matlab软件Polymath软件使用方法输入需要求解的方程输入其中参数的计算方法及参数值输入自变量的初始值、终值和计算步长或者函数值的初始猜值在工具栏中选定求解方法、输出格式（图形或表格）运行程序，获得计算结果应用举例线性方程组求解输入窗口运行结果方程数运行按钮给定一组数据（xi,yi）,i=1,2,…,m，做拟合直线p(x)=a+bx,均方误差为：（1）

Q(a,b)的极小值需满足：

线性方程组求解可用消元法或克莱姆方法解出方程：线性方程组求解下表为实验测得的某一物性和温度之间的关系数据，表中x为温度数据，y为物性数据。请用线性函数拟合温度和物性之间的关系。x131516212223252930313640y111011121213131214161713x42556062647072100130y142214212124172334练习一解：设拟合直线

，并计算得下表：编号xyxyx212345…21Σ1315162122…1309561110111212…34344143150176252264…442018913121100121144144…115661640将数据代入法方程组（1-12）中，得到：

解方程得：a=8.2084,b=0.1795。拟合直线为：

线性方程组求解应用举例非线性方程组求解输入方程输入初值运行结果练习二用PLOYMATH求解下列非线性方程组：农场主的财产有500只兔子（x）和200只狐狸（y）。用POLYMATH画一条最多到500d的期间内作为时间函数的狐狸和兔子浓度的曲线。微分方程组求解微分方程组求解练习二微分方程组求解原程序运行报告图形输出表格输出问题的提出化工设计及化工模拟计算中，有大量的物性参数及各种设备参数。实验测量得到的常常是一组离散数据序列（xi,yi）图1-1所示为“噪声”图1-2所示为无法同时满足某特定的函数图1-1含有噪声的数据图1-2无法同时满足某特定函数的数据序列数据拟合与分析问题的提出确定在其他条件不变的情况下，转化率y和温度T的具体关系，现拟用两种模型去拟合实验数据，两种模型分别是：

(1-2)

(1-3)

拟合的标准

向量Q与Y之间的误差或距离有以下几种定义方法：（1）用各点误差绝对值的和表示（2）用各点误差按绝对值的最大值表示（3）用各点误差的平方和表示(均方误差)(1-4)

(1-5)

(1-6)

拟合的标准由于计算均方误差的最小值的原则容易实现而被广泛采用。按均方误差达到极小构造拟合曲线的方法称为最小二乘法。拟合的标准实例如果以直线拟合p=a+bt,即拟合函数是一条直线。通过计算均方误差Q(a,b)最小值而确定直线方程（见图1-3）图1-3DME饱和蒸汽压和温度之间的线性拟合拟合得到得直线方程为：方差为0.0050081。

(1-8)

(1-7)

拟合的标准实例如果采用二次拟合，通过计算下述均方误差：拟合得二次方程为：（1-9）（1-10）方差为0.