MATLAB语言与控制系统仿真参考的答案第3章.docVIP

3.5 MATLAB绘图实训 3.5.1 实训目的 1.学会MATLAB绘图的基本知识； 2.掌握MATLAB子图绘制、图形注释、图形编辑等基本方法； 3.学会通过MATLAB绘图解决一些实际问题； 4.练习二维、三维绘图的多种绘图方式，了解图形的修饰方法； 5.学会制作简单的MATLAB动画。 图3-46 炮弹发射示意图 3.5.2 实训内容 炮弹发射问题 （1）炮弹发射的基础知识 炮弹以角度射出的行程是时间的函数，可以分解为水平距离和垂直距离。 %水平方向的行程； %垂直方向的行程； 其中，是初速度； 是重力加速度，为9.8； 是时间。 （2）炮弹发射程序举例： 分析以下程序以及图3-47各个图形的实际意义。 a=pi/4; v0=300; g=9.8; t=0:0.01:50; x=t*v0*cos(a); y=t*v0*sin(a)-0.5*g*t.^2; subplot(221);plot(t,x);grid;title(‘时间-水平位移曲线); subplot(222);plot(t,y);grid;title(‘时间-垂直位移曲线); subplot(223);plot(x,y);grid;title(‘水平位移-垂直位移曲线); subplot(224);plot(y,x);grid;title(‘垂直位移-水平位移曲线); 图3-47 角发射曲线 （3）编程解决炮弹发射问题 ①假设在水平地面上以垂直于水平面的角度向上发射炮弹，即发射角，假设初速度分别为[310，290，270]m/s,试绘制时间-垂直位移曲线，编程求取最高射程； 绘图要求： 标题设为“炮弹垂直发射问题”； 在图上通过添加文本的方式表明初速度； 在x轴标注“时间”； 在y轴上标注“垂直距离”； 添加网格线； 将310m/s的曲线改为线粗为2的红色实线； 将290m/s的曲线改为线粗为3的绿色点划线； 将270m/s的曲线改为线粗为2的蓝色长点划线； a=pi/2; v1=310; g=9.8; t=0:0.01:50; x1=t*v1*cos(a); y1=t*v1*sin(a)-0.5*g*t.^2; plot(t,y1);grid; title(炮弹垂直发射问题); xlabel(时间); ylabel(垂直距离); hold on; v2=290; x2=t*v2*cos(a); y2=t*v2*sin(a)-0.5*g*t.^2; plot(t,y2); v3=270; x3=t*v3*cos(a); y3=t*v3*sin(a)-0.5*g*t.^2; plot(t,y3); zgsc=[max(y1); max(y2); max(y3)] %三次发射的最高射程 运行结果如下： zgsc = 1.0e+003 * 4.9031 4.2908 3.7194 最高射程分别为:4903.1米,4290.8米,3719.4米。 ②假设距地面3750米的高空（火炮正上方）有一个半径为2000米的圆形漂浮物，火炮初始速度为280m/s, 若发射角度为90度，问能击中吗？若发射角度为80度，能击中吗？ 绘图要求： 在同一幅图上绘制和两条射程曲线； 把曲线设为红色粗度为2的点划线形式； 通过gtext()命令在图上添加文本“高空漂浮物体击打”； 通过insert‘legend’添加图例，分别注明“90度角发射垂直位移”、“80度角发射垂直位移”；“80度角水平位移”； 在的曲线上通过添加箭头方式画一箭头，通过添加文本方式注明最高射程数据。 a1=pi/2; a2=80*pi/180; v1=280; v2=280; g=9.8; t=0:0.01:50; x1=t*v1*cos(a1); x2=t*v2*cos(a2); y1=t*v1*sin(a1)-0.5*g*t.^2; y2=t*v2*sin(a2)-0.5*g*t.^2; plot(t,y1); hold on; plot(t,y2); plot(t,x2) grid; gtext(高空漂浮物体击打); zgsc=[max(y1); max(y2)] zgsc = 1.0e+003 * 4.0000 3.8794 通过图中可以看出： 90度角发射时最高射程为4000米，因为4000＞3750，因而能够击中高空漂浮物。 80度角发射时，炮弹至最高点时，水平位移接近1000米时，垂直位移尚未达到3750米，因此80度角发射不能击中漂浮物。 ③以初速度2