标题: 数学建模算法总结(二) [打印本页] 作者: 杨利霞 时间: 2019-3-24 11:16 标题: 数学建模算法总结(二) 数学建模算法总结(二) ~, f4 z) j% D5 v- z" i7 x% E
第十一章 方差分析 c' [+ ?) W# j* ?! f$ t' I f2 |. b. c! x( Q. Z- B
为了使生产过程稳定,达到优质、高产,需要对影响产品质量的因素进行分析,找出有显著影响的那些因素,除了从机理方面进行研究外,常常要作许多试验,对结果作分析、比较,寻求规律。用数理统计分析试验结果、鉴别各因素对结果影响程度的方法称为方差分析(Analysis Of Variance),记作 ANOVA。 & X) e5 B6 D) U" _' f+ U! E" c人们关心的试验结果称为 指标,试验中需要考察、可以控制的条件称为 因素或因子,因素所处的状态称为 水平。上面提到的灯泡寿命问题是单因素试验,小麦产量问题是双因素试验。处理这些试验结果的统计方法就称为单因素方差分析和双因素方差分析。$ L5 d% i: P+ X' v0 A- F4 Y( D9 z
§1 单因素方差分析 ' y% H/ J2 h/ [$ |. f7 H9 T( _: f只考虑一个因素 A 对所关心的指标的影响, A 取几个水平,在每个水平上作若干个试验,试验过程中除 A 外其它影响指标的因素都保持不变(只有随机因素存在),我们的任务是从试验结果推断,因素 A 对指标有无显著影响,即当 A 取不同水平时指标有无显著差别。; i# U/ h; B5 ?* M
A 取某个水平下的指标视为随机变量,判断 A 取不同水平时指标有无显著差别,相当于检验若干总体的均值是否相等。0 {8 b# R% j1 X# g0 ^: {
5 O1 R l8 w' F5 E' p: v) h§2 双因素方差分析3 p6 X" S9 |( e# H6 l
如果要考虑两个因素 B A, 对指标的影响, B A, 各划分几个水平,对每一个水平组+ o/ k) z9 ]" i2 W9 }
合作若干次试验,对所得数据进行方差分析,检验两因素是否分别对指标有显著影响,; U5 h1 I3 M2 j# l1 w
或者还要进一步检验两因素是否对指标有显著的交互影响。5 }# V$ f2 f ~" v/ P4 `: G
1 Q% b, w; D8 D§3 正交试验设计与方差分析1 @$ b1 b* P- h1 w* _& k
由于因素较少时,我们可以对不同因素的所有可能的水平组合做试验,这叫做全面试验。当因素较多时,虽然理论上仍可采用前面的方法进行全面试验后再做相应的方差分析,但是在实际中有时会遇到试验次数太多的问题。如果考虑更多的因素及水平,则全面试验的次数可能会大得惊人。因此在实际应用中,对于多因素做全面试验是不现实的。于是我们考虑是否可以选择其中一部分组合进行试验,这就要用到试验设计方法选择合理的试验方案,使得试验次数不多,但也能得到比较满意的结果。 9 s, a1 u) q" ]6 y& M+ g ! d( x9 ^5 p. }; @$ R' r6 I6 u # m( r2 s/ }/ e8 b0 d X; q9 X3 ]3 v% V3 h. {2 F" ^" V
第十二章 回归分析 ! n4 w9 g8 Y, x- `, k8 C7 f3 _8 o$ H& U& Y8 m; j! c# l
曲线拟合问题的特点是,根据得到的若干有关变量的一组数据,寻找因变量与(一个或几个)自变量之间的一个函数,使这个函数对那组数据拟合得最好。通常,函数的形式可以由经验、先验知识或对数据的直观观察决定,要作的工作是由数据用最小二乘法计算函数中的待定系数。从计算的角度看,问题似乎已经完全解决了,还有进一步研究的必要吗?) Y6 n, u' @% k( u4 k
从数理统计的观点看,这里涉及的都是随机变量,我们根据一个样本计算出的那些系数,只是它们的一个(点)估计,应该对它们作区间估计或假设检验,如果置信区间太大,甚至包含了零点,那么系数的估计值是没有多大意义的。另外也可以用方差分析方法对模型的误差进行分析,对拟合的优劣给出评价。简单地说,回归分析就是对拟合问题作的统计分析。 - Q* |$ a9 M, z4 T& r具体地说,回归分析在一组数据的基础上研究这样几个问题:! T% i1 `! L, p- _3 H5 R
(i)建立因变量 y 与自变量 x1,x2,……,xm之间的回归模型(经验公式);" W4 y3 C, Y+ ^. u0 d+ ?/ C
(ii)对回归模型的可信度进行检验;8 f. r3 T, t S3 M
(iii)判断每个自变量xi=(i=1,2,……,m)对 y 的影响是否显著;( h5 \) u& ^& O3 s9 ^4 c
(iv)诊断回归模型是否适合这组数据;2 [/ B2 ?3 k2 p
(v)利用回归模型对 y 进行预报或控制。. a+ j5 O$ M) a' j! T# m' [% E; F Q
8 H$ v0 R" P% D / P/ |$ Y0 i# R K6 s8 C* b+ T- k3 A1 h. K- L4 k* N# K2 `5 m( |
第十三章 微分方程建模, _5 z$ N7 l" S2 u
: y7 ^7 i6 m' r+ |0 j& m7 G4 l6 ~
微分方程建模是数学建模的重要方法,因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题,大体上可以按以, \ B, s& l" d6 K
下几步: 2 s4 c- R' l/ L4 R: H! m7 l/ N1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。+ T% z7 z! @2 r* Z ~
2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。9 E. h2 i% O* f& l
3. 运用这些规律列出方程和定解条件。+ \6 J; b& r8 Q3 T6 \& ?$ h$ c2 h
列方程常见的方法有: $ f' o& Q% H. w+ @' O(i)按规律直接列方程 ) Q% y* v# j. z# j/ k5 I在数学、力学、物理、化学等学科中许多自然现象所满足的规律已为人们所熟悉,并直接由微分方程所描述。如牛顿第二定律、放射性物质的放射性规律等。我们常利用这些规律对某些实际问题列出微分方程。 . @3 W/ d+ E. K3 X, M, N(ii)微元分析法与任意区域上取积分的方法 7 u5 a: p: ^1 ~* E4 t自然界中也有许多现象所满足的规律是通过变量的微元之间的关系式来表达的。对于这类问题,我们不能直接列出自变量和未知函数及其变化率之间的关系式,而是通过微元分析法,利用已知的规律建立一些变量(自变量与未知函数)的微元之间的关系式,然后再通过取极限的方法得到微分方程,或等价地通过任意区域上取积分的方法来建立微分方程。 / u, y8 S0 U" N" }(iii)模拟近似法 9 _1 x; f0 c; r. l0 w h在生物、经济等学科中,许多现象所满足的规律并不很清楚而且相当复杂,因而需要根据实际资料或大量的实验数据,提出各种假设。在一定的假设下,给出实际现象所满足的规律,然后利用适当的数学方法列出微分方程。在实际的微分方程建模过程中,也往往是上述方法的综合应用。不论应用哪种方法,通常要根据实际情况,作出一定的假设与简化,并要把模型的理论或计算结果与实际情况进行对照验证,以修改模型使之更准确地描述实际问题并进而达到预测预报的目的。3 a9 @( W9 ^1 Y
1 K+ }) @8 |/ w! ?
2 X" [' c+ w( Y
* j. P1 d) V7 X8 A; k第十四章 稳定状态模型 - o' W) b- x8 k4 Q* l q/ q$ O/ c4 J# J; b
虽然动态过程的变化规律一般要用微分方程建立的动态模型来描述,但是对于某些实际问题,建模的主要目的并不是要寻求动态过程每个瞬时的性态,而是研究某种意义下稳定状态的特征,特别是当时间充分长以后动态过程的变化趋势。譬如在什么情况下描述过程的变量会越来越接近某些确定的数值,在什么情况下又会越来越远离这些数值而导致过程不稳定。为了分析这种稳定与不稳定的规律常常不需要求解微分方程,而可以利用微分方程稳定性理论,直接研究平衡状态的稳定性就行了。 " z6 ?9 ^8 {" y& h2 f ; `, h( i' O8 U- V. J) U$ v- Y1 s. l1 ~
) u# Q9 f8 ^- l6 V7 }
第十五章 常微分方程的解法 0 C; a! A' z$ U - ^) T. t: d) p7 I; S建立微分方程只是解决问题的第一步,通常需要求出方程的解来说明实际现象,并加以检验。如果能得到解析形式的解固然是便于分析和应用的,但是我们知道,只有线性常系数微分方程,并且自由项是某些特殊类型的函数时,才可以肯定得到这样的解,而绝大多数变系数方程、非线性方程都是所谓“解不出来”的,即使看起来非常简单的方程,于是对于用微分方程解决实际问题来说,数值解法就是一个十分重要的手段。 1 c* r+ |7 Q- ?; j% c 1 i2 `. t( B1 r1 P( @$ L; k欧拉(Euler)方法、龙格—库塔(Runge—Kutta)方法、线性多步法、一阶微分方程组与高阶微分方程的数值解法/ K! P1 ]) [4 F$ S' N