Digital Signal Processing (ENGT5111)

Digital signal processing
%PART I
%problem 1
clear all
close all
clc
sfs = 98.5; %
Sampling frequency
fcuts = [0.5 1.0 45 46]; %
Frequency Vector
mags = [0 1 0]; %
Magnitude (Defines Passbands & Stopbands)
devs = [0.05 0.01 0.05]; %
Allowable Deviations
[n,Wn,beta,ftype] = kaiserord(fcuts,mags,devs,sfs);
n = n + rem(n,2);
hh = fir1(n,Wn,ftype,kaiser(n+1,beta),'scale');
figure(1)
freqz(hh, 1, 2^14, sfs)
%problem 2
[y,fs]=audioread('flyin_high.wav');
sound(y,fs);
figure(2)
plot(y)
grid on
title('Flyin_high Sound Signal')
xlabel('Samples')
ylabel('Continuous time sound signal')
x=y(10000:20000);
figure(3)
plot(x)
grid on
title('Flyin_high Sound Signal')
xlabel('Samples')
ylabel('Continuous time sound signal')
%Sampling the speech spectrogram
figure(4)
subplot(2,1,1)
specgram(y)
subplot(2,1,2)
specgram(x,256,fs) %outputs the line spectra
1

2