- 1

- 0

Hey everyone,

So what I'm trying to do is to find the outgoing signal in function of the time, when

the transfer function of the filter is known and the incoming signal is a periodic piecewise function.

The incoming signal I(t)= 0 for -T/2<t<-T/4

[tex]=\cos(\frac{t\pi}{T}) [/tex] for -T/4<=t<T/4

= 0 fot T/4<=t<T

And the transfer function I'm using is

[tex] S(\omega)=\sqrt{10}(\frac{i\omega/\omega_2}{1+i\omega/\omega_2})^2, [/tex]

where [tex]\omega_2=9301.85 rad/s [/tex].

I have a solution for the outgoing signal U(t) but I'm not sure I'm doing it right.

This is my source code:

T=1*10^(-3);w2=9301.82; t=-0.005:1/100000:0.005;

y=incoming(t,T); %Sampling my incoming signal

Y=fft(y); %taking the fast fourier transform

N=length(Y);

freq=[0:N-1];

U=Y.*filter(freq,w2) %where filter is the transfer function

u=ifft(U);

plot(t,ifftshift(u));

So what I'm trying to do is to find the outgoing signal in function of the time, when

the transfer function of the filter is known and the incoming signal is a periodic piecewise function.

The incoming signal I(t)= 0 for -T/2<t<-T/4

[tex]=\cos(\frac{t\pi}{T}) [/tex] for -T/4<=t<T/4

= 0 fot T/4<=t<T

And the transfer function I'm using is

[tex] S(\omega)=\sqrt{10}(\frac{i\omega/\omega_2}{1+i\omega/\omega_2})^2, [/tex]

where [tex]\omega_2=9301.85 rad/s [/tex].

I have a solution for the outgoing signal U(t) but I'm not sure I'm doing it right.

This is my source code:

T=1*10^(-3);w2=9301.82; t=-0.005:1/100000:0.005;

y=incoming(t,T); %Sampling my incoming signal

Y=fft(y); %taking the fast fourier transform

N=length(Y);

freq=[0:N-1];

U=Y.*filter(freq,w2) %where filter is the transfer function

u=ifft(U);

plot(t,ifftshift(u));

Last edited: