[理学]数学建模讲座 2.ppt
《[理学]数学建模讲座 2.ppt》由会员分享,可在线阅读,更多相关《[理学]数学建模讲座 2.ppt(131页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、理学理学数学建模讲座数学建模讲座 2 建模中常用的软件建模中常用的软件 Matlab简介简介 Matlab使用使用 Matlab画图画图 Matlab工具箱工具箱提纲建模中常用的软件 Lindo&Lingo LINGO是Linear Interactive and General Optimizer的缩写,即“交互式的线性和通用优化求解器”,由美国LINDO系统公司(Lindo System Inc.)推出的,可以用于求解非线性规划,也可以用于一些线性和非线性方程组的求解等, SPSS SPSS(Statistical Product and Service Solutions),“统计产品与
2、服务解决方案”软件,是一款优秀的统计分析软件。 MATLAB 美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵运算、科学数据可视化以及非线性动态系统的建模和仿真等强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言的编辑模式。 Mathematica、Maple 优秀的数学计算软件 Matlab简介 MATLAB源于MATrix LABoratory一词,原意为矩阵实验室,一开始它是一种专门用于矩阵数值计算的软件。之后
3、发展了数据可视化功能,符号运算功能,图形处理功能以及目前具有的数十个工具箱。 Matlab可作为线性代数, 数理统计, 自动控制, 数字信号处理, 模拟与数字通信, 时间序列分析, 动态系统仿真, 神经网络等课程的学习工具 MATLAB语言主要特点 1、功能强大 Matlab可以很方便地处理线性代数中的矩阵计算, 方程组的求解, 微积分运算, 多项式运算, 偏微分方程求解, 统计与优化等问题.并具有完善的图形可视化功能。Matlab简介2、语言简单 MATLAB语言以向量和矩阵为基本的数据单元, 包括流程控制语句(顺序, 选择, 循环, 条件, 转移和暂停等), 大量的运算符, 丰富的函数,
4、多种数据结构, 输入输出以及面向对象编程. 这些既可以满足简单问题的计算, 也适合于开发复杂的大型程序. Matlab简介3、 扩充能力强, 可开发性强 MATLAB是一个开放的系统, 用户可以方便地看见函数的源程序, 也可以很方便地开发自己的程序, 甚至可以创建自己的工具箱, 以解决本领域内常见的计算问题. 也可以方便地与FORTRAN, C等语言接口, 以充分利用各种资源. Matlab简介 4. 编程容易, 效率高 从形式上看, MATLAB程序文件是一个纯文本文件, 扩展名为m. 用任意字处理软件都可以对它进行编写和修改, 因此程序易调试, 人机交互性强. Matlab简介MATLAB
5、使用环境MATLAB使用环境MATLAB使用环境MATLAB使用运算与常用函数1、关系运算符:小于 :大于 =:大于等于=:不等于2、逻辑运算符&:与运算 |:或运算 :非运算3、for循环、while循环4、break语句:跳出上述循环5、定义函数: function y=new_fun(x)MATLAB使用运算与常用函数MATLAB使用运算与常用函数常用的基本数学函数(1)abs(x):纯量的绝对值或向量的长度sqrt(x):开平方round(x):四舍五入至最近整数fix(x):无论正负,舍去小数至最近整数floor(x):地板函数,即舍去正小数至最近整数ceil(x):天花板函数,即加
6、入正小数至最近整数rat(x):将实数x化为分数表示rats(x):将实数x化为多项分数展开sign(x):符号函数 (Signum function)。MATLAB使用运算与常用函数sin(x):正弦函数cos(x):馀弦函数tan(x):正切函数asin(x):反正弦函数acos(x):反馀弦函数atan(x):反正切函数atan2(x,y):四象限的反正切函数sinh(x):超越正弦函数cosh(x):超越馀弦函数tanh(x):超越正切函数asinh(x):反超越正弦函数acosh(x):反超越馀弦函数atanh(x):反超越正切函数min(x): 向量x的元素的最小值max(x):
7、向量x的元素的最大值mean(x): 向量x的元素的平均值median(x): 向量x的元素的中位数std(x): 向量x的元素的标准差diff(x): 向量x的相邻元素的差sort(x): 对向量x的元素进行排序length(x): 向量x的元素个数norm(x): 向量x的欧氏长度sum(x): 向量x的元素总和prod(x): 向量x的元素总乘积cumsum(x): 向量x的累计元素总和cumprod(x): 向量x的累计元素总乘积dot(x, y): 向量x和y的内 积常用的基本数学函数(2)MATLAB使用运算与常用函数矩阵函数det(A): 行列式计算 A/ : 矩阵的转置inv(
8、A): 矩阵的逆 orth(A): 正交化poly(A): 特征多项式 rank(A): 矩阵的秩trace(A): 矩阵的迹 zeros(m,n): m行n列的零矩阵ones(m,n): m行n列的全1矩阵eys(n): n阶单位矩阵d=eig(A), v,d=eig(A): 特征值与特征向量rand(m,n): m行n列均匀分布随机数矩阵randn(m,n): m行n列正态分布随机数矩阵MATLAB使用图形功能二维图形功能plot: x轴和y轴均为线性刻度注:1、图形颜色:y(黄) k(黑) w(白) b(蓝) g(绿) r(红)2、线条形状:-(实线) -.(点虚线) -(虚线)3、数轴
9、范围调整: axis(xmin,xmax,ymin,ymax)3、图形的各种注解与处理:x轴注解 :xlabel(Input Value)y轴注解 :ylabel(Function Value);图形标题 :title(Two Trigonometric Functions)图形注解 :legend(y = sin(x),y = cos(x)显示格线:grid on4、画出数个小图形于同一个视窗之中: subplotMATLAB使用图形功能bar 长条图 errorbar 图形加上误差范围fplot 较精确的函数图形 polar 极座标图hist 累计图 rose 极座标累计图stairs 阶
10、梯图 stem 针状图fill 实心图feather 羽毛图compass 罗盘图 quiver 向量场图更多的二维作图命令matlab 画图(一)plot 最基本的二维图形指令plot的功能: plot命令自动打开一个图形窗口Figure 用直线连接相邻两数据点来绘制图形 根据图形坐标大小自动缩扩坐标轴,将数据标尺及单位标注自动加到两个坐标轴上,可自定坐标轴,可把x, y 轴用对数坐标表示matlab 画图二维作图 如果已经存在一个图形窗口,plot命令则清除当前图形,绘制新图形 可单窗口单曲线绘图;可单窗口多曲线绘图;可单窗口多曲线分图绘图;可多窗口绘图 可任意设定曲线颜色和线型 可给图形
11、加坐标网线和图形加注功能plot的调用格式 plot(x) plot(x,y) plot(x1,y1,x2,y2) plot(x,y,s) S的标准设定值如下: 字母 颜色 标点 线型 y 黄色 点线 m 粉红 圈线 c 亮蓝 线 r 大红 字线 g 绿色 实线 b 蓝色 星形线 w 白色 : 虚线 k 黑色 (-) 点划线 v 下三角 上三角 右三角 square 正方形 diamond 菱形 pentagram 五角星 hexagram 六角星1. 单窗口单曲线绘图figure,plot(x) 01234567-1-0 .8-0 .6-0 .4-0 .200 .20 .40 .60 .81
12、2. 单窗口多曲线绘图例2:t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y,t,y1,t,y2)yy1y23. 单窗口多曲线分图绘图subplot 子图分割命令调用格式: subplot(m,n,p) 按从左至右, 从上至下排列 行列绘图序号subplot(1,3,1); plot(t,y1)subplot(1,3,2); plot(t,y2)subplot(1,3,3); plot(t,y3)051 0-1-0 .8-0 .6-0 .4-0 .200 .20 .40 .60 .81051 0-1-0 .8-0 .6
13、-0 .4-0 .200 .20 .40 .60 .81051 0-1-0 .8-0 .6-0 .4-0 .200 .20 .40 .60 .81subplot(3,1,1);plot(t,y)subplot(3,1,2);plot(t,y3)subplot(3,1,3);plot(t,y2)01234567-10101234567-10101234567-1014. 多窗口绘图 figure(n) 创建窗口函数,n为窗 口顺序号。t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y) 自动出现第一个窗口figure(2)
14、plot(t,y1) 在第二窗口绘图figure(3)plot(t,y2) 在第三窗口绘图0 1 2 3 4 5 6 7- 1- 0 .8- 0 .6- 0 .4- 0 .200 .20 .40 .60 .810 1 2 3 4 5 6 7- 1- 0 .8- 0 .6- 0 .4- 0 .200 .20 .40 .60 .810 1 2 3 4 5 6 7- 1- 0 .8- 0 .6- 0 .4- 0 .200 .20 .40 .60 .815.可任意设置颜色与线型例 4(例2):plot(t,y,r-,t,y1,g:,t,y2,b*)01234567-1-0 .8-0 .6-0 .4-0
15、 .200 .20 .40 .60 .816.图形加注功能将标题、坐标轴标记、网格线及文字注释加注到图形上,这些函数为: title 给图形加标题 xlable 给x轴加标注 ylable 给y轴加标注 text 在图形指定位置加标注 gtext 将标注加到图形任意位置 grid on(off) 打开、关闭坐标网格线 legend 添加图例 axis 控制坐标轴的刻度例:t=0:0.1:10y1=sin(t);y2=cos(t);plot(t,y1,r,t,y2,b-);x=1.7*pi;1.6*pi;y=-0.3;0.8;s=sin(t);cos(t);text(x,y,s);title(正
16、弦和余弦曲线);legend(正弦,余弦)xlabel(时间t),ylabel(正弦、余弦)gridaxis square7.fplot 绘制函数图函数fplot的调用格式:fplot(fun,lims) 绘制函数fun在x区间 lims=xmin xmax的函数图。fplot(fun,lims,corline) 以指定线形绘图。x,y=fplot(fun,lims) 只返回绘图点的值,而不绘图。用plot(x,y)来绘图。fplot(sin(x),tan(x),cos(x),2*pi*-1 1 -1 1)fplot(humps,0 1,rp)-6-4-20246-6-4-2024600.10
17、.20.30.40.50.60.70.80.9101020304050607080901008.ezplot 符号函数的简易绘图函数ezplot的调用格式:ezplot(f) 这里f为包含单个符号变量x的符号表达式,在x轴的默认范围 -2*pi 2*pi内绘制f(x)的函数图ezplot(f,xmin,xmax) 给定区间ezplot(f,xmin,xmax,figure(n) 指定绘图窗口绘图。ezplot(sin(x)-6-4-20246-1-0.500.51xsin(x)(二) fill 基本二维绘图函数 fill的功能: 绘制二维多边形并填充颜色例:x=1 2 3 4 5;y=4 1
18、5 1 4;fill(x,y,r) (三)特殊二维绘图函数 bar 绘制直方图 polar 绘制极坐标图 hist 绘制统计直方图 stairs 绘制阶梯图 stem 绘制火柴杆图 rose 绘制统计扇形图 comet 绘制彗星曲线 errorbar 绘制误差棒图 compass 复数向量图(罗盘图) feather 复数向量投影图(羽毛图) quiver 向量场图 area 区域图 pie 饼图 convhull 凸壳图 scatter 离散点图例,绘制阶梯曲线x=0:pi/20:2*pi;y=sin(x);stairs(x,y)01234567-1-0.8-0.6-0.4-0.200.20
19、.40.60.81例:阶梯绘图h2=1 1;1 -1;h4=h2 h2;h2 -h2;h8=h4 h4;h4 -h4;t=1:8;subplot(8,1,1);stairs(t,h8(1,:);axis(off)subplot(8,1,2);stairs(t,h8(2,:);axis(off)subplot(8,1,3);stairs(t,h8(3,:);axis(off)subplot(8,1,4);stairs(t,h8(4,:);axis(off)subplot(8,1,5);stairs(t,h8(5,:);axis(off)subplot(8,1,6);stairs(t,h8(6,:
20、);axis(off)subplot(8,1,7);stairs(t,h8(7,:);axis(off)subplot(8,1,8);stairs(t,h8(8,:);axis(off)h2=1 1;1 -1;h4=h2 h2;h2 -h2;h8=h4 h4;h4 -h4;t=1:8;for i=1:8subplot(8,1,i);stairs(t,h8(i,:)axis(off)end例:绘制极坐标绘图t=0:2*pi/90:2*pi;y=cos(4*t);polar(t,y) 0.2 0.4 0.6 0.8 13021060240902701203001503301800例:绘制火柴杆绘图
21、 t=0:0.2:2*pi; y=cos(t); stem(y)05101520253035-1-0.8-0.6-0.4-0.200.20.40.60.81例:绘制直方图 t=0:0.2:2*pi; y=cos(t); bar(y)05101520253035-1-0.8-0.6-0.4-0.200.20.40.60.81例:绘制彗星曲线图 t= -pi:pi/500:pi; y=tan(sin(t)-sin(tan(t); comet(t,y)-3-2-10123-2.5-2-1.5-1-0.500.511.522.5x=magic(6);area(x)11.522.533.544.555.
22、56020406080100120 x=1 2 3 4 5 6 7;y=0 0 0 1 0 0 0;pie(x,y)4%7%11%14%18%21%25%pie(x,y,North,South,East,West,middle,fa,white)NorthSouthEastWestmiddlefawhiteload seamountscatter(x,y,50,z)210.8210.9211211.1211.2211.3211.4211.5211.6211.7211.8-48.45-48.4-48.35-48.3-48.25-48.2-48.15-48.1-48.05-48-47.9500.1
23、0.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.70.80.91a=rand(200,1);b=rand(200,1);c=rand(200,1);scatter(a,b,100,c,p)二、三维绘图三维绘图的主要功能: 绘制三维线图 绘制等高线图 绘制伪彩色图 绘制三维网线图 绘制三维曲面图、柱面图和球面图 绘制三维多面体并填充颜色(一)三维线图 plot3 基本的三维图形指令调用格式:plot3(x,y,z) x,y,z是长度相同的向量plot3(X,Y,Z) X,Y,Z是维数相同的矩阵plot3(x,y,z,s) 带开关量plot3(x1,y1
24、,z1,s1, x2,y2,z2,s2, ) 二维图形的所有基本特性对三维图形全都适用。定义三维坐标轴 大小 axis(xmin xmax ymin ymax zmin zmax ) grid on(off) 绘制三维网格 text(x,y,z,string) 三维图形标注 子图和多窗口也可以用到三维图形中例:绘制三维线图t=0:pi/50:10*pi;plot3(t,sin(t),cos(t),r:)010203040-1-0.500.51-1-0.500.51(二) 三维饼图pie3(4 3 6 8 9)27%30%20%13%10%(三)三维多边形 fill3 = fill 三维多边形的
25、绘制和填色与二维多边形完全相同 调用格式: fill3(x,y,z,s) 与二维相同例: 用随机顶点坐标画出5个粉色的三角形,并用黄色的表示顶点y1=rand(3,5);y2=rand(3,5);y3=rand(3,5);fill3(y1,y2,y3,m);hold on;plot3(y1,y2,y3,yo)(四)三维网格图 mesh 三维网线绘图函数 调用格式: mesh(z) z为nm的矩阵,x与y 坐标为元素的下标 mesh(x,y,z) x,y,z分别为三维空 间的坐标位置矩阵的三维网线图z=rand(6);z=round(z);mesh(z)0246024600 .20 .40 .6
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 理学 理学数学建模讲座 数学 建模 讲座
限制150内