2-脚本文件、函数及数组.pptxVIP

MATLAB 程序设计MATLAB程序设计第二讲 脚本文件、函数及数组张 兵Email: zhangbing_end@163.com第页教学内容脚本文件函数和自定义函数MATLAB数组第页脚本文件（1/6）什么是脚本文件例 计算圆柱表面积和体积r = 5输入命令法h=10无法保存，不易修改脚本文件：扩展名为“.m”的文本文件，存储了一系列计算变量、命令、表达式。第页脚本文件（2/6）创建和保存新建脚本文件保存按钮第页脚本文件（3/6）块注释以%{开头，%}结尾注释以%开头代码单元以%%开头第页编写代码添加注释代码单元脚本文件（4/6）运行和调用执行按钮第页脚本文件（5/6）文件扩展名是.m通过调用文件名来执行文件名必须满足MATLAB的命名规则一般的脚本文件结构第页脚本文件（6/6）例子 第页函数完成一定功能的子程序使用方式如：sin(a) – 计算以弧度为单位的角度a的正弦值 linspace(a,b,c) – 生成一个以a为起始值、b为结束值、元素个数为c的向量函数名(参数1，参数2，…)变量=函数名()cosaexp[60]sind{60}a=1+sind(60)第页编写自定义函数(1/4)匿名函数定义方式：函数名＝@(输入变量1,输入变量2,…)表达式cylinder_volume=@(r,h)pi*r.^2.*h;circle_area = @(r) pi*r.^2;完成上述定义后，采用如下方法使用：%计算半径2、高度3的圆柱体体积，将结果存储在vol中vol = cylinder_volume(2,3);area = circle_area(5); %圆面积第页编写自定义函数(2/4)一般函数定义方式，在m文件中采用如下方式定义：function [输出值1,输出值2,…] = 函数名(输入参数1,输入参数2,…) 函数执行体end名称要一致tri_area.mfunction [area] = tri_area(a,b,c) s=(a+b+c)/2; area = sqrt(s.*(s-a).*(s-b).*(s-c));end调用方式：tri_area(3,4,5);第页编写自定义函数(2/4)add.mfunction [] = add(a,b,c) c = a + b;end答案：c = 0命令窗口：参数的作用域仅为函数内：在函数体内对输入参数的修改不影响调用时的值。 a = 1; b = 2; c = 0; add(a,b,c);问题：c = ???第页输入参数的作用域编写自定义函数(3/4)function […] = funname(…) ... fun_1(…); ...endfunction […] = fun_1(…) … fun_2(…); …endfunction […] = fun_2(…) …end...主函数funname.m子函数子函数第页函数脚本文件一般格式编写自定义函数(4/4)全局变量 global varname定义全局变量，并初始化为空数组定义方法： 在各个程序、函数中共享数据作用： clear global varname从全局工作空间中删除清除方法： clear varname从当前工作空间中删除 isempty(varname)如何判断全局变量是否已初始化：演示：Test_2_1第页再论MATLAB编程（1/2）MATLAB程序MATLAB软件操作系统计算机硬件MATLAB程序的执行层次结构第页再论MATLAB编程（2/2） 处理器MATLAB操作系统 存储器R=2;H=5;V=pi*R^2*H;ALU: pi*^2*RH第页MATLAB数组（1/4）什么是数组a=int32(2);R=2;A=[1,2,3,4,5,6,7,8];B=[1,2,3,4; 5,6,7,8];第页MATLAB数组（2/4）一个4维数组第页MATLAB数组（3/4）数组种类[]空数组标量，只包含一个元素一维数组二维数组三维数组更高维数组数组（array）……第页MATLAB数组（4/4）稀疏矩阵10×10二维数组，非零元素17个100000×100000二维数组非零元素5×100000第页数组的性质(1/5)元素：数组中的单个数据项数组a：索引第四个元素：a( 4 )元素值为13元素值和索引是一一对应的第页数组的性质(2/5)确定数组大小：size函数用法：n = size(a)输入：a为数组输出：n为一维数组，每个元素表示a每一维的大小第页数组的性质(3/5)用法：n = numel(a)输入：a为数组输出：n为数组a的所有元素个数 a = [1,2,3;4,5,6]; n = numel(a)n = 6第页确定数值元素个数：nu