t=[0 : 1/Fs : L/Fs - 1/Fs].';
d=(1/2048)*9.81*data1(:,ax);
doff=repmat(-mean(d),[L 1]);
FigHandle = figure('Position', [100, 100, 1280, 720]);
title(['Time domain of recorded vibration, axis = ' num2str(ax)],'FontSize',12)
xlabel('Time (s)','FontSize',14);
ylabel('Acceleration (m/s^2)','FontSize',14);
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f,2*abs(Y(1:NFFT/2+1)));
title(['FFT of recorded vibration, axis = ' num2str(ax)],'FontSize',12)
xlabel('Frequency (Hz)','FontSize',14)
ylabel('Acceleration (m/s^2)','FontSize',14)
w = window(@chebwin,nfft);
[B,f,t] = spectrogram(x,w,[],nfft,fs);
contour3(t,f,abs(B),128);
xlabel('Time (s)','FontSize',14);
ylabel('Frequency (Hz)','FontSize',14);
zlabel('Spectrogram amplitude','FontSize',14);
title(['STFT, Chebyshev Window size = ' num2str(nfft) ', axis = ' num2str(ax)],'FontSize',12);
xlabel('Time (s)','FontSize',14);
ylabel('Frequency (Hz)','FontSize',14);
axis([0 max(t) 0 0.5*max(f)]);
title(['STFT, Chebyshev Window size = ' num2str(nfft) ', axis = ' num2str(ax)],'FontSize',12);
Error using pwelch
Expected x to be finite.
Error in signal.internal.spectral.welchparse>parse_inputs (line 111)
validateattributes(x2,{'single','double'}, {'finite','nonnan'},'pwelch','x')
Error in signal.internal.spectral.welchparse (line 31)
parse_inputs(x1,esttype,varargin{:});
Error in pspectrogram (line 30)
[xw,nx,~,yw,ny,win,~,~,noverlap,~,~,options] = signal.internal.spectral.welchparse(x,esttype,inpArgs{:});
Error in spectrogram (line 191)
[varargout{1:nargout}] = pspectrogram({x},'spect',inpArgs{:});
Error in shorttimefouriertransform (line 97)
[B,f,t] = spectrogram(x,w,[],nfft,fs);