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

 

 

 

总结

以上是生活随笔为你收集整理的matlab 多次求解偏微分方程 ode45的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。