用matlab解二阶微分方程
2个回答
展开全部
function Myexamp1
clear,clc,
tspan = [0 20]; % t 的取值范围
y0=[0;1]; % 初始值 y(t0) 和 y'(t0)
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[t,y] = ode45(@fun,tspan,y0 ,options);
plot(t,y(:,1),'-',t,y(:,2),'-.')
title('y''(t)= -3 cos(2t) +2sin(t)+ t-3.8')
function dy = fun(t,y)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = -3-cos(2*t) + 2*sin(t)+t-3.8;
clear,clc,
tspan = [0 20]; % t 的取值范围
y0=[0;1]; % 初始值 y(t0) 和 y'(t0)
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[t,y] = ode45(@fun,tspan,y0 ,options);
plot(t,y(:,1),'-',t,y(:,2),'-.')
title('y''(t)= -3 cos(2t) +2sin(t)+ t-3.8')
function dy = fun(t,y)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = -3-cos(2*t) + 2*sin(t)+t-3.8;
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
