A-weighting Filter with Matlab
- 53 files
- 979 downloads
The present code is a Matlab function that provides an A-weighting of a given signal via spectral processing.
An example is given in order to clarify the usage of the function. For convenience, the input and output arguments are given in the beginning of the function.
The code is based on the theory described in:
[1] Design Response of Weighting Networks for Acoustical Measurements ANSI S1.42-2001. Washington, American National Standards Institute, 2001.
[2] Electroacoustics Sound Level Meters Part 1: Specifications IEC 61672-1:2002. Geneva, International Electrotechnical Commission, 2002.
Cite As
Hristo Zhivomirov (2019). A-weighting Filter with Matlab (https://www.mathworks.com/matlabcentral/fileexchange/46819-a-weighting-filter-with-matlab), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (8)
Updates
| 1.5.0.0 | Change of the title of the submission. |
|
| 1.4.0.0 | A new version of the code has been uploaded. |
|
| 1.3.0.0 | A new version of the code has been uploaded. |
|
| 1.2.0.0 | A new version of the code has been uploaded. |
|
| 1.1.0.0 | Typos correction has been made. |
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Hristo Zhivomirov (view profile)
Hi, Ifigenia!
In order to resolve the issue you can filter the signal (time series) x[n] as y = filterA(x, fs); and then you can find the dB(A) values as Y = 20*log10(y/2e-5);. If you need some additional help, do not hesitate to contact me.
All best,
Hristo
Ifigenia Aslanidou (view profile)
Hello!Can anyone help me calculating the SPL in Db(A)?
I have the Sound Pressure signal in Pa and I need it dB(a).
Using this function, and using for x the Sound pressure Level I computed , I get a result of NaN.
I appreciate any help :)
Hristo Zhivomirov (view profile)
Hi Jonas! Thank you for your comment! Everything is alright with the function now! You may test and share your opinion about it!
Best regards,
Hristo
Jonas Stehli (view profile)
Hi. I think this implementation is wrong. When plotting 20*log10(A) over frequency, this does not at all match the blue curve in https://upload.wikimedia.org/wikipedia/en/c/c2/Lindos3.svg.
I think you are off by a square, i.e., a factor 2 in logarithmic domain. If you'd take the square root of the weighting curve before multiplying with the signal spectrum, you'd be more or less correct.
Tom Ransegnola (view profile)
Very helpful! Thank you!
micpro (view profile)
Hi,
Does anyone help me related with dbA filter please?
Hristo Zhivomirov (view profile)
Hi Steve! I can't resolve the problem at this stage, since I haven't enough information. Please, send me the figure, it will be a good starting point.
Steve (view profile)
I made a comparison of the resulting spectrum using this time domain filter and the resulting spectrum using a direct implementation of the frequency domain definition of A-weighting as per the ANSI Standard S1.42-2001. Agreement was less than satisfactory.
My sampling rate was 25.6 kHz. Maybe this is too low. Let me know if you have any suggestions. I can send a figure showing the comparison if you'd like.