Find the integral of a function at an arbitrary location
Course:- MATLAB Programming
Reference No.:- EM13465

Assignment Help
Expertsmind Rated 4.9 / 5 based on 47215 reviews.
Review Site
Assignment Help >> MATLAB Programming

1. We want to find the integral of a function at an arbitrary location x from the origin.  Thus,

2403_find the integral of a function at an arbitrary location.png

where I(x=0) is the value of the integral for all times less than 0.  (Essentially, I(x=0) is the unknown constant of integration or the initial condition.)  From the class lecture on the trapezoidal rule for numerical integration, it can be seen that this function can be approximated by

52_find the integral of a function at an arbitrary location1.png

Write a Matlab function to perform numerical integration of a set of evenly spaced data points using the trapezoidal rule.  Your function should accept two vectors as inputs, x and f.  The first vector (x) contains the independent variable data (the points at which the values of the function are known.  The second vector (f) should contain the values of the function at the points provided in the first vector.  Your function should return the integral of f with respect to x, as a function of x.

2.  The mass of the spring-mass-damper system shown below was initially displaced by 1 meter to the right and then released. The subsequent acceleration of the mass was measured using an accelerometer; this data is contained in the file smd_data.txt posted on the class web site. The file contains ascii data which provides the acceleration of the mass as a function of time; the first column of data in the file is time in seconds, the second column is the acceleration in meters/sec2. We wish to determine the maximum compressive force experienced by the spring during the time range provided in the data.  The spring force is equal to the spring rate times the displacement,  Fspring = kx.

1788_find the integral of a function at an arbitrary location2.png

Write a script file which will determine the displacement of the mass as a function of time from the acceleration data.  Trapezoidal integration will provide sufficient accuracy to infer the position and velocity of the seat as a function of time.  Your program should plot the position of the mass as a function of time.  You can use this plot to estimate the maximum negative displacement of the mass; the maximum compressive force is this displacement times the spring rate (k = 100 N/m).

 3.  We wish to numerically perform the following integration: 

1175_find the integral of a function at an arbitrary location3.png

(a) Use Matlab's quad function to perform the above integral.

(b) Analytically evaluate the integral (e.g. do the integration by hand).

Put your comment

Ask Question & Get Answers from Experts
Browse some more (MATLAB Programming) Materials
Create a namespace titled global Type that contains four members. The members should be two named constants, n and rate. It should also contain a variable count and the func
Write a MATLAB function called newtons Method that accepts two inputs and returns one output. The inputs should be a vector representation of a polynomial and an initial gue
Write a mileage per Gallons and miles per dollar calculator program using two functions, one to calculate MPG and one to calculate MPD. The application will prompt the user
Two arrays are considered to be equal when the result of one array subtracted from the other is an array of zeros. In the script file, create the MATLAB commands to prove th
Write mass balance equations for the remaining 4 reactors. Put the known values in the 5 equations and convert each equation into a form where left hand side contains all th
Write all of the above sets of equations in the matrix format and identify the sets that you cannot solve by using an iterative method such as Jacobi and/or Gauss-Seidel met
Implement specific MATLAB functions to encode the categorical attribute values to numeric values. For example, column 1 consists of only categorical values "a" and "b", and
Create a loop for each parameter in the magnitude vector and You can do this with another for loop, but that is a bad way to do it. MATLAB is a matrix based programming langu