您好,欢迎访问三七文档
当前位置:首页 > 商业/管理/HR > 管理学资料 > MATLAB程序设计与应用(第二版)部分实验参考答案
1MATLAB程序设计与应用(第二版)实验参考答案实验一MATLAB运算基础第一题(1)z1=2*sin(85*pi/180)/(1+exp(2))(2)x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))(3)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)(4)t=0:0.5:2.5;z4=t.^2.*(t=0&t1)+(t.^2-1).*(t=1&t2)+(t.^2-2*t+1).*(t=2&t3)第二题A=[1234-4;34787;3657];B=[13-1;203;3-27];A+6*BA-B+eye(size(A))A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]第三题A=[12345;678910;1112131415;1617181920;2122232425]B=[3016;17-69;023-4;970;41311]C=A*BF=size(C)D=C(F(1)-2:F(1),F(2)-1:F(2))whos第四题(1):A=100:999;B=rem(A,21);C=length(find(B==0))(2):A='lsdhKSDLKklsdkl';k=find(A='A'&A='Z');2A(k)=[]实验二MATLAB矩阵分析与处理第一题E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([2,3]);A=[E,R;O,S];A^2B=[E,(R+R*S);O,S^2]第二题H=hilb(5)P=pascal(5)Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)第三题:A=fix(10*rand(5))H=det(A)Trace=trace(A)Rank=rank(A)Norm=norm(A)第四题:A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)数学意义略第五题方法一:(1):A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];b=[0.95,0.67,0.52]';x=inv(A)*b(2):B=[0.95,0.67,0.53]';x=inv(A)*B(3):cond(A)方法二:A=hilb(4)A(:,1)=[]A(4,:)=[]B=[0.95,0.67,0.52]';3X=inv(A)*BB1=[0.95,0.67,0.53]';X1=inv(A)*B1N=cond(B)N1=cond(B1)Na=cond(A)%矩阵A为病态矩阵第六题A=[1,4,9;16,25,36;49,64,81]B=sqrtm(A)C=sqrt(A)%sqrtm函数是以矩阵为单位进行计算,sqrt函数是以矩阵中的元素进行计算实验三选择程序结构设计第一题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[];%建立存放所有y值的矩阵forx0=xifx00&x0~=-3y=[y,x0*x0+x0-6];elseifx0=0&x05&x0~=2&x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y第一题程序二x=[-5,-3,1,2,2.5,3,5];y=[];fora=1:7ifx(a)0&x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseifx(a)=0&x(a)5&x(a)~=2&x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend第二题程序一:x=input('请输入一个百分制成绩:');ifx100|x0disp('您输入的成绩不是百分制成绩,请重新输入。');elseifx=100&x=904disp('A');elseifx=89&x=80disp('B');elseifx=79&x=70disp('C');elseifx=69&x60disp('D');elsedisp('E');endend第二题程序二:s=input('请输入一个成绩(0分到100分之间):');%s用于存放成绩while1%判断输入成绩的合理性ifs0|s100disp('输入的成绩需在0到100之间,请重新输入:')s=input('请输入一个成绩(0分到100分之间):');elsebreak;endendswitchfix(s/10)%对成绩做出等级判断case{9,10}disp('A')case8disp('B')case7disp('C')case6disp('D')otherwisedisp('E')end第三题n=input('请输入员工工号:');h=input('该员工工作时数是:');ifh120x=(h-120)*84*(1+0.15)+120*84;elseifh60x=h*84-700;elsex=h*84;enddisp([num2str(n),'号员工','的应发工资为',num2str(x)]);5第四题(还可以用switch语句实现)a=fix(10+(99-10)*rand(1,2))%产生两个随机整数x=a(1);y=a(2);t=input('请输入运算符号:','s');ift=='+'z=x+y;elseift=='-'z=x-y;elseift=='*'z=x*y;elseift=='/'z=x/y;enddisp([num2str(x),t,num2str(y),'=',num2str(z)])%输出运算结果第五题a=rand(5,6)%产生5x6的随机矩阵n=input('请输入您要输出矩阵的第几行:');ifn5disp('超出了矩阵的行数,矩阵的最后一行为:')a(5,:)elsedisp(['矩阵的第',num2str(n),'行为:'])a(n,:)end实验四循环结构程序设计第一题程序一s=0;n=input('n=?');fori=1:ns=s+1/i/i;endPI=sqrt(6*s)pi第一题程序二n=input('n=?');a=1:n;b=1./a.^2;PI=sqrt(6*sum(b))pi第二题y=0;n=1;6while(y3)y=y+1/(2*n-1);n=n+1;endy=y-1/(2*(n-1)-1)n=n-2第三题a=input('a=?');b=input('b=?');Xn=1;Xn1=a/(b+Xn);n=0;whileabs(Xn1-Xn)1e-5Xn=Xn1;Xn1=a/(b+Xn);n=n+1;ifn==500break;endendnXn1r1=(-b+sqrt(b*b+4*a))/2r2=(-b-sqrt(b*b+4*a))/2第四题fori=1:100ifi==1f(i)=1;elseifi==2f(i)=0;elseifi==3f(i)=1;elsef(i)=f(i-1)-2*f(i-2)+f(i-3);endendmax(f)min(f)sum(f)length(find(f0))length(find(f==0))length(find(f0))第五题:s=0;n=0;7fori=2:49b=i*(i+1)-1;m=fix(sqrt(b));forj=2:mifrem(b,j)==0breakendendifj==mn=n+1;s=s+b;endendns实验五函数文件第一题functiony=mat1(x)%建立函数文件mat1.my=[exp(x),log(x),sin(x),cos(x)];%在命令窗口调用上述函数文件:y=mat1(1+i)第二题程序一function[a,b,N,M]=shiyanwu2(m,n,t)A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];B=[0,9.8*m,0,9.8*n];C=inv(A)*B';a=C(1);b=C(2);N=C(3);M=C(4);在命令窗口调用该函数文件:m1=input('m1=');m2=input('m2=');theta=input('theta=');[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)第二题程序二functionX=mat2(m1,m2,t)g=9.8;A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0;m1*sin(t*pi/180),0,cos(t*pi/180),0;0,m2,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];B=[0;m1*g;0;m2*g];X=inv(A)*B;8在命令窗口调用该函数文件:X=mat2(1,1,60)第三题functionflag=mat3(x)flag=1;fori=2:sqrt(x)ifrem(x,i)==0flag=0;break;endend在命令窗口调用该函数文件:fori=10:99j=10*rem(i,10)+fix(i/10);ifmat3(i)&mat3(j)disp(i)endend第四题functiony=fx(x)y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);在命令窗口调用该函数文件:y=fx(2)a=[1,2;3,4];y=fx(a)第五题(1)functionf1=mat5(n)f1=n+10*log(n*n+5);在命令窗口中调用该函数文件:y=mat5(40)/(mat5(30)+mat5(20))(2)方法一functionf2=mat6(n)f2=0;fori=1:nf2=f2+i*(i+1);end在命令窗口中调用该函数文件如:y=mat6(40)/(mat6(30)+mat6(20))(2)方法二functionf2=mat7(n)i=1:n;m=i.*(i+1);f2=sum(m);9end在命令窗口中调用该函数文件如:y=mat7(40)/(mat7(30)+mat7(20))实验六高层绘图操作第一题:x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);plot(x,y)第二题:(1)x=linspace(-2*pi,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'b-',x,y2,'r:',x,y3,'y--');text(4,16,'\leftarrowy1=x^2');text(6*pi/4,-1,'\downarrowy2=cos(2*x)');text(-1.5*pi,-2.25*pi*pi,'\uparrowy3=y1*y2');(2)x=linspace(-2*pi,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);%分区plot(x,y1);title('y1=x^2');%设置标题subplot(1,3,2);plot(x,y2);title('y2=cos(2*x)');subplot(1,3,3);plot(x,y3);title('y3=x^2*cos(2*x)');(3)x=linspace(-2*pi,2*pi,20);y1=x.^2;subplot(2,2,1);%分区bar(x,y1);title('y1=x^2的条形图');%设置标题subplot(2,2,2);stairs(x,y1);titl
本文标题:MATLAB程序设计与应用(第二版)部分实验参考答案
链接地址:https://www.777doc.com/doc-2882209 .html