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...