matlab 部分答案,matlab第三章课后部分答案.doc

matlab第三章课后部分答案.doc

习题三3-2从键盘输入一个三位整数，将它反向输出。如输入639，输出为936程序如下：m=( 请输入一个三位整数： );m1=fix(m/100);%求m的百位整数m2=rem(fix(m/10),10);%求m的十位数字m3=rem(m,10);%求m的个位数字m=m3*100+m2*10+m1%反向输出m3-3输入一个百分制成绩，要求输出成绩等级A,B,C,D,E。其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。要求：(1)分别用if语句和switch语句实现。(2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。程序如下：(1)if语句c=( 请输入成绩： );ifc>=90fori=2:20ifn>当n等于100时：1.635当n等于1000时：1.6439当n等于10000时：1.6448(2)1-1/3+1/5-1/7+···(=pi/4)程序如下：clearalln=1:100;b=1./(4.*n-3)-1./(4.*n-1);c=sum(b);disp([ 当n等于100时： ,num2str(c)])n=1:1000;b=1./(4.*n-3)-1./(4.*n-1);c=sum(b);disp([ 当n等于1000时： ,num2str(c)])n=1:10000;b=1./(4.*n-3)-1./(4.*n-1);c=sum(b);disp([ 当n等于10000时： ,num2str(c)])结果：>>当n等于100时：0.78415当n等于1000时：0.78527当n等于10000时：0.78539(3)1/4+1/16+1/64+···+1/4^n+···(=1/3)程序如下：clearalln=1:100;b=1./(4.^n);c=sum(b);disp([ 当n等于100时： ,num2str(c)])n=1:1000;b=1./(4.^n);c=sum(b);disp([ 当n等于1000时： ,num2str(c)])n=1:10000;b=1./(4.^n);c=sum(b);disp([ 当n等于10000时： ,num2str(c)])结果：>>当n等于100时：0.33333当n等于1000时：0.33333当n等于10000时：0.33333(4)···((2n)(2n)/(2n-1)(2n+1))···(=pi/2)程序如下：clearalln=1:100;b=((2*n).*(2*n))./((2*n-1).*(2*n+1));c=prod(b);disp([ 当n等于100时： ,num2str(c)])n=1:1000;b=((2*n).*(2*n))./((2*n-1).*(2*n+1));c=prod(b);disp([ 当n等于1000时： ,num2str(c)])n=1:10000;b=((2*n).*(2*n))./((2*n-1).*(2*n+1));c=prod(b);disp([ 当n等于10000时： ,num2str(c)])结果：>>当n等于100时：1.5669当n等于1000时：1.5704当n等于10000时：1.57083-9若两个连续自然数的乘积减1是素数，则称这两个连续自然数是亲密数对，该素数是亲密素数。例如，2*3-1=5，由于5是素数，所以2和3是亲密数对，5是亲密素数。求[2,50]区间内：(1)亲密数对的对数。(2)与上述亲密数对对应的所有亲密素数之和。程序如下：clearallp=0;l=0;H=[];=2:49a=m;b=m+1;c=a*b-1;fork=1:cifrem(c,k)==0l=l+1;endendifl==2disp([ 亲密数对是 , ( ,num2str(a), , ,num2str(b), ) ])p=p+1;H=[H,c];endl=0;endsum=sum(H);disp([ 亲密数对的个数是： ,num2str(p)])disp([ 亲密数对的和是： ,num2str(sum)])3-10根据y=1+1/3+1/5+···+1/(2n-1)，求：(1)y3y=0;while(y<3)x=1/(2*n-1);n=1+n;y=y+x;enddisp([ y<3时n的最大值是： ,num2str(n-2)])disp([ 相应的y值是： ,num2str(y-x)])3-12编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在命令文件中调用该函数。程序如下：mul函数%用于求两个矩阵的乘积和点乘function[x,y]=mul(a,b)x=a*b;y=a.*b;main%主函数a=( 请输入一个矩阵： );b=( 请再输入一个矩阵：(注意：两矩阵要可以相乘) );[x,y]=mull(a,b);disp( 这两个矩阵的乘积为： );xdisp( 这两个矩阵的点乘为： );y3-14先用函数的递归调用定义一个函数文件求∑I，然后调用该函数文件求∑k+∑k^2+∑1/k程序如下：定义一个求和函数：functionsum=add(n,m)ifn<1sum=0;elsesum=n^m+add(n-1,m);end主函数：clearally=add(100,1)+add(50,2)+add(10,-1);disp(y)