1樓:匿名使用者
該方程為2階微分方程,可以化為兩個1階方程構成的方程組來求解。
方法為:設回x1=y,x2=y'
先建立myode.m檔案用來描述該系統
%%% 以下答為myode.m的內容
function dx=myode(t,x)dx=zeros(2,1);
dx=[x(2);
(5120-4*x(1)-0.05*(x(2)^2-0.2*x(2)))/160;];
%%%%%%%%%%%%%%%%%%%%%%%%%%以下為畫出y關於t的影象的測試程式test.m的內容tspan = [0 5]; % solve from t=1 to t=5
ic = [-160,0]; %
[t x] = ode23('myode',tspan,ic); % solve ode
plot(t,x(:,1))
xlabel('t')
ylabel('y(t)')
以下為執行結果:
matlab微分方程求解 用ode23函式
2樓:不好好改就革
請給初始值,o(∩_∩)o謝謝!!
比如我給初始值:初始值:
θ(0)=0.5; θ'(0)=0;
下面用y代表θ,然後,函式程式:
function yprime=verderpol(t,y)global g;
global l;
g=9.8;
l=10;
yprime=[y(2);-g/l*sin(y(1))];
end主程式:
clc; clear all
global g;
global l;
g=9.8;
l=10;
y0=[0.5;0];%初始值:y(0)=0.5;y'(0)=0;
hold on
[t,y]=ode23('verderpol',[0,3*pi],y0);
plot(t,y(:,1),'r','linewidth',2);
plot(t,y(:,2),'g','linewidth',2);
legend('θ-t圖','θ''-t圖')grid on
hold on
畫出的影象:
希望我的回答對你有所幫助,滿意敬請採納喲,o(∩_∩)o謝謝!
3樓:匿名使用者
給的條件不夠,只好自己假設了
l = 1m, g=9.8 θ0=π/12%存為fun.m
function xdot=fun(t,theta)g=9.8; l=1;
xdot=[theta(2);-g*sin(theta(1))/l];
%存為main.m
theta0=[pi/12,0];
[t, theta]=ode23('fun',[0,10],theta0);
plot(t,theta(:,1));
4樓:匿名使用者
g=10;
l=1;
fun=@(t,x) [x(2),;-g/l*sin(x(1))];
ts=[0,20];
x0=[pi/20,0];
[t,x]=ode23(fun,ts,x0);
t=t';
theta=x(:,1)';
plot(t,theta);
xlabel('t');
ylabel('\theta');
請用euler和matlab函式ode23求解下列微分方程
5樓:宇宙無敵小陽
1--尤拉法
t=0;
h=0.02;
y(1)=10;
for i=2:100000000
y(i)=y(i-1)+0.02*h*(1-0.001*y(i-1))*y(i-1);
if i*h>=20
break;
endend
2-----ode23
先建立function dy=dfun(t,y)dy=0.02*(1-0.001*y)*y;
再編寫[t,y]=ode23(@dfun,[0 120],10);
plot(t,y)
微積分 求下列微分方程的通解,求微分方程通解,要詳細步驟
a dy dx 2xy 0 dy dx 2xy dy y 2x dx ln y x 2 c y c.e x 2 b dy dx xy 2x dy dx x y 2 dy y 2 xdx ln y 2 1 2 x 2 c y 2 ce 1 2 x 2 y 2 ce 1 2 x 2 a dy dx 2x...
高數,求微分方程。要過程,求大佬解下啊
r 1 2 1 r1 1 i,r2 1 i 齊次方程的通解為 y e x c1 cosx c2 sinx 因為非齊次項為4e x cosx 所以設非齊次方程的特解為 y xe x mcosx nsinx y e x mcosx nsinx xe x mcosx nsinx xe x ncosx ms...
怎麼求下列微分方程滿足所給初始條件的特解
1 dy dx 2 2x 2 y 2 ydy 2 2x dx 兩邊積分 2 y ln2 2 2x ln2 1 2 c2 y 2 2x 1 c 令x 0 1 1 2 c,c 1 2 所以2 y 2 2x 1 1 2 2 y 1 2 2x 1 2 y ytanx secx 因為 ye f x e f x...