order analysis of the electric motor

I have to do a order analysis for the electric motor. when i am doing the analysis there are lot of errors showing in my code. please help me to cmplete this analysis. my code is given below and i have attached a file.

NFFT = 2000; % then df = 1 Hz

OVERLAP = 0.75;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% load signal

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[rpm, fs]=audioread('09_52_00.wav');

dt = 1/fs;

[samples,channels] = size(rpm);

% time vector

t = (0:samples-1)*dt;

vib = vib - mean(vib); % Remove the DC component

subplot(2,1,1)

plot(t,rpm) % Plot the motor rotational speed

xlabel('Time (s)')

ylabel('Engine Speed (RPM)')

title('Engine speed')

subplot(2,1,2)

plot(t,vib) % Plot the vibration signal

title('Accelerometer Vibration Data')

rpmfreqmap(vib,fs,rpm)

rpmfreqmap(vib,fs,rpm,1)

rpmordermap(vib,fs,rpm,0.005)

[map,mapOrder,mapRPM,mapTime] = rpmordermap(vib,fs,rpm,0.005);

figure

orderspectrum(map,mapOrder)

[spec,specOrder] = orderspectrum(map,mapOrder);

[~,peakOrders] = findpeaks(spec,specOrder,'SortStr','descend','NPeaks',2);

peakOrders = round(peakOrders,3)

ordertrack(map,mapOrder,mapRPM,mapTime,peakOrders)

orderWaveforms = orderwaveform(vib,fs,rpm,peakOrders);

helperPlotOrderWaveforms(t,orderWaveforms,vib)

mainRotorOrder = mainRotorEngineRatio;

tailRotorOrder = tailRotorEngineRatio;

ratioMain = peakOrders/mainRotorOrder

ratioTail = peakOrders/tailRotorOrder

[Mathieu NOE]

hello again ! (seen this wav file already ! )

FYI , the rpm cannot be loaded from your wav file (that contains 3 vibration signals and no RPM)

so this line is wrong :

[rpm, fs]=audioread('09_52_00.wav');

this is supposed to be :

[vib, fs]=audioread('09_52_00.wav');

and from there we don't have the rpm data (where do you record it ? )