二元二階微分方程組求數值解的matlab程式

時間 2021-08-14 05:14:58

1樓:花開勿敗的雨季

方程或者初始條件有問題,請認真核實。

原因:如果按照現在的初始條件,則第一個微分方程各項為0,方程不起作用。而且就演算法而言,二階(最高階)導數的係數為0,會導致數值積分演算法無法啟動。

參考**大致如下:12

3456

78910

1112

1314

1516

17m=1;

g=9.8;

k=100;

l0=1.1;

x0=[0.1 0 0 0];

% 定義x1=l, x2=l', x3=a, x4=a'

dx=#(t,x)[x(2); (m*x(1)*x(4)^2-m*g*cos(x(3))-k*(x(1)-l0))/m; ...

x(4); (-2*m*x(1)*x(2)*x(4)+m*g*x(1)*sin(x(3)))/(m*x(1)^2)];

[t,x]=ode15s(dx,[0 10],x0);

lstr = ''', '\alpha', '\alpha''' };

for i=1:length(lstr)

subplot(2,2,i)

plot(t, x(:,i));

xlabel('time')

ylabel( lstr )end

2樓:匿名使用者

function dz=yueh(t,z) %儲存為 yueh.m 檔案

dz=[z(2);

10-z(1)-0.001*z(3);

z(4);

-0.0001*z(1)-0.001*z(3)];

執行一下:

[t,z]=ode45('yueh',[0 1],[0 0 0 0]);

% 感覺初始值不完整啊,這裡假設x(0)=θ(0)=x'(0)=θ'(0)=0

plot(z(:,1),z(:,3)) % 這裡是 x-θ 影象

3樓:光之痕跡

微分方程初值問題,用ode45可以解

龍格庫塔求解二階微分方程組的matlab程式設計

4樓:匿名使用者

matlab求解x''+0.7x'+0.8x'|x'|+25.6x-25.6x³=0二階微分方程組的方法,可以按下列步驟進行:

1、建立自定義函式func()

function f = func(t,x)

%x''+0.7x'+0.8x'|x'|+25.6x-25.6x³=0

f(1)=x(2);

f(2)=25.6*x(1)^3-25.6*x(1)-0.8*x(2)*abs(x(2))-0.7*x(2);

f=f(:);

2、建立龍格庫塔演算法函式runge_kutta()

呼叫格式:[t,x] = runge_kutta(@(t,x)func(t,x),x0,h,a,b);

3、然後根據x和x'資料,繪製出x(t)、x′(t)的圖形。

plot(x(:,1),x(:,2))

怎麼用matlab求解二階微分方程組

如何用matlab求解二階微分方程,以及程式例項

5樓:用著追她

1、首先看一下matlab求解方程

的方法,首先指明所解方程的變數,然後指明方專程,未知數和限制條件,屬最後求解方程。

2、在matlab命令列視窗中輸入syms x [x,params,conds]=solve(sin(x)==1,'returnconditions', true) ,按回車鍵可以得到方程解。

3、轉換一下,可以看到sin(x)=1方程的解是如下圖所示 。

4、也可以求解下面的一個方程。

5、輸入syms a b c y x;[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])。

6、按回車鍵可以得到方程解。

6樓:匿名使用者

^1、對

copy於求數值解的微分方程,你可以用ode45()函式求解。如求下列微分方程

func。m   %自定義微分方程的函式

function z = func(t,y)z =[y(2);(1-y(1)^2)*y(2)-y(1)];

main。m  %主程式

clear all;close all;clcy0 = [0.25;0];

h = 0.1;

a = 0;

b = 20;

[t1 y1] = ode45(@func,y0,h,a,b)2、對於求解析解的微分方程,你可以用dsolve()函式求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解,可以下列步驟計算

>> syms y(x)

>>dy = diff(y);d2y = diff(y, 2);

>>dsolve(x*d2y+x*(dy)^2-dy==0,'x')

7樓:閻輝門曼雲

^dsolve('d2y=x','x')

ans=

x^3/6

+c2*x+c3

d2y的意思就是復y的二階微分項制

不明白你的問題什麼意思,要輸入的話直接定義符號變數輸入syms

d2xx

d2x=x;

如果是矩陣,那就用矩陣表示

求二階微分方程的通解,高等數學,二階微分方程,求通解,需要詳細步驟,謝謝

2y y y 3e x,先求齊次方程通解。令2t 2 t 1 0,解得t 1或1 2即齊次解為y a e x b e 1 2x 其中a,b r 再求1個特解即可。令y c e x,則2c c c 3,即c 3 2故問題的解為3 2 e x a e x b e x 2 其中a,b r 北極灬寒冰 可以...

關於微分方程已知某二階線性微分方程三解為X,X 1,eX X求其通解及原方程

通解寫成這樣的形式是線性方程解的結構決定的。高等數學 同濟大學第五版p296 若原方程為y p x y q x y f x 則對應的齊次方程y p x y q x y 0 y1 x,y2 x 1,y3 ex x為原方程的特解,則可知y1 y3,y2 y3是對應齊次方程的兩個線性無關的特解,於是 的通...

可降階的二階微分方程和二階常係數線性微分方程的區別

倔強的水蘿蔔 可降階的就是把y 換成y,算出y後再積分!實際上就是一階的! 不要說話 可降階的二階微分方程 1,y f x 型的微分方程 此類方程特點是 方程右端僅含有自變數x,只需積分兩次便可得到方程的通解。2,y f x,y 型的微分方程 此類方程特點是 方程右端不顯含未知函式y。作變數代換y ...