endof-

Chapter 8 The Fourier Transform11.We could also have used numerical integration with the integralthereby avoiding the integration by parts procedure. Shown below are the functionfourierxfm2 which is saved as an .m file and the program execution using this function.function y2=fourierxfm2(x)x=x+(x==0)*eps;y2=(sin(x)./x).^2;and after this file is saved, we execute the statement below observing that the limits of integrationare from 0 to π .value2=quad('fourierxfm2',0,pi)value2 =1.4182fplot('abs(2.*exp(−j.*w)*(sin(w)/w))',[0 4*pi 0 2])∫π02xsin------------x 2 dx21.510.500 2 4 6 8 10 128−60Signals and Systems with MATLAB ® Computing and Simulink ® Modeling, Fourth EditionCopyright © Orchard Publications