《MATLAB论文.docx》由会员分享，可在线阅读，更多相关《MATLAB论文.docx（16页珍藏版）》请在第壹文秘上搜索。
1、MAT1.AB在数学建模中的作用姓名：冯文俊【摘要】学号：201964100129通过对实际问题的抽象和简化，引入一些数学符号、变量和参数,运用某些规律,用数学语言和数学方法建立变量、参数间的内在联系，得出一个数学结构，该数学结构是实现的一个近似刻画，称之为数，学模型。建立和求解数学模型的全过程就是数学建模，它包括模型的建立、求解、分析、检验循环来回的全过程，MAT1.AB语言正是处理此类问题的很好工具，既能进行数值求解，乂能绘制有关曲线，特别便利好用。1.数学建模的概念定义：数学模型，就是运用科学抽象法，把困难的探讨对象转化为数学问题，经合理简化后，建立起揭示探讨对象定量的规律性的数学关系式
2、（或方程式）。这既是数学方法中最关键的一步，也是最困难的一步；数学建模是指对现实世界的一特定对象，为了某特定目的,做出一些重要的简化和假设，运用适当的数学工具得到一个数学结构，用它来说明特定现象的现实性态，预料对象的将来状况，供应处理对象的优化决策和限制，设计满意某种须要的产品等。2.数学建模的一般步骤第一步：依据探讨对象的特点，确定探讨对象属哪类自然事物或自然现象，从而确定运用何种数学方法及建立何种数学模型。即首先确定对象及应当运用的数学模型的类别归属问题，是屈于“必定”类，还是“随机”类：是“突变”类，还是“模糊”类。其次步：确定几个基本量和基本的科学概念，用以反映探讨对象的状态。这须要依
3、据已有的科学理论或假说及试验信息资料的分析确定。例如在力学系统的探讨中，首先确定的摹本物理量是质主(m)、速度(v)、加速度（八）、时间(t)、位矢(r)等。必需留意确定的基本量不能过多，否则未知数过多，难以简化成可能数学模型，因此必需诜择出实质性、关键性物理量才行。第三步：抓住主要冲突进行科学抽象。现实探讨对象是困难的，多种因素混在一起，因此，必需变困难的探讨对象为简洁和志向化的探讨对象，做到这一点相当困难，关键是分清主次。如何分清主次只能详细问题详细分析，但也有两条基本原则：一是所建数学模型肯定是可能的，至少可给出近似解：二是近似解的误差不能超过实际问题所允许的误差范围。第四步：对简化后的
4、基本量进行标定，给出它们的科学内涵。即标明哪些是常量，哪些是已知量，哪些是待求量，哪些是矢量，哪些是标量，这些量的物理含义是什么？第五步：按数学模型求出结果。第六步：验证数学模型。验证时可依据状况对模型进行修正，使其符合程度更高，当然这以求原模型及实际状况基本相符为原则。筒而言之I)明确问题进行合理的假设；一般模型假设遵从以下原则：目的性原则:从原型中抽象出及建模目的有关的因素，简化掉无关的因素或关系不大的因素。简明性原则:所给的假设条件要简洁,精确,有利于构造模型。真实性原则：设条款要符合情理，简化带来的误差应满足实际问题所允许的范围内。全面性原则：在对事物原型本身作出的假设的同时，还要给出
5、原型所处的环境条件。3）构造模型4）模型求解5）模型的检验及修正3 .数学建模的一般方法1）机理分析法从基本物理定律以及系统的结构数据来推导出模型。比例分析法一建立变量之间函数关系的最基本最常用的方法。代数方法一求解离散问题（离散的数据、符号、图形）的主要方法。逻辑方法一是数学理论探讨的重要方法，对社会学和经济学等域的实际问题，在决策，对策等学科中得到广泛应用。常微分方程一解决两个变量之间的改变规律，关键是建立”瞬时改变率”的表达式。偏微分方程一解决因变量及两个以上自变量之间的改变规律。2)数据分析法从大量的观测数据利用统计方法建立数学模型。回来分析法一用于对函数f(x)的一组观测值(xi,f
6、i)i=l,2,n,确定函数的表达式，由于处理的是静态的独立数据，故称为数理统计方法。时序分析法一处理的是动态的相关数据，又称为过程统计方法。回来分析法一用于对函数f(x)的一组观测值(xi,fi)i=l,2,on,确定函数的表达式，由于处理的是静态的独立数据，故称为数理统计方法。时序分析法一处理的是动态的相关数据，乂称为过程统计方法。3)仿真和其他方法计算机仿真(模拟)一实质上是统计估计方法,等效于抽样试验。1.离散系统仿真一有一组状态变量。2.连续系统仿真一有解析表达式或系统结构图。因子试验法一在系统上作局部试验,再依据试验结果进行不断分析修改,求得所需的模型结构。人工现实法一基于对系统过
7、去行为的了解和对将来希望达到的目标,并考虑到系统有关因素的可能改变,人为地组成一个系统。4 .应用MAT1.AB进行数学建模探讨数学是在实际应用的需求中产生的，我们把遇到的实际问题进行分析，发觉其中的可以用数学语言来描述的关系或规律，把这个实际问题转化成一个数学问题，建立了数学模型！。但数学模型迫切须要一个便利、快捷且功能强大的工具去实现并解决，特殊是随着科技的进步，人们在解决问题的时候经常要用到很多比较困难的数学学问和大量的数据计算，这无疑加大了人们解决问题的难度，也要耗费更长的时间。而VAT1.BA正是在数学计算和大量数据处理方面具备其它软件所不具备的优势，且操作简洁，运算速度快，所以应用
8、MAT1.BA进行数学建模也就大大提高了人们的效率。而且MAT1.BA还有很强的绘图功能,这就可以使得模型图象化，使得探讨人员对建模成果的优劣一目了然,简洁进行修正及改进。5 .简洁数模问题中MAT1.AB的应用例1：某公司2000年各月生产总值(单位：万元)分别为：22608895562391014815623,试显示折线图表示该厂生产总值的改变状况。在MAT1.AB吩咐窗口中键入：p=22,60,88,95,56,23,9,10,14,81,56,23;plot(p)；可得：如图：RlFigure1问2:依据表1数据作一个多子图。要求：第一个图各类网井产油量及年份曲线图，其次个图为0105
9、年各类网井产油量的对比直方图，第三个图为03年各类网井产油量的饼图，第四个图为七五井和卜五井产油量年份的双座标图。给如卜数据：各类井网在近几年的产油量(万吨)年份：201920192019200020192019201920192019MAT1.AB语言来对此例题做以下解析：figure(,position*,50f50,800,650)%在图形窗口左卜角建立横纵坐标都为50的，宽度800,高度650的窗11=2019:2019;%产生行向量tlt2=2019:2019;%产生行向量t2yl=500.6442.4428.6370.1343.1;%一个行矩阵yiy2=354.7318.0280.
10、7246.6229.0;与建立一个行矩阵y2y3=197.4297412.8547.0579.8547.5527.0492.3437.0；%y3y4=72.3218.2297.1416.1508.7居建立一个行矩阵y4SUbPlOt(2,2,1)%2*2个区中的1号区plot(t2,yl,t2,y2,thy3,t2,y4);%绘制二维图像title(各类网井产油量及年份曲线图)N标题为各类网井产油量及年份曲线图Iegend(七五井，八五井，九五井，十五井，1)；%列出图标12=2019:2019;%产生行向量t2y2=500.6354.7579.872.3442.4318.0547.5218.
11、2428.6280.7527.0297.1370.1246.6492.3437.0343.1229.0437.0508.4;与建立一个行矩阵y2SUbPlOt(2,2,2)%2*2个区中的2号区bar(t2,y2)title(,005年各类网井产油量的对比直方图);%标题为005年各类网井产油量的对比直方图IegendC七五井八五井，九五井，十五井，D；%列出图t3=343.1229.0437.0508.7产生行向量t3SUbPIOt(2,2,3)%2*2个区中的3号区pie(t3):title(103年各类网井产油量的饼图)；$标题03年各类网井产油量的饼图IegendC七五井，八五井，九五
12、井，十五井，1)；%列出图标t4=2019:1:2019;%产生行向量t4yl=354.7318.0280.7246.6229.0;y2=72.3218.2297.1416.1508.7;舟建立一个行矩阵ylsubplot(2,2,4);%2*2个区中的4号区plotyy(t4,yl,U,y2);先绘制二维图像title(,七五井和十五井产油量及年份的双座标图)港标题七五井和十五井产油量及年份的双座标图IegendC七五井，十五井)；先列出图标得如1.线性系统的模型建立G(三)=Y(三)U(三)=(6s3+12s2+6+10)(s4+2s3+3s2+s+l)可以将其用下列MAT1.AB语句来实
13、现：num=612610：den=l2311：printsys(num,den)执行结果为num/den=6s3+12s2+6s+10s4+2s3+3s-2+s+12.若给定的系统的传递函数为：G(三)=4(s+2)(s2+6s+6)(s(s+l)3(s3+3s2+2s+5)解：num=4*conv(12,166)num=4327248den=conv(l0,conv(l1,conv(l1,1325)den=159121250关于从状态空间形式到零极点形式的转换MT1.B函数ss2zp()的调用格式为：Z,P,K=ss2zp(A,B.C,I),iu)其中，A,B,CD为状态空间形式的个系数矩阵
14、：iu为输入的代号。对于单变量系统iu=l,对于多变量系统iu表示要求的输入序号；返回量矩阵P储存传递函数的极点；而零点储存在矩阵Z中，Z的列数等于输出y的维数，每列对应一个输出：对应增益则在列向量K中。如下示例：对于以给出状态空间表达式，是依据以上函数求取系统的传递函数矩阵。X(t)=O01J-3/2-2-l/2;-30-4x(t)+l1;-1-1;-1-3u(t)Y(t)=l00:010x(t)此系统可由以下MAT1.AB语句唯一地表示出来A=001:-3/2-2-l/2;-30-4B=11：-11；-1-3C=100;010,D=zeros(2,2)试验过程如下A=001:-3/2-2-l2j-30-4A=1.0000-1.5000-2.0000-0.5000-3.00000-4.0000B=l1;-1-1:-1-3B=C=lO0,D=zeros(2,2)再运行:Z1,P1,Kl=ss2zp(,B,C,D,1),Z2,P2,K2=ss2zp(,B,C,D12)窗口中出现Zl=-3.OOOO-2.OOOO-2.OOOO-3.OOOO-1.0000-3.OOOO-2.OOOO-1.0000P2再看一个将此多变量系统求其传递函数矩阵的例子还是上一个例子的系统解利用以下MAT1.AB语句numl,denl=ss2tf(A,B,C,D,1),num2,den2=ss2tf(A,
