三本数学实验课件3.pptVIP

M 文件 大家已熟悉了 MATLAB 命令窗口的一些操作,它简单,得结果方便.但这个窗口也有不便之处:当需要几十乃至更多的语句去得到一个结果时,命令窗口的操作将很不方便.一是每一次换行（回车）要执行语句,使得窗口很乱,二是当一个命令因为错误而执行了,与它相关的命令都要重新执行一遍,这很麻烦.克服这一缺陷的方法是编制M 文件，这需在Editor窗口进行. Editor窗口的打开:在命令窗口中点 File→New → M-file或直接点击上面的白色标志,就打开了Editor窗口.如果过去已制作过M 文件,则通过File-Open打开那个文件也打开了Editor窗口. 在Editor窗口,所有的语句当编辑完成后才统一的执行,方便应用. M文件的起名:文件的名字不要用中文、数字开头,要用英文字母开头,而且要保存在默认的地址中,不要保存在个人的文件夹中. 对文件的内容可以做注释,注释用”%”即可,在”%”后的内容将不被执行. 执行方法,点击标志Debug → run或快捷标志或在命令窗口中输入文件名并按回车,而结果则在命令窗口显示. 例如 X=1.2 % self variable Y=23.5, % self variable Z=xy-3x+y/5-23 % result 练习： 已知x=4.5,y=67.2,制作文件求 z=x?y+x?y-1.4x2+2.8y3-123 的值并给予注释 控制流 当对大量有规律的数据处理时，需要用控制流。它们由几个结构组成，分别介绍如下 if-else-end结构，其一般结构有三种 if 条件表达式 语句体 end 语句的执行过程是：首先检验条件表达式，满足了就执行语句体，不满足就不执行。 2。 if 条件表达式 语句体 else 语句体 end 3. if 表达式 语句体 elseif 语句体…… else 语句体 end 执行过程是：遇到第一个满足条件的执行相应的语句体，否则执行else后面的语句，而后跳出循环。 例: 输入任意的x值，求函数值y，其中 程序如下 clc x=input(x=) if x=0 y=x+1, else y=x-1, end 例: 输入任意的x值，求函数值，其中 程序如下 clc x=input(x=) if x0 y=pi-1, elseif x3 y=2*exp(x)-1, else y=sin(2*x)+3^4.5-6 end 练习 编制文件对任意的自变量求函数值，函数如下 for 循环，结构如下 for x=一维数组 内容 end 执行规则是:对数组中的每个元素，内容依次执行一次。 如：对n=100分别求 clc n=100; he=0; ji=1; for k=1:n he=he+1/k; ji=ji*k; end he,ji 练习，分别求级数的前1000项和 while循环，结构如下 while 表达式 内容 end 执行规则是:当表达式成立时对内容执行，接着再检查表达式、执行，直到表达式不成立才跳出循环。 例: 3X+1问题（世界难题） 对任意正整数n， 若n为偶数，则用n/2得另一个正整数；若n为奇数，则用（3n+1）/2得另一个正整数。如此下去，对任意正整数最后必然得到1，如3→5 → 8 → 4 → 2 → 1。试编程验证。 clc n=3; while n(end)1 if mod(n(end),2)==0 n=[n n(end)/2]; else n=[n (3*n(end)+1)/2]; end end n 另：switch-case语句的用法类似于if-else语句，但一般针对符号变量，这里省略。此外, break语句的作用是中断循环，continue也类似，return用的少，易出问题，此处略去. 又一个例子:用二分法找出方程2x3-4x2+3x-8=0在（-10，10）内的一个根，精确到10-4. 分析: 连续函数f(x)在[a,b]上若满足, f(a)f(b)0，则f(x)=0在 (a,b)内必有一根。二分法的原理是求出f((a+b)/2)，若 f(a) f((a+b)/2) 0，则这个根在（a, ((a+b)/2 ）内，反之，在（ ((a+b)/2 ,b）。反复下去，直到满足精度后取区间中点作为根即可，程序如下 clc x=[-10,10] f1=2*x(1)