数模培训—matlab1

1、学建模培训数学建模培训陈华陈华 讲师讲师234命令窗口命令窗口当前工当前工作目录作目录当前工当前工作空间作空间输入命令的输入命令的历史记录历史记录命令命令提示符提示符5q Matlab 区分大小写区分大小写，它的，它的命令和函数全是小写命令和函数全是小写的的 q 一行可以输入几个命令，用分号一行可以输入几个命令，用分号“;” 或逗号或逗号“,”隔开隔开 q 续行符：续行符： “”（三个点），如果语句很长，可用续行符（三个点），如果语句很长，可用续行符将一个语句写成多行。将一个语句写成多行。续行符的前最好留一个空格。续行符的前最好留一个空格。q 注释符：注释符： “%”

2、，其后面的内容为注释，对，其后面的内容为注释，对 Matlab 的计的计算不产生任何影响算不产生任何影响 6对于线性系统有Ax=b789列优先存贮列优先存贮10三维数组的示意图三维数组的示意图1112q 数据输入数据输入 (1) 直接输入：直接输入： a=1 2 ; 3, 4 (2) 冒号冒号 “:” 运算符运算符：初值初值 : : 步长步长 : : 终值终值 同一行中各元素之间用同一行中各元素之间用“空格空格”或或“,” (英文状态下英文状态下) 分开；分开；行与行之间用行与行之间用 “;” 或或“回车回车”分开分开 a=1:5， b=0:pi/4:pi (3) 由向量或小矩阵生成：由向量或

3、小矩阵生成：X=a ; b(5) 交互式输入：交互式输入： input n=input(Please input n: ) (4) 由由 数据文件生成数据文件生成13n赋值语句的一般形式：var = expressionn例如：nx=1;ny=2; nvar=40*x;nvar2=var/5;narray=1 2 3 4;14151617一些生成特殊矩阵的函数一些生成特殊矩阵的函数18%数值数值12%字符字符1、2ASCII码码192021 pi 圆周率圆周率 的近似值的近似值3.1415926 inf或或Inf 表示正无穷大表示正无穷大,定义为定义为1/0 NaN 非数，它产生于非数，它产生

4、于0 ，0/0， / 等运算等运算 i，j 虚数单位虚数单位 ans 对于未赋值运算结果，自动赋给对于未赋值运算结果，自动赋给ans eps 容差变量，定义为容差变量，定义为1.0到最近浮点数的距离到最近浮点数的距离.22q 矩阵元素的提取矩阵元素的提取 q(1) (1) 单个元素：单个元素：A ( 2, 3 )A ( 2, 3 )q(2) (2) 整行或整列：整行或整列：A ( 2, : )A ( 2, : )， A ( :, 3 )A ( :, 3 )q(3) (3) 子矩阵：子矩阵：A ( 2:5, 4:8 )A ( 2:5, 4:8 )， A ( 1,3, 2,4 ) A ( 1,3,

5、 2,4 ) ， A A ( 3,2, 2,4 )( 3,2, 2,4 )q(4) (4) 删除矩阵的行列删除矩阵的行列: :A =， A ( 3, : )=， A ( :, 2,4 )= 23q 矩阵的转置：矩阵的转置：共轭共轭“”，非共轭，非共轭“.”q 矩阵的翻转和旋转：矩阵的翻转和旋转： fliplr、flipud、rot90q 矩阵元素重组：矩阵元素重组： reshape(A, m, n)q 查看矩阵的大小：查看矩阵的大小：size(A)、size(A, 1)、size(A, 2)将将 A 排成一个排成一个m n 的矩阵，满足的矩阵，满足 m n = #Aq 数组运算（点运算）：数组

6、运算（点运算）：.*，./，.，.q 矩阵算术：矩阵算术： +，-，*，/，参与参与 “+, -, .+, -, .* *, ./ , ., ./ , .” 运算的对象必须具有相同的形运算的对象必须具有相同的形状状24252627282930313233343536利用记事本创建m文件3738392.6 matlabmatlab联机帮助命令联机帮助命令 help 功能提供功能提供matlab大部分主题的在线帮助信息大部分主题的在线帮助信息help 显示显示help 主题一览表主题一览表help plotxyz 显示有关三维做图指令帮助信息显示有关三维做图指令帮助信息help 显示特殊字符与符号

7、帮助信息显示特殊字符与符号帮助信息help help 显示显示help 的帮助信息的帮助信息虽然虽然help可以随时提供帮助，但必须知道准确的函数可以随时提供帮助，但必须知道准确的函数名称。当不能确定函数名称时，名称。当不能确定函数名称时，help就无能为力了就无能为力了。40 Lookfor函数函数 它可提供通过一般的关键词，它可提供通过一般的关键词，搜索出一组与之相关的命令搜索出一组与之相关的命令lookfor fourier 寻找含有傅立叶变换的相关指寻找含有傅立叶变换的相关指令令2.6 matlabmatlab联机帮助命令联机帮助命令help、lookfor 两个指令构成了两个指令构成

