MATLAB程序设计与应用(第二版)第4章课后题答案.doc

1. a=input(请输入一个4位数:); while (a1000|a9999) a=input(输入错误，请重新输入一个4位数:); end b=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(d,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. 逻辑表达式法： a=input(请输入a: ); b=input(请输入b: ); c=input(请输入c: ); x=0.5:1:5.5; x1=(x=0.5x1.5); x2=(x=1.5x3.5); x3=(x=3.5x=5.5); y1=a.*(x.^2)+b.*x+c; y2=a*(sin(b)^c)+x; y3=log(abs(b+c./x)); y=y1.*x1+y1.*x2+y3.*x3; disp(y) if语句法： a=input(请输入a: ); b=input(请输入b: ); c=input(请输入c: ); for x=0.5:1:5.5 if x=0.5 x1.5 y=a.*(x.^2)+b.*x+c elseif x=1.5 x3.5 y=a*(sin(b)^c)+x elseif x=3.5 x5.5 y=log(abs(b+c./x)) end end switch语句法： a=input(请输入a: ); b=input(请输入b: ); c=input(请输入c: ); for x=0.5:1:5.5 switch 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,8,9,10} y=log(abs(b+c./x)); end disp(y) end 3. x=fix(rand(1,20)*89)+10; x1=mean(x); n=find(rem(x,2)==0 xx1); disp([小于平均数的偶数是：,num2str(x(n))]); 4.(1) A=input(请输入20个数的一个行向量：); a=A(1); b=A(1); for m=A if a=m a=m; elseif b=m b=m; end end disp([最小数是：,num2str(a)]) disp([最大数是：,num2str(b)]) (2) A=input(请输入20个数的一个行向量：); maxval=max(A) minval=min(A) 5. s=0; for a=0:63 c=2^a; s=s+c; end disp([2的0次方到63次方的和是：,num2str(s)]) k=0:63 n=2.^k s=sum(n) 6.（1） sum1=0; for n=1:100 x=(-1)^(n+1)*(1/n); sum1=sum1+x; end disp([当n取100时： sum=,num2str(sum1)]) sum2=0; for n=1:1000 x=(-1)^(n+1)*(1/n); sum2=sum2+x; end disp([当n取1000时： sum=,num2str(sum2)]) sum3=0; for n=1:10000 x=(-1)^(n+1)*(1/n); sum3=sum3+x; end disp([当n取10000时：sum=,num2str(sum3)]) （2） sum1=0; n1=0; for n=1:2:100 x=(-1)^n1*(1/n); sum1=sum1+x; n1=n1+1; end disp([当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; end disp([当n取1000时： sum=,num2str(sum2)]) sum3=0; n3=0; for n=1:2:1000