matlab实验程序分析和总结.docxVIP

实验一 T-1: z1=2*sin(85*pi/180)/(1+exp(2)); x=[2 1+2*i;-.45 5]; %x=[2,1+2i;-0.45,5]; z2=1/2*log(x+sqrt(1+x^2)); 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); t=0:0.5:2.5; z4=(t=0t1).*(t.^2)+(t=1t2).*(t.^2-1)+(t=2t3) .*(t.^2-2*t+1); T-2: A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; c1=A+6*B; c2=A-B+eye(3); c3=A*B; c4=A.*B; c5=A^3; c6=A.^3; c7=A/B; c8=B\A; c9=[A,B]; c10=[A([1,3],:);B^2]; T-3: A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25]; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; C=A*B; D=C(3:end,2:end); whos; T-4-1: m=100:999; n=find(mod(m,21)==0); length(n) T-4-2: ch=ABC123d4e56Fg9; k=find(ch=Ach=Z); ch(k)=[]; ch T-1： E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([2,3]);%S=diag(diag(R)) A=[E,R;O,S];  实验二 B=A^2 C=[E,(R+R*S);O,S^2] %B==C T-2： H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) T-3： % A=fix(10*rand(5)) A=rand(5) H=det(A) Trace=trace(A) Rank=rank(A) Norm=norm(A) T-4： A=[-29,6,18;20,5,12;-8,8,5] [V,D]=eig(A) A*V V*D T-5a： %(1): A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; B1=[0.95,0.67,0.52]; x1=inv(A)*B1 %(2): B2=[0.95,0.67,0.53]; x2=inv(A)*B2 B3=[0.95,0.67,0.54]; x3=inv(A)*B3 %(3): cond(A)%矩阵A 为病态矩阵T-5a： %(1): A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; B1=[0.95,0.67,0.52]; x1=inv(A)*B1 %(2): B2=[0.95,0.67,0.53]; x2=inv(A)*B2 % C=A; % C(9)=1/7; % x3=inv(C)*B2 %(3): cond(A) %矩阵A 为病态矩阵 T-5b： A=hilb(4) A(:,1)=[] A(4,:)=[] B=[0.95,0.67,0.52]; X=inv(A)*B B1=[0.95,0.67,0.53]; X1=inv(A)*B1 N=cond(B) N1=cond(B1) Na=cond(A) %矩阵A 为病态矩阵T-6： A=[1,4,9;16,25,36;49,64,81] B=sqrtm(A) C=sqrt(A) %sqrtm 函数是以矩阵为单位进行计算，sqrt 函数是以矩阵中的元素进行计算实验三 T-1： x=input(请输入x 的值) if x0x~=-3 y=x^2+x-6; elseif x=0x5x~=2x~=3 y=x^2-5*x+6; else y=x^2-x-1; end y T-1a： x=[-5.0 -3.0 1.0 2.0 2.5 3.0 5.0]; y=[]; for i=x if i0i~=-3 y=[y,i^2+i-6]; elseif i=0i5i~=2i~=3 y=[y,i^2-5*i+6]; else y=[y,i^2-i-1]; end end y T-2a： x=input(请输入成绩) if x=90x=100 disp(等级为A) elseif x=80x=89 disp(等级为B) elseif x=70x=79 disp(等级为C) elseif x=60x=69 disp(等级为D) el