matlab数米粒实验报告,文化学习上机实验中遇到的一些Matlab函数.ppt

文化学习上机实验中遇到的一些Matlab函数.ppt

上机实验中遇到的一些Matlab函数,(1)* .* . (2)abc linspace(a,b,c) (3)log log10 (4)exp (5)length size (6)round ceil floor (7)inv mod,3.1 Matlab 的程序类型 3.2 程序流程控制 3.3 关系与逻辑运算符 3.4 函数 3.5 绘图,第三章 MATLAB 编程基础,3.1 Matlab程序类型,1、 脚本文件 以.m格式进行存取，包含一连串的MATLAB指令和必要的注解，所有的变量都存在Workspace中，没有输入参数，也不返回参数。 2、 函数文件(help function) 以.m格式存在，函数接受输入参数，然后执行并输出结果。,3.2 程序流程控制,1、 If, else, elseif 2、 For 3、While 4、Switch,1、 If, else, elseif,TYPE 1if expression statements end Example, if moda,2 0 aa3; end,TYPE 2 if expression statements else statements end ExampleifsumA,210 A1,1A1,1B1,1 A2,12*A2,1 else A1,1A1,1A1,1; A2,10.5 end,TYPE 3 if expression statements elseif expression statements else statements end,If语句可以嵌套 if expression 1 statements else if expression 2 statements else statements end end,步长缺省值为1，可以在正实数或负实数范围内任意指定。对于正数，循环变量的值大于终止值时，循环结束；对于负数，循环变量的值小于终止值时，循环结束。循环结构可以嵌套使用。,2、for循环语句 基本格式 for 循环变量起始值步长终止值 循环体 end,例help for,3、while循环语句 基本格式 while 表达式 循环体 end,若表达式为真，则执行循环体的内容，执行后再判断表达式 是否为真，若不为真，则跳出循环体，向下继续执行。,While循环和for循环的区别在于，while循环结构的循环体被执行的次数不是确定的，而for结构中循环体的执行次数是确定的。,,4、Switch语句 switch expression can be a scalar or a string case test1 command set 1 case test 2 command set 2 ...................................... otherwise command set last end,3.3 关系和逻辑运算符,关系运算符 假设有A1 2 -1 -5 B0 2 3 1 大于 AB ans1 0 0 0 A1 ans0 1 0 0 大于等于 AB ans1 1 0 0 等于 AB ans0 1 0 0 ; A1 ans1 0 0 0 不等于 AB ans1 0 1 1 ; A1 ans0 1 1 1, 非 A0 0 1 0 10,注意在处理逻辑运算时，运算元只有两个值即0和1，所以如果指定的数为0，MATLAB认为其为0，而任何数不等于0，则认为是1。,设有A5 -4 0 -0.5 B0 1 0 9, k2;exactpi2/6 while abs1-ss/exact1e-4 ssss1/k2; kk1; end,3.4 函数的编写,格式function 返回变量列表函数名输入变量列表 function outarg1, outarg2, fnameinarg1, inarg2,,求n function kfactorialn if nargin1,error 输入变量个数错误，只能输入 个数量;end if nargout1,error 输出变量个数过多;end if absn-floorn epsn0 判定n是否为非负整数 error n应该为非负数; end if n1 kk*factorialn-1; else ifn0 | n1 k1; end end,3.4 画图,axis，plot， subplot， figure， hold on, hold off,,ploty y 4 2 3 4 6 plotx,y t 00.14*pi y sint plott,y plotx1,y1,x2,y2,.. t1 00.12*pi t2 00.14 plott1,sint1,t2,sqrtt2 plotx,y,option 如 t1 00.42*pi t2 10.1;4; plott1,sint1, ,t2;cost2, g;,semilogx 横坐标为对数 semilogy 纵坐标为对数 loglog 横纵坐标均为对数 实例演示程序 t0.10.12*pi; ysint; semilogxt,y grid on,图形窗口的分割 Subplotm,n,I 分成m行n列，并选定I 实例演示程序 subplot2,1,1 t0.10.12*pi ysint semilogxt,y grid on subplot2,1,2 x10.015; yexpx plot x,y,线的颜色、线的类型、标记类型 The line colors have yellow, magenta, cyan, red, green, blue, white and black. The line styles have solid line - , dotted line , dot-dash line -., dash line - - and no line . The markers have point ., circleo, x-markx, plus, star*, squares, diamondd, trianglev, up-triangle, pentagramp, etc.,图形标注,坐标加标注 图形加标题 xlabel， ylabel， zlabel title, text,坐标轴 axis xmin xmax ymin ymax; hold on； hold off,实例演示程序 t00.14*pi； ysint； plott,y xlabelx轴(04pi)； ylabely轴 title正弦波,Fontsize, 12,图例 为图形加图例便于观察分析 legendstring1,string2,string3, 实例演示 t00.14*pi plott,sint, ,t,cost, * legendsin wave, cos wave, None,三维图例 mesh , surf , plot3,xx1linspace-3,3,20 xx2linspace-3,13,17 x1,x2meshgridxx1,xx2; zx1.43*x1.2x2.2-2*x1-2*x2-2*x1.2.*x26 surfx1,x2,z,,,,,,,,,,,,Q1(米3/秒),Q2(米3/秒),截面积A(米2),,,,,,,水位高度H(米),仿真实例水位系统,,,,,,,,,,,,,系统建模,实体,水池、管道、阀门,属性,参数信息A,输入信息Q1 Q2,输出(状态)信息H,活动,由于Q1， Q2的变化，H在不断地变化。,,A Q1 Q2,H,数学模型,不难看出，在上面的方程中有微分运算，这种运算在计算机中是无法实现的，可以将其以采样时间T近似地转化为差分方程,仿真建模,通过方程可以看出如果知道初始时刻的水位高度和每个时刻的流入流出水量，就能够求出任意时刻的水位的高度。,程序设计,假设,在Matlab上设计程序sw.m,仿真结果,