- 综合排序
- 最热优先
- 最新优先
- 来自专栏javascript趣味编程
导热控制偏微分方程的推导Edition2
导热控制偏微分方程的推导： 板书： ? 直角坐标系导热控制偏微分方程的推导 ? 直角坐标系导热控制偏微分方程的推导及柱坐标系的补充 视频中的草稿： ? 几个《传热学》相关的小程序总结如下，可在微信中点击体验： 有限元三角单元网格自动剖分 Delaunay三角化初体验 （理论戳这） Contour等值线绘制 （理论戳这） 2D非稳态温度场有限元分析 1D稳态导热温度场求解 （源码戳这) 1D非稳态导热温度场求解程序 （源码戳这) 2D稳态导热温度场求解 （源码戳这) 普朗克黑体单色辐射力 《传热学》相关小程序演示动画如下（其中下图1D非稳态导热计算发散 《(计算)流体力学》中的几个小程序，可在微信中点击体验： Blasius偏微分方程求解速度边界层 （理论这里） 理想流体在管道中的有势流动 （源码戳这） 涡量-流函数法求解顶驱方腔流动
97120发布于 2019-12-05 - 来自专栏算法之名
微分方程整理
微分方程 例：求一曲线方程，使其满足过点(1,2)，且其上任意一点处的切线斜率为其横坐标的2倍。 设曲线方程为y=f(x)，有\({dy\over dx}=2x\)，这样的方程与代数方程的不同在于包含了未知函数对自变量的导数，为一阶的微分方程。 由该式可得 dy=2xdx 两端积分，有 \(y=x^2+C\) 将(1,2)代入，可得 C=1 故该曲线方程为 \(y=x^2+1\) 其中\(y=x^2+C\)为微分方程\({dy\over dx} 基本概念 微分方程分为常微分方程和偏微分方程，之前的示例就为常微分方程，偏微分方程例如 \({∂^2u\over ∂x^2}+{∂^2u\over ∂y^2}=0\) 的多元函数的方程。 这里又把微分方程分为一阶微分方程和高阶微分方程。 从线性和非线性的角度，又可以把微分方程分为线性方程和非线性方程。
60410编辑于 2023-11-24 - 来自专栏全栈程序员必看
matlab微分方程组_matlab求微分方程特解
主要内容：matlab参数识别应用，主要适用于微分方程、微分方程组参数识别、simulink模型参数识别，领域不限。 1 使用matlab识别微分方程参数以及微分方程组(多个微分方程)参数 2 使用matlab调用simulink并识别simulink模型的参数(m函数与simulink交互) 内容为本人在学习过程中总结的知识 当然也可以用xdata1=sim(‘pm’,tspan,options,ut); clc; clear; global k1 global k2 global xdata; k1=3; k2=-4; kc =[k1,k2] tspan=[0,3]; load_system(‘pm’) opt=simget(‘pm’); x=min(tspan):opt.FixedStep:max(tspan); xdata ,tspan); plot(x’,yest,’or’,x’,ydata,’b*’) pm.jpg (13.7 KB, 下载次数: 10) 2014-6-23 23:14 上传 待识别模型,k1,k2参数
73510编辑于 2022-10-04 - 来自专栏CSDN小华
数学建模--微分方程
求解微分方程：对于能够求得解析解的微分方程，可以直接求解；对于复杂的微分方程，则需要利用数值方法进行近似求解。 dydt = -2 * y # 以 dy/dt = -2y为例 return dydt # 设置初始条件 y0 = 1 # 设置时间点 t = np.linspace (0, 5, 100 grid.copy() for i in range(1, N-1): for j in range(1, N-1): new_grid[i, j] += h**2 常微分方程（ODE）与偏微分方程（PDE）在数学建模中的优缺点分别是什么？ 在数学建模中，常微分方程（ODE）和偏微分方程（PDE）各有其优缺点。 在进行微分方程模型求解时，选择最有效的数值方法取决于微分方程的类型和复杂性。
1.8K10编辑于 2024-10-16 - 来自专栏往期博文
数学建模暑期集训5：matlab求解常微分方程偏微分方程
1.Matlab求常微分方程的数值解 1.1非刚性常微分方程的数值解法： 功能函数：ode45，ode23，ode113 例：用RK方法（四阶龙格—库塔方法）求解方程 f=-2y+2x^2+2*x 2.Matlab求常微分方程的解析解 2.1求常微分方程的通解 syms x y diff_equ='x^2+y+(x-2*y)*Dy=0' dsolve(diff_equ,'x') 注：'x ’代表x为自变量，D代表求导 2.2求常微分方程的初边值问题 syms x y diff_equ='D3y-D2y=x' dsolve(diff_equ,'y(1)=8,Dy(1)=7,D2y(2)= 4','x') 2.3求常微分方程组 equ1='D2f+3*g=sin(x)'; equ2='Dg+Df=cos(x)'; [general_f,general_g]=dsolve(equ1,equ2 ,t); subplot(2,2,2),pdesurf(p,t,u) subplot(2,2,3),pdesurf(p,t,u-exact') 4.Matlab pdetool工具箱求解偏微分方程 对于一般的区域
1.7K20编辑于 2022-06-14 - 来自专栏灰灰的数学与机械世界
考研（大学）数学 微分方程（1）
微分方程（2） 第三节 可降阶的高阶微分方程 ---- 3.1形如 y^{n}=f(x) 的方程 方法：对方程直接左右两边进行不定积分，重复 n 次。 解：根据 yy^{'}-y^2=1\Leftrightarrow 2yy^{'}-2y^2=2 ，即 (y^{2})^{'}-2y^2=2 ，则 \displaystyle y^2=(\int2e^{\ int-2dx}+C){e^{-\int -2dx}}=(C-e^{-2x})e^{2x}=Ce^{2x}-1 ，根据 y(0)=0 ，解出 C=1 ，所以 y^2=e^{2x}-1 ---- 10. }} ，由题意 y(0)=1 ,解得 C_{2}=1 ，所以特解 \displaystyle y=e^{\frac{x^2}{2}} . ---- 基础题二 7.求微分方程 (2x+3)y^{''}=4y 前两题是关于常微分方程的特殊方法，一个是凑微分，另外一个利用导数的除法公式；化成常见的方程，例如一阶齐次线性微分方程和一阶非齐次线性微分方程，再利用初始条件，得出解；后面两题是关于缺 x 型和 y 型的真正解法
69140编辑于 2022-11-14 - 来自专栏AI科技评论
求解微分方程，用seq2seq就够了，性能远超 Mathematica、Matlab
论文地址：https://arxiv.org/abs/1912.01412 这篇论文提出了一种新的基于seq2seq的方法来求解符号数学问题，例如函数积分、一阶常微分方程、二阶常微分方程等复杂问题。 由于seq2seq模型的特点，作者所提方法能够对同一个公式得出不止一个的运算结果，例如如下的微分方程 ? 该模型能够反馈这么多的结果： ? 可以验证一下，这些结果都是正确的，至多差一个常数 c。 学过高等数学的我们都有过求积分和解微分方程的痛苦经历，对计算机软件来讲，求解这些问题事实上也同样困难。 基于这种思路，作者首先提出了将数学表达式转换为seq2seq表示形式的方法，并用多种策略生成了用于监督学习的数据集（积分、一阶和二阶微分方程），然后将seq2seq模型用于这些数据集，便得出了比最新计算机代数程序 理想情况下，我们应该生成问题空间的代表性样本，即随机生成要积分的函数和要求解的微分方程。但我们知道，并不是所有的函数都能够积分（例如f=exp(x^2)和f=log(log(x))）。
1.5K10发布于 2019-12-19 - 来自专栏全栈程序员必看
常微分方程初值问题数值解法MATLAB(泛函微分方程)
Matlab 解常微分方程的初值问题 题目：Matlab 解常微分方程的初值问题 设计目的： 1、熟练掌握Matlab的基本编程方法，及其编程风格。 2、熟练掌握Matlab常用函数的使用。 设计内容： 已知一个三阶微分方程:，利用matlab软件求这个三阶微分方程在初值 下的解。 原三阶微分方程可化为： 令 则原三阶微分方程可化为微分方程组 在初值 下的解。 dy(2) =y(3); dy(3) = 2*(1-y(1)^2)*y(3)-y(1)*y(2); %调用函数ode45求解，时间区间为[0,10] [t,Y] = ode45(@rigid,[0 10 参考文献： [1] 张圣勤编 MATLAB7.0 机械工业出版社 [2]周义仓 靳祯 秦军林编 常微分方程极其应用 科学出版社 [3]韩明 王家宝 李林编 数学实验(matlab版) 同济大学出版社 [ 4]汪晓银 皱庭荣编 数学软件与数学实验 科学出版社 2、把这个三阶微分方程化为形如 的标准形式 1、已知一个三阶微分方程 3、编写函数文件rigid.m 4、调用函数文件rigid.m，利用ode45
1.3K20编辑于 2022-07-28 - 来自专栏灰灰的数学与机械世界
考研（大学）数学 微分方程（3）
微分方程（3） 第四节 高阶微分方程 ---- 4.1 高阶齐次线性微分方程 4.1.1 高阶齐次微分方程的基本概念 1.n阶齐次线性微分方程的定义 例如 y^{n}+a_{1}(x)y^{n-1}+\ dotsb+a_{n-1}(x)y^{'}+a_{n}(x)y=0 \qquad (1) 称为n阶齐次线性微分方程 2.n阶非齐次线性微分方程的定义 例如 y^{n}+a_{1}(x)y^{n-1}+\ \varphi_{2}(x)+\dotsb+k_{n}\varphi_{n}(x)+\varphi_{0}(x) 是 (2) 的一个解 ---- 4.1.3 高阶常系数微分方程 1.二阶常系数齐次微分方程的解法 +2=0 ,解得特征值为 \lambda_{1,2}=1\pm i ,则原方程的通解为 y=e^{x}(C_{1}\sin x+C_{2}\cos x) ---- 2.二阶常系数非齐次线性微分方程的特解 lambda_{2}=-2 ,所以对应的齐次线性微分方程的通解为： y=(C_{1}+C_{2}x)e^{-2x} ( C_{1},C_{2} 均是常数），由于在方程的右边含有参数，故对 a 进行讨论
74810编辑于 2022-11-14 - 来自专栏司六米希
【高等数学】【8】微分方程
【高等数学】【8】微分方程 1. 微分方程的基本概念 1.1 微分方程 2.可分离变量的微分方程 3. 齐次方程 4.一阶线性微分方程 4.1 线性方程 5. 可降阶的高阶微分方程 6. 高阶线性微分方程 7. 常系数齐次线性微分方程 8. 常系数非齐次线性微分方程 1. 微分方程的基本概念 1.1 微分方程 2.可分离变量的微分方程 例子 3. 齐次方程 4.一阶线性微分方程 4.1 线性方程 5. 可降阶的高阶微分方程 6. 高阶线性微分方程 7. 常系数齐次线性微分方程 8. 常系数非齐次线性微分方程
55220编辑于 2022-11-15 - 来自专栏全栈程序员必看
求微分方程的特解matlab_二阶微分方程求解
求解微分方程 desolve函数 实例1 实例2 实例3 实例4 求解有条件的微分方程 微分方程显示隐式解 未找到显式解决方案时查找隐式解决方案 求微分方程级数解 为具有不同单边限制的函数指定初始条件 （特解） 练习题 desolve函数 S = dsolve(eqn)求解微分方程eqn，其中eqn是符号方程。 使用diff和==来表示微分方程。例如，diff(y,x) == y表示方程dy / dx = y。通过指定 eqn为这些方程的向量来求解微分方程组。 ,{\mathrm{e}}^{-\sqrt{a}\,t}+C_{2}\,{\mathrm{e}}^{\sqrt{a}\,t} C1e−a t+C2ea t 求解有条件的微分方程 d =-1附近微分方程的级数解，请将“ExpansionPoint”设置为 -1。
1.4K10编辑于 2022-11-17 - 来自专栏全栈程序员必看
matlab求解微分方程组(matlab解微分方程的数值解)
如何用matlab来求解简单的微分方程？举例来说明吧。 求解三阶常微分方程。我们知道，求解高阶常微分方程可以化为求解一阶常微分方程组。 编写函数eq3.m： %解常微分方程 3*y'''+5*y''+6*sin(t)*y=cost function ydot = eq3(t,y) ydot=[y(2);y(3);(cos(t)-5*y 如图： 二阶常微分方程 编写函数eq2.m function ydot= eq2(t,y) ydot=[y(2);-3-cos(2*t) + 2*sin(t)+t-3.8]; 主函数 clc y(0)=0 y'(0)=1 y''(0)=-1 求无初始条件的微分方程的解析通解各项 clc clear syms x y diff_equ='x^2+y+(x-2*y)*Dy=0'; dsolve (diff_equ,'x') %求无初始条件的微分方程的解析通解各项 求线性系统的解析解并画相图 clc,clear equ1='Dx1 - x2 = 0'; equ2='Dx2 + x1 + 2*
2.5K30编辑于 2022-08-01 - 来自专栏数据结构和算法
Scipy 中级教程——积分和微分方程
演示了如何使用 Scipy 进行定积分： import numpy as np from scipy import integrate # 定义被积函数 def func(x): return x**2 2. 微分方程求解 Scipy 提供了 odeint 函数用于求解常微分方程组。 ) plt.title('简单的一阶微分方程求解') plt.show() 在这个例子中，model 函数定义了一阶微分方程 dy/dt = -y。 更复杂的微分方程 如果需要求解更复杂的微分方程组，可以通过定义更复杂的 model 函数和初始条件，然后使用 odeint 函数进行求解。 总结 Scipy 提供了强大的积分和微分方程求解工具，方便科学计算和工程应用。通过这篇博客的介绍，你可以更好地理解和使用 Scipy 中的积分和微分方程求解功能。
88610编辑于 2024-01-10 - 来自专栏巴山学长
带你用matlab轻松搞定微分方程
联系一些未知函数的一组微分方程称为微分方程组。微分方程中出现的未知函数的导数的最高阶称为微分方程的阶。 有些微分方程比较简单可直接通过积分求解。例如一阶常系数线性常微分方程： ? x) + diff(y(x), x, x) == 0 S = C4*exp(-x)*cos(3^(1/2)*x) + C5*exp(-x)*sin(3^(1/2)*x) 演示了两个比较简单的微分方程用符号解微分方程的方法解出通解 考虑一阶常微分方程组初值问题： ? 其中y=(y1,y2,...,ym)T，f=(f1,f2,...,fm)T，y0=(y10,y20,... 代码如下： %第一种解微分方程的方法 syms y(x) a b eqn1 = diff(y,x,1)==-2*y(x)+2*x^2+2*x; eqn2 = y(0)==1; y1= dsolve([eqn1 ,eqn2]); %第二种解微分方程的方法 y2=dsolve('Dy=-2*y+2*x^2+2*x','y(0)=1','x'); %第三种解微分方程的方法 warning off feature jit
1.9K31发布于 2020-11-03 - 来自专栏Python编程 pyqt matplotlib
差分法求解微分方程
差分法求解微分方程的示例： import numpy as np from numpy import exp from sympy import symbols, solve, Eq from math P[0] + 4 * P[1] - P[2]) / (2 * dx) D2p_dx2[0] = (P[0] - 2 * P[1] + P[2]) / (dx * dx) # 中间节点 for i in range(1, n-1): Dp_dx[i] = (P[i+1] - P[i-1]) / (2 * dx) D2p_dx2[i] = (P[i+1] - 2*P[i] + P[i-1 ]) / (dx * dx) # 右边界 Dp_dx[n-1] = (3*P[n-1] - 4*P[n-2] + P[n-3]) / (2 * dx) D2p_dx2[n-1] = (P[n-1] - range(1, n-1): # 内部满足微分方程 eq = Eq(D2p_dx2[i] + Dp_dx[i] - 2 * P[i], 0) eqs.append(eq) # print
38310编辑于 2024-04-25 - 来自专栏YoungGy
微分方程与欧拉法
微分方程概述 微分方程在各个领域应用颇多。 形如 [图片] 的微分方程表示了系统的变化信息， 如果在加上初始条件(x0,y0)，那么就可以求出系统整体随时间变化的信息。 可以说，正是微分方程将物理世界模型化。 y的方向场与积分曲线') hold off 微分方程的解析解法 微分方程的解析解法通常是将x,y分别移到等式的一边。 下面以 [图片] 为例，移项后 [图片] 所以有 [图片] 进而有 [图片] 最后解得： [图片] 其实， [图片] 就是根据微分方程y′=y在(0,1)(0,1)的初始条件下确定的 =1','x') %输出为： exp(x) 微分方程的数值解法 欧拉法 欧拉法的核心是，设定步长为h，然后已知y′和(x0,y0)，根据下面方法迭代： [图片]
1.4K50发布于 2018-01-03 - 来自专栏全栈程序员必看
matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究
固定时滞的微分方程：满足下面的形式，也就是微分方程右边包含时滞部分，且时滞为常数。 使用dde23函数求解： 问题： (1)微分方程定义：多了一个时滞部分 创建myddefun.m文件，文件里的内容如下： function dy = myddefun(t,y,Z) dy=[ Z(1,1); Z(1,1)+Z(2,2); y(2) ]; 说明：其中Z表示时滞部分：即y(t-T)，T是时滞。 Z(i,j)表示y(i)(t-T)，即y(i)的时滞形式；j表示T选取第j个时滞值 (2)外部调用方程，输入参数求解 lags=[1,0.2]; history=[1;1;1]; tspan=[0,5] ; sol = dde23(@myddefun,lags,history,tspan) 其中lags表示时滞数组，如上面输入了2个时滞，Z(1,1)表示选用第一个时滞T=1，Z(2,2)表示选取第二个时滞
1.1K40编辑于 2022-11-10 - 来自专栏全栈程序员必看
小程序：matlab解微分方程「建议收藏」
* y(6)); W = sqrt(y(7) * y(7) + y(8) * y(8) + y(9) * y(9)); f = zeros(9, 1); f(1) = y(4); f(2)
49930编辑于 2022-07-31 matlab求解时滞微分方程
在MATLAB中求解时滞微分方程（DDE），可以使用dde23函数。这个函数是为了解决具有常数时滞的微分方程而设计的。 使用dde23求解时滞微分方程的基本步骤： 定义时滞微分方程：首先，需要定义时滞微分方程的函数。这个函数应该接受当前时间t，当前状态y，以及在时滞点的状态Z作为输入，并返回导数dydt。 调用dde23函数：使用dde23函数求解时滞微分方程。需要提供时滞微分方程的函数、时滞、初始条件函数、以及求解的时间范围。 绘制结果：求解完成后，可以使用plot函数绘制结果。 下面的例子，展示如何使用dde23求解一个具有常数时滞的微分方程： % 定义时滞微分方程 function dydt = ddefun(t, y, Z) dydt = y - Z; end % 定义时滞 tau = 0.5; % 定义初始条件 function y = history(t) y = 1; end % 求解时滞微分方程 sol = dde23(@ddefun,
29810编辑于 2025-07-24- 来自专栏大白技术控的技术自留地
使用Maxima求解常微分方程~
使用Maxima求解常微分方程~ 含带导数符号或带微分符号的未知函数的方程称为微分方程。 如果在微分方程中未知函数是一个变元的函数，这样的微分方程称为常微分方程。 (%i5) eq3:'diff(y,x,3)+y=0; sol3:ode2(eq3,y,x); ? 上面的例子用了ode2函数来求解常微分方程。 sol1 中的%c 和 sol2 中的 %k1 %k2 是任意常数。 ode2函数只能求解一阶和二阶常微分方程，第三个例子给出的是一个三阶常微分方程，无法求解，因此输出 false。 2 初值问题 函数ic1 (solution, xval, yval)和ic2 (solution, xval, yval, dval)分别用来解一阶和二阶微分方程的初值问题，其中solution是用 3 边值问题 函数bc2 (solution, xval_1, yval_1, xval_2, yval_2)用来求解二阶微分方程的边值问题， 其中solution是ode2解得的通解，xval_1
2.1K20发布于 2019-03-05
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号