MATLAB课后试验答案

第四章作业 4.1、 a=( 请输入一个 4 位数: ); while (a9999) a=( 输入错误，请重新输入一个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(c,10); e=rem(e,10); g=b;b=d;d=g; g=c;c=e;e=g; a=1000*d+100*e+10*b+c; disp([ 加密后： ,num2str(a)]) 4.2、 a=( 请输入 a: ); b=( 请输入 b: ); c=( 请输入 c: ); x=0.5:1:5.5; x1=(x=0.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) 4.3、 x=fix(rand(1,20)*89)+10; x1=fix(sum(x)/20); disp([ 平均数是： ,num2str(x1)]) m=(rem(x,2)==0 elseif b=m b=m; end end disp([ 最小数是： ,num2str(a)]) disp([ 最大数是： ,num2str(b)]) 4.5、 s=0; a=0; for b=1:64 c=2^a; a=a+1; s=s+c; end disp([ 2 的 0 次方到 63 次方的和是： ,num2str(s)]) 4.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)]) 4.8、clear all a=( 请输入一个矩阵： ); b=( 请再输入一个矩阵：(注意：两矩阵要可以相乘) ); [f1,f2]=juzhenji(a,b); 5.1 x=-2*pi:pi/1000:2*pi; y=100./(1+x.^2); plot(x,y) x=-2*pi:pi/1000:2*pi; y=(1/2*pi)*exp(-0.5*x.^2); plot(x,y) t=-2*pi:pi/1000:2*pi; x=sin(t); y=cos(t); plot(x,y) t=-2*pi:pi/1000:2*pi; x=t.^2; y=5*t.^3; plot(x,y, r: ) 实验 3 1、clear m=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]; for x=m if x=0 disp([ 您本月工资是： ,num2str(a)]) elseif y5|n0 b=a(5:5,:); disp([ 输入错误，程序输出最后一行的元素： ,num2str(b)]) else c=a(n:n,:); disp([ 输出结果是： ,num2str(c)]) end 实验 4 1、n=(“请输入 n 值：”)； s1=0; for k=1:n; x=1/(k^2); s1=s1+x; pi=sqrt(6*s1); 2、clear n=1; y=0; while (y3) x=1/(2*n-1); n=1+n; y=y+x; end disp([ y3 时 n 的最大值是： ,num2str(n-2)]) disp([ 相应的的 y 值是： ,num2str(y-x)]) 5、clear p=0; l=0; H=[]; for m=2:49 a=m; b=m+1; c=a*b-1; for k=1:c if rem(c,k)==0 l=l+1; end end if l==2 disp([ 亲密数对是 , ( ,num2str(a), , ,num2str(b), ) ]) p=p+1; H=[H,c]; end l=0; end sum=sum(H); disp([ 亲密数对的个数是： ,num2str(p)]) disp([ 亲密数对的和是： ,num2str(sum)]) 实验 6 1、clear all x=linspace(0,2*pi,101); y=(0.5+3*sin(x)/(1+x.^2)).*cos(x); plot(x,y);title( y=[0.5+3sin(x)/(1+x^{2})]cos(x)的图像如下： ) 2/1 clear all x=linspace(-2*pi,2*pi,500); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1, b- ,x,y2, r-- ,x,y3, k: ); 2//3 clear all x=linspace(-2*pi,2*pi,500); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(3,3,1); bar(x,y1);%条形图 subplot(3,3,2); stairs(x,y1);%阶梯图 subplot(3,3,3); stem(x,y1);%杆图 subplot(3,3,4); bar(x,y2); subplot(3,3,5); stairs(x,y2); subplot(3,3,6); stem(x,y2); subplot(3,3,7); bar(x,y3); subplot(3,3,8); stairs(x,y3); subplot(3,3,9); stem(x,y3); 3、clear all t=linspace(-5,5,200); H=[]; for x=t if x0 y1=0.5*log(x+sqrt(1+x.^2)); end H=[H,y1]; end plot(t,H, r: ); 4、clear all a=( 请输入一个数 a： ); b=( 请输入一个数 b： ); n=( 请输入一个数 n： ); theta=linspace(0,2*pi,500); rho=a*sin(b+n.*theta); polar(theta,rho); 5、clear all x1=linspace(-5,5,21); y1=linspace(0,10,31); [x,y]=meshgrid(x1,y1); z=cos(x).*cos(y).*exp(-