数学实验M文件.pptVIP

* 三 M文件设计 命令文件(脚本式) 变量、表达式、运算符 程序设计中的流程控制 函数文件的编写 ? ? ? ? 1、命令文件(脚本式M文件) 无输入参数，不返回参数 操作工作空间中的数据 对于需要多次执行的一系列命令特别有用 2、函数文件(函数式M文件) 可输入参数，也可返回参数 函数的变量为局部变量，也可设外部变量。 扩充MATLAB的函数库 M文件分为命令文件和函数文件两种 命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息 M文件分为命令文件和函数文件两种 ①在编辑窗口中编写; ②保存并对文件命名； ③命令窗口键入文件名运行； ④观察运行结果； ······························· 例2.1 将地球模型取为半径为R=6400(km)的球体,计算高度为1万公里的地球同步卫星对地面的覆盖率 在编辑窗口编写文件 planet.m 在命令窗口键入 planet 结果显示: ans = 30.4878 R=6400; S0=4*pi*R*R; d=10000; S=2*pi*R*R*d/(R+d); S/S0*100 + 加 - 减 ^ 矩阵幂 .^ 点幂 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 1 .算术运算符 2.关系运算符 ? 小于 ? = 小于等于 ? 大于 ? = 大于等于 ? == 等于 ? ~= 不等于 3.逻辑运算符 与 (and) | 或 (or) ~ 非 (not) MATLAB 程序中的流程控制 条件控制; 循环控制; 错误控制;终止运行控制 条件控制——有选择地运行程序块 1. if / elseif 语句 分析下列程序，说明运行结果及原因 clear a=100;b=10; if ab disp(ab) %显示字符串ab else disp(ab) %显示字符串ab end 例2 .2 判润年程序 ————程序设计中的流程控制———— year=input(input year:=); n1=year/4; n2=year/100; n3=year/400; if n1==fix(n1)n2~=fix(n2) disp(是润年) elseif n1==fix(n1)n3==fix(n3) disp(是润年) else disp(不是润年) end 润年条件有二 ①能被4整除，但不能被100整除； ②能被4整除，又能被400整除。 1. for循环 语法： for 变量 =初值：步长：终值 可执行语句 end 求裴波拉奇数列的第 n 项( n 2) n:=10, f(10)=55 n:=44, f(44)=701408733 循环语句 n=input(input n:=); f(1)=1;f(2)=1; for k=3:n f(k)=f(k-1)+f(k-2); end f(n) 2. while循环 while循环将循环体中的语句循环执行不定次数。 语法： while 条件表达式 循环体语句 end 条件表达式一般由变量、数字、逻辑运算、关系运算符和一般运算符组成，以判断循环的进行和停止； 只要表达式的值(逻辑值)结果为正确(非0)，循环继续；直到表达式值为0，循环停止。 角谷猜想. n=input(‘input n=’); %输入数据 while n~=1 r=rem(n,2); %求n/2的余数 if r ==0 n=n/2 %第一种操作 else n=3*n+1 %第二种操作 end end 例2.3 对任一自然数n，按如下法则