北科大MATLAB考试用参考代码.doc

北科大MATLAB考试都是书后题稍微改动，本人再次将所有老师留过的作业里的代码整理在这，希望给后来者一点帮助哟P27.1用不同数据格式显示自然底数e的值 （见书p6）formatshort,exp(1)ans=2.7183formatlong,exp(1)ans=2.718281828459046formatshorte,exp(1)ans=2.7183e+000formatlonge,exp(1)ans=2.718281828459046e+000formatshortg,exp(1)ans=2.7183formatlongg,exp(1)ans=2.71828182845905formatrat,exp(1)ans=1457/536formathex,exp(1)ans=4005bf0a8b14576aformat+,exp(1)ans=+formatbank,exp(1)ans=2.72formatcompact,exp(1)ans=2.72formatloose,exp(1)ans=2.72P27.2矩阵A=1 2 3；4 5 6；7 8 9 ，B=4 6 8;5 5 6;3 2 2A=1,2,3;4,5,6;7,8,9A=123456789B=4,6,8;5,5,6;3,2,2B=468556322A*Bans=23222659617495100122A.*Bans=41224202536211618P27.3已知矩阵A=5 2;9 1B=1 2;9 2A=5,2;9,1A=5291B=1,2;9,2B=1 292ABans=1000A=Bans=0110A(A=B)&(A(A=B)&(AB)ans=0000P79.1绘制的图像，要求用蓝色的星号画图clf;x=0:pi/50:4*pi;y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3);plot(x,y1,b*,x,y2,r-.,x,y3,r-.),gridonlegend(y1=exp(x/3).*sin(3*x),y2=+-exp(x/3)p79.3在同一图形窗口画图要求使用指令 gtext，axis，legend，title，xlabel，ylabelclf; x1=-pi:pi/50:pi; x2=pi:pi/50:4*pi; x3=1:0.1:8; y1=x1.*cos(x1); y2=x2.*tan(x2.(-1).*sin(x2.3); y3=exp(x3.(-1).*sin(x3); subplot(2,2,1);plot(x1,y1,m.);grid on;title(y=x*cosx);xlabel(x00);ylabel(y00);gtext(y=x*cosx);legend(y=x*cosx);subplot(2,2,2);plot(x2,y2,r*);grid on;title(y=x*tan(1/x)*sin(x3)gtext(y=x*tan(1/x)*sin(x3);legend(y=x*tan(1/x)*sin(x3);subplot(2,2,3);plot(x3,y3,bp);grid on;title(y=e(1/x3)*sinx);xlabel(x00);ylabel(y00);gtext(y=e(1/x3)*sinx);legend(y=e(1/x3)*sinx)p79.5绘制圆锥曲线的图像并加各种标注t=0:pi/50:20*pi;x=t.*cos(t*pi/6);y=t.*sin(t*pi/6);z=2*t;plot3(x,y,z)（标注见上一题）P79.7用mesh和surf命令 绘制三维曲面的图像t=-1:0.1:1;x,y=meshgrid(t);z=x2+3*y2;subplot(1,2,1),mesh(x,y,z),colormap(bone),light(position,20,20,5);subplot(1,2,2),surf(x,y,z),colormap(cool)p79.8绘制由函数=1形成的立体图xx,yy,zz=sphere(40);x=xx*2;y=yy*3;z=zz*4;subplot(2,2,1),surf(x,y,z);subplot(2,2,2),surf(x,y,z);view(0,90)subplot(2,2,3),surf(x,y,z);view(90,0)subplot(2,2,4),surf(x,y,z);view(0,0)p79.9画三维曲面z=与平面x=3的交线t=-2:0.1:2;x,y=meshgrid(t);z1=5-x.2-y.2;subplot(1,3,1),mesh(x,y,z1),title(z1=5-x.2-y.2);z2=3*ones(size(x);subplot(1,3,2),mesh(x,y,z2),title(z=3);r0=abs(z1-z2)=1;zz=r0.*z2;yy=r0.*y;xx=r0.*x;subplot(1,3,3);subplot(1,3,3),plot3(xx(r0=0),yy(r0=0),zz(r0=0),.),title(交线)p114.12随机输入一个6阶方阵a=rand(6);adet(a)rank(a)rref(a)p114.14求矩阵A=2 1 1;1 2 1;1 1 2的特征多项式、特征值、特征向量a=2 1 1;1 2 1;1 1 2;poly(a)d=eig(a)V,D=eig(a)P115.21.(1)求线性方程组通解A=1 1 2 -1;-1 1 3 0;2 -3 4 -1;rref(A)分析过程为原方程组等价于x1=0.56x4x2=0.2x4x3=0.12x4取x4=1则x1=0.56,x2=0.2,x3=0.12基础解系为=0.56 0.2 0.12 1所以通解为x1 x2 x3 x4=0.56 0.2 0.12 1P115.21.(2)B=1-1-110;1-11-31;1-1-23-1/2;rref(B)分析过程为原方程组等价于X1=x2+x4+0.5X3=2x4+0.5其中一个特解为k=0.5 0 0.5 0对应的其次线性方程组为X1=x2+x4X3=2x4基础解系为1=1 1 0 0 2=1 0 2 1所以通解为x=c1*1+c2*2+kP130.8在钢线含碳量对于电阻的效应的研究中，得到以下数据x=0.10 0.30 0.4 0.55 0.70 0.80 0.95;y=15 18 19 21 22.6 23.8 26;p1=polyfit(x,y,1);p3=polyfit(x,y,3);p5=polyfit(x,y,5);disp(一次拟合函数),f1=poly2str(p1,x)disp(三次拟合函数),f1=poly2str(p3,x)disp(五次拟合函数),f1=poly2str(p5,x)x1=0.10:0.01:0.95;y1=polyval(p1,x1);y3=polyval(p1,x1);y5=polyval(p1,x1);plot(x1,y1,rp,x1,y3,-,x1,y5);legend(一次拟合,三次拟合,五次拟合)p130.10在某种添加剂的不同浓度下对铝合金进行抗拉强度试验x=10:5:30;y=25.2 29.8 31.2 31.7 29.4;x1=10:1:30;y1=interp1(x,y,x1,nearest);y2=interp1(x,y,x1,linear);y3=interp1(x,y,x1,spline);y4=interp1(x,y,x1,cubic);subplot(2,2,1),plot(x,y,rp,x1,y1),title(nearest);subplot(2,2,2),plot(x,y,rp,x1,y2),title(linear);subplot(2,2,3),plot(x,y,rp,x1,y3),title(spline);subplot(2,2,4),plot(x,y,rp,x1,y4),title(cubic);a=interp1(x,y,18,spline)b=interp1(x,y,26,cubic)p130.12用不同方法对z=x2/16-y2/9在（-3,3）上的二维插值效果进行比较x,y=meshgrid(-3:.5:3);z=x.2/16-y.2/9;x1,y1=meshgrid(-3:0.1:3);z1=x1.2/16-y1.2/9;figure(1)subplot(1,2,1),mesh(x,y,z),title(数据点)subplot(1,2,2),mesh(x1,y1,z1),title(原函数图象)xi,yi=meshgrid(-3:0.125:3);zi1=interp2(x,y,z,xi,yi,*nexarest);zi2=interp2(x,y,z,xi,yi,*linear);zi3=interp2(x,y,z,xi,yi,*spline);zi4=interp2(x,y,z,xi,yi,*cubic);figure(2)subplot(2,2,1),mesh(xi,yi,zi1),title(最近点插值)subplot(2,2,2),mesh(xi,yi,zi2),title(线性插值)subplot(2,2,3),mesh(xi,yi,zi3),title(样条插值)subplot(2,2,4),mesh(xi,yi,zi4),title(立方插值)p168.24求解微分方程y=xsinx/cosydsolve(Dy=x*sin(x)/cos(y),x)p168.27用数值方法其求下列微分方程，用不同颜色和线形将y和y画在同一个图形窗口里原微分方程为即写成矩阵形式为令u=3cost令xdot= 红字是分析过程，需手写录入第一步先编第一个m文件function xdot=exf(t,x)u=1-2*t;xdot=0,1;1,-t*x+0 1*u;第二步编写第二个m文件clf;t0=0;tf=pi;x0t=0.1;0.2;t,x=ode23(exf,t0,tf,x0t);y=x(:,1);Dy=x(:,2);plot(t,y,r-,t,Dy,b*);legend(y,Dy);xlabel(t轴)解析解为dsolve(D2y+y=1-t2 /(2*pi),y(0)=0,Dy(0)=0,t) P190.15描绘一下 数组的频数直方图load A.txthist(A,10)（题目所用数据）6.829.633.635.736.945.254.865.843.453.863.769.970.779.597.9139.4157.0p190.16求样本均值、标准差、中位数、极差和方差load B.txtM=mean(B),std(B),median(B),range(B),var(B)P167.17(2)h=pi/10;x=0:h:1;y=x.*sin(x)./(1+(cos(x).2);formatlongt=length(x);s1=sum(y(1:(t-1)*hs2=sum(y(2:t)*hs3=tra