资源描述：
《matlab与工程计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《MATLAB与工程计算》实验题目1.在一个图形窗口中分别画两幅横向排列的子图：和。(基本要求：用不同颜色绘制两条曲线，并使用text等命令实现对图形的文字说明)答案：x=linspace(-6*pi,6*pi,100);y=x.^2.*cos(x);z=sin(x)./x;subplot(1,2,1);plot(x,y,'k-');title('x.^2.*cos(x)');text(0,0,'y=x.^2.*cos(x)');subplot(1,2,2);plot(x,z,'r-');title('sin(x)./x');
2、text(pi/2,2/pi,'sin(x)./x')2.分别用数值积分和符号积分命令求。答案：（1）数值积分结果：>>g=inline('log(x+(x.^2+1).^(1/2))');i=quadl(g,0,1)i=0.4672（2）符号积分结果：>>x=sym('x');f=log(x+(x.^2+1).^(1/2));int(f,0,1)ans=log(2^(1/2)+1)-2^(1/2)+13.分别通过数值计算和符号计算的方法求解非线性方程的根：答案：（1）数值计算functionf=f(x);f=x.^2+10*
3、(x-1).*x.^(1/2)+14*x+1;（先建立函数文件）>>fzero('f',1)ans=0.1716（2）符号计算>>f=sym('x^2+10*(x-1)*x^(1/2)+14*x+1');x=solve(f)x=3-2*2^(1/2)33-8*17^(1/2)4.设A(0,0)为一导弹发射点，发现位于B(0,100)处一架敌机沿水平方向逃离（如图），随即发射一枚导弹予以打击，现已知导弹时刻对准敌机，且速率为飞机速率的两倍（设飞机速度为1）。试编程模拟导弹打击敌机的动态过程。如果敌机飞行60单位距离之外即逃出我方
4、空域，那么，要想在我方空域内击落敌机，则导弹的速度至少应提高到敌机速度的多少倍？(要求：给出模型(或主要公式)及matlab程序)答案：(1),holdonaxis([03000120])grid;A=[0,0];B=[0,100];d=norm(A-B);k=0;v=1;dt=1;K=input('K=');whilek<1000k=k+1;B=B+[v*dt,0];plot(A(1),A(2),'r.','markersize',15);plot(B(1),B(2),'b.','markersize',15);e=B-A;
5、d=norm(e);fprintf('k=%.0fB(%.0f,100)A(%.2f,%.2f)d=%.2f',k,B(1),A(1),A(2),d)ifd<=0.5breakende=e/d;A=A+K*v*dt*e;pause(0.01)end(2),利用二分法对K取不同的值，经过多次试验发现K=2.17时，k=60B(60,100)A(59.91,100.02)d=0.09;K=2.175时，k=60B(60,100)A(60.16,100.00)d=0.16，所以在误差允许的范围内，K>2.17即可在60个单位距离
6、内实施打击。附加题：计算椭圆的周长，使计算结果至少具有五位有效数字。答案：>>formatshorteh=1.e-6;x=0:h:2;y=sqrt(1-0.25*x.*x);dx=h;dy=diff(y);dl=sqrt(dx.^2+dy.^2);l=4*sum(dl)l=9.6884e+000
