How can I plot the integration result?

2 次查看（过去 30 天）
Muhammad Hadyan Utoro
评论： Star Strider 2021-6-21
I used loop to do an integration to get velocity and position, but when I want to plot them, it showed error:
Error using plot
Vectors must be the same length.
Error in Droptest (line 30)
plot(time, V, 'm--')
Here is script I used:
%% Deformation Velocity
for i = 1:length(time)-1;
j = 1:length(xaccel)-1;
V0 = 0;
V = V0 + (xaccel(j+1)+xaccel(j)/2)*(time(i+1)-time(i));
end
I can't figure it out how to make them the same length and using the iterative at the same time.
I appreciate your help

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

采纳的回答

Star Strider
Star Strider 2021-6-21
Ran in:
If you want to plot the cumulative acceleration to get velocity and position, use the cumtrapz function:
t = linspace(0, 10); % Create Data
xaccel = exp(-0.25*t) .* sin(2*pi*t/10); % Create Data
xveloc = cumtrapz(t, xaccel); % Integrate
xposit = cumtrapz(t, xveloc); % Integrate
figure
plot(t, xaccel)
hold on
plot(t, xveloc)
plot(t, xposit)
hold off
grid
legend({'Acceleration','Velocity','Position'}, 'Location','best')
.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by