matlab第三章习题解答

3-1>> N=[1 2 3 4 5]N =1 2 3 4 5(1)>> N1=2*NN1 =2 4 6 8 10(2)>> N2=N/2N2 =0.5000 1.0000 1.5000 2.0000 2.5000 (3)>> N3=1./NN3 =1.0000 0.5000 0.3333 0.2500 0.2000 (4)>> N4=1./(N.*N)N4 =1.0000 0.2500 0.1111 0.0625 0.0400 3-2m=input('请输入一个三位整数');i=fix(m/100);j=mod(fix(m/10),10);k=mod(m,10);n=100*k+10*j+i;disp('反向输出的数为');disp(n);3-3用if来实现score=input('请输入成绩');if score>100|score<0disp('成绩必须>=0且<=100');elseif score>=90&score<=100disp('A');elseif score>=80&score<=89disp('B');elseif score>=70&score<=79disp('C');elseif score>=60&score<=69disp('D');elsedisp('E');end用switch来实现score=input('请输入成绩');switch scorecase num2cell(90:100)disp('A');case num2cell(80:89)disp('B');case num2cell(70:79)disp('C');case num2cell(60:69)disp('D');case num2cell(0:59)%case{0,1,2,3,4,5}也可以 disp('E');otherwisedisp('成绩必须>=0且<=100');end注意:以下程序是错的score=input('请输入成绩');switch fix(score/10)case {9,10}disp('A');case {8}disp('B');case {7}disp('C');case {6}disp('D');case num2cell(0:5)%case{0,1,2,3,4,5}也可以disp('E');otherwisedisp('成绩必须>=0且<=100');end因为在这个程序中101到109都将对应被当作A3-4clear;A=input('请输入一个5行六列的矩阵');n=input('请输入矩阵的行数n');if n>5B=A(5,:);disp('行数大于5了');disp(B);elseB=A(n,:);disp(B);end3-5方法1:x=90*rand(1,20)+10y=fix(x)average_num=mean(y)index1=find(rem(y,2)==0&y>average_num)result=y(index1)方法2:sum=0;for i=1:20p(i)=fix(90*rand+10);%用rand函数产生20个随机数endpfor i=1:20sum=sum+p(i);%累加数组p的所有的元素endaverage_num=sum/20%数组p的所有的元素的平均数for i=1:20if(p(i)<average_num&rem(p(i),2)==0)disp(p(i))endend方法3:sum=0;for i=1:20p(i)=fix(90*rand+10);%用rand函数产生20个随机数endpfor i=1:20sum=sum+p(i);%累加数组p的所有的元素endaverage_num=sum/20%数组p的所有的元素的平均数b=find(p<average_num&rem(p,2)==0)%依次访问p中的每一个元素，如果它满足条件小于average_num %且是偶数，那么记录下它的位置p(b)3-6用循环来实现x=input('请输入20个数');n=20;maxnum=x(1);minnum=x(1);for i=1:20if x(i)>maxnummaxnum=x(i);endif x(i)<minnumminnum=x(i);enddisp('最大的数为');disp(maxnum);disp('最小的数为');disp(minnum);合法的输入输出的例子>> hw3_6a请输入20个数[ 1 2 3 4 5 6 7 8 9 10 20 19 18 17 16 15 14 13 12 11]最大的数为20最小的数为1用max和min来实现x=input('请输入20个数');n=20;maxnum=max(x);minnum=min(x);disp('最大的数为');disp(maxnum);disp('最小的数为');disp(minnum);3-7解:用顺序结构a=-3.0:0.1:3.0;y=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)用循环结构a=-3.0:0.1:3.0;m=size(a);n=m(2);%矩阵a的列数for i=1:ny(i)=(exp(0.3*a(i))-exp(-0.3*a(i)))/2*sin(a(i)+0.3)+log((0.3+a(i))/2); enddisp(y(i));end3-8(1)用循环语句实现sum1=0;n=10000;%n=100,n=1000只需修改代码中n的值得for i=1:nsum1=sum1+1/(i*i);enddisp(sum1);用向量运算n=10000;%n=100,n=1000只需修改代码中n的值得a=1:n;b=a.*a;c=1./b;sum1=sum(c);disp(sum1);直接写如下的代码也对n=10000;%n=100,n=1000只需修改代码中n的值得a=1:n;sum1=sum(1./(a.*a));disp(sum1);(2) 用循环语句实现sum1=0;n=10000;for i=1:nsum1=sum1+(-1)^(i+1)/(2*i-1);enddisp(sum1);用向量运算sum1=0;a=1:n;b=(-1).^(a+1);c=2*a-1;d=b.*(1./c);sum1=sum(d);disp(sum1);(3) 用循环语句实现sum1=0;n=10000;for i=1:nsum1=sum1+1/(4^i);enddisp(sum1);用向量运算sum1=0;n=10000;a=1:n;b=4.^(a);c=1./b;sum1=sum(c);disp(sum1);(4)用循环语句实现sum1=1;n=10000;for i=1:nsum1=sum1*((2*i)*(2*i)/((2*i-1)*(2*i+1)));enddisp(sum1);用向量运算(注意：因为这里是连乘,不是用sum而是用prod,其用法见教材133页)n=10000;a=1:n;c=2*a-1;d=2*a+1;e=b.*b./(c.*d);prod1=prod(e);disp(prod1);3-9clear;counter=0;sum=0;for i=2:49u=i*(i+1)-1;if (is_prime(u)==1)counter=counter+1;sum=sum+u;endendcountersum其中is_prime函数对应的m文件的内容为function f=is_prime(m)f=1;u=fix(sqrt(m));for(i=2:u)if(rem(m,i)==0)f=0;break;endendend3-10y=0;i=1;while y<3y=y+1/(2*i-1);endn=i-2;y=y-1/(2*(i-1)-1);disp(n);disp(y);3-11(1)function result=f3_11a(n)result=n+10*log(n.*n+5);y=f3_11a(40)/(f3_11a(30)+f3_11a(20))(2)function result=f3_11b(n)result=0;for i=1:nresult=result+i*(i+1);endy=f3_11b(40)/(f3_11b (30)+f3_11b (20))3-12function [chengji,diancheng]=matrixprod(A,B) chengji=A*B;diancheng=A.*B;例如:a=[1 2;3 4];b=[1 0 ;0 1]调用函数>> [c,d]=matrixprod(a,b)c =1 23 4d =0 4。