常用数学软件教程035第3章matlab使用基础5节编程

第3 使用基编CC++语则都比较简单它是一种解释性编程语言修饰和调试也很方便前面学习和使M存盘后将和软件中已有程序一样可以反复调用。进入和退出编辑调试窗3-编辑调试窗口中的快捷File（文件、Edit（编辑、Text（正文、Debug（调试）Hekp（帮助）等。每个主菜单133-14中，同时列出了相应子菜单及其功能。3-按 相应的菜 功主菜 子菜New/M-M-文件Edit（编辑粘贴剪下或的内ShowSet/Clear/Clearall M-文凡在编辑调试窗口中用M-文件，扩展名为“.mM-文件分为两类：M-指令文（scriptfile，也叫文件和M-函数文（functionfile1.M-指令文M-指令文件就是在编辑调试窗口中用语言编写的一连串指令集%comment（ment为去除除变量指令“clear3.5.1M-AAa的y=f(a)=a2+sin3a-e-a。解M-按钮（或按Ctrl+s）存盘，文件取名为A中各元素ay=f(x)=a^2+(sina)^3-exp(-clearcif指令擦去原有曲线A=input('输入数值、向量或矩阵自变量A=A的内容栏内填写“ex67（（>>ex67（debug→Ru>>然后，再键入（输入[3*pi52;710>>A=[3*pi52;710A f1 ※注2.M-函数文函数文件是在M-文件，每个函数文件定义了sin、log等，函数文件的结构有下述特点。的第一行必须以function开头，然后用函数的通用形式y=f(x)规定输入变量x、输入函y和函f，xy均不限于一个变量，这里确保的函数名f就是存盘时的M-函数※注Error（’message’）——显示出错信息，终止程序运行；Warning（’message’）例 编写M-函数文件，由它求矩阵A中各元素a的函数值：y=f(a)=a2+sin3a-a解M-在编辑调试窗口中键入程在编辑调试窗口中键入程（或按Ctrl+s）存盘,文件取名为A中各元素a%y=f(a)=a^2+(sina)^3-exp(-a)functionf1=ex68(A)3.5.1ex683.5.1f1=ex68(A)中，f1为输出的函数值，ex68为M-函数文件名，A为函数这时若欲求A 2中各元素a为自变量的函数y=f(a)=a2+sin3a-e-a取值 5 >>y=ex68([3*pi52;710y >>A=[3*pi52;710y 关系和逻辑运它只能取10，其中1表示真（T，0表示假（F。两个同维矩阵间的关系运算规定为它布尔矩阵（01的矩阵。3-153-<> 3 例 已知a1 6,a25。求a3a1a2（注意：=是赋值号） 9解>>>>a1=[123;456;78a3 的“异或（xor”逻辑运算。在逻辑运算中，非零元素的逻辑量为1，表示“真（或T）;零0（3-16中。3-符 意 中允许数与矩阵之间进行逻辑运算，规则与关系运算相同，是数与矩阵各个3-17中。3-逻辑量及其运 B1100110真值（逻辑量 00101101100110011000 9例 求数组矩阵a4 0的“非，a4和0的“异或 8解>>在命令窗口键>>a4=[10-509;3-2060;0057键a5 a6 组成计算机程序的一系列指令语句可以分成两类：运算语句和控制顺序语句。编程的主要任务是安排调整好控制顺序语句，使运算语句的运行顺序合理，运算循环结构语即循环使用这同一组指令，直到满足某种条件（循环结束条件。为此，中设有两3- while逻辑关系式

功能特 规定了执行“语句组”的条1（真）时，往下执行语句组；0（假while…end循环，执行后面的程序。例

n

210

mm m解

n

1

在编辑调试窗口中键入程在编辑调试窗口中键入程（或lt3_5_5Debug→Run（F5）运whilea<1e4a=a+n^2;n=n+1;%语句组end,m=n-1,a%结束循环，显示m，am在指令窗口中a如果是

3-for循环次数表达式功能特 规定了执行“语句组”的次

3-例 创建一个范德矩 1 3 解在编辑调试窗口中键入程（或ex366 解在编辑调试窗口中键入程（或ex366 n=length(a);%规定执行语句组的次数fork=0:n-1 （>>输入 德矩阵的基本向量然后，再键入（输入[452>>a=[452D1 分支结构（条件转移）复杂的计算中，有时需要根据表达式满足的条件，确定下一步进行哪一项工作，中设有分支（条件）结构语句。

3-if逻辑关系式

3-例 按钮（或按Ctrl+s）存盘，文件取名为forx=0:0.3:4按钮（或按Ctrl+s）存盘，文件取名为forx=0:0.3:4y=x^2-2*sin(x);%计算函数值语句组ify<0disp([x,y])%显示x和yy=x^2-2sinx') （>>xy=x^2-----

3-if

B”

3-例 nn>100“y>100y100按钮（按钮（或按Ctrl+s）存盘，文件取名为ifn^2-n^3*sin(n)<=100（>>>>然后，再键入（>>nans分层条件语句叫递阶选择分支语句322329.

3-ifA1elseif2规定了执行“语句组2A23

3-例3.5.9 由键盘输入任一自变量时判断函数y=x4-19x3+32是负数偶数还是其他。解先判断函数y的正负，再判断其奇偶性，据此可在编辑调试窗口中键入下列程序按钮（或按Ctrl+s）存盘，文件取名为clear,x=input('输入自变量x=')%由键盘输入自变量xify<0%判断函数y是否为负数elseifrem(y,2)==0%判断y是否为偶数>>（>>然后，再键入（>>x3a3.5.1036个头、100按钮（或按Ctrl+s）存盘，文件取名为n=1;%设鸡起码有1只，可以多设while1%while后的逻辑值为1时将无限循环下去ifrem(q-n*2,4)==0&(n+(q-2*n)/4)==p按钮（或按Ctrl+s）存盘，文件取名为n=1;%设鸡起码有1只，可以多设while1%while后的逻辑值为1时将无限循环下去ifrem(q-n*2,4)==0&(n+(q-2*n)/4)==pbreak（>>>>>>lt3_5_10q=然后，再键入（>>lt3_5_10q=然后，再键入（>>nans最后的显示语句disp可以用sprintf代替，换成“sprintf('总共有鸡兔%d只，它们共有%d只脚，则其中有%d只鸡，%d只兔.',[p,q,n,p-n])这样可将全部内容都显示在同一行按钮（或按Ctrl+s）存盘，文件取名为n=1;%设鸡起码有1只，可以多设while1%while后的逻辑值为1时将无限循环下去ifrem(q-n*2,4)==0&(n+(q-2*n)/4)==pbreak（>>>>lt3_5_10_0p=36q=然后，键入（lt3_5_10_0p=36q=最后，再键入（>>ans22只鸡，14只兔※注smsrtindent，则自动调节成缩进格式。※※注lt3_5_10_0p=37lt3_5_10_0p=37q=101作>>Errorin==>lt3_5_10_0at3-233-30。是逐级向下进行的，而是由switch后的“表达式结果”一次确定的，表达式结果可以是数，casecase后面的语句组。

3-switchcasecase（或字符串）case后的数值i（或字符串）一致时，就Ai 组 3-例3.5.11 数（不分正负，还是其他数（0视为偶数。按钮（按钮（或按Ctrl+s）存盘，文件取名为clear,x=input('x=');%由键盘输入自变量xy=x^4-19*x^3+x+32%计算函数值switchmod(y,2)%得出y被2除的余数case1%余数为1时，y是奇数case0%余数为0时，y是偶数（>>>>然后，再键入（>>y-两类M-文件的转两类例3.5.12 编制一个M-指令文件，列出由键盘输入的正整数a到b之间的所有素数，并将它改编成M-函数文件。forforifrem(m,n)==0%检验m可否被n整除ifm==n%选取只能被自己整除的数按钮（或按Ctrl+s）存盘，文件取名为>>（>>然后，键入（>>ex3612>>ex3612最后，再键入（>>在命令窗口键入b0c2357c=ex78(a0,b0)后一行中的c，就使这个M-指令文件变成M-函数文件，函数文件在调用时必须要有输入参※补充forforifrem(m,n)==0%检验m可否被n整除ifm==n%选取只能被自己整除的数在编辑调试窗口中键入程（或按Ctrl+s）存盘,文件取名为function3.5.1ex783.5.1c=ex78(a0,b0)，c为输出的函数值，ex78M-函数文件名，a0,b0为>>在命令窗口键入命c 57编程中的一些控编写M-文件时，常常需要设置一些查看、检验。修改和人机方面的控制指令，它键盘输入指令a=input('Pleaseinputanumbera，然后继续往后运行。>>a=input('Pleaseinputa在命令窗口键入命Pleaseinputaa2a=input（’Pleaseinputanumber何数据（数值或字符串）a。暂停运行指令pause（nn秒后继续运行。人机切换指令return，回车后，程序才继续运行。程序显示指令M-文件时，屏幕上并不现实运行程序的内容，若在运行程序前编入“echo在命令窗口键入命>>在命令窗口键入命>>>>echoex78???Inputargument"a0"is>>forforifm==n forn=2:mif %mnifm==n forif %mnif %mnifm==n forif %mnifm==n%选取只能被自己整除的数forn=2:mifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mnifm==n forif %mnifm==n%选取只能被自己整除的数forn=2:mifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mnif %mnif %mnifm==n forif %mnifm==n forn=2:mif %mnif %mnifm==n forn=2:mif %mnifm==n%选取只能被自己整除的数forn=2:mifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mnifm==n forif %mnifm==n%选取只能被自己整除的数forn=2:mifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mnifm==n forif %mnifm==n forn=2:mif %mnif %mnifm==n forn=2:mif %mnifm==n%选取只能被自己整除的数forn=2:mifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mnifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除ifrem(m,n)==0%mn整除i