MATLAB错误：产品：不合格的参数（OP1是131072x1，OP2 IS 128000x1）

Question

我正在尝试在MATLAB代码中乘以两个向量。但是在这一行中，

fHatNoisy  = indices.*fHat 

它说了错误，

 error: product: non-conformant arguments (op1 is 131072x1, op2 is 128000x1)

尽管我尝试添加操作员'。我仍然会遇到错误..

知道为什么会发生此错误？

代码总数：

[samples,samplingRate] = audioread('car.wav')
dt=1/samplingRate
sampleLength=length(samples)
duration = sampleLength/samplingRate
time = 0:dt:duration #time vector
n=length(time)

%disp("sampling Rate"+Fs)

fHat = fft(samples, n)
%psd = periodogram (fHat)
PSD = periodogram (fHat)
indices = PSD<0.01 #Find all freqs with large power
PSDClean1 = PSD.*indices #find out all others(While there are 0s, PSDs which are there will be 0 from indices vector)
fHatNoisy  = fHat.*indices #Zero out small Fourier coeffs. in Y

更新：这是因为两个阵列，FHAT和索引不相等，这就是我怀疑发生错误的原因。试图弄清楚为什么FFT（）返回与PSD = fhat * conj（fhat）/n

更多有关PSD，N，其长度等的信息：

[答案]更新2，计算DFT/FFT的PSD：

pkg load signal;

[samples,samplingRate] = audioread('car.wav')

dt=1/samplingRate

sampleLength=length(samples)

n = length(samples);          % number of samples
frequencyRange = (0:n-1)*(samplingRate/n);     % frequency range
fHat = fft(samples)
power = abs(fHat).^2/n;    % power of the DFT

plot(frequencyRange,power)
xlabel('Frequency')
ylabel('Power')