師兄和我討論了一個問題，就是在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的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。