matlab 多次求解偏微分方程 ode45

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

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

 

其中，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