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

 

 

 

©️2020 CSDN 皮肤主题: 游动-白 设计师: 上身试试 返回首页
实付0元
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值