优质解答
x1=ones(1,2);
x2=[x1,zeros(1,6)];
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk;
magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k,magXk);
xlabel('k');ylabel('Xtilde(k)');
title('DFS of SQ.wave :L=2,N=8')%该信号频谱图
不好意思写掉了一点
x1=ones(1,2);
x2=[x1,zeros(1,6)];
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
k1=[-N/2:N/2];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk; magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k1,magXk);
xlabel('k1');ylabel('Xtilde(k)');
title('DFS of SQ.wave :L=2,N=8')%该信号频谱图
你再试一试看看,我运行了是正确的