1樓: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畫乙個微分方程組的圖,具體要求如下
2樓:匿名使用者
題主給出的微分方程組缺初始條件。為了解題方便,自行給出其初始條件x(0)=1,y(0)=0.5。
用matlab求解微分方程組可以用ode()函式求得其數值解。求解方法:
1、首先,自定義微分方程組函式odefun(t,x),其主要內容
dy(1)=x(1)*(1-x(1))-x(2)+1/x(1);
dy(2)=x(2)*(-1+x(1));
2、再用ode()函式求得其數值解。使用方法
[t,x]=ode45(@odefun,[0 20],y0);
3、用plot()函
數繪製t—x'(t)、y'(t)和x'(t)—y'(t)相平面圖
計算結果
請問在matlab中如何畫出此二階微分方程組的y-x的解?
3樓:二中混過
^^令x1 = x, x2 = x', x3=y, x4=y'
得zhi
x1' = x2;
x2' = -x2*sqrt(x2^dao2+x4^2);
x3' = x4;
x4' = -x4*sqrt(x2^2+x4^2) - 9.8;
程式為版
f = @權(t, x)[x(2);
-x(2)*sqrt(x(2).^2+x(4).^2);
x(4);
-x(4)*sqrt(x(2).^2+x(4).^2)-9.8];
[t, y ]= ode45(f, [0,5], [0;35.355;0;35.355]);
plot(t,y)
legend('x1', 'x2', 'x3', 'x4')
matlab問題:使用ode45求如下微分方程的數值解,繪製y-t圖
4樓:匿名使用者
用ode45()求微分方程的格式為
[t,y] =ode45(odefun,tspan,y0)odefun——微分方程自定義函式;tspan——t的區間;y0——y,dy的初始值。
用plot()繪出微分方程的解,即y(t)函式圖形。
對於本提問,可以按下列**實施。
下圖是精確解與數值解的圖形比較。
二元函式的全微分求積,高數二元函式的全微分求積
看圖,來ab段的方程為y 0 將y 0代入 源積分後,對於dy來說,由於y是常數,dy就是0,因此這個積分為0,不用計算 對於dx這個積分來說,由於前面乘了個y,因此y 0代入後結果也為0,所以ab段的積分為0.高數 二元函式的全微分求積 類似於積分上限函式,這裡需要利用二元函式的全微分求積,先證明...
解微分方程mxkx,如何用二階微分方程來解mx kx 0,聽說用這個很簡單,但是不知道怎麼用,請高手詳細解答,謝謝
又是你提問來了麼?呵呵,你還在對那個簡諧運動的問題耿耿於懷吧!你應該說明白x 是位移x對時間t的二階變化率,也就是加速度。否則會引起象樓上這樣的誤解。我已經說過,這個方程的通解是x t c1 cos w t c2 sin w t 其中c1,c2為任意常數,而w等於k m的平方根。也可以用三角公式整合...
如何用matlab求解微分方程並畫圖
1 找到關於求解常微分方程的習題。2 這裡用matlab求解,主要用到的就是 dsolve 具體的格式如下,但是不僅僅侷限下面兩種,其中eq就是代表方程式,而con為初始條件。3 這裡的 x 關於x的函式。知道語法就可以進行計算了。4 在matlab中輸入如下,dsolve d2y 2 dy exp...