Hi i am stuck with a program in MATLAB to find the time domain impulse response of a system from its continuous transfer function in the frequency domain.(adsbygoogle = window.adsbygoogle || []).push({});

Here is the program-

delt=1.5625e-9; %definition of delta t(sampling time).To be taken sufficiently small depending upon the time constant of the original signal.

delomg=1/(N*delt); % Kept sufficiently small to prevent repetition of curves in the plot due to aliasing .

komg=[-100*pi/(delt):delomg:100*pi/(delt)]; %matrix containing the omegas over which the summation is done.Should be sufficiently large to cover all the frequencies present the original signal.

s=sqrt(-1)*komg; %definition of s=jw

transfer_func_image= 1 -(s.*r1*(c1+c2) +1)./(s.^2.*r1.*r2.*c1.*c2 + s.*(r1*c1 + r2*c2+ r1*c2) +1);

h_image_dirac=zeros(1,N); %initialisation of the h_image matrix for better speed

for n=1:N % first sample at t=0

t=(n-1)*delt;

sumint_image = sum(transfer_func_image.*exp(s.*t).*delomg);

h_image_dirac(n)=(1./(2*pi)).*sumint_image;

end

This program defines a transfer function in terms of s.

The last part of the program starting from for n= 1:N , then integrates the transfer function by multiplying the transfer function vector with exp(j*omega*t) and summing the vector obtained after multiplying it with delomg.This duplicates the formula of laplace transform for the special case where the real part of s =0 and is thus effectively inverse fourier transform.

The problem is the program is not giving the right result.If i find the value of the transfer function above at 2Meg, it is of the order of 10^-4

but if i use this program to obtain the impulse response and then convolve that with a sine wave at 2meg the answer is of the order of 100 and not 10^-4.

Where am i going wrong?This program works fine for another transfer function which i have tested.

please help.

thanks a lot.Please contact me if the problem is not clear.any help wud be appreciated.It is really important.

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Homework Help: MATLAB user defined inverse fourier transform

**Physics Forums | Science Articles, Homework Help, Discussion**