What changes are required in the code

3 次查看（过去 30 天）

显示更早的评论

Manav Divekar 2021-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1567478-what-changes-are-required-in-the-code

关闭： Cris LaPierre 2021-10-20

G2_adult.txt

在 MATLAB Online 中打开

clc
clear all
close all
C = readtable('G2_adult.txt', 'NumHeaderLines',8);
%Convert to SI Units
accG2test = C{:,1}.*9.81;
forceG2test = C{:,2}.*0.00981;
dispG2test = C{:,3}./1000;
% time = dispG2test./accG2test;
dispG2test2 = smoothdata(dispG2test, 'sgolay', 100);
% Valleys Of Smoothed Data
[Vlys,vlocs] = findpeaks(-dispG2test2, 'MinPeakProminence',0.5);   
Vlys = [-dispG2test2(1); Vlys];
vlocs = [1;vlocs];
% Sampling Interval
Ts = 1;                                                             
t = 1:height(C);
figure
plot(t, dispG2test2)
hold on
plot(t(vlocs), -Vlys, '+r')
hold off
for k = 1:numel(vlocs)-1
	% Index Range For This Segment
    idx{k} = vlocs(k):vlocs(k+1);                                   
		% Time Vector
    tc{k} = t(idx{k});                                              
		% Signal Segment
    sc{k} = forceG2test(idx{k});                                    
end
 % Maximum Segment Length
tcmax = max(cellfun(@numel,tc));                                   
% Preallocate
ensb = zeros(tcmax,numel(idx));                                     
for k = 1:numel(idx)
	% Create 'Ensemble'
    ensb(1:numel(idx{k}),k) = sc{k};                                
end
figure
 % Plot Ensemble
plot(1:tcmax, ensb) 
title ('C:Diflection Vs Time (Cycles)');
xlabel('Time[sec]');
ylabel('Diflection[m]');