【Matlab】解个微分方程
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=alv˙f=−m0.1+5∗vf+0.25∗vf2D˙=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】解个微分方程相关推荐
- 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 ...
- matlab ode45 二阶微分,matlab关于ode45解二阶微分方程的困惑
matlab关于ode45解二阶微分方程的困惑 matlab关于ode45解二阶微分方程的困惑 一个二阶微分方程: y''+y'+y=sin(t) 初始条件为y(0)=5,y'(0)=6. 过程: 先 ...
- matlab 向前欧拉公式,向前欧拉公式在Matlab解微分方程初值解的问题
向前欧拉公式在Matlab解微分方程初值解的问题0 fuqilin1202013.07.04浏览527次分享举报 用向前欧拉公式(10.8)求解初值问题,dy/dx=-3x+8x-7,y(0)=1,分 ...
- matlab求微分方程的初值问题,Matlab 解常微分方程的初值问题.doc
Matlab 解常微分方程的初值问题 题目:Matlab 解常微分方程的初值问题 设计目的: 1.熟练掌握Matlab的基本编程方法,及其编程风格. 2.熟练掌握Matlab常用函数的使用. 3.与本 ...
- matlab解无解析解微分方程组,数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解...
<数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解>由会员分享,可在线阅读,更多相关<数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解(12页 ...
- matlab求偏微分方程程序,用MATLAB解偏微分方程.pdf
用MATLAB解偏微分方程.pdf 年 月 阴 山 学 刊 第 卷 第 期 丫叫 加 用 解偏微分方程 田 兵 包头师范学院 学报编辑部 , 内蒙古 包头 摘 要 讨论 了以 中偏徽分方程工具箱的用法 ...
- Matlab符号数学(Symbolic Math with MATLAB)MATLAB解方程
创建符号变量 用符号而不是数值表示数字. 存储数字符号 % 你可以把一个数存储为数值格式,就像我们经常做的那样 x = 3; % 你也可以把这个数存储为符号格式,使用sym()函数 xsym = sy ...
- 如何用matlab解异或方程,Matlab-6:解非线性方程组newton迭代法
函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\ ...
- 用MATLAB求解非线性微分方程
总结一下MATLAB中求解微分方程的思路和步骤.固然,网上很多关于此类的技术型文章,但往往一看下来发现,文章中的友情链接比文章字数还多,要了解这一篇文章,你要先了解那个:要了解那个,你又要了解那个那个 ...
- matlab 解函数方程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
最新文章
- HTML5column属性布局页脚,利用column多列属性调整页面文字列布局
- Jupyter Notebook 使用小记
- JDBC批量操作批量增加批量修改
- [swift]-类的构造函数
- char*赋值给std::string是深拷贝
- matplotlib——直方图
- 硬盘分区工具找不到硬盘该怎么办
- ZABBIX3.04设置短信报警
- 网络安全之密码安全基础
- 深聊性能测试,从入门到放弃之:如何对IO进行性能调优
- Tagged Pointer分析
- uni-app 页面生命周期
- 刀柄锥度气动量仪检测
- 自考深圳大学计算机专业难吗,深圳大学全日制自考本科有用吗?计算机与软件学院自考办通过率...
- 基于STM32的温度遥测系统
- 强化学习和 Q-Learning
- NUMECA系列: FINE / Open 10.1 最新版64位CFD流体集成环境
- 新手小白可做的22个Python迷你项目(附源码),建议收藏练习
- 各省、地级市-2020春运期间高德地图人口迁徙数据-excel、shp格式
- 【Oracle】触发器
热门文章
- 记一次事务报错问题 Transaction synchronization is not active
- BZOJ - 2244 拦截导弹 (dp,CDQ分治+树状数组优化)
- 必须要调整心态,积极起来,不能再偷懒
- UA MATH564 概率论VI 数理统计基础4 t分布
- 求正负数交互序列的和的MFC版
- VC++ 添加UNICODE和_UNICODE预处理选项
- 【洛谷p1058】立体图(已完结)
- AES CBC模式下的CBC bit flipping Attack
- 【STM32H7教程】第4章 STM32H7工程模板建立(MDK5)
- J20170521-ts