《MATLAB程序设计与应用(第二版)实验参考答案》由会员分享，可在线阅读，更多相关《MATLAB程序设计与应用(第二版)实验参考答案（8页珍藏版）》请在金锄头文库上搜索。
1、MATLAB程序设计与应用第二版实验参考答案%实验一 MATLAB运算根底%第一题%1z1=2*sin(85*pi/180)/(1+exp(2)%2x=2,5;z2=0.5*log(x+sqrt(1+x.2)%(3);z3=(exp(0.3*a)-exp(-0.3*a)/2.*sin(a+0.3)+log(0.3+a)/2)%(4);z4=t.2.*(t=0&t=1&t=2&t=A&A=Z);A(k)=%实验二 MATLAB矩阵分析与处理%第一题E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag(2,3);A=E,R;O,S;A2B=E,(R+R*S);O,S2%
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,5V,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):3;x=inv(A)*B%(3):cond(A)%第五题方法二:A=hilb(4)A(
3、:,1)=A(4,:)=B=0.95,0.67,0.52;X=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,81B=sqrtm(A)C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算%实验三 选择程序结构设计%第一题程序一x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0;y=; %建立存放所有y值的矩阵for x0=x if x0=0&x05&x0=2&x0=3 y=y,x
4、0*x0-5*x0+6; else y=y,x0*x0-x0-1; endendx %输出所有xy %输出所有y%第一题程序二x=-5,-3,1,2,2.5,3,5;y=;for a=1:7if x(a)=0&x(a)100|x0disp(您输入的成绩不是百分制成绩，请重新输入。);elseif x=90disp(A);elseif x=80disp(B);elseif x=70disp(C);elseif x60disp(D);elsedisp(E);endend%第二题程序二：s=input(请输入一个成绩0分到100分之间：); %s用于存放成绩while 1 %判断输入成绩的合理性 i
5、f s100 disp(输入的成绩需在0到100之间，请重新输入：) s=input(请输入一个成绩0分到100分之间：); else break; endendswitch fix(s/10) %对成绩做出等级判断 case 9,10 disp(A) case 8 disp(B) case 7 disp(C) case 6 disp(D) otherwise disp(E)end%第三题n=input(请输入员工工号：);h=input(该员工工作时数是：);if h120x=(h-120)*84*(1+0.15)+120*84;elseif h5 disp(超出了矩阵的行数，矩阵的最后一行
6、为：) a(5,:)else disp(矩阵的第,num2str(n),行为：) a(n,:)end%实验四 循环结构程序设计%第一题程序一s=0;n=input(n=?);for i=1:n s=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;while(y1e-5 Xn=Xn1; Xn1=a/(b+Xn); n=n+1; if n=500 break; endendnXn1r1=(-b+sqrt(b*b+4*a)/2r2=(-b-sqrt(b*b+4*a)/2
7、%第四题for i=1:100if i=1f(i)=1;elseif i=2f(i)=0;elseif i=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;for i=2:49b=i*(i+1)-1;m=fix(sqrt(b);for j=2:mif rem(b,j)=0breakendendif j=mn=n+1;s=s+b;endendns%实验五 函数文件%第一题y=exp(x),log(x)
8、,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=
9、input(theta=);a1,a2,N1,N2=shiyanwu2(m1,m2,theta)%第二题程序二function X=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;%在命令窗口调用该函数文件：X=mat2(1,1,60)%第三题function flag=mat3(x)flag=1;for i=2:sqrt(x)if rem(x,i)=0 flag=0; br
