(完整版)安徽师范大学《计算方法》实验课程序+答案(数学专业)六

1、实验六第3章 MATLAB程序设计第二节程序控制结构四、其他流程控制语句break语句continue 语句return语句try语句echo命令keyboard 命令pause命令1. break 和 continue 语句与循环结构相关的语句还有 break语句和continue语句.它们一般 与if语句配合使用.break语句用于终止循环的执行.当在循环体内执行到该语句时， 程序将跳出循环，继续执行循环语句后的下一语句.continue语句控制跳过循环体中的某些语句.当在循环体内执行到 该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环 . 例15求100, 200之间第一个能被

2、21整除的整数.for n=100:200if rem(n,21)=0 continueendbreakendnfor n=100:200if rem(n,21)=0breakendend n 12. try语句语句格式为：try语句组1catch语句组2endtry语句先试探性执行语句组1,如果语句组1在执行过程中出现错误，则将错误信息赋给保留的lasterr变量，并转去执行语句组 2. 这种试探性执行语句是其他高级语言所没有的.例16矩阵乘法运算要求两矩阵的维数相容，否则会出错.先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘.function C=rc(A,B)%A,B 为矩阵tryC

3、=A*B;catchC=A.*B;endlasterr%显示出错原因3. 请求键盘输入命令 keyboard暂停运行程序并处于等待键盘输入状态(便于 调试及程序运行中 修改变量)，键入“return'序继2运行.如function y=f(x) keyboardy=(5*x-6)/(x-3)*(x-3)+2);运行如下：>> y=f(4)K>> returny =4.6667>> fK>> x=4;K>> return ans =4.66674. echo命令(在命令窗口，程序调用前使用)一般情况下，M文件执行时，文件中的命令

4、不会显示在命令窗口 .echo命令可使文件命令在执行时可见.这对程序的调试及演示很有 用.对命令文件和函数文件，echo的作用稍微有些不同.对命令文件，其格式如下：echo on打开命令文件的回应命令echo off关闭回应命令对函数文件，其格式如下：echo file on使指定的file文件的命令在执行中被显示出来echo file off关闭指定文件的命令在执行中的回应习题:一、阅读理解下列程序，并在计算机上运行.1. qjsh.mfunction s=qjsh(a,x) % a 为多项式系数向量(按降幂)n=length(a);s=a(1);for k= 1:ns=s*x+a(k);e

5、nd(以x2 2x 1 测试，并计算3x4 7x3 9x2 23的值)2. wanshu.m function wanshu(n) %n 为正整数A=;for m=1:ns=0;for k=1:m/2if rem(m,k)=0 s=s+k;endendif m=sA=A,m;endenddisp(A)(对n 10测试，并对不同的n计算)3. dazhe.mprice=input(' 请输入商品价格:');switch fix(price/100)case 0,1rate=0;case 2,3,4rate=3/100;case num2cell(5:9) rate=5/100;c

6、ase num2cell(10:24) rate=8/100;case num2cell(25:49) rate=10/100;otherwiserate=14/100;endprice=price*(1-rate)（对各种商品价格进行计算）4. num.m% Classic “3n+1” problem from number theory. while 1n=input( Enter n, negative quits. ) if n<=0, break, endwhile n>1if rem(n,2)=0n=n/2elsen=3*n+1endendend5. dd.mfunc

7、tion k, y=dd(x,e) % x 为迭代初值，e 为精度y(1)=x ; er=1;k=1;w=7/25;while er>ex=w*x+(1-w)*f(x);k=k+1;y(k)=x;er= abs(y(k)-y(k-1);endfunction y=f(x) % 为子函数y=1/(1+x);(课本P266,对不同的初值及精度运行)二、阅读理解下列程序，并在计算机上调试修改运行.1. jc.mfunction s=jc(n)%n 为正整数s=n;for k=n:1s=(s+1)*k;2. la.m(以课本P32例2、P33例3测试)function y=la(X,Y,x)% X为节点向量，Y为节点相应的函数值向量，x为插值点 n=size(X);y=0;for i=1:nk=1;for j=1:nif j=ik=k*(x-X(j)/(X(i)-X(j); endy=y+k*Y(i); end end三、编程并在计算机上调试修改运行.1 . 编写构造Newton 插值多项式求值程序. 要求：输出差商表及计算结果；以