数学实验 matlab第三次作业(修订版).docVIP

作业1 猜数游戏。首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。用户最多可以猜7次。 M=ceil(rand*100) for i=1:7 n=input(Input a number between 1 to 100:); if n==M disp(You Won);break; elseif nM disp(High); elseif nM disp(Low); end i=i+1; end M = 82 Input a number between 1 to 100:23 Low Input a number between 1 to 100:89 High Input a number between 1 to 100:45 Low Input a number between 1 to 100:37 Low Input a number between 1 to 100:89 High Input a number between 1 to 100:76 Low Input a number between 1 to 100:82 You Won 作业2 用筛选法求某自然数范围内的全部素数。 素数是大于1，且除了1和它本身以外，不能被其他任何整数所整除的整数。用筛选法求素数的基本思想是：要找出2~m之间的全部素数，首先在2~m中划去2的倍数(不包括2)，然后划去3的倍数(不包括3)，由于4已被划去，再找5的倍数 (不包括5)，…，直到再划去不超过的数的倍数，剩下的数都是素数。 m=input(Input a number:); for i=2:m for p=2:i if rem(i,p)==0 break; else p=p+1; end if p==i disp(i); end end i=i+1; end Input a number:12 3 5 7 11 作业3 求函数f(x)在[a,b]上的定积分，其几何意义就是求曲线y=f(x)与直线x=a，x=b，y=0所围成的曲边梯形的面积。为了求得曲边梯形面积，先将积分区间[a,b]分成n等分，每个区间的宽度为h=(b-a)/n，对应地将曲边梯形分成n等分，每个小部分即是一个小曲边梯形。近似求出每个小曲边梯形面积，然后将n个小曲边梯形的面积加起来，就得到总面积，即定积分的近似值。近似地求每个小曲边梯形的面积，常用的方法有：矩形法、梯形法以及辛普生法等。 function dfx f=inline(x*x,x); a=input(输入下限：); b=input(输入上限：); n=10000; h=(b-a)/n; s=0; t=0; for i=0:n t=(f(a+h*i)+f(a+h*(i+1)))*h/2; s=s+t; end R=vpa(s,4); disp(R); 函数下限:1 函数上限:2 2.334 作业4 Fibonacci数列定义如下： f1=1 f2=1 fn=fn-1+fn-2 (n2) 求Fibonacci数列的第20项。 n=input(Enter a number :); f1=1; f2=1; i=2; while(in) f1=f1+f2; f2=f1+f2; i=i+2; end disp(f2); Enter a number :20 6765 作业5 根据矩阵指数的幂级数展开式求矩阵指数。 A=input(Enter a matrix: ); n=input(input a number); t=1; s=eye(size(A)) for i=1:n t=t*i; p=(A^i)/t; s=s+p; end disp(s); B=expm(A); disp(B); Enter a matrix:[1 2 3;4 5 6;7 8 9] input a number100 s = 1 0 0 0 1 0 0 0 1 1.0e+06 * 1.1189 1.3748 1.6307 2.5339 3.1