Matlab程序设计与应用第二版刘卫国

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性能更好);elseifabdisp(希尔伯特矩阵H性能更好);elsedisp(两个矩阵性能相同);endT3:a=1:25;A=reshape(a,5,5)disp(行列式的值：);disp(det(A);disp(矩阵的秩：);disp(rank(A);disp(矩阵的迹：);disp(tr

2、ace(A);disp(矩阵的范数：);disp(norm(A);T4:A=-29,6,18;20,5,12;-8,8,5V,D=eig(A)T5:B=0.95,0.67,0.52X1=ABB(3)=0.53X2=ABdisp(系数矩阵A的条件数：);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.0fori=1:length(x);if(x(i)=0

3、)&(x(i)5)&(x(i)=2)&(x(i)=3);y(i)=x(i2-5*x(i)+6;elsey(i)=x(i)A2-x(i)-1;endendyT2-if:s=input(pleaseenterthescore:);while(s100)disp(thescoreisnotreasonable);s=input(pleaseenterthescore:);endif(s=90&s=80&s=70&s=60&s70);disp(D);elsedisp(E);endT2-switch:s=input(pleaseenterthescore:);while(s100)disp(thesco

4、reisnotreasonable);s=input(pleaseenterthescore:);endswitchfix(s/10)case9,10disp(A);case8disp(B);case7disp(C);case6disp(D);case0,1,2,3,4,5disp(E);endT3:t=input(请输入工时time=);ift120w=120*84+(t-120)*1.15*84;elseift60w=t*84-700;elsew=84*t;enddisp(应发工资为：)；disp(w);T4:a=10+floor(rand(1)*89)b=10+floor(rand(1)

5、*89)s=input(请输入+或-或*或/,s);while(s=+&s=-&s=*&s=/)disp(输入的符号错误，请重新输入);s=input(请输入+或-或*或/,s);endswitchscase+c=a+b;case-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:n=input(请输入n:);m=0;fori=1:n;m=m+1/iA2;endpii=sqrt(6*m)T1-2:n=input(请输入n:);m=1./(1:n).

6、A2;s=sum(m);pii=sqrt(6*s)T2:n=1000;y=0;fori=1:n;y=y+1/(2*i-1);ify=3disp(最大n值为：);i-1disp(对应的y值为：);y-1/(2*i-1)break;endendT3:formatlonga=input(请输入a:);b=input(请输入b:);x(1)=a/(1.0+b);i=1;x(2)=a/(x(i)+b);while(abs(x(i+1)-x(i)0.00001&i500)i=i+1;x(i+1)=a/(x(i)+b);endr1=(-b+sqrt(bA2+4*a)/2;r2=(-b-sqrt(bA2+4*

7、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(r2,8);formatshortT4:f(1)=1;f(2)=0;f(3)=1;i=4;while(i0)fn=length(find(f0)f0=length(find(f=0)T5:j=0;s=0;fori=3:50m=i*(i-1)-1;ifisprime(m)=1ms=s+m;j=j+1;endenddisp(亲密数的对数：,num2str(j);disp

8、(所有亲密素数之和：,num2str(s);实验五：fn1:functionf=fn1(n)f=n+10*log(nJ2+5);fn2:functionf=fn2(n)f=0;fori=1:nf=f+i*(i+1);endfx:functionf=fx(x)f=1./(x-2)八2+0.1)+1./(x-3)八4+0.01);T1:x=input(请输入一个复数：);fe=exp(x)flg=log(x)fsin=sin(x)fcos=cos(x)T2:m1=input(请输入ml:);m2=input(请输入m2:);t=input(请输入theta:);A=m1*cos(t),-m1,-s

9、in(t),0;m1*sin(t),0,.cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1B=0,m1*9.8,0,m2*9.8disp(a1,a2,N1,N2=);disp(AB);T3:j=1;n=1;k=1;fo门=10:99%挑选出1099中所有的素数ifisprime(i)=1m(j)=i;j=j+1;endendfort=1:length(m);s(n)=10*rem(m(t),10)+fix(m(t)/10);%挑选出的素数进行位置交换ifisprime(s(n)=1%M挑选出交换绝对素数r(k)=m(t);k=k+1;endn=n+1;enddisp(

10、所有两位绝对素数：,num2str(r)T5:y1=fn1(40)/(fn1(30)+fn1(20)%函数调用y2=fn2(40)/(fn2(30)+fn2(20)实验六：T1:x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+xd2).*cos(x);plot(x,y);T2:x=linspace(0,2*pi,100);y1=x.A2;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)

11、;plot(x,y2,y-);subplot(2,2,4);plot(x,y3,b-.);T2-3:x=linspace(0,2*pi,10);y1=x.A2;y2=cos(2*x);y3=y1.*y2;subplot(3,4,1);bar(x,y1,r);axis(0,7,0,40);subplot(3,4,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

12、);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);subplot(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,

13、4,12);fill(x,y3,g);axis(0,7,-30,50);T3:x=linspace(-5,5,100);y=;forx0=x;ifx0=0;y=y,(x0+sqrt(pi)/exp(2);elseendendsubplot(1,2,1);plot(x,y);fori=1:length(x)ifx(i)0.5);b=length(k);disp(大于0.5的随机数个数占总数的百分比：);sprintf(%f%,b/300)T2:P=45+50*rand(100,5)disp(最高分及序号：)；Y,U=max(P,1)disp(最低分及序号：);Y1,U1=min(P,1)disp

14、(每门课平均分：);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,descend)%descend降序T3:h=6:2:18;t1=18.0,20.0,22.0,25.0,30.0,28.0,24.0;t2=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:

15、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,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)AA=polyvalm(P,A)实验九：T1:f

16、=inline(det(xxA2xA3;12*x3*xA2;026*x);i=1;x=1;whilex=3.01g(i)=f(x);i=i+1;x=x+0.01;enddx=diff(g)/0.01;dx(1)dx(101)dx(length(g)-1)T2:f1=inline(sqrt(cos(t.A2)+4*sin(2*t).A2+1)f2=inline(log(1+x)./(1+x.A2)I1=quad(f1,0,2*pi)I2=quad(f2,0,2*pi)T3: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=l

17、u(A);x2=U(LB)Q,R=qr(A);x3=R(QB)fun:functionF=fun(X);x=X(1);y=X(2);z=X(3);F(1)=sin(x)+yA2+log(z)-7;F(2)=3*x+2Ay-zA3+1;F(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)实验十：T1:x=sym(6);y=sym(5);z=(x+1)/(sqrt(3+x)-sqrt(y)T2:1symsxy;t=sym(5135);factor(xA4-yA4)fa

18、ctor(t)T3:symsbeta1beta2x;simplify(sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2)simplify(4*xA2+8*x+3)/(2*x+1)T4:symsabcdefghk;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)B1*B;tril(B)det(B)T5:symsxtayf1=(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=aAx,tA3;t*cos(x),log(x);diff(A,x,1)&(4)diff(A,t,2)diff(diff(A,x,1),t,1)f5=(xA2-2*x)*exp(-xA2