Matlab求解微分方程組

我們採用ode方法：

（1）求解普通微分方程組：使用ode45方法

1. 建立乙個函式檔案eq2.m，在函式檔案中描述這個解的微分方程組：

%eq2.m檔案

%描述微分方程組
function dy=eq2(t,y) 
%說明微分變數是二維的，令y(1)=x,y(2)=y
dy=zeros(2,1); 
%微分方程組
dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2);
dy(2)=5*(1-y(2))/sqrt((1-y(1))^2+(t-y(2))^2);
end

[t,y]=ode45(@eq2,[0,2],[0,0]);

[t,y]=ode45(@eq2,[t1,t2],[y1(0),y2(0)]);

3.繪製響應曲線

plot(t,y(:,1),'-',t,y(:,2),'*')

微分方程 微分方程 高階微分方程組理論

附 微分方程的分類 常微分方程和偏微分方程。1 常微分方程 ode 是指微分方程的自變數只有乙個的方程。最簡單的常微分方程，未知數是乙個實數或是複數的函式，但未知數也可能是乙個向量函式或是矩陣函式，後者可對應乙個由常微分方程組成的系統。2 偏微分方程 pde 是指微分方程的自變數有兩個或以上，且方程...

matlab求解時滯微分方程

matlab求解時滯微分方程，dde23呼叫格式 sol dde23 ddefun,lags,history,tspan ddefun函式控制代碼，求解微分方程y f t,y t y t 1 y t k 必須寫成下面形式 dydt ddefun t,y,z 其中t對應當前時間t，y為列向量，近似於y...

matlab練習程式（常微分方程組向量場）

過去有畫過常微分方程的向量場，通過向量場能夠很形象的看出方程解的狀態。這裡用matlab也實現一下，同時對三維情況也做了乙個實現。繪製的方法就是計算方程在二維或三維某個點的方向，然後把方向歸一化，畫出歸一化的向量即可。二維微分方程組如下 三維微分方程組如下 matlab 如下 二維情況 clear ...