MATLAB期末上机试题带答案.pdf

MATLAB 期末上机考试试题带答案版 姓名：学号：成绩： 1.请实现下图： 50100150200250 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel(x) ylabel(y) title(y=sin(x) 2.请实现下图： x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,k-) grid on xlabel(x) ylabel(y) title(sin(x) legend(y=sin(x) y2=cos(x); subplot(2,2,2) plot(x,y2,r-) grid on xlabel(x) ylabel(y) title(cos(x) legend(y=cos(x) y3=tan(x); subplot(2,2,3) plot(x,y3,k-) grid on xlabel(x) ylabel(y) title(tan(x) legend(y=tan(x) y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel(x) ylabel(y) title(cot(x) legend(y=cot(x) 3.解方程组： a=3 2 1;1 -1 3;2 4 -4； b=7;6;-2 ； x=ab 4.请实现下图： 02468101214 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 x y 自 己 的 姓 名 sin(x) sin(2*x) x=linspace(0,4*pi,1000); y1=sin(x); y2=sin(2*x); plot(x,y1,-,x,y2,b*) grid on xlabel(x);ylabel(y);title(耿蒙蒙耿蒙蒙) legend(sin(x),sin(2*x) 5.请在 x，y 在（-2，2）内的 z=xexp (-x2-y2) 绘制网格图 x,y=meshgrid(-2:0.1:2); z=x.*exp (-x.2-y.2); mesh(x,y,z) 6.请实现 peaks 函数： -3 -2 -1 0 1 2 3 -2 0 2 -5 0 5 x Peaks y x,y=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z) surf(x,y,z) shading flat axis(-3 3 -3 3 -8 8) xlabel(x);ylabel(y);title(Peaks) 7.请在 x=0,2,y=-0.5*pi,7.5*pi,绘制光栅的振幅为 0.4 的三维正弦光栅。 x=0:0.1:2;y=-0.5*pi:0.01*pi:7.5*pi; x,y=meshgrid(x,y) z=sin(y); mesh(x,y,z) surf(x,y,z) shading flat axis(0 2-0.5*pi 7.5*pi-6 6) 8.请用 ezplot 函数绘制 sin（x2） ，x 的区间为 0 到 8*pi。 ezplot(x,sin(x2),0,8*pi) 9.样本点; x=00.25*pi0.5*pi0.75*pipi1.25*pi1.5*pi1.75*pi2*pi; y=00.5*2.0.510.5*2.0.50-0.5*2.0.5-1-0.5*2.0.50; 对样本点进行 spline 插值。并将样本点和插值后的数据进行绘图。 x0=00.25*pi0.5*pi0.75*pipi1.25*pi1.5*pi1.75*pi2*pi; y0=00.5*2.0.510.5*2.0.50-0.5*2.0.5-1-0.5*2.0.50; plot(x0,y0,o) hold on x=0:0.01:2*pi; y=interp1(x0,y0,x,spline) plot(x,y) 10.请实现下图：利用多项式求根方法求解 x3-x2-3=0。 y=1 -1 0 -3; r=roots(y) 11.A= 2 i + 5j + 7 k B = 8 i + 4j + 6 k 求C=AB A=2 5 7; B=8 4 6; C=conv(A,B) 12.A= 2 i + 5j + 7 k B = 8 i + 4j + 6 k 求 C=AB A=2 5 7; B=8 4 6; C=A*B 13. 用不同标度在同一坐标内绘制曲线 y1=e-0.3xcos(2x)及曲线 y2=10e-1.5x。 x=0:pi/180:2*pi; y1=exp(-0.3*x).*cos(2*x);y2=10*exp(-1.5*x); plotyy(x,y1,x,y2) 14.请实现下图： 200400600800100012001400 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 x=linspace(0,8*pi,1400); y=sin(x); plot(x,y) area(y,0) 15. n=dblquad(exp(-(x.2)/2).*sin(x.2+y),-1,1,-2,2) 16.请实现:t=0:0.1:10*pi,x=tcos(t),y=tsin(t),z=t,三维曲线。 ezplot3(t.*cos(t),t.*sin(t),t,0,10*pi) 17. a=1 6 11 6; r=roots(a) poly(r) 18. syms x f=5*x3+6*x2+3*x+9; diff(f,x,1) 19. 已知样本点 x=-2.8 -1 0.2 2.1 5.2 6.8; y=3.1 4.6 2.3 1.2 2.3 -1.1;求其三次拟合，并绘出样 本点和拟合图像。 x=-2.8 -1 0.2 2.1 5.2 6.8; y=3.1 4.6 2.3 1.2 2.3 -1.1; plot(x,y,o) hold on p=polyfit(x,y,3) x0=-3:0.01:7; y0=polyval(p,x0) plot(x0,y0,r-) 20.构建内联函数 y=sin(x)exp(x2);并求出 x=1 4 2 5 8的 y 值。 y=inline(sin(x).*exp(x.2);xi= 1 4 25 8; yi=y(xi) 21.请实现从距离地面 20 米高处，以水平速度 5m/s 跳下的实际运动轨迹。 解： 2 1 2 hgt ， 0 xv t，得 2 2 0 2 g hx v ，即 2 0.2hx （0x10m） x=0:0.1:10;h=-0.2*x.2; plot(x,h,-) grid on xlabel(x/m);ylabel(h/m);title() 22.请绘出斜抛运动的实际轨迹。 初速度为 10m/s，与地面的夹角为 300。 解： 2 0 1 2 y hgtv t , 0x xv t,得 2 11 153 hxx x=0:0.1:13;h=-15(-1)*x.2-3(-1/2)*x; plot(x,h) grid on xlabel(x/m);ylabel(h/m);title() 23.请求出df(x)/dx=ax3+x2-bx-c syms x a b c f=a*x3+x2-b*x-c; diff(f,x,1) 24. x,y=meshgrid(-3:0.1:3); z=1./(x+1).2+(y+1).2+1)-1./(x-1).2+(y-1).2+1); mesh(x,y,z) 25. x=-10:0.01:10 subplot(1,2,1) plot(x,sin(2*x).*cos(3*x) xlabel(x);ylabel(y);title(sin(2x)cos(3x) subplot(1,2,2) plot(x,0.4*x) xlabel(x);ylabel(y);title(0.4x) 26. x=0:0.01:25; y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.8; y2=1.6*cos(3*x)+sin(x); plot(x,y1,b-,x,y2,r-) legend(y1=2.6exp(-0.5x)cos(0.6x)+0.8,y2=1.6cos(3x)+sin(x) grid on 27. y=int(sin(x)+2,x,0,pi/6) 28. solve(sin(x)+tan(x)+1=0,x) 29. syms x y=dsolve(Dy=(x+y)*(x-y),x) 30. 解：插值法 拟合法 31. 请用三种方法求解 sin(x)在 0pi 之间的积分。 1：a=quad(sin(x),0,pi) 2：x=linspace(0,pi,1000);y=sin(x); a=trapz(x,y) 3：y=int(sin(x),x,0,pi) 32. x,y=meshgrid(-2:0.1:2);z=x.2.*exp(-x.2-y.2); mesh(x,y,z) surf(x,y,z) 33. 22 yx xez ，当 x 和 y 的取值范围均为-2 到 2 时，用建立子窗口的方法在同一个图形窗 口中绘制出网线图、表面图和去网格效果的表面图。 x,y=meshgrid(-2:0.1:2); z=x.*exp(-x.2-y.2); subplot(1,3,1) mesh(x,y,z) subplot(1,3,2) surf(x,y,z) subplot(1,3,3) surf(x,y,z) shading flat 34. 有一组测量数据满足 -at ey，t 的变化范围为 010，用不同的线型和标记点画出 a=0.1、 a=0.2 和 a=0.5 三种情况下的曲线。 t=0:1:10;y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t); plot(t,y1,b-o,t,y2,r-*,t,y3,g-) 35. 有一正弦衰减数据 y=sin(x).*exp(-x/10)，其中 x=0:pi/5:4*pi，用三次样条法进行插值。 x=0:pi/180:4*pi;y=sin(x).*exp(-x/10);xi=0:pi/5:4*pi; yi=interp1(x,y,xi,cubic) plot(x,y,b-,xi,yi,ro) 36. 求解多项式x3-7x2+2x+40的根。 a=1 -7 2 40; r=roots(a) 37. 对于BAX ，如果 753 467 294 A， 28 26 37 B，求解 X。 a=4 9 2;7 6 4;3 5 7;b=37;26;28; x=ab 38. 请建立隐函数，y=x2exp(x),并求出 x=7 8 9时的函数值。并将函数绘图 y=inline(x.2.*exp(x); xi=7 8 9; yi=y(xi) 39. 在-10,10,-10,10的范围内会三维图 22 22 s