Why my calculation result is not matrix result？

The result of my matla calculation is a numerical value instead of a matrix result, so I can't draw the relevant graph(figure3)：

——The output of AweightedSPL and SPL is a numerical rather than a matrix result,So I can't plot the "frequency-AweightedSPL "spectrum. Input audio has been uploaded

can someone help me see where the code can go wrong?How should I change it?

thanks~

[x,fs] = audioread('123.wav'); 

L = length(x);
y=fft(x); 
f  =fs/2.*linspace(0,1,L/2+1);
y = y(1:L/2+1);
plot(f,abs(y),'b'); 

A= filterA(f);
yweighted = A'.*y; 
figure(1);
plot(f,abs(yweighted),'r');

xweighted = ifft(yweighted);
xweighted = real(xweighted); 

tempdata = xweighted;
tempL = length(tempdata);
temp1 = fft(tempdata)/tempL*2;
figure(2);
plot(f,abs(temp1),'g');


linelevel =1.8;
preampout = abs(xweighted.*1.8);
preampin = preampout;
pressure = 1000.*preampin./22.78;
pref = 0.00002;
pressure = sqrt(mean(pressure.^2));
AweightedSPL  =20.*log10(pressure/pref);
figure(3);
plot(f,AweightedSPL,'b');grid

[Matt J]

We can't run your code without the input file.

[Stephen23]

"can someone help me see where the code can go wrong?"

https://www.mathworks.com/help/matlab/matlab_prog/array-vs-matrix-operations.html

[ZHU z]

thank you @Stephen

I'm going to study your link

[Jeffrey Clark]

These three lines of yours produce a scaler (except in the case of a matrix, mean returns a scaler), while f is a vector, so the last plot will fail. Perhaps you meant to use the mean to normalize each element of the [initial] vector pressure?

pref = 0.00002;
pressure = sqrt(mean(pressure.^2));
AweightedSPL  =20.*log10(pressure/pref);

[Jan]

What is filterA?

A= filterA(f);

[ZHU z]

thank you @Jan

This is a function of A frequency weighting~

[ZHU z]

thank you @Matt J

Sorry for the long wait, the input file has been uploaded~

[ZHU z]

thank you, @Jeffrey Clark

I will go check it