实验一 快速傅里叶变换及其应用(7)

xlabel('t/T'); ylabel('xa(n)'); xk1=abs(fft(xa1)); title('p=8 q=8') subplot(5,2,6); stem(n,xk1) xlabel('k'); ylabel('Xa(k)'); title('p=8 q=8')

p=13;q=8; %p=13;q=8; xa1=exp(-((n-p).^2)/q); subplot(5,2,7); plot(n,xa1,'-*'); xlabel('t/T'); ylabel('xa(n)'); xk1=abs(fft(xa1)); title('p=13 q=8')

xa(n)

0 xa(n)

xa(n)

xa(n)

xa(n)

stem(n,xk1) xlabel('k'); ylabel('Xa(k)'); title('p=13 q=8')

p=14;q=8; %p=14;q=8; xa1=exp(-((n-p).^2)/q); subplot(5,2,9); plot(n,xa1,'-*'); xlabel('t/T'); ylabel('xa(n)'); title('p=14 q=8') xk1=abs(fft(xa1)); subplot(5,2,10); stem(n,xk1) xlabel('k'); ylabel('Xa(k)'); title('p=14 q=8’)

15

0 0

p=8 q=2

5 10 5 10 15

t/T

k

t/T k p=8 q=8

t/T k

t/T k

t/T

5

k

10 15

n=0:1:15;

a=0.1;f=0.0625;

xb1=exp(-(a*n)).*sin(2*pi*f*n) subplot(3,2,1);