Math 447 Project 5

Project B. Audio Compression

Reference: T. Sauer, 2012, Numerical Analysis | Click here for more

Jiajing Guan



Main Goal:

To explore the process of audion compression.



Methods/Equations:

Reality Check 11
1) Create modified discrete cosine transformation (MDCT) matrix
2) Quantize and dequantize
3) (Optional) Apply windowing function
4) Apply inverse of MDCT



Question 1

Steps:
1) Set bits to be 4 and window size to be 32.
2) Input signal of frequency \(64Hz\).
3) Compare the original and decoded signal and the difference (RMSE) between them.
4) Try different even and odd frequencies and compare the result.

Matlab code: here

Result:

\(64Hz\):

Original Audio
Output Audio
Fig1. Original and decoded signal
Fig1. Original and decoded signal of 64 Hz
Fig2. Original and decoded signal zoomin
Fig2. Original and decoded signal of 64 Hz zoomin
Fig3. RMSE between original and decoded signal
Fig3. RMSE between original and decoded signal

\(65Hz\):

Original Audio
Output Audio
Fig4. Original and decoded signal
Fig4. Original and decoded signal of 65 Hz
Fig5. Original and decoded signal zoomin
Fig5. Original and decoded signal of 65 Hz zoomin
Fig6. RMSE between original and decoded signal
Fig6. RMSE between original and decoded signal

\(128Hz\):

Original Audio
Output Audio
Fig7. Original and decoded signal
Fig7. Original and decoded signal of 128 Hz
Fig8. Original and decoded signal zoomin
Fig8. Original and decoded signal of 128 Hz zoomin
Fig9. RMSE between original and decoded signal
Fig9. RMSE between original and decoded signal

\(129Hz\):

Original Audio
Output Audio
Fig10. Original and decoded signal
Fig10. Original and decoded signal of 129 Hz
Fig11. Original and decoded signal zoomin
Fig11. Original and decoded signal of 129 Hz zoomin
Fig12. RMSE between original and decoded signal
Fig12. RMSE between original and decoded signal


Question 2


Steps:
1) Adjust code to apply the windowing funciton.
2) Set bits to be 4 and window size to be 32.
3) Input signal of frequency \(128Hz\).
4) Compare the original and decoded signal and the difference (RMSE) between them.
5) Try different even and odd frequencies and compare the result.

Matlab code: here

Result:

\(128Hz\):

Original Audio
Output Audio
Fig13. Original and decoded signal
Fig13. Original and decoded signal of 128 Hz
Fig14. Original and decoded signal zoomin
Fig14. Original and decoded signal of 128 Hz zoomin
Fig15. RMSE between original and decoded signal
Fig15. RMSE between original and decoded signal

\(129Hz\):

Original Audio
Output Audio
Fig16. Original and decoded signal
Fig16. Original and decoded signal of 129 Hz
Fig17. Original and decoded signal zoomin
Fig17. Original and decoded signal of 129 Hz zoomin
Fig18. RMSE between original and decoded signal
Fig18. RMSE between original and decoded signal


Question 3

Solution



Question 4

Steps:
1) Set bits to be 4 and window size to be 32.
2) Input a chord signal 261.63Hz, 329.63Hz, 392.00Hz and 523.25Hz (C)
3) Compare the original and decoded signal and the difference (RMSE) between them.
4) Try different bits and compare the result.

Matlab code: here

\(4 bits\):

Original Audio
Output Audio
Fig19. Original and decoded signal
Fig19. Original and decoded signal of C chord with 4 bits
Fig20. Original and decoded signal zoomin
Fig20. Original and decoded signal of C chord with 4 bits zoomin
Fig21. RMSE between original and decoded signal
Fig21. RMSE between original and decoded signal

\(16 bits\):

Original Audio
Output Audio
Fig22. Original and decoded signal
Fig22. Original and decoded signal of C chord with 16 bits
Fig23. Original and decoded signal zoomin
Fig23. Original and decoded signal of C chord with 16 bits zoomin
Fig24. RMSE between original and decoded signal
Fig24. RMSE between original and decoded signal

\(32 bits\):

Original Audio
Output Audio
Fig25. Original and decoded signal
Fig25. Original and decoded signal of C chord with 32 bits
Fig26. Original and decoded signal zoomin
Fig26. Original and decoded signal of C chord with 32 bits zoomin
Fig27. RMSE between original and decoded signal
Fig27. RMSE between original and decoded signal


Question 5

Steps:
1) Input an audio signal.
3) Compare the original and decoded signal and the difference (RMSE) between them.
4) Try different bits with and without the windowing function, then compare the result.

Without Windowing:

\(4 bits\):

Original Audio
Output Audio
Fig28. Original and decoded signal
Fig28. Original and decoded signal of the audio file with 4 bits
Fig29. Original and decoded signal zoomin
Fig29. Original and decoded signal of the audio file with 4 bits zoomin
Fig30. RMSE between original and decoded signal
Fig30. RMSE between original and decoded signal

\(8 bits\):

Original Audio
Output Audio
Fig31. Original and decoded signal
Fig31. Original and decoded signal of the audio file with 8 bits
Fig32. Original and decoded signal zoomin
Fig32. Original and decoded signal of the audio file with 8 bits zoomin
Fig33. RMSE between original and decoded signal
Fig33. RMSE between original and decoded signal

\(16 bits\):

Original Audio
Output Audio
Fig34. Original and decoded signal
Fig34. Original and decoded signal of the audio file with 16 bits
Fig35. Original and decoded signal zoomin
Fig35. Original and decoded signal of the audio file with 16 bits zoomin
Fig36. RMSE between original and decoded signal
Fig36. RMSE between original and decoded signal

With Windowing:

\(4 bits\):

Original Audio
Output Audio
Fig28. Original and decoded signal
Fig28. Original and decoded signal of the audio file with 4 bits
Fig29. Original and decoded signal zoomin
Fig29. Original and decoded signal of the audio file with 4 bits zoomin
Fig30. RMSE between original and decoded signal
Fig30. RMSE between original and decoded signal

\(8 bits\):

Original Audio
Output Audio
Fig31. Original and decoded signal
Fig31. Original and decoded signal of the audio file with 8 bits
Fig32. Original and decoded signal zoomin
Fig32. Original and decoded signal of the audio file with 8 bits zoomin
Fig33. RMSE between original and decoded signal
Fig33. RMSE between original and decoded signal

\(16 bits\):

Original Audio
Output Audio
Fig34. Original and decoded signal
Fig34. Original and decoded signal of the audio file with 16 bits
Fig35. Original and decoded signal zoomin
Fig35. Original and decoded signal of the audio file with 16 bits zoomin
Fig36. RMSE between original and decoded signal
Fig36. RMSE between original and decoded signal