实验4(延伸) MATLAB与微分方程
由于 实验4 的导弹飞机追击问题需要用到,微分方程的概念,但MATLAB更偏向于离散型概念理解,直接应用并理解比较困难,所以这里单独将微分方程概念抽出来。本文将着重阐述微分方程的概念以及常见应用。关于解微分方程的 MATLAB 命令介绍在本文最下面。
一、微分方程的基本概念
% 初值条件:微分方程过点(1,2)
% 方程建立:任一点(x,y)处的切线的斜率为2x
clear;close all;clc;
syms f(x);
eqn=diff(f,x)==2*x;
cond=f(1)==2;
dsolve(eqn,cond);
disp(ans);% 输出 1 + x^2
% 2. 方程建立:加速度:(ds/dt)/dt=-0.4,本体是利用微分方程的特性求解
% 利用初值条件求s特解,能不能直接在dsolve中加上两个特解去解啊?能
clear;close all;clc;
syms s(t);
eqna=diff(s,t,2)==-0.4; % 基于要求的通解,先建立一个微分方程;但注意这是二阶,一般需要两个初值条件
v=diff(s,t); % 再基于要求的通解,建立一个相关函数,这里是v,目的是给第二个初值条件创建函数
cond=[s(0)==0,