师兄和我讨论了一个问题,就是在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