基于MATLAB的PCM编码

% Removing all variables functions and MEX-files from memory leaving the
% workspace empty.
clear all


% Deleting all figures whose handles are not hidden.
close all


% Deleting all figures including those with hidden handles.
close all hidden


% Clearing all input and output from the Command Window display giving us a clean screen.
clc


% Opening the file ‘TEOTH.mp3‘ in the read access mode.
fid = fopen （‘TEOTH.mp3‘‘r‘）;


% Generating the input signal ‘m（t）‘ by reading the binary data in 16 bit
% integer format from the specified file and writing it into a matrix
% ‘m（t）‘.
m = fread （fid‘int16‘）;


% Defining the count for efficiency.
count = 8500;


% Calculating maximum value of the input signal ‘m（t）‘.
Mp = max （m）


% Setting number of bits in a symbol.
bits = 8;


% Defining the number of levels of uniform quantization.
levels = 2^bits;


% Calculating the step size of the quantization.
step_size = （2*Mp）/levels


% Setting the sampling frequency.
% because the audio signal has a maximum frequency of 4K and according to
% Nyquist criteria we get the following sampling frequency.
Fs = 8000;


% Setting the sampling instant.
Ts = 1;


% Setting the number of samples to be used.
No_Samples = （2*Fs）+Ts;


% Define the time vector for the calculations.
time = [1:Fs/64];


% Calculating the bit rate.
bit_rate = 8000*bits;