1樓：愚人談娛樂
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)的相平面圖');
2樓：匿名使用者
>> clear
>> syms a b c d e;
>> y=dsolve('a*d2y+b*dy+c*y=0','y(0)=d','dy(0)=e')
y =(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2))
>> a=1;b=1;c=1;d=1;e=1;%若常數已知
>> t=1;%desolve中沒指定自變數x，這裡預設為t，帶入任意一個x值
>> y=eval(y)
y =1.1932 + 0.0000i
這是解方程的方法，畫圖也差不多，我在命令窗寫的改麻煩，下面是我在editor寫的。
syms a b c d e;
y=dsolve('a*d2y+b*dy+c*y=0','y(0)=d','dy(0)=e','x');
x=200:0.1:400;%取值範圍，步長
a=1;b=1;c=1;d=1;e=1;%帶入引數
y=eval(y);
plot(x,y)
你補充問題的方程解得也是一個空的，並報錯，或許就是無解的，我也很疑問。
3樓：匿名使用者
這麼簡單的常微分方程問題高數應該有教怎麼解吧，去看看高數吧，常微分方程部分有關λ的部分，又是其次方程，兩步就出來了。
第一次補充回答也是最後一次補充回答：
開始以為abcde是常數，沒想到不是。
dsolve後面的函式應該是一次的微分方程，你應該將函式標準化處理。令y'=y1，帶到原方程中去形成兩個一次的微分方程，注意初始條件也要變一下y'(0)=e變成y1（0）=e，然後再用dsolve函式求解，係數為非常數，可能這個函式解不了，可能需要手動程式設計，建議你去找本常微分的書看看，高數上可能沒有講過這種方程的求解。剛翻了下書，發現這問題還是比較煩的，你自己去研究吧 o(∩_∩)o~
另外用一個函式多看看幫助，matlab幫助寫的很好的
4樓：善冰夏侯
doslove函式可以解決
5樓：匿名使用者
用dsolve函式來解.
有關matlab解二階微分方程
6樓：匿名使用者
你的問題錯就錯在，y(1)=x; y(2)=dxdt; 這兩句上，而且是多餘的。去掉後，用ode45（）求解可以執行得到其數值解。
用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)的相平面圖');
matlab 求解二階微分方程並畫出y與t函式曲線
8樓：匿名使用者
用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)
9樓：匿名使用者
你好，求解這個方程的一般步驟用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);
二階微分方程求通解，高等數學，二階微分方程，求通解，需要詳細步驟，謝謝
求微分方程 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 ...
求二階微分方程的通解，高等數學，二階微分方程，求通解，需要詳細步驟，謝謝
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 北極灬寒冰 可以...
可降階的二階微分方程和二階常係數線性微分方程的區別
倔強的水蘿蔔 可降階的就是把y 換成y，算出y後再積分！實際上就是一階的！ 不要說話 可降階的二階微分方程 1,y f x 型的微分方程 此類方程特點是 方程右端僅含有自變數x，只需積分兩次便可得到方程的通解。2,y f x,y 型的微分方程 此類方程特點是 方程右端不顯含未知函式y。作變數代換y ...