数学建模微分方程导弹问题matlab求解,Matlab建模—导弹追踪问题
数学建模期末复习,撰写博客做总结之用,主要侧重于算例的模型建立与部分代码的实现,其中不足之处望读者多多指正。
Matlab微分求解工具箱使用
求数值解
dsolve(‘方程1’,‘方程2’,…,‘方程n’,‘初始条件’,‘自变量’)
求解析解
[t,x]=solver(’f’,ts,x0,options)
详细的使用说明可参考博文Matlab微分方程求解
导弹问题
问题介绍
设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰发射导弹,导弹头始终对准乙舰.如果乙舰以最大的速度v0v_0v0(常数)沿平行于y轴的直线行驶,导弹的速度是5v0v_0v0,求导弹运行的曲线方程.乙舰行驶多远时,导弹将它击中?
模型建立与求解
1. 解析法
模型建立:设t时刻导弹的位置为P(x(t),y(t)),乙舰位于Q(1,V0tV_0tV0t),由题(导弹头始终对准乙舰)可得,导弹在运动弧线p点处的切线y′=v0t−y1−xy^{\prime}=\frac{v_{0} t-y}{1-x}y′=1−xv0t−y即 v0t=(1−x)y′+y(1)v_{0} t=(1-x) y^{\prime}+y(1)v0t=(1−x)y′+y(1)
又由题意(导弹速度是舰的5倍)故有:∫0x1+y′2dx=5v0t(2)\int_0^x {\sqrt {1 + y{'^2}} } {\rm{d}}x = 5{v_0}t (2)∫0x1+y′2dx=5v0t(2)
有(1)(2)式可得微分微分方程为:(1−x)y′′=151+y′2(3)(1-x) y^{\prime \prime}=\frac{1}{5} \sqrt{1+y^{\prime 2}}(3)(1−x)y′′=511+y′2(3)
初值方程为:y(0)=0y(0) = 0y(0)=0,y′(0)=0y'(0) = 0y′(0)=0,利用Matlab求解可得:y=−58(1−x)45+512(1−x)65+524y = - \frac{5}{8}{(1 - x)^{\frac{4}{5}}} + \frac{5}{{12}}{(1 - x)^{\frac{6}{5}}} + \frac{5}{{24}}y=−85(1−x)54+125(1−x)56+245
2. 数值解法
将上述方程化为一阶方程,不妨设y1=y,y2=y′y_1=y,y_2=y'y1=y,y2=y′,代入(3)中得到{y1′=y2y2′=151+y12/(1−x)\left\{\begin{array}{l}y_{1}^{\prime}=y_{2} \\ y_{2}^{\prime}=\frac{1}{5} \sqrt{1+y_{1}^{2}} /(1-x)\end{array}\right.{y1′=y2y2′=511+y12/(1−x)
利用Matlab求解:
%定义函数
function dy=eq1(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);
end
调用
x0=0;
xf=0.9999;
[x,y]=ode15s('eq1',[x0 xf],[0 0]);
plot(x,y(:,1),'-')
hold on
y=0:0.01:2;
plot(1,y,'*')
复制代码
由图可得,导弹大致在点(1,0.2)处击中乙舰。
3. 参数方程
上述结题假设不变,在此基础再将导弹的速度设为www,则在P点处有:(dxdt)2+(dydt)2=w2(1){({{{\rm{d}}x} \over {{\rm{d}}t}})^2} + {({{{\rm{d}}y} \over {{\rm{d}}t}})^2} = {w^2}(1)(dtdx)2+(dtdy)2=w2(1)
同时,导弹始终对准乙舰,即导弹速度向量与乙舰的位置向量平行,故有:
\frac{\mathrm{d} x}{\mathrm{d} t} \\ \frac{\mathrm{d} y}{\mathrm{d} t} \end{array}\right)=\lambda\left(\begin{array}{l} X-x \\ Y-y \end{array}\right), \quad \lambda>0(2)$$ 将 $$\left\{\begin{array}{l} \frac{\mathrm{d} x}{\mathrm{d} t}=\frac{w}{\sqrt{(X-x)^{2}+(Y-y)^{2}}}(X-x) \\ \frac{\mathrm{d} y}{\mathrm{d} t}=\frac{w}{\sqrt{(X-x)^{2}+(Y-y)^{2}}}(Y-y) \end{array}\right.(3)$$ 不妨假设乙舰速度为1,则w=5,X=1,Y=t,此时导弹的参数轨迹为: $$\left\{\begin{array}{l} \frac{\mathrm{d} x}{\mathrm{d} t}=\frac{5}{\sqrt{(1-x)^{2}+(t-y)^{2}}}(1-x) \\ \frac{\mathrm{d} y}{\mathrm{d} t}=\frac{5}{\sqrt{(1-x)^{2}+(t-y)^{2}}}(t-y) \\ x(0)=0, y(0)=0 \end{array}\right.$$ matlab求解: ```php function dy=eq2(t,y) dy=zeros(2,1); dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2); dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2); end %调用 [t,y]=ode45('eq2',[0 2],[0 0]); Y=0:0.01:2; plot(1,Y,'-') hold on plot(y(:,1),y(:,2),'*') ``` ![在这里插入图片描述](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/c12063ad33dd4333b64448e9ae4d3075~tplv-k3u1fbpfcp-zoom-1.image) 由图可得,导弹最终大致在(1,0.2)处命中目标。这里还以尝试用二分法的思想进一步得到更精确的答案,……tf=0.21是可得求解图像为: ![在这里插入图片描述](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/b4f9a6123c24433383fcd64a18b1adfa~tplv-k3u1fbpfcp-zoom-1.image) 据此可得更为精确的答案。
数学建模微分方程导弹问题matlab求解,Matlab建模—导弹追踪问题相关推荐
- 偏微分方程matlab求解,Matlab求解微分方程(2)——偏微分方程的求解
从写完上一篇常微分方程的求解到现在已经很长时间了,这周也一直忙于报到的各种事宜,无暇坐下来写些东西,趁着这个周末,终于完成了这个姊妹篇. 对于偏微分方程的求解,Matlab提供了两种工具.第一种是pd ...
- evalin matlab,求解MATLAB问题
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 这个是代码: syms x y z [deltafd]=xlsread('附件2.xls','FDOA','B2:B290'); [x1]=xlsread ...
- 线性规划 - 用单纯形法解决整数规划问题 - (Matlab、Lingo建模)
现实生活中,比如机器的台数,参与工作的人数,可调动的车辆数,这些数据都是整数.因此对于变量中包含整数.或者完全是整数的规划问题,我们称之为整数规划.在解决整数规划常用的算法便是单纯形法. 课题名称:任 ...
- 数学建模微分方程导弹问题matlab求解,数学建模之微分方程(符实现例题和MATLAB源码)...
微分方程的基本概念 微分方程:一般的,凡表示未知函数.未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程. 微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的 ...
- 数学建模专栏 | 第七篇:MATLAB连续模型求解方法
连续模型是指模型是连续函数的一类模型总称,具体建模方法主要是微分方程建模.微分方程建模是数学建模的重要方法,因为许多实际问题的数学描述将导致求解微分方程的定解问题.把形形色色的实际问题化成微分方程的定 ...
- 数学建模之微分方程(符实现例题和MATLAB源码)
微分方程的基本概念 微分方程:一般的,凡表示未知函数.未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程. 微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的 ...
- MATLAB数学建模-规划模型总结| MATLAB求解
目录 1 线性规划问题(LP) 风格1 风格2 2 非线性规划 3 动态规划 A星算法 基于dijkstra的概率路线图 4 多目标规划 帕累托最优 支配(Dominace) 不可支配解集 帕累托最优 ...
- matlab求解微分方程6,牛津大学出版社数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解:经典爱情语录大全...
漳州理工职业学院-酒会礼仪 注:上机作业文件夹以自己的班级姓名学号命名,文件夹包括如下上机报告和Matlab程序. 上机报告模板如下: 佛山科学技术学院 上 机 报 告 课程名称 数学应用软件 上机项 ...
- matlab求解最优化问题(数学建模)
matlab求解最优化问题(数学建模) 1.线性规划 matlab中线性规划优化计算方法和实例 在matlab中用于线性规划优化计算的是linprog()函数. 公式:[x,fval,exitflag ...
最新文章
- Python使用PIL将png图片转化为jpg图片
- 面试干货——年底干货大放送,你准备好了吗?
- select * 映射错误_高性能IO模型分析-浅析Select、Poll、Epoll机制(三)
- 抖音回应李小璐PGone视频曝光:草稿视频不会上传到后台
- windows如何安装MySql(包含一些安装时问题的解决)
- [51单片机] 串口通讯 简单通信
- modelsim安装_Modelsim10.7安装教程
- jmeter 安装与配置
- 将PDF论文的公式截图后转成Word可编辑公式
- IANA已注册的TCP/UDP/SCTP/DCCP传输协议端口及服务名称
- php支付宝接口开发流程,对接支付宝支付接口开发详细步骤(证书签名方式)
- 弗洛伊德:本我、自我和超我
- 大学生学C语言用什么笔记本电脑,有哪些适合大学生用的笔记本电脑
- 楼氏电子推出Raspberry Pi 开发工具包,为新物联网应用和行业语音整合提供支持
- 如花搞笑图片集锦(转贴)
- JVM中的monitorenter和monitorexit
- 再生核希尔伯特空间(RKHS)
- 高通狂吹新GPU:赶超桌面显卡
- 【厚积薄发系列】Python项目总结1—后端常驻程序的基本要求
- HTML期末学生作业 HTML+CSS+JavaScript仿猫眼电影在线网站 Hbuilder网页制作
热门文章
- 武汉云控系统,马云,王健林都惊艳了
- 为什么Lisp如此先进,却永远成为不了编程主流语言?深度解析Lisp的优势与劣势
- 为什么用vue,它解决了什么问题,如何使用它?
- 混淆矩阵--在图像精度的应用(mIou指标)
- php程序员的出路,php程序员有前途吗
- Leetcode 14.最长公共前缀(Longest Common Prefix)
- php调用天气预报接口,PHP调用百度天气接口API实现查询实时天气
- 计算机网络定义记不住,记不住是怎么回事?——记忆是个神奇的小东西
- 逻辑回归(Logistic Regression)
- 三、Bugku----手机热点------流量分析题------obex