matlab問題使用ode45求如下微分方程的數值解繪

2021-03-27 08:01:36 字數 2567 閱讀 7134

1樓:匿名使用者

用ode45()求微分方程的格式為

[t,y] =ode45(odefun,tspan,y0)odefun——微分方程自定義函式;tspan——t的區間;y0——y,dy的初始值。

用plot()繪出微分方程的解,即y(t)函式圖形。

對於本提問,可以按下列**實施。

下圖是精確解與數值解的圖形比較。

matlab中用ode45求微分方程的數值解問題

2樓:匿名使用者

步長在起始時間和結束時間中間設定,如想設定步長為0.1,可以用

[t,y]=ode45('function1',[0 0.1 3000],[2 0]);

即在0和3000之間版加上步長,如果不設定步權長,預設為1。

用matlab解微分方程數值解並作圖

3樓:高等數學答案

先定義個函

bai數:

function f=fffff(t,y)f=[y(2);cos(y(1))];

然後du

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

plot(t,y(:,1),'r',t,y(:,2),'b')即可。我給你設了初值zhi【0,0】,你可以自己定dao義。紅線版為y,藍線為y'.

你先把權上面的函式fffff寫成m檔案,儲存到你的matlab工作目錄裡,然後再在命令視窗輸入:

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

plot(t,y(:,1),'r',t,y(:,2),'b')我作過,能做出影象的。

如何用matlab解二元微分方程組,分別做出x-t,y-t的座標圖? 10

4樓:king手舞足蹈我

% 求解微分方程(時間copy

範圍0-2)

dx=inline('[(2-3*x(1)+x(2)).*x(1); (4-x(2)+2*x(1)).*x(2)]','t','x');

[t,x]=ode45(dx,[0 2],[0;3])y=x(:,2);

x=x(:,1);

% 繪製x-t, y-t曲線

% 注意:由於初始條件x(0)=0,計算出導數dx/dt=0,求出的x一直為0

figure(1)

subplot(2,1,1)

plot(t,x)

xlabel('time (sec)')

ylabel('x')

subplot(2,1,2)

plot(t,y)

xlabel('time (sec)')

ylabel('y')

% 繪製x-y曲線

% 由於x一直為0,x-y圖為一條豎線

figure(2)

plot(x,y)

xlabel('x')

ylabel('y')

在matlab中用ode45命令求解微分方程組,畫隨時間變化的曲線。

5樓:匿名使用者

function test()

p0=0.6;q0=0.7;

a1=0.07; b1=0.04; s1=0.3; k1=0.5; m1=1800; w1=0.5; c1=3;

a2=0.04; b2=0.02; s2=0.4; k2=0.4; m2=1500; w2=0.4; c2=2;

u1 = (a1 - b1) * s2 * k2 * m2 + w1;

u2 = (a2 - b2) * s1 * k1 * m1 + w2;

[t, y] = ode45(@func, [0, 5], [p0; q0], , u1, -c1, u2, -c2);

plot(t, y);

function f=func(t, y, ua, ub, va, vb)

f = [y(1) * ( 1 - y(1)) * (y(2) * ua + ub); y(2) * (1 - y(2)) * (y(1) * va + vb)];

matlab ode45解二階常微分方程怎麼把y算出來

6樓:匿名使用者

你說的y是什麼概念?

程式中有繪圖啊,那條藍線就是y。

matlab使用ode45解微分方程組的基本問題,下面有一個簡單的微分方程組求**

7樓:匿名使用者

^建立myfun.m檔案,把以下**

複製進去

function dx=myfun(t,x)dx=zeros(2,1);

dx(1)=4*x(1)+x(2).^2;

dx(2)=8*x(1)+2*x(2).^2+5*x(2);

執行以內下**

t=[0 0.1];

x0=[1 1];

[tt,xx]=ode45('myfun',t,x0);

plot(tt,xx(:,1),'-',tt,xx(:,2),'-.')%畫出兩者影象容

怎麼用ode45求解常微分方程,matlabode45求解二階常微分方程

上述微分方程組比較複雜,用dsolve求解解析解非常耗時,且不一定可行。在實際應用中,對內於複雜的微分容方程 組 很多時候解析解是很難得到的,因此經常以求解數值解來代替解析解。本節介紹一個求解微分方程數值解的函式 ode45。ode45採用四階和五階runge kutta單步演算法,用變步長求解器求...

微分方程問題,見下圖,高數。求微分方程的通解。題目見下圖。

y 2 y f 0,x y t dt 1 等式兩邊求導 2yy y y 2 y 2 y 0 2yy y y 2 y 3 0 同除以y 2 y y 3 y 2 2y 2 y 0 設y e t y dy dt dt dx e t t y e t t 2 e t t t 2e t t e t t 3e 3...

微分方程的通解問題,微分方程的通解怎麼求

這和微分方程的階數有關,微分方程是幾階的,通解中就有幾個常數項,直觀上是好理解的,微分方程就是還有變數導數的方程,解方程的過程和不定積分類似 y x可看做是最簡單的微分方程 導數最高是幾階,就要積分幾次,而不定積分每進行一次,表示式中就多出一個常數c。例如你說的那兩個微分方程,dy 2xdx是一階微...