Compare filter with that given by matlab function bilinear

Assignment Help MATLAB Programming
Reference no: EM13869640

Assignment: COMPUTATIONAL AND EXPERIMENTAL METHODS

Experimental Methods- Signal Filtering

Several sets of data have been generated which contain ~4000 data points, sampled at 100 Hz and are made up of the superposition of several sine waves at different frequencies and of different magnitude. You can find your specific file and damping ratio to use via the VLE. This should be used in a filter as discussed below. The data files can be downloaded from the assignments section.

After copying to your user area and unzipping the files, the appropriate file can be loaded into a Matlab program by entering
load filename

The data to be analyzed will be contained in a variable y and is of the form y = f(t)

1) Using Matlab, read in the data file and, for the first 100 data points, produce a graph of amplitude against time.

2) Fourier transform this data (noting that there are 4000 data points) and produce a power spectrum, clearly stating the units on each axis. Hence, for each sine wave determine the constants A and W in the equation y = A sin(Wt).

3) Manually design a second order digital filter, with the damping ratio given in the table, to attenuate the highest frequency component found as much as possible while changing the amplitude of the other waves as little as possible. Clearly justify any choice of filter parameters. Manually this filter using bilinear transformation.

4) Compare your filter with that given by the Matlab function ‘bilinear' for the same conditions. The syntax for this filter is

[n,d] =bilinear(num ,[ den1 den2 den3] ,sample _frequency )

Where num is the numerator of the second order filter and den1, den2 and den3 are the three terms in the denominator. The numerator and denominator frequencies should be in rad/sec while the sample frequency should be in Hz, to match the input data.

5) Using the Matlab function ‘filter', use the second order filter to remove the highest frequency component from the original signal. Produce a graph of amplitude against time and a periodogram (power spectrum) for the filtered data and original data.

6) Manually calculate the amplitude ratio of the two highest frequencies, i.e. the frequency that should have been removed and the next highest one that should have been ‘passed' by the filter) and compare the values with those suggested by the periodogram. Suggest reasons for any differences. What would have been the effect of (a) using a first order filter, and (b) changing the damping ratio.

You may find it helpful to work through the file MatlabExercise_FrequencyAnalysis_LowPassFiltering.pdf before attempting the assignment

Reference no: EM13869640

Questions Cloud

Identify key characteristics of the facility : Identify key characteristics of the facility
Identify the criteria used in ranking top global companies : Find one such ranking system and identify the criteria used in ranking top global companies. Although some of these rankings require subscriptions.
Calculate the financial ratios : Calculate the financial ratios
What is the standard deviation of return of stock : What is the standard deviation of return of stock
Compare filter with that given by matlab function bilinear : Compare your filter with that given by the Matlab function ‘bilinear' for the same conditions. Manually design a second order digital filter, with the damping ratio given in the table.
What is juniors weighted average cost of capital : What is Juniors Weighted Average Cost of Capital
Difference between calling a bond and a bond refunding : difference between calling a bond and a bond refunding.
Thesis statement is poverty a choice : Write 700 to 1000 words about - Thesis statement Is poverty a choice
What is overgrazing? : What is overgrazing? What effect does this have on rangelands?

Reviews

Write a Review

 

MATLAB Programming Questions & Answers

  Matlab inline functions

Consider g(t)=e^(-2t) cos?(πt)(u(t)-u(t-2.5)), use Matlab inline functions for the following questions

  Mobile robot control and potential field

Write a MATLAB code to compute forces due to each obstacle and goal. Compute total force on vehicle at point(x,y).

  Consider a cstr with a feed stream

Consider a CSTR with a feed stream containing only A at a concentration of cA and the reactions A?B?C taking place in the CSTR. Both reactions are first order in the reactant

  Calculate an approximation to the energy consumed

Calculate an approximation to the energy consumed over the period and calculate an approximation to the energy consumed over the period t = a to t = b.

  Design simple interface hardware

Evaluate the capability in selecting and configuring a microcontroller unit designing simple interface hardware and writing subroutines to operate those interfaces

  An image histogram for a color image

goal is to build an image histogram for a color image based on its 6-bit color code.

  Vb and matlab to calculate the cumulative days

Create a function TotalDaysMonth(intYear, intMonth, strIndOrCum) in both VB and Matlab to calculate the cumulative days or number of days in a given month, depending on inputs, by employing array vector in VB or a vector matrix in Matlab.

  How much will an object weigh on this point

The value of acceleration due to gravity (g) on a point 10,000 kilometers above sea level is about 1.49 meters/second2. How much will an object, which weighs 98 newtons on the surface of Earth, weigh on this point

  Find the probability

Find the probability that more than 25 requests are received in the one minute interval by applying the central limit theorem. (Hint: apply the de Moivre-Laplace approximation)

  Calculates the required input torque and the magnitudes

Calculates the required input torque (T) and the magnitudes of the pin reaction forces (A, B, and C, along both x and y directions) for two full revolutions (0 ≤ θ ≤ 4Π) of the crank.

  Estimation of the rayleigh limit

Propose a method that would allow estimation of the Rayleigh limit experimentally - Rayleigh frequency resolution limit

  Using dimensional analysis nondimensionalize the period t

the period of a pendulum t is assumed correctly or incorrectly to depend on the mass m the length of the pendulum l the

Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd