Matlab computer science

1. CLEANING A REAL-WORLD SIGNAL (35 PTS) You’ve just recorded a biosignal, but unfortunately someone was running a motor in the lab throwing off all sorts of interference such that your signal is corrupted with a massive amount of noise as shown in Fig 1. This is really important and you don’t have the opportunity to rerecord it, so you’re going to have to clean it up and figure out what it is from what you have.
2. Display and review the signal (5 pts, 1 pts each)
3. Load the data file Lab8Data. This will load a structure that includes the time each sample was recorded as well as the amplitude.
4. Plot this raw signal as a subplot with appropriate title and labels

iii.Since this is real data and you have the recording times, find the times between each sample

1. Plot the histogram of the time differences between samples as a second subplot.
2. Compute the average time between samples. Does this match with your histogram?
3. Analyzing the Frequency Content (5 pts, 2.5 pts each)

You know there’s noise, but you can’t do much about it until you know what frequency it’s at.Take a look bycompleting the following steps:

1. Compute the Fourier transform of the corrupted signal using your estimated dt from part A.
2. Plot the Fourier transform with appropriate title and labels. Make a note of the larger spikes – these arethe dominant frequencies of the signal.
3. Display Spectrogram (10 pts)

The time domain analysis of the corrupted signal doesn’t show much. You should see something in the frequencydomain. Create a spectrogram using the following steps to see if there are any time varying components to thesignal.

1. Create a vector of frequencies from 0 to half the sampling frequency (2 pts).
2. Compute the spectrogram using a window width that is a power of 2 (16, 32, 64) and an overlap of 80%of the window width. You may have to try a few window widths to find the best one (2 pts).

iii.Plot the spectrogram with appropriate title and labels. Can you identify the frequency of the noise? Doesit match what you saw in the Fourier transform plot? (6 pts)

1. Filter the Signal (15 pts)

OK, now you’ve identified the frequency of the noise, it’s time to clean it up. Thing is, you’re not sure how much cleaning you should do. You don’t want to over-filter it, but you also need to be able to identify the signal you recorded. Try a range of filter orders by following these steps.

1. You’ll be looking over five different filter orders. For comparison, plot the raw signal in the first of 6 subplots. (3 pts)
   1. For close examination, set your x axis limits to 0.5 to 5 s.
2. Use a loop to (10 pts):
   1. Compute the parameters for the appropriate type of Butterworth filter using a single cutoff frequency for all iterations. The exact frequency is up to you, but it should be away from the “legitimate” frequencies of your signal, but well enough before your noise so that it’s well attenuated.
   2. Iterate the order of the filter from 1 to 5.
   3. Plot the filtered signal with appropriate title (including the order in the title) and axis labels as separate subplots.
      1. For close examination, set your x axis limits to 0.5 to 5 s.

• Identify the signal. Use the disp command to send a message to the workspace with you guess as to what the biosignal recorded (hint: it’s regular and follows the frequency shown when someone is at rest). (2 pts)

2. MAKE CASE ORANGE (6 PTS)

It’s Fall and just like the leaves are changing, so is Case. Download the logo from Canvas and use Matlab to revise it to be orange instead of blue. Hint: check what kind of image it is and use the information in one matrix to change the values of the others to set the blend of colors.

3. FIND THE HIDDEN MESSAGE (4 PTS)

An important message has been recorded in the file Lab8Q3MSG.wav on Canvas, but someone botched the sampling rate. Play the signal and find the best sampling frequency to play it back at in order to determine the message. Use the disp function to write this message to the Command Window.