%各函數N=16的譜分析
n=0:15;
x1=[ones(1,4) zeros(1,12)];
x2=cos(pi*n/4);
x3=sin(pi*n/8);
f1= abs (kqsfft(x1));
f2= abs (kqsfft(x2));
f3= abs (kqsfft(x3));
subplot(3,1,1);
stem(n,f1)
xlabel('n');ylabel('幅度');title('x1譜分析幅度譜');
subplot(3,1,2);
stem(n,f2)
xlabel('n');ylabel('幅度');title('x2譜分析幅度譜');
subplot(3,1,3);
stem(n,f3)
xlabel('n');ylabel('幅度');title('x3譜分析幅度譜');