Analyze EEG data in frequency domain. EEG data is given in eeg_data.mat (data is taken from https://physionet.org/content/hmc-sleep-staging/1.1/ ). This file contains raw EEG data from different sleep stages (data1, , data5) and the sampling rate (fs). You will hand in a report and all files related to your code (functions, m-file, data etc.). Solve the following problems and add the results to your report. The report should include every plot requested in the problems. Make sure to add titles to your plots.

1) Use Fourier transform to obtain the frequency spectrum for each signal given in the data file. Comment on the active frequency bands and discuss possible sleep stages for each signal. You may use this article for information about sleep stages.

2) Design an FIR and an IIR filter for each EEG frequency band (alpha, beta etc.). Explain your reasoning for choosing a specific type of filter (Butterworth, Chebyshev etc. for IIR; type I, type II etc. for FIR). Plot the frequency and phase response of your filters (use Hz for x-axis) and comment on your results. (You may use built-in MATLAB functions for your filters)

3) Use the filters you have obtained in 2nd part to filter the signals to separate each signal to EEG frequency bands. Plot each band in time domain (there should be 2 plots for each band, IIR filtered and FIR filtered). Show the frequency spectrum of the filtered signals. Is there a difference between FIR or IIR filtered signals? Your MATLAB code should give the following outputs when run (separate from the results requested for the report);

1) Frequency response for data3

2) FIR filter frequency and phase response for theta band

3) IIR filter frequency and phase response for alpha band

4) Alpha filtered signal for data1 in time domain (both IIR and FIR filtered)

5) Beta filtered signal for data4 in frequency domain (both IIR and FIR filtered)

neccesary data is given in drive link

https://drive.google.com/file/d/1CH9U32XynCmjhe-SUuwHLcM1W6k11X1B/view?usp=sharing