Write a matlab function, MATLAB Programming

Assignment Help:

There are many approaches to numerically estimating the derivative of the function.  The relationship:

1605_Write a Matlab function.png

 
is called a forward difference, since the estimate of the derivative at the point x is generated from the value of the function at x and the value of the function at x+?x (a point forward from x).  A backward difference formula estimates the derivative at the point x is from the value of the function at x and the value of the function at x-?x (a point backward from x) as follows:

2479_Write a Matlab function1.png

 
A central difference formula uses values of the function on either side of the point x in order to estimate the value of the function at x, as follows:

2074_Write a Matlab function2.png

 
The central difference formula has the advantage that it is more accurate, for a given ?x, than the forward or backward difference formulas.

In this problem, we will use a combination of the above techniques to estimate the derivative of a function from a set of measured data points.  We will use the following algorithm to estimate the derivative of the function: 
 
where fk is the value of the function at the kth data point, k f ′  is the estimate of the derivative of the function at the kth data point, N is the total number of data points, and  x ?  is the spacing between data points.  Note that we are using a forward difference formula for the first data point, a backward difference formula for the last data point, and a central difference formula for all other data points.

25_Write a Matlab function3.png

(a) Write a Matlab function which implements the above algorithm to estimate the derivative of a function described by a set of data points.  Your function should accept two inputs:

  • the first input is an N-element vector containing the measured values of the function (the fks)
  • the second input is the spacing between data points ( x ? )

Your function should return an N-element vector containing the estimated derivative of the function.  The size of the vectors should be determined by your function.  Your function does not need to do any error checking on the inputs.

(b) Write a script file which uses the function you created in part (a) to calculate and plot the derivative of the function  )x sin( )x(f π 2 = , for  1 0 ≤ ≤ x .  Your script file should do the following:

i.  Create data points for f(x) and use your function to estimate its derivative with  x ? =0.1
ii.  Create data points for f(x) and use your function to estimate its derivative with  x ? =0.01
iii.  Calculate the exact value of the derivative of f(x) for 1 0 ≤ ≤ x at increments of  x ? =0.01
iv.  Plot the above three sets of data on a single figure.  Include a legend to indicate which line corresponds to which case, and label the vertical and horizontal axes appropriately.


Related Discussions:- Write a matlab function

GAUSS ELIMNATION, 1;write a matlab of geuss elimination with scaling and pi...

1;write a matlab of geuss elimination with scaling and pivoting under consideration?

Illustration of variable number of input arguments, Illustration of Variabl...

Illustration of Variable number of input arguments: In this situation, it was supposed that the radius will always be passed to the function. The function header can hence be

Illustration of script documentation, Illustration of script Documentation:...

Illustration of script Documentation: For illustration, the script to compute the area of a circle could be modified to have comments: The first comment at the startin

Calculate the signal-to-noise ratio, Use the MATLAB randn function to gener...

Use the MATLAB randn function to generate 1000 points for x. Generate the output of the unknown system with the ?lter function and b=[1232 1] and a=[1]. Normalise the ?lter output

Elastic stiffness of a shaft segment, The natural frequencies of a fixed-fi...

The natural frequencies of a fixed-fixed uniform shaft made of homogeneous material undergoing longitudinal vibration is given by rad/s, where L is the length of the bar, G is th

Determine standard deviation, Find f(t) for each F(s): Problem: ...

Find f(t) for each F(s): Problem: Based on an examination given to a large class in which the maximum score is 100 points, assuming that 20 grades taken at random f

Printing vectors and matrices, Printing Vectors and Matrices: For vect...

Printing Vectors and Matrices: For vector, if the conversion character and newline character are in the format string, it will print in a column in spite of of whether the vec

Solution by using pdepe function, Solution by using pdepe function fu...

Solution by using pdepe function functionpdex1 m = 0; x = linspace(0,1,100); t = linspace(0,0.2,10);   sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t); % Ext

Illustration of modifying elements, Illustration of Modifying Elements: ...

Illustration of Modifying Elements: Illustration, the fifth element in the vector newvec is 9 >> newvec(5) ans = 9 The subset of a vector, that would be a vector it

Critical path method, can you please help me with matlab coding for CMP? I ...

can you please help me with matlab coding for CMP? I am new to matlab and hence need help

Write Your Message!

Captcha
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