matlab 多次求解偏微分方程 ode45

师兄和我讨论了一个问题，就是在matlab中求解偏微分方程，

其中，偏微分方程中有的常数是一直变化的，要求很多次，而不是一个固定的常数求一次就行了。

$\frac{dA1}{dx} = k * A2 * exp(-i\delta x)$

$\frac{dA2}{dx} = k * A1 * exp(-i\delta x)$

其中，A1和A2是要求解的因变量，x是自变量，其他为常数

求解微分方程如下：

%定义函数如下
function  c = abss(x,A)
k=1;
% global k;
deta = 1;
c(1,1) =k*A(2)*exp(-1i*deta*x);
c(2,1) =k*A(1)*exp(1i*deta*x);
% c = c';
end%注意 函数abss返回的值   必须返回列矢量
%所以c为2行1列，如果为1行2列，则 c = c'; %求解
xx = [1 2];
yy = [1,0]';
d,e = ode45(@abss,xx,yy);

当，我们想让k的值的变化的，求解不同k值下的解，则设定k为全局变量。

而不能在函数内加入k使得k成为变量。

function  c = abss(x,A)
%k=1;
global k;
deta = 1;
c(1,1) =k*A(2)*exp(-1i*deta*x);
c(2,1) =k*A(1)*exp(1i*deta*x);
% c = c';
end%求解
xx = [1 2];
yy = [1,0]';
for k = 1:10;[d{k},e{k}] = ode45(@abss,xx,yy);
end

matlab 多次求解偏微分方程 ode45相关推荐

matlab采取有限差分法求解偏微分方程
原理如上求解以下方程 function u=ExplicitFDM_Heat(M,N)if(~exist('M','var'))M = 50; endif(~exist('N','var'))N = ...
matlab中PDE工具箱如何使用,使用PDE工具箱求解偏微分方程
在科学技术各领域中,有很多问题都可以归结为偏微分方程问题.在物理专业的力学.热学.电学.光学.近代物理课程中都可遇见偏微分方程. 偏微分方程,再加上边界条件.初始条件构成的数学模型,只有在很特殊情况下 ...
matlab 离散偏微分,有限差分法求解偏微分方程MATLAB.doc
有限差分法求解偏微分方程MATLAB.doc 南京理工大学课程考核论文课程名称: 高等数值分析论文题目: 有限差分法求解偏微分方程姓名: 罗晨学号: 115104000545 成绩: ...
matlab pdepe函数边界,[求助]用pdepe求解偏微分方程时遇到分段边界条件的问题
[求助]用pdepe求解偏微分方程时遇到分段边界条件的问题用Matlab中的pdepe函数求解偏微分方程时,左边界条件是一个分段函数.左边界开始的时候是1,在t时刻后,左边界突变为0.在写左边界的时 ...
matlab 文件名分离_MATLAB偏微分方程
4.8.2 偏微分方程在自然科学的很多领域内,都会遇到微分方程初值问题,特别是偏微分方程,它的定解问题是描述自然界及科学现象的最重要的工具.可以说,几乎自然界和各种现象都可以通过微分方程(特别是偏 ...
10种基于MATLAB的方程组求解方法
线性方程组的求解包括直接法和迭代法,其中迭代法包括传统的高斯消元法,最速下降法,牛顿法,雅克比迭代法,共轭梯度法,以及智能启发式算法求解法和神经网络学习算法,传统算法可以相互组合改进,智能仿生启发式算 ...
偏微分方程数值解法python_基于python求解偏微分方程的有限差分法资料
基于python求解偏微分方程的有限差分法资料 Computer Era No. 11 2016 0 引言在数学中, 偏微分方程是包含多变量和它们的偏导数在内的微分方程.偏微分方程通常被用来求解 ...
数学建模——差分算法(求解偏微分方程)
差分算法(求解偏微分方程) 差分算法是数学建模比赛中的一种十分常见的代码,在2018A题和2020A中均用到一维热传导模型,模型的求解用的就是差分算法,具体如何解可以自己去查看相关论文. 定义差分方 ...
matlab微分方程求法,matlab微分方程的求解的方法ppt课件
<matlab微分方程的求解的方法ppt课件>由会员分享,可在线阅读,更多相关<matlab微分方程的求解的方法ppt课件(44页珍藏版)>请在人人文库网上搜索. 1.定义:含 ...

matlab 多次求解偏微分方程 ode45

matlab 多次求解偏微分方程 ode45相关推荐

最新文章

热门文章