已阅读5页,还剩5页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数学模型脚本l MATLAB入门 启动 MATLAB后,就进入 MATLAB命令窗口(Command Window)或称工作空间(Workspace),见图01。这里约定 (l)所有在命令窗口输入的命令都用 Courier New字体,并以开头,请读者注意为系统提示符,不要以为是输入字符。 (2)显示结果用 Times New Roman字体。 (3)号后面的文字用于注释,并不参与运算。图0-1 1 简单的运算 看一个简单的例子,计算 , 其中r =2用MATLAB计算如下: r2; 分号“;”使此运算结果不显示 v= 4/3*pi*r3 没有分号,显示结果v 33.5103 系统直接显示结果,就像计算器那样注:(l) MATLAB命令书写格式灵活,可多命令写一行,也可一个命令写多行。同一行命令用逗号或分号分开,若命令很长,一行不够,可用三点“”续行。 (2) MATLAB允许使用箭头键来编辑前面的命令。试一试你就明白了。 2变量和数据 (1)变量类型 MATLAB使用变量最常用的有数值数组(double array)和字符串(char array)两类。所有数值变量以双精度(double)方式存储,不区分整数、实数、复数等,变量类型和数组大小也无需预先定义,例如a=1+2*ia l.00002.0000i 复数 al=This is a stringalThis is a string 字符串A1 2;1 3A= 1 2 1 3 二维数组,即矩阵 (2)常量和特殊变量 MATLAB中一些常用的常量见表0-1。 表0-1常量一览表ans是系统本身一个特殊变量名,若运算结果没有赋予任何变量,系统就自动将其赋予ans。 (3)变量查询和清除 MATLAB命令窗口中用户变量可用who或whos查询,键入 whoYour variables are:A a a1r vwhos Name Size Bytes Class A 2x2 32 double array a lxl 16 double array(complex) al lx16 32 char array r lxl 8 double arrayv lxl8 double arrayGrand total is 23 elements using 96 bytes以上信息也可从工具栏中的工作空间测览器(Workspace Browser)观察到。变量的值可通过键人变量名得到,例如 a 显示 a的值a 1000020000i若命令窗口中有些变量不再使用,可使用命令clear来清除。 clear a A 清除变量 a和 Aa? Undefined function or variable a 说明 a已清除ala1This is a string al未清除clear 清除工作空间所有变量(慎用!)注意 clear与菜单 EditClear Command Window的区别。后者作用是将稿纸(窗口显示)擦干净。 (4)变量命名规则 MATLAB变量名总以字母开头,由字符、数字和下划线组成。有效字符长度为31个,且区分大小写,如a与A表示不同变量。用户在定义变量时,要尽量防止它与系统的常量名、特殊变量名、函数名等冲突,系统不会告诉你冲突发生,造成的后果是系统有些原来的功能暂不能使用。当这些变量被清除或MATLAB重新启动后,这些功能得以恢复。 (5)数据显示格式 MATLAB缺省的数据显示格式为:当结果为整数,就作为整数显示;当结果是实数,以小数点后4位的精度显示。若结果的有效数字不在这一范围,以科学计数法显示(如 1e-6表示10-6)。数据显示格式可通过命令format改变。需要指出的是,显示格式的改变不会影响数据的实际值,所以不会影响计数精度。其计数精度约为16位有效数字。cpic 3.1416 format rational;c c = 355/113 最接近的有理数之一format long;c c = 3.14159265358979 小数点后 14位format;cc 3.1416 恢复MATLAB还允许使用fprintf格式化输出,其使用方法与C语言基本一致。fprintf( 20.6f,c) 3.141593 (6)数据保存和调用 当我们退出MATLAB时,命令窗口中变量不复存在。为了保留变量值,我们可在退出之前用命令save将变量连同它的仅用二进制方式存储在数据文件中(详见save的帮助信息)。例如 A=1,2;1,3; save 所有变量和数据写入数据文件matlab.mat 好了,现在退出MATLAB,再次启动,用 A ? Undefined function or variable A 工作空间是空的 load 从 matlab.mat调出A A = 1 2 -1 3 说明矩阵A已存在于命令窗口 菜单FileSave Workspace As的功能与save命令等价。需要注意的是,save只能保存变量和数据,不能保存命令(保存命令行须通过程序设计M文件方式实现)。 为了与其他应用程序交换数据,可能需要用ASCll码方式或格式化文件来传递数据。save和load提供了写和读ASCll码数据文件的选项(详见save和load的帮助信息)。MATLAB还允许使用C语言读写命令fprintf,fscanf,fopen,fread等来传递格式化数据文件,其使用格式与C语言基本一致。l 矩阵和数组运算 MATLAB基本数据单元是无需指定维数的矩阵,标量可看作 IXI矩阵,n维行向量或列向量分别可看作IXn或nXI矩阵。 1矩阵的输入 输入矩阵最基本的方法是直接输入矩阵的元素,用中括号 表示矩阵,同行元素间用空格或逗号分隔,不同行间用分号或回车分隔,例如 clear;al,2,3;4,5,6;7,8,9 a= 123 456 789或 a1 2 3 这种方式特别适用于大型矩阵 4 5 67 8 0 a 123 456 780对于等差数列构造的向量,可用冒号运算生成 b=0:pi/10:pi 初值:增量:终值b =Columns 1 through 7 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850Columns 8 through 11 2.1991 2.5133 2.8274 3.1416b2 =l:10 增量缺省值为1b2 = 1 2 3 4 5 6 7 8 9 10size(a),length(b) 矩阵的阶数,向量的长度ans 3 3ans 11矩阵元素双下标编址按通常方式,单下标编址按列排序,一律从1开始编址(注意下标不能用0)。利用下标可获得矩阵元素或修改矩阵,例如a(3,2),a(6)ans 8ans 8 a(3, 4)1 将第3行第4列元素设为 1a 1230 4560 7801可以从大矩阵中抽取一个小矩阵, c = b(l:4)c 0 0.3142 0.6283 0.9425d = a(1 3,:) 取a的第1,3行,所有列d = 1230 7801小矩阵也可拼成大矩阵, ec;de 0 0.3142 0.6283 0.9425 1.0000 2.0000 3.0000 0 7.0000 8.0000 0 1.0000一些特殊矩阵可用内部函数产生,如 ones( m, n)生成元素全为1的m行n列矩阵, zeros(m,n)生成元素全为0的m行n列矩阵等。 矩阵的输入还有其他方式,如用文件产生,用剪贴板写入,从外部数据文件读入等。 2矩阵运算符 矩阵有多种运算。矩阵运算符见表02。 表02矩阵运算符 clear; A =1 -1;0 2; B0 1; 1 -1; A*Bans -12 2 -2AB 即 A-1Bans 0.5000 0.5000 0.5000 -0.5000 A/B 即AB-1ans 01 20A100ans 101 99 100 102 3.数组运算符(点运算符)点运算符.* 、/、.表示矩阵对应元素的运算,在作图、编写函数时经常使用。A.*B 注意与A*B的区别ans 0-1 0-2A.Bans 0-1.0000 Inf-0.5000A.2ans 11 041./A等价于ones(size(A)./Aans = 1.0000-1.0000 Inf O.50O0 4关系与逻辑运算 MATLAB有6个关系运算和3个逻辑运算符,都是对于元素的操作: 小于; = 不等于; 大于;|或; =大于等于;非。 =等于;在MATLAB中,“真”用1表示,“假”用0表示,而逻辑运算中,所有非零元素作为1处理。 ABans 10想一想,什么意思? 01ABans 01 01 5.数学函数 矩阵的数学函数也是按元素的运算,使用通常的函数符号,如 sin(A), asin(A), tan(A),eexp(A),sqrt(A),log(A)等,例如 sqrt(A)ans 1.0000 01.0000i 0 1.4142l 图 形以下是有关作图的常用函数与命令函数与命令注 释函数与命令注 释plotfplotplot3meshgridmeshsurfcontourcontour3titlexlabelylabelzlabel基本二维图形;一元函数图象;空间曲线;网格数据生成;网面图;曲面图;等高线图;3维等高线图;标题;X轴说明;y轴说明;Z轴说明;clabeltextgridlegendholdaxxsVllWsubplotfigureclfclose等高线高度标志;文本;格栅;图例;图形保持;定制坐标轴;改变视点;子图;新图形窗口;清除图形;关闭图形窗口。1曲线图plot(x,y)作出以数据(x(i), y(i)为节点的折线图,其中 x,y为同长度的 向量;plot(x1,y1,x2,y2,)作出多组数据折线图;fplot(fun,a,b)作出函数 fun在区间a,b上的函数图;plot3(x,y,z)空间曲线图,其中x,y,z为同长度的向量。 plot(1 4 2 5,3 5 1 4) 依次连接(l,3),(4,5),(2,1),(5,4)四点(如图0-2)图0-2 plot(1 4 2 5,3 5 1 4,1:3,2 10 3) 注意颜色图形可直接保存为M文件,也可利用图形窗口菜单 EditCopy figure将其作为图片剪贴到Word或其他应用程序中。 图形的线型,标记,颜色均可设定。常用的见表0-3。表03图形元素的设定例1一元函数图和(结果见图0-3) fplot(x.3-x-1,-1,2); hold on;hold on在作下一幅图时保留已有图象x-1:0.2:2; yabs(x).2.*sin(5*x);plot(x,y,:ro);hold off; hold off 释放hold on图0-3 曲线图2曲面图x,y= meshgrid(xa,ya) 当 xa,ya分别为 m维和 n维行向量,得到 x和 y均为 n 行m列矩阵。meshqrid常用于生成x-y平面上的网格数据;mesh(x,y,z)绘制网面图,是最基本的曲面图形命令,其中x、y、z是同阶矩阵,表示曲面三维数据;mesh(xa, ya,z) xa,ya分别为 m维和 n维向量,z为 n行m列矩阵。等价于先x,ymeshqrld(xa,ya)再mesh(x,y,z);surf(x, y, z)绘制曲面图,与mesh用法类似;contour(x, y, z)绘制等高线图,与mesh用法类似,可指定 z的范围;contour3(x, y, z)绘制三维等高线图,与mesh用法类似,可指定z的范围。例如 xa1:3, yal:4;x,ymeshgrid(xa,ya); z=xy;mesh(x,y,z)x,y,z 1 2 3 1 1 1 2 3 4 1 2 3
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 供应链沟通协调策略
- 公共关系与社交礼仪学习通超星期末考试答案章节答案2024年
- 舞蹈工作室兼职教练协议
- 渔业渔政管理机构五险协议书
- 家庭实+用护理技术学习通超星期末考试答案章节答案2024年
- python程序设计(电信和通信)学习通超星期末考试答案章节答案2024年
- 纺织品行业合同保管规定
- 自然辩证法学习通超星期末考试答案章节答案2024年
- 水毁修复施工合作协议
- 城市照明工程框架合同样本
- 矿热炉(电炉)运行有渣冶炼熔炼特性、电极插入深度解析与控制方法
- 妊娠晚期促子宫颈成熟与引产指南
- 基金委托募集合作协议
- GB/T 4942-2021旋转电机整体结构的防护等级(IP代码)分级
- 马工程《刑法学(下册)》教学课件 第20章 侵犯公民人身权利、民主权利罪
- GB/T 3820-1997纺织品和纺织制品厚度的测定
- GB/T 3215-1982炼厂、化工及石油化工流程用离心泵通用技术条件
- GB/T 2885.5-2008矿用窄轨车辆第5部分:平板车
- GB/T 25021-2010轨道检查车
- GB/T 2492-2003普通磨具交付砂轮允许的不平衡量测量
- 医院运行与医疗业务指标数据统计收集管理规定
评论
0/150
提交评论