ode45 or dsolve

6 次查看（过去 30 天）
felik kamwera
felik kamwera 2019-6-24
评论： felik kamwera 2019-6-24
whats wrong with my code.. its from a book but cant run
tspan=[0 4];
y0=[0.02;0];
[t,y]=ode45('unforced1',tspan,y0);
plot(t,y(:,1));
grid on
xlabel(‘time’)
ylabel(‘Displacement’)
title(‘Displacement Vs Time’)
hold on;
its function is function yp = unforced1(t,y)
yp = [y(2);(-((c/m)*y(2))-((k/m)*y(1)))];
  4 个评论
显示 2更早的评论
Chinmay Anand
Chinmay Anand 2019-6-24
What are c and m defined in your function ? Have you assigned any values to them ?
felik kamwera
felik kamwera 2019-6-24
should i give them values in the function definition or call?

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

madhan ravi
madhan ravi 2019-6-24
tspan=[0 4];
y0=[0.02;0];
m=5;
k=1000;
c=1000;
[t,y]=ode45(@(t,y)unforced1(t,y,c,m,k),tspan,y0); % function call
plot(t,y(:,1));
grid on
xlabel('time')
ylabel('Displacement')
title('Displacement Vs Time')
%function definition
function yp = unforced1(t,y,c,m,k) % if your version is prior to 2016b save this in a separate file name unforced1.m
yp = [y(2);
(-((c/m)*y(2))-((k/m)*y(1)))];
end
  1 个评论
felik kamwera
felik kamwera 2019-6-24
thank u very much...i have seen where i was wrong..again thank u

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by