8、了matlab语言相语言相当完善的在线帮助查询系统当完善的在线帮助查询系统 41q 三角函数三角函数42q 基本数学函数基本数学函数43q 取整函数取整函数q 矩阵相关矩阵相关函数函数44q 矩阵的转置：矩阵的转置：共轭共轭“”，非共轭，非共轭“.”q 矩阵的翻转和旋转：矩阵的翻转和旋转： fliplr、flipud、rot90q 矩阵元素重组：矩阵元素重组： reshape(A, m, n)q 查看矩阵的大小：查看矩阵的大小：size(A)、size(A, 1)、size(A, 2)将将 A 排成一个排成一个m n 的矩阵，满足的矩阵，满足 m n = #Aq 数组运算（数组运算（点运算点运

9、算）：）：.*，./，.，.q 矩阵算术矩阵算术： +，-，*，/，参与参与 “+, -, .+, -, .* *, ./ , ., ./ , .” 运算的对象必须具有相同的形运算的对象必须具有相同的形状状45464748492.7 二维绘图二维绘图-50plot(t,y,red,t,y1,gre,t,y2,blu)5152t=0:pi/100:2*pi;y=sin(t); y1=sin(t+0.25); y2=sin(t+0.5);plot(t,y,r-,t,y1,g:,t,y2,b*)title(不同的正弦曲线);legend(sin(t),y1=sin(t+0.25),y2=sin(t+

10、0.5)xlabel(角度t),ylabel(y值)53540246810-1-0.8-0.6-0.4-0.200.20.40.60.81sin(t)cos(t)正 弦和 余 弦曲 线时 间 t正 弦、 余 弦正 弦余 弦55565701234567-1-0.8-0.6-0.4-0.200.20.40.60.815805101520253035404550-8-6-4-2024681059051 0-1-0 .8-0 .6-0 .4-0 .200 .20 .40 .60 .81051 0-1-0 .8-0 .6-0 .4-0 .200 .20 .40 .60 .81051 0-1-0 .8-0

11、 .6-0 .4-0 .200 .20 .40 .60 .81行列 绘图序号6001234567-10101234567-10101234567-101610 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 .81620 1 2 3 4 5 6 7- 1- 0 . 8- 0 . 6- 0 . 4-

12、 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 . 8163fplot(sin(x),2*pi*-1 1) fplot(sin(x),2*pi*-1:0.01:1)64-6-4-20246-6-4-2024665例如: x2+3y2=5ezplot绘制隐函数ezplot(x2+3*y2-5) %绘制隐函数图形绘制隐函数图形ax

13、is(-4,4,-4,4) %设置坐标轴范围设置坐标轴范围ezplot(f) f 为包含单个符号变量x的符号表达式， 在x轴的默认范围-2*pi 2*pi内绘制f(x)的函数图66假设只想绘制 之间的图形，可以按照如下操作(/4, ),( 1,3)xy ezplot(x2+3*y2-5,-pi/4,pi,-1,3) %绘制隐函数图形绘制隐函数图形axis(-4,4,-4,4) %设置坐标轴范围设置坐标轴范围ezplot(f,xmin,xmax) 给定区间ezplot(f,xmin,xmax,figure(n) 指定绘图窗口绘图。67例如: 绘制参数方程对应的图形 sin3sin3 sin ,(

14、0, )xtcostytt tezplot(sin(3*t)*cos(t), sin(3*t)*sin(t),0,pi) %绘制隐函数图形绘制隐函数图形axis(-4,4,-4,4) %设置坐标轴范围设置坐标轴范围6869707172x=0:pi/20:2*pi;y=sin(x);stem(x,y)73x=0:pi/20:2*pi;y=sin(x);area(x,y)74x=0:pi/20:2*pi;y=sin(x);scatter(x,y)757677787980例:参数方程t=0:.1:2*pi;x=t.3.*sin(3*t).*exp(-t);y=t.3.*cos(3*t).*exp(-

15、t);z=t.2;figure(1)plot3(x,y,z)grid onfigure(2)stem3(x,y,z)figure(1)figure(2)81例:参数方程sin( ),cos( ),xtytztt=0:pi/50:2*pi;x=sin(t);y=cos(t);z=t;figure(1)plot3(x,y,z)grid onfigure(2)stem3(x,y,z)figure(3)grid offfill3(x,y,z,g)8227%30%20%13%10%8384850246024600.20.40.60.818687坐标矩阵坐标向量纵坐标矩阵-505-505010203040

16、508889调和色0.5 0.5 05 灰色0.5 0 0 暗红色1 0.62 0.4 铜色0.49 1 0.8 浅绿0.49 1 0.83 宝石兰9091929394-202-202-505xyPeaks95-202-202-505xyP eaks9697989901020300102030-10-5051010001020300102030-10-5051010102040600204060-10-5051010202040600204060-10-505101031040204002040-10010az=-37.5,el=300204002040-10010az=-15,el=6002040-10-50510az=-9