demo1.m

这是一个 York大学开发的FDTD代码。很好用的工具

disp('Demonstrate Pulse constructed from Fourier series')

disp('Generate time and frequency steps')
Ts=Tp/(N)
fs=1/Ts
f=linspace(0,fs*(N-1)/N,N);
disp('Generate time steps for one period (t)');
t=linspace(0,Tp*(N-1)/N,N);



disp('Generate time sequence (x) and plot it against time (t)')
x=prect(t,Tw,Tp);
title('Periodic Rectangular pulse train')
plot(t,x)
input('Any key to continue')

disp('Generate & plot Fourier coefficients for pulse')
k=-K:K;
X = Tw * sinc(k*(Tw)/Tp);
title('Fourier Coeficients')
plot(k,X,'*')
input('Any key to continue')

disp('Generate time sequence (xr) from coefficients and plot it against time (t)')
xr=idft(t/Tp,X,k);
title('Original and Reconstructed pulse')
plot(t,x,t,real(xr))