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是一階微...