1樓：二中混過
^^令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求解二階微分方程組
如何用matlab求解二階微分方程，以及程式例項
2樓：用著追她
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、按回車鍵可以得到方程解。
3樓：匿名使用者
^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')
4樓：閻輝門曼雲
^dsolve('d2y=x','x')
ans=
x^3/6
+c2*x+c3
d2y的意思就是復y的二階微分項制
不明白你的問題什麼意思，要輸入的話直接定義符號變數輸入syms
d2xx
d2x=x;
如果是矩陣，那就用矩陣表示
用matlab求解二階微分方程 100
5樓：匿名使用者
用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中畫出影象，看系統辨識出的模型，和自己原來
syms s gs sym 5.6068 1 4.5965 s 1 4.0785 s 1 0.035903 s u 360,365,370,375,380,390,400,405,410 s u ds subs gs plot s,ds,g 上面是第一個問題，下面是第二個 試求下列函式的部分分式式 ...
matlab中如何畫出冪函式指數為分數時比如yx
當x從正方向趨 近於bai0時 即0.0000000.y趨近於正du無窮大 zhi dao當x從負方向趨近於0時 內即 0.0000000.y趨近於負無窮容大 當x從趨近於正無窮大時 即9999999.y趨近於0 當x從趨近於負無窮大時 即 9999999.y趨近於0 類似於反比例函式y 1 x,過...
matlab如何在圓形影象中畫出經過圓心的任意剖面圖
把這些點的所有座標存入陣列xy，假設圓心為 x0,y0 半徑r,節點數為版n，xy的第一列為x座標，權第二列為y座標。如下 x0 1 y0 1 r 2 n 100 xy zeros 100,2 h 2 pi n 1 for i 1 n xy i,1 x0 r cos i 1 h xy i,2 y0 ...