matlab求解二階微分方程並畫出y與t函式曲線

時間 2021-05-05 14:36:46

1樓:匿名使用者

用matlab 求解二階微分方程並畫出y與t函式曲線,可以用dsolve函式和plot函式來實現。

給出的二階微分方程存在著幾個疑點:1、表示式中的z是否是筆誤,還是y。如是z,其表示式是什麼?2、求解二階微分方程的數值解,必須有兩個初始條件,本題只有乙個。

下列**中,z按y來處理,增加初始條件y'(0)=0

>> m=0.00267;k=0.5;b=0.000001;a=0.000004;u=90;p=-19.5216;a=0.001256;

>> syms y(t)

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

>> y=dsolve(m*d2y+(k+1/b)*y==a/b*u-p*a,y(0)==0,dy(0)==0);

>> t=0:0.1:20;

>> y=eval(y);

>> plot(t,y)

2樓:匿名使用者

你好,求解這個方程的一般步驟用ode45。

**如下

clc;clear

m = 0.00267;    alpha = 4e-6;

beta = 1e-6;    u = 90;

k = 0.5;        p = 19.5216;

a = 0.001256;

dy=@(t, x, m, alpha, beta, u, k, p, a)[x(2);...

m^(-1)*(-(k+1/beta)*x(1)+alpha*u/beta-p*a)];

[t,x] = ode45(@(t,x)dy(t, x, m, alpha, beta, u, k, p, a),[0,1e-3],[0,0]);

subplot(1,2,1)

plot(t,x(:,1))

xlabel('$$t$$','interpreter','latex','fontsize',16);

ylabel('$$y$$','interpreter','latex','fontsize',16)

subplot(1,2,2)

plot(t,x(:,2))

xlabel('$$t$$','interpreter','latex','fontsize',16);

ylabel('$$y''$$','interpreter','latex','fontsize',16);

請問在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求解二階微分方程,以及程式例項

4樓:用著追她

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、按回車鍵可以得到方程解。

5樓:匿名使用者

^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')

6樓:閻輝門曼雲

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

ans=

x^3/6

+c2*x+c3

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

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

d2xx

d2x=x;

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

用matlab求解二階微分方程 100

7樓:匿名使用者

用matlab的dsolve函式計算該二bai階微分方程的解析du值是有一定的zhi難度。只能用

dao數值分析的方法來求版解其數權值解。當然求其數值解,必須還得已知a、θ的具體數值。該問題可以用matlab的ode45函式求得。求解方法如下:

1、建立微分方程的自定義函式,odefun(x,y)2、當a=10、θ=π/6時,執行下列**theta=pi/6;

[x,y]= ode45(@odefun,[0,0.18],[0.5,theta])

plot(x,y),grid on

legend('y (x)','y』(x)')xlabel('x'),ylabel('y (x),y』(x)')figure(2)

plot(y(:,1),y(:,2)),grid onxlabel('y (x)'),ylabel('y』(x)')title('y (x)—y』(x)的相平面圖');

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

求微分方程 y 2y y 5e x 的通解 解 齊次方程 y 2y y 0的特徵方程 r 2r 1 r 1 0的根r r 1 因此齊次方程的 通解為 y e x c c x 因為原方程右邊的5e x 中的指數所含 1正好是特徵方程的重根,因此要設特解為 y ax e x y 2axe x ax e ...

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

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 北極灬寒冰 可以...

二階線性齊次微分方程通解求法

朋秀愛薩棋 解求特徵方程r 2 p x r q x 0解出兩個特徵根r1,r2 若r1 r2且r1,r2為實數,則y c1 e r1 x c2 e r2 x 若r1 r2且r1,r2為實數,則y c1 xc2 e r1 x 若r1,r2即a bi為複數,則y e ax c1 cosbx c2 sin...