Matlab程序設(shè)計與應(yīng)用第二版劉衛(wèi)國

1、實驗一：T1：%第一小題z1=2*sin(85*pi/180)/(1+exp(2)%第二小題x=2,1+2i;-0.45,5;z2%第三小題a=-3.0:0.1:3.0;0.3)+log(0.3+a)/2)%第四題t=0:0.5:2.5z4=(t=0&t=1&t=2&t=A&abdisp(帕薩卡矩陣 P 性能更好 );elseif abdisp(希爾伯特矩陣 H 性能更好 );elsedisp(兩個矩陣性能相同 );endT3:a=1:25;A=reshape(a,5,5)disp(行列式的值： );disp(det(A);disp(矩陣的秩： );disp(rank(A);disp(矩陣的跡

2、： );disp(trace(A);disp(矩陣的范數(shù)： );disp(norm(A);T4:A=-29,6,18;20,5,12;-8,8,5V,D=eig(A)4/27T5:B=0.95,0.67,0.52X1=ABB(3)=0.53X2=ABdisp(系數(shù)矩陣 A 的條件數(shù)： );disp(cond(A);T6:a=1:25;A=reshape(a,5,5)disp(B1=sqrtm(A);B1=sqrtm(A)disp(B2=sqrt(A);B2=sqrt(A)disp(B1*B1);B1*B1disp(B2.*B2);B2.*B2實驗三：T1:x=-5.0:2:5.0for i=1

3、:length(x);5/27if (x(i)=0)&(x(i)5)&(x(i)=2)&(x(i)=3);y(i)=x(i)2-5*x(i)+6;elsey(i)=x(i)2-x(i)-1;endendyT2-if:s=input(please enter the score:);while (s100)disp(the score is not reasonable);s=input(please enter the score:);endif (s=90&s=80&s=70&s=60&s70);6/27disp(D);elsedisp(E);endT2-switch:s=input(ple

4、ase enter the score:);while (s100)disp(the score is not reasonable);s=input(please enter the score:);endswitch fix(s/10)case9,10disp(A);case8disp(B);case7disp(C);case6disp(D);case0,1,2,3,4,5disp(E);end7/27T3:t=input( 請輸入工時 time=);if t120w=120*84+(t-120)*1.15*84;elseif t60w=t*84-700;elsew=84*t;enddis

5、p(應(yīng)發(fā)工資為： );disp(w);T4:a=10+floor(rand(1)*89)b=10+floor(rand(1)*89)s=input( 請輸入 +或-或* 或/,s);while(s=+&s=-&s=*&s=/)disp(輸入的符號錯誤，請重新輸入);s=input(請輸入 +或-或* 或/,s);endswitch scase+c=a+b;8/27case-c=a-b;case*c=a*b;case/c=a/b;endcT5:A=rand(5,6)n=input( 請輸入 n:);while (n5B=A(5,:)elseB=A(n,:)End實驗四：T1:9/27n=inpu

6、t( 請輸入 n:);m=0;for i=1:n;m=m+1/i2;endpii=sqrt(6*m)T1-2:n=input( 請輸入 n:);m=1./(1:n).2;s=sum(m);pii=sqrt(6*s)T2:n=1000;y=0;for i=1:n;y=y+1/(2*i-1);if y=3disp(最大 n 值為： );i-1disp(對應(yīng)的 y 值為： );y-1/(2*i-1)break;10/27endendT3:format longa=input( 請輸入 a:);b=input( 請輸入 b:);x(1)=a/(1.0+b);i=1;x(2)=a/(x(i)+b);wh

7、ile (abs(x(i+1)-x(i)0.00001&i500)i=i+1;x(i+1)=a/(x(i)+b);endr1=(-b+sqrt(b2+4*a)/2;r2=(-b-sqrt(b2+4*a)/2;disp(x,num2str(i+1),=,num2str(x(i+1),8);disp(x,num2str(i),=,num2str(x(i),8);disp(r1=,num2str(r1,8),r2=,num2str(r 2,8);format shortT4:f(1)=1;f(2)=0;f(3)=1;i=4;while(i0)fn=length(find(f0)f0=length(f

8、ind(f=0)T5:j=0;s=0;for i=3:50m=i*(i-1)-1;if isprime(m)=1ms=s+m;j=j+1;endenddisp(親密數(shù)的對數(shù)： ,num2str(j);disp(所有親密素數(shù)之和： ,num2str(s);實驗五：fn1:12/27function f=fn1(n)f=n+10*log(n.2+5);fn2:function f=fn2(n)f=0;for i=1:nf=f+i*(i+1);endfx:function f=fx(x)f=1./(x-2).2+0.1)+1./(x-3).4+0.01);T1:x=input( 請輸入一個復(fù)數(shù)： )

9、;fe=exp(x)flg=log(x)fsin=sin(x)fcos=cos(x)T2:m1=input( 請輸入 m1:);m2=input( 請輸入 m2:);t=input( 請輸入 theta:);A=m1*cos(t),-m1,-sin(t),0;m1*sin(t),0,.13/27cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1B=0,m1*9.8,0,m2*9.8 disp(a1,a2,N1,N2=);disp(AB);T3:j=1;n=1;k=1;for i=10:99%挑選出 1099 中所有的素數(shù)if isprime(i)=1m(j)=i;j=j

10、+1;endendfor t=1:length(m);s(n)=10*rem(m(t),10)+fix(m(t)/10);%挑選出的素數(shù)進行位置交換if isprime(s(n)=1%再挑選出交換絕對素數(shù)r(k)=m(t);k=k+1;endn=n+1;end14/27disp(所有兩位絕對素數(shù)： ,num2str(r)T5:y1=fn1(40)/(fn1(30)+fn1(20)%函數(shù)調(diào)用y2=fn2(40)/(fn2(30)+fn2(20)實驗六：T1:x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.2).*cos(x);plot(x,y);T2:x

11、=linspace(0,2*pi,100);y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,r:,x,y2,y-,x,y3,b-.);subplot(2,2,2);plot(x,y1,r:);subplot(2,2,3);plot(x,y2,y-);subplot(2,2,4);plot(x,y3,b-.);15/27T2-3:x=linspace(0,2*pi,10);y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(3,4,1);bar(x,y1,r);axis(0,7,0,40);subplot(3,4,

12、2);stairs(x,y1,y);axis(0,7,0,40);subplot(3,4,3);stem(x,y1,b);axis(0,7,0,40);subplot(3,4,4);fill(x,y1,g);axis(0,7,0,40);subplot(3,4,5);bar(x,y2,r);axis(-1,7,-1,1);subplot(3,4,6);stairs(x,y2,y);axis(-1,7,-1,1);subplot(3,4,7);stem(x,y2,b);axis(-1,7,-1,1);subplot(3,4,8);fill(x,y2,g);axis(-1,7,-1,1);16/2

13、7subplot(3,4,9);bar(x,y3,r);axis(0,7,-30,50);subplot(3,4,10);stairs(x,y3,y);axis(0,7,-30,50);subplot(3,4,11);stem(x,y3,b);axis(0,7,-30,50);subplot(3,4,12);fill(x,y3,g);axis(0,7,-30,50);T3:x=linspace(-5,5,100);y=;for x0=x;if x0=0;y=y,(x0+sqrt(pi)/exp(2);elseendendsubplot(1,2,1);plot(x,y);for i=1:leng

14、th(x)if x(i)0.5);b=length(k);disp(大于 0.5 的隨機數(shù)個數(shù)占總數(shù)的百分比：);19/27sprintf(%f%, b/300)T2:P=45+50*rand(100,5)disp(最高分及序號： );Y,U=max(P,1)disp(最低分及序號： );Y1,U1=min(P,1)disp(每門課平均分： );mean(P,1)disp(每門課的標準方差： );std(P,0,1)s=sum(P,2);disp(總分最高分： );Y2,U2=max(s,1)disp(總分最低分： );Y3,U3=min(s,1)zcj,xsxh=sort(s,1,desce

15、nd)%descend 降序T3:h=6:2:18;t1=18.0,20.0,22.0,25.0,30.0,28.0,24.0;20/27t2=15.0,19.0,24.0,28.0,34.0,32.0,30.0;t=6.5:2:17.5;T1=interp1(h,t1,t,spline)T2=interp1(h,t2,t,spline)T4:x=1:10:101;y=log10(x);p=polyfit(x,y,5)x1=1:0.01:101;y1=log10(x1);y2=polyval(p,x1);plot(x1,y1,r-,x1,y2,b-.);T5:P1=1,2,4,0,5;P2=0

16、,0,0,1,2;P3=0,0,1,2,3;P4=conv(P2,P3);sp1=length(P1);sp2=length(P4);P=zeros(1,sp2-sp1),P1+P4x=roots(P)A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5;Ax=polyval(P,A)21/27AA=polyvalm(P,A)實驗九：T1:f=inline(det(x x2 x3;1 2*x 3*x2;0 26*x);i=1;x=1;while x=3.01g(i)=f(x);i=i+1;x=x+0.01;enddx=diff(g)/0.01;dx(1)dx(101)dx(leng

17、th(g)-1)T2:f1=inline(sqrt(cos(t.2)+4*sin(2*t).2+1)f2=inline(log(1+x)./(1+x.2)I1=quad(f1,0,2*pi)I2=quad(f2,0,2*pi)22/27T3:A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;B=-4,13,1,11;x1=ABL,U=lu(A);x2=U(LB)Q,R=qr(A);x3=R(QB)fun:function F=fun(X);x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y2+log(z)-7;F(2)=3*x+2y-z3+1;F(

18、3)=x+y+z-5;T5:f=inline(3*x+sin(x)-exp(x);fzero(f,1.5)X=fsolve(fun,11,optimset(Display,off)實驗十：23/27T1:x=sym(6);y=sym(5);z=(x+1)/(sqrt(3+x)-sqrt(y)T2:1syms x y;t=sym(5135);factor(x4-y4)factor(t)T3:syms beta1 beta2 x;simplify(sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2)simplify(4*x2+8*x+3)/(2*x+1)T4:sy

19、ms a b c d e f g h k;p1=0,1,0;1,0,0;0,0,1;p2=1,0,0;0,1,0;1,0,1;A=a,b,c;d,e,f;g,h,k;B=p1*p2*AB1=inv(B)24/27B1*B;tril(B)det(B)T5:syms x t a yf1=(x*(exp(sin(x)+1)-2*(exp(tan(x)-1)/sin(x)3;limit(f1,x,0)%(1)f2=(sqrt(pi)-sqrt(acos(x)/sqrt(x+1);limit(f2,x,-1,right)%(2)f3=(1-cos(2*x)/x;diff(f3,x,1)%(3)diff(f3,x,2)A=ax,t3;t*cos(x),log(x);diff(A,x,1)&(4)diff(A,t,2)diff(diff(A,x,1),t,1)f5=(