验证频率采样定理
clc
clear all;
close all;
a=0.9;
L=10;
n1=-L:L;
xn=a.^abs(n1);
% n1=-10:10;
% xn=(0.9).^abs(n1);
xn(12)
figure
stem(n1,xn,'.');
axis([-11,11,0,1.2]);
title('x(n)');
xlabel('n1');
ylabel('x(n)');
% N=30; sample
for k=0:29
Ck(k+1)=1;
for n2=1:10
Ck(k+1)=Ck(k+1)+2*xn(n2+1+10)*cos(2*pi*k*n2/30);
end
end
%attention:xn(n2+1+10)is not begin from xn(1),it's from xn(12)
x30n=ifft(Ck,30);
figure
k=0:29;
stem(k,x30n,'.');
axis([0,30,0,1.2]);
% N=15; sample
for k=0:14
Ck(k+1)=1;
for n2=1:10
Ck(k+1)=Ck(k+1)+2*xn(n2+1+10)*cos(2*pi*k*n2/15);
end
end
%attention:xn(n2+1+10)is not begin from xn(1),it's from xn(1