The VIBES Toolbox has comprehensive plot functionality for the following classes:
vibes.TimeBlocks
vibes.TimeSeries
vibes.FreqBlocks
vibes.FRFMatrix
This tutorial focuses on the plotting capabilities of the vibes.FRFMatrix and vibes.FreqBlocks classes.
All plot settings are available by using the Property, Value syntax. By default, the plot methods of the classes in the above configure the axes according to the units of the Channels and/or RefChannels.
Property
Value
Channels
RefChannels
An object of class vibes.FRFMatrix has lots of predefined plot settings for easy and convenient visualization of the data. This ranges from logarithmic, phase, dB(A), Nth-octave spectra, sum-levels, and many more.
% Load a FRFMatrix object of which FRFs have to be visualized. Y = vibes.load('VIBES,Examples,Datasets','mat','Benchmark','YA.FRFMatrix.mat');
Create a figure and plot a transfer function. By default, an amplitude plot is shown with semilog Y-scale (base 10). One can select specific FRFs by inserting the Channel/RefChannel indices as input arguments, in this example the 1st Channel and 3rd RefChannel, in a frequency range of 0-1000 Hz.
Channel/RefChannel
Channel
RefChannel
vibes.figure('FRF'); Y.plot(1,3,[0 1000]);
Plotting a transfer function
The first two arguments, hence Channel and RefChannel, accept an array of class vibes.Channel as well. This can be especially useful for smart selection of the (Ref-)Channels using the vibes.Channel/find and/or vibes.Channel/select methods on the properties of the Channel objects.
vibes.Channel
(Ref-)Channels
vibes.Channel/find
vibes.Channel/select
Here’s an example of partial and exact channel matching using the logical == in the latter case:
partial
exact
ch_i = Y.Channels.select('Name','Sensor 1'); idx_j = Y.RefChannels.find('Name','==','Impact 1'); Y.plot(ch_i,idx_j,[0 1000]);
Plotting multiple channels of an FRF
Choose plot settings with property and value arguments. The plot method accepts, amongst others, the following properties:
property
value
PlotType
Style
2b-
LineSpec
LineWidth,Color,LineStyle
Label
Smoothing
Y.plot(1,3,[],'PlotType','phase','Style','2b-','Label','FRF-Plot');
Changing plot settings
The plot function also accepts a structure S with properties as its fields:
plot
S = struct(); S.PlotType = 'log'; S.Label = 'Linear smoothing'; S.Style = '2b-'; S.Smoothing = 10; clf Y.plot(1,3,[100 3000], S);
Changing plot settings using a structure
The short-hand options for PlotType are:
log,maglog,loglog,phase,magphase,area,...
db,magdb,logdb,dba,magdba,logdba,... oct,oct3,oct12,...
octdb,oct3db,oct12db,... octdba,oct3dba,oct12dba,...
sum,sumlog,sumdb,suma,sumloga,sumdba;
You can add linear (i.e. moving average) smoothing when Smoothing is a positive scalar:
S.Smoothing = 10;
Logarithmic smoothing when Smoothing is a scalar fraction < 1. This is interpreted as N-th octave smoothing, N being the denominator of the octave fraction, for instance:
S.Smoothing = 1/3; S.Label = '3rd-octave smoothing'; S.Style = '2g--'; hold on Y.plot(1,3,[100 3000],S);
Linear and logarithmic smoothing
Nth-octave plots are available as well; note that these are not energy-normalized (the energy-normalized variant is chosen by setting the option Density to PSD):
Density
PSD
hold on Y.plot(1,3,[100 3000], 'plottype','oct3','Style','2y--','Label','3rd-octave plot'); xlim([100 3000])
Plotting octaves
All of the functionality mentioned above can also be accessed through the vibes.graph.Spectrum class. A vibes.graph.Spectrum object is created automatically for every plot you make. You can access it by right-clicking on a line in your spectrum plot. You can change settings from the context menu such as the plot type, plot style, and smoothing.
vibes.graph.Spectrum
You can also access the spectrum object from your workspace. Edit some properties of the spectrum plot and notice how the figure updates accordingly.
h = Y.plot(); h.Style = '2r--';
Editing properties of the spectrum plot
The same possibilities are available for plotting of data from vibes.FreqBlocks objects. This will be illustrated using the data from a vehicle measurement.
% Load the data FB = vibes.load('VIBES,Examples,Datasets','mat','i3','i3Measurement.FreqBlocks.mat'); close all
When plotting measurement data, it is often useful to define two structures: one for the plot options and one for the plot style.
plotOpts = struct(); plotOpts.PlotType = 'oct'; plotOpts.XLim = [100 2000]; % Define two structures for the styling. st1 = struct(); st1.Style = '1b-'; st2 = struct(); st2.Style = '2k:'; vibes.figure('Spectra'); clf % Plot the acceleration channels using one style FB.plot(1:2,3,plotOpts,st1); hold on % Plot the microphone channel with the same options but a different style FB.plot(3,3,plotOpts,st2);
Plotting measured data
In order to visualize the spectra of all blocks at once, a waterfall plot can be created. This is done through the plot3d method. The [] argument means all blocks should be included in the plot. The plot3d method returns a vibes.graph.WaterfallPlot object, which can be used in a similar way to the vibes.graph.Spectrum objects discussed before.
plot3
[]
plot3d
vibes.graph.WaterfallPlot
vibes.figure('Waterfall diagram'); clf % Create a waterfall plot wf = FB.plot3d(1,[],[100 3000]); % Set limits of the values (YLim) wf.YLim = [60 100];
A waterfall plot
Waterfall plots are not limited to have time on their Z-axis. In fact, any operational data can be used if it is stored in the blocks of the dataset. This will be illustrated by plotting the measured spectra against an estimate of the car’s velocity. This will make a so-called Campbell diagram.
% Load the velocity data from the examples folder Z = vibes.load('VIBES,Examples,Datasets','mat','i3','i3Measurement.ZData.mat'); % Add the operational data to the Z data of the vibes.Block array. FB.Blocks.setZData('Speed (km/h)',Z.v); % Update the waterfall plot, only select the run-up section wf.Columns = {'tc','<',10.5}; wf.ZScale = 'Speed';
A Campbell diagram
For vibes.FreqBlocks objects, an additional class of plots is available: partial plots. These plots help to visualize and compare the contributions of individual datasets to a reference. This is useful for, for example, a Transfer Path Analysis comparison. The use cases will be illustrated with some data from a simulated operational measurement.
% Load the measurement data and FRF of the system. YAB = vibes.load('VIBES,Examples,Datasets','mat','Benchmark','YABm.FRFMatrix.mat'); TS = vibes.load('VIBES,Examples,Datasets','mat','Benchmark','uAB.TimeSeries.mat'); % Convert the simulated data to the frequency domain. TB = TS.toTimeBlocks(); u = TB.toFreqBlocks('hann',YAB.Freq); % Specify the measurement channels and force reference channels. ch_u4 = {'Grouping',1}; ch_f2 = {'Grouping',1}; % Calculate excitation forces for the measured operational condition f2 = YAB.subs(ch_u4,ch_f2) \ u; % Calculate the partial contribution of each force channel. Notice that % this operation returns an array of FreqBlocks objects. u_partials = YAB .* f2; % Select which channels, blocks and frequencies to plot. i = u_partials(1).Channels.select('Label','Sensor 3 +X'); b = {'StartTime','><',[2 3]}; frq = [400 3000];
Show the partial contribution of each force channel by calling the plotPartials method on the entire array of FreqBlocks objects. Request the vibes.graph.PartialPlot object as output.
plotPartials
FreqBlocks
vibes.graph.PartialPlot
vibes.figure('Partial contributions'); clf pp = u_partials.plotPartials(i,b,frq); % Add some white space before the next plot with the 'addSeparator' method. pp.addSeparator(); % Use the 'add' method to add the original measurement as a reference. % Specify 'Reference' as its label. pp.add(u,'Reference')
Plotting partial contributions and reference measurement
Calling the plotPartials method on a single FreqBlocks object will plot the spectrum data per channel. In this example, we will use this functionality to compare the force contributions calculated in the previous section.
% Plot partials for all channels in the dataset by specifying '[]' for the % channels. clf(); pp = f2.plotPartials([],b,frq); % Add a small, 0.1 width separator after the third channel to make a % distinction between forces and moments. pp.addSeparator(0.1,3)
Comparing force contributions
The last class of plots that will be discussed in this tutorial is the bar plot. A bar plot can be created by calling the plotBars method on an array of vibes.FreqBlocks objects, a single vibes.FreqBlocks object or a vibes.FRFMatrix object. It allows you to select one or more metrics (plot types), that are calculated over the specified frequency range and shown in a bar graph. Currently, the supported metrics are: min, mean, max and sum. Note that sum calculates the energetic sum.
plotBars
min
mean
max
sum
% Open a new, full screen figure and create two subplots. vibes.figure('Comparison','Units','Normalized','OuterPosition',[0.25 0.25 0.5 0.5]); ax1 = subplot(1,2,1); % Copy the previously created partial plot to the new axes. pp.duplicate('Parent',ax1); % Create a bar graph in the second subplot. In this case, we will show the % maximum force over the previously defined frequency range. Set the % 'Horizontal' property to true, so that the channel axis is aligned with % that of the partial plot. ax2 = subplot(1,2,2); bp = f2.plotBars([],b,frq,'PlotType','max','Parent',ax2,'Horizontal',true); % Add the same separator as for the partial plot to completely align the % two. bp.addSeparator(0.1,3)
Bar plots
Plot a dashed line crossing the x-axis or y-axis at a specified value using the dashedline function. This function is especially handy when used in combination with the vibes.TimeSeries and the values in the vibes.Block, to visualize the selection of time-data.
dashedline
vibes.Block
% Activate the left subplot axes(ax1); % Plot a dashed line at two points on the frequency axis. dashedline('x', [440 1178],'color','k')
Plotting dashed lines on the frequency axis
Contact us for more VIBES
Contact our support team or call us on +31 85 822 50 49
The supported file type is .txt. Note: The maximum file size is 100kb. If your file is larger than the set limit, please send us an email at support@vibestechnology.com with the file in attachment.
Subscribe to our newsletter
Privacy statement. By submitting this form you agree to having your details registered in our database. At your request, we'll remove your details from our database immediately *