v˙l=alv˙f=−0.1+5∗vf+0.25∗vf2mD˙=vl−vf\begin{aligned} & \dot{v}_l = a_l \\ & \dot{v}_f = -\frac{0.1 + 5*v_f + 0.25 * v_f^2}{m} \\ & \dot{D} = v_l - v_f \\ \end{aligned}​v˙l​=al​v˙f​=−m0.1+5∗vf​+0.25∗vf2​​D˙=vl​−vf​​

%% Solve Differential Equations
% Author: Zhao-Jichao
% Date: 2021-07-18
clear
clc%% Initial Parameters
% Variable initial values
vl0 = 1;
vf0 = 2;
D0  = 3;%% Solve ODE
% Time parameters
tBegin = 0;
tFinal = 10;
in = [vl0 vf0 D0]';
[t, out] = ode45(@odeFun, [tBegin, tFinal], in);%% Draw Results
subplot(2,2,1)
plot(t,out(:,1), t,out(:,2), t,out(:,3), 'linewidth',1.5);
legend('v_l', 'v_f', 'D');
grid on;subplot(2,2,2)
plot(t,out(:,1), 'linewidth',1.5); hold on ;
legend('v_l');
grid on;subplot(2,2,3)
plot(t,out(:,2), 'linewidth',1.5); hold on ;
legend('v_f');
grid on;subplot(2,2,4)
plot(t,out(:,3), 'linewidth',1.5); hold on ;
legend('D');
grid on;%% Ode SubFunction
function out = odeFun(~, in)vl = in(1,:);vf = in(2,:);D  = in(3,:);al = 1;m  = 1;% Set up differential relationsdvl = al;dvf = -(0.1 + 5*vf + 0.25*vf^2) / (m);dD  = vl - vf;out = [dvl; dvf; dD];
end

of

【Matlab】解个微分方程相关推荐

  1. matlab坡度计算公式,matlab解二阶微分方程怎么用matlab来解呢?x的定义域是(0,l),i是道路横坡坡度,w是降雨强度,k是路面横向...

    共回答了16个问题采纳率:81.3% 你的方程即为: y'*y' + y''*y' +w/k =0 s = dsolve('Dy*D2y + Dy^2 + w/k ', 'y(l) = a', 'Dy ...

  2. matlab ode45 二阶微分,matlab关于ode45解二阶微分方程的困惑

    matlab关于ode45解二阶微分方程的困惑 matlab关于ode45解二阶微分方程的困惑 一个二阶微分方程: y''+y'+y=sin(t) 初始条件为y(0)=5,y'(0)=6. 过程: 先 ...

  3. matlab 向前欧拉公式,向前欧拉公式在Matlab解微分方程初值解的问题

    向前欧拉公式在Matlab解微分方程初值解的问题0 fuqilin1202013.07.04浏览527次分享举报 用向前欧拉公式(10.8)求解初值问题,dy/dx=-3x+8x-7,y(0)=1,分 ...

  4. matlab求微分方程的初值问题,Matlab 解常微分方程的初值问题.doc

    Matlab 解常微分方程的初值问题 题目:Matlab 解常微分方程的初值问题 设计目的: 1.熟练掌握Matlab的基本编程方法,及其编程风格. 2.熟练掌握Matlab常用函数的使用. 3.与本 ...

  5. matlab解无解析解微分方程组,数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解...

    <数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解>由会员分享,可在线阅读,更多相关<数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解(12页 ...

  6. matlab求偏微分方程程序,用MATLAB解偏微分方程.pdf

    用MATLAB解偏微分方程.pdf 年 月 阴 山 学 刊 第 卷 第 期 丫叫 加 用 解偏微分方程 田 兵 包头师范学院 学报编辑部 , 内蒙古 包头 摘 要 讨论 了以 中偏徽分方程工具箱的用法 ...

  7. Matlab符号数学(Symbolic Math with MATLAB)MATLAB解方程

    创建符号变量 用符号而不是数值表示数字. 存储数字符号 % 你可以把一个数存储为数值格式,就像我们经常做的那样 x = 3; % 你也可以把这个数存储为符号格式,使用sym()函数 xsym = sy ...

  8. 如何用matlab解异或方程,Matlab-6:解非线性方程组newton迭代法

    函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\ ...

  9. 用MATLAB求解非线性微分方程

    总结一下MATLAB中求解微分方程的思路和步骤.固然,网上很多关于此类的技术型文章,但往往一看下来发现,文章中的友情链接比文章字数还多,要了解这一篇文章,你要先了解那个:要了解那个,你又要了解那个那个 ...

  10. matlab 解函数方程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值

    MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1  线性方程组求解 7.2  非线性方程数值求解 7.3  常微分方程初值问题的数值解法 7. ...

最新文章

  1. HTML5column属性布局页脚,利用column多列属性调整页面文字列布局
  2. Jupyter Notebook 使用小记
  3. JDBC批量操作批量增加批量修改
  4. [swift]-类的构造函数
  5. char*赋值给std::string是深拷贝
  6. matplotlib——直方图
  7. 硬盘分区工具找不到硬盘该怎么办
  8. ZABBIX3.04设置短信报警
  9. 网络安全之密码安全基础
  10. 深聊性能测试,从入门到放弃之:如何对IO进行性能调优
  11. Tagged Pointer分析
  12. uni-app 页面生命周期
  13. 刀柄锥度气动量仪检测
  14. 自考深圳大学计算机专业难吗,深圳大学全日制自考本科有用吗?计算机与软件学院自考办通过率...
  15. 基于STM32的温度遥测系统
  16. 强化学习和 Q-Learning
  17. NUMECA系列: FINE / Open 10.1 最新版64位CFD流体集成环境
  18. 新手小白可做的22个Python迷你项目(附源码),建议收藏练习
  19. 各省、地级市-2020春运期间高德地图人口迁徙数据-excel、shp格式
  20. 【Oracle】触发器

热门文章

  1. 记一次事务报错问题 Transaction synchronization is not active
  2. BZOJ - 2244 拦截导弹 (dp,CDQ分治+树状数组优化)
  3. 必须要调整心态,积极起来,不能再偷懒
  4. UA MATH564 概率论VI 数理统计基础4 t分布
  5. 求正负数交互序列的和的MFC版
  6. VC++ 添加UNICODE和_UNICODE预处理选项
  7. 【洛谷p1058】立体图(已完结)
  8. AES CBC模式下的CBC bit flipping Attack
  9. 【STM32H7教程】第4章 STM32H7工程模板建立(MDK5)
  10. J20170521-ts