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
Some systems show a marked change in behavior as parameters change. For example, a damped harmonic oscillator with exhibit over damped, under damped, or critically damped be
how do you solve this problem and why does MATLAB keep telling that you need more input arguments - Create a function that returns true if the input is a numeric array and fal
Analyse the bar using the Galerkin finite element method. The coding solution for linear h-refinement and its manual solution for a limited number of linear elements areon B
the script that will graphically show the range of frequencies for which the ouput amplitude is less than 70% of the input amplitude. The original problem statement
Write a Matlab script using the "randn" command to generate a vector of 10,000 random numbers that are normally distributed around zero, with a variance of 1. Use the "hist"
create an animation of a square according to the directions below. Create a square with a side length of 3 with the square centered at (5,5). Translate this square along a cos
Create two vectors of length 10 of your choice, and label them x and z. Show, element by element, where x and z are equal. Show, element by element, where x is greater than z.
MECH 3010: Elementary Numerical Methods and Programming - Calculate the approximate percent relative error in θ1 at t = 4 s between the current step size and the next largest