Plot a second order equation and plot two equations on the same grafic

Question

Igor 2012-2-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/28444-plot-a-second-order-equation-and-plot-two-equations-on-the-same-grafic

编辑： Paul 2013-10-8

采纳的回答： Kevin Holst

Hi everyone,

1. I can't plot this function in a grafic, I think the code in the editor is not right. The second order equation is:

(a) dx^2/dt^2 (t) = -xg0*w0^2*sin(wf*t)-w0^2*x-2*M*w0*beta*dx/dt

xg0, w0, M, beta are constants.

2. Plot two solutions in the same grafic.

I already solve other system with ode45, now i want to plot that second order equation with one solution of this sistem.

It is a enginneering problem that consist in a structure with and without a vibration control system. One equation (a) give me the response without that control and the system of equations give me the response with that control... I need to compare that in the same grafic.

Thank's everyone.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Kevin Holst 2012-2-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/28444-plot-a-second-order-equation-and-plot-two-equations-on-the-same-grafic#answer_36746

在 MATLAB Online 中打开

You can plot multiple sets of data on one figure multiple ways. The first way would be:

plot(x1,y1,x2,y2,...)

The next most common would involve using hold:

plot(x1,y1)
hold on
plot(x2,y2)

11 个评论
显示 9更早的评论隐藏 9更早的评论

Igor 2012-2-8

I can plot all the solutions off the system s, funtion of time (t) with the comand plot(T,Y(:,1)) (for example). What I can't do is define the equation:

dx^2/dt^2 (t) = -xg0*w0^2*sin(wf*t)-w0^2*x-2*M*w0*beta*dx/dt

on the EDITOR and than plot it. After define that funtion, in the end I want to plot the two functions in the same grafic.

I send you my EDITOR:

M=0.577;

mv=1000;

A=1.157*10^-5;

xg0=0.05;

w0=5.17;

wf=5.3;

beta=2/100;

alfa=0.7;

L=0.734;

g=9.81;

delta=0.5;

%Define Equation1 -I CAN'T DEFINE THIS

dx^2/dt^2 (t) = -xg0*w0^2*sin(wf*t)-w0^2*x-2*M*w0*beta*dx/dt ???

%Define System - OK!

s=@(t,y) [y(3); y(4); (-(M+mv*A*L)*xg0*w0^2*sin(wf*t)-M*w0^2*y(1)-2*M*w0*beta*y(3)-mv*A*alfa*L*((mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)))/(M+mv*A*L);(mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)];

[T,Y]= ode45(s,[0 30],[0 0 0 0]);

%PLOT TWO SOLUTIONS - I CAN'T DO THIS

plot (T,Y(:,1), dx^2/dt^2 (t)) ?????

Igor 2012-2-8

It doesn't work =\. I think I explain something not right.

1. System of equation

s=@(t,y) [y(3); y(4); (-(M+mv*A*L)*xg0*w0^2*sin(wf*t)-M*w0^2*y(1)-2*M*w0*beta*y(3)-mv*A*alfa*L*((mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)))/(M+mv*A*L);(mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)];

[T,Y]= ode45(s,[0 30],[0 0 0 0]);

- This is ok, I can plot the 4 solutions of the system with plot(T,Y(:,1)), solution 1 for example.

2. I define that second order equation and solve with ode45

r=@(t,y) [y(3); y(4); -xg0*w0^2*sin(wf*t)-w0^2*y(1)-2*w0*beta*y(3); 0];

[T,Y]= ode45(r,[0 30],[0 0 0 0]);

3. If the 2. is right, now I want to plot two solution on the same grafic. Isolate I can plot with the command plot(T,Y(:,1)).. But I don't now the command to plot the two. Already try:

plot(T,Y1(:,1),T,Y2(:,1)) , and other things that doesn't work.

I really appreciate the help you given to me, really ;)

Kevin Holst 2012-2-8

Also, don't forget to click accept answer under my name up there.

Igor 2012-2-10

I have some other problem with my dissertation, sorry stop answering. I am a finalist of civil engineering.

All you say it works, thank's to you I have the program correct. Thank you.

I have another question. Now I have two functions in the same graphic. I need to calculate each maximum to calculate the effectiveness of the control system. What is the code for each function?

I send to a reduce EDITOR:

M=0.577;

mv=1000;

A=1.157*10^-5;

xg0=0.05;

w0=5.17;

wf=5.17;

beta=2/100;

alfa=0.7;

L=0.734;

g=9.81;

delta=0.5;

%Equation 1

r=@(t,y) [y(3); y(4); -xg0*w0^2*sin(wf*t)-w0^2*y(1)-2*w0*beta*y(3); 0];

[Tr,Yr]= ode45(r,[0 50],[0 0 0 0]);

%System of equations

s=@(t,y) [y(3); y(4); (-(M+mv*A*L)*xg0*w0^2*sin(wf*t)-M*w0^2*y(1)-2*M*w0*beta*y(3)-mv*A*alfa*L*((mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)))/(M+mv*A*L);(mv*A*alfa*L*((M+mv*A*L)*xg0*w0^2*sin(wf*t)+M*w0^2*y(1)+2*M*w0*beta*y(3))+(M+mv*A*L)*(-mv*A*alfa*L*xg0*w0^2*sin(wf*t)-2*mv*A*g*y(2)-(1/2)*mv*A*delta*abs(y(4))*y(4)))/(mv*A*L*(M+mv*A*L)-(mv*A*L)^2)];

[Ts,Ys]= ode45(s,[0 50],[0 0 0 0]);

%Graphic

plot(Ts,Ys(:,1),Tr,Yr(:,1))

% Effectiveness of the system ????

Tank's.

请先登录，再进行评论。

Plot a second order equation and plot two equations on the same grafic

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

11 个评论
显示 9更早的评论隐藏 9更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Plot a second order equation and plot two equations on the same grafic

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

11 个评论 显示 9更早的评论隐藏 9更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

11 个评论
显示 9更早的评论隐藏 9更早的评论