1、MATLAB 程序设计与应用刘卫国高等教育出版社第4章课后答案1.a=input(请输入一个4位数:);while (a9999)a=input(输入错误,请重新输入一个4位数:);endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*b+100*c+10*d+e;disp(加密后:,num2str(a)2
2、.逻辑表达式法:a=input(请输入a: );b=input(请输入b: );c=input(请输入c: );x=0.5:1:5.5;x1=(x=0.5&x=1.5&x=3.5&x=0.5 & x=1.5 & x=3.5 & x5.5 y=log(abs(b+c./x)endendswitch语句法:a=input(请输入a: );b=input(请输入b: );c=input(请输入c: );for x=0.5:1:5.5switch floor(x/0.5) case 1,2 y=a.*(x.2)+b.*x+c; case 3,4,5,6 y=a*(sin(b)c)+x; case 7,
3、8,9,10 y=log(abs(b+c./x); enddisp(y)end3.x=fix(rand(1,20)*89)+10;x1=mean(x);n=find(rem(x,2)=0 & x=m a=m; elseif b=m b=m; endenddisp(最小数是:,num2str(a)disp(最大数是:,num2str(b) (2)A=input(请输入20个数的一个行向量:);maxval=max(A)minval=min(A)5.s=0;for a=0:63 c=2a; s=s+c;enddisp(2的0次方到63次方的和是:,num2str(s) k=0:63n=2.ks=s
4、um(n)6.(1)sum1=0;for n=1:100 x=(-1)(n+1)*(1/n); sum1=sum1+x;enddisp(当n取100时: sum=,num2str(sum1) sum2=0;for n=1:1000 x=(-1)(n+1)*(1/n); sum2=sum2+x;enddisp(当n取1000时: sum=,num2str(sum2) sum3=0;for n=1:10000x=(-1)(n+1)*(1/n);sum3=sum3+x;enddisp(当n取10000时:sum=,num2str(sum3)(2)sum1=0;n1=0;for n=1:2:100 x
5、=(-1)n1*(1/n); sum1=sum1+x; n1=n1+1;enddisp(当n取100时: sum=,num2str(sum1) sum2=0;n2=0;for n=1:2:1000 x=(-1)n2*(1/n); sum2=sum2+x; n2=n2+1;enddisp(当n取1000时: sum=,num2str(sum2) sum3=0;n3=0;for n=1:2:10000x=(-1)n3*(1/n);sum3=sum3+x;n3=n3+1;enddisp(当n取10000时:sum=,num2str(sum3)(3)sum1=0;for n=1:100 x=1/(4n
6、); sum1=sum1+x;enddisp(当n取100时: sum=,num2str(sum1) sum2=0;for n=1:1000 x=1/(4n); sum2=sum2+x;enddisp(当n取1000时: sum=,num2str(sum2) sum3=0;for n=1:10000 x=1/(4n); sum3=sum3+x;enddisp(当n取10000时:sum=,num2str(sum3)(4)sum1=1;for n=1:100 x=4*n*n/(2*n-1)/(2*n+1); sum1=sum1*x;enddisp(当n取100时: sum=,num2str(su
7、m1) sum2=1;for n=1:1000 x=4*n*n/(2*n-1)/(2*n+1); sum2=sum2*x;enddisp(当n取1000时: sum=,num2str(sum2) sum3=1;for n=1:10000 x=4*n*n/(2*n-1)/(2*n+1); sum3=sum3*x;enddisp(当n取10000时:sum=,num2str(sum3)7.函数文件function f=fibnacci(n)if n=1 | n=2 f=1;else f=fibnacci(n-1)+fibnacci(n-2);end命令文件:shulie=;for k=1:n sh
8、ulie=shulie fibnacci(k);endshulie8.function f1,f2=juzhenji(x1,x2)f1=x1*x2;f2=x1.*x2;命令文件:clear alla=input(请输入一个矩阵:);b=input(请再输入一个矩阵:(注意:两矩阵要可以相乘);f1,f2=juzhenji(a,b);disp(f1)disp(f2)9.function sum=qiuhe(n,m)if n=1sum=0;elsesum=nm+qiuhe(n-1,m);end命令文件:clear ally=qiuhe(100,1)+qiuhe(50,2)+qiuhe(10,-1);disp(y)10.s=0;a=12,13,14;15,16,17;18,19,20;21,22,23;for k=a for j=1:4 if rem(k(j),2)=0 s=s+k(j); end endendss = 108(2)global xx=1:2:5;y=2:2:6;sub(y);xyfunction fun=sub(z)global xz=3*x;x=x+z;x = 4 12 20y = 2 4 6
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1