数学建模微分方程导弹问题matlab求解,数学建模之微分方程(符实现例题和MATLAB源码)...

本文介绍了微分方程的基本概念,并通过MATLAB求解了导弹追踪问题的解析解和数值解。文章提供了解析解的MATLAB代码,展示了解析解的过程,同时对比了数值解的方法,包括欧拉法和龙格-库塔法。最终,通过两种数值解法验证了导弹击中目标的坐标,结果接近解析解。
摘要由CSDN通过智能技术生成

微分方程的基本概念

微分方程:一般的,凡表示未知函数、未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程。

微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的阶。

微分方程的解:使得微分方程成立的函数称为微分方程的解。不含任意常数的解称为微分方程的特解。若微分方程的解中所含的相互独立的任意常数的个数与微分方程的阶数相等,称这个解为微分方程的通解。

微分方程的通解:如果微分方程的解中含有任意常数,且任意常数的个数与微分方程的阶数相同,这样的解叫做微分方程的通解。

微分方程的解析解

一般只有比较简单的微分方程才能求出解析解,后面做建模一般都是求数值解。

这里我们讲的基本都是直接用MATLAB程序算的

MATLAB求解命令为

dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’)

记号:在表达微分方程时,用字母D表示求微分,D2、D3等表示求高阶微分.任何D后所跟的字母为因变量,自变量可以指定或由系统规则选定为确省.

NzA36j.png

例一:求 d u d t = 1 + u 2 \frac{du}{dt}=1+u^2dtdu​=1+u2通解

dsolve('Du=1+u^2','t')

tan(C1 + t)%通解

1i%虚数+i、-i数值解

-1i

例二:

{ d 2 y d x 2 + 4 d y d x + 29 y = 0 y ( 0 ) = 0 y ′ ( 0 ) = 15 \begin{cases} \frac{d^2y}{dx^2}+4\frac{dy}{dx}+29y=0 \\ y(0)=0\\ y'(0)=15 \end{cases}⎩⎪⎨⎪⎧​dx2d2y​+4dxdy​+29y=0y(0)=0y′(0)=15​

y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')

y =

3*sin(5*x)*exp(-2*x)

BfUrqu.png

[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z', 't');

x=simplify(x) % 将x化简

y=simplify(y)

z=simplify(z)

x =

C1*exp(-t) + C3*exp(2*t)

y =

exp(-2*t)*(C2 + C1*

  • 2
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值