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')%该信号频谱图

你再试一试看看，我运行了是正确的

◆◆

评论读取中....

请登录后再发表评论!

◆◆

修改失败，请稍后尝试

信号分析与处理第一次上机作业

例：矩形脉冲周期信号频谱分析的MATLAB实现。

①矩形脉冲周期信号的时域波形；

②矩形脉冲周期信号的频谱：实部和虚部；

③矩形脉冲周期信号的频谱：幅值和相位。

MATLAB程序：

t=-10:0.01:10;

y=0.5*(square(0.4*pi*(t+0.5),20)+1);

plot(t,y);grid;axis([-10,10,-0.1,1.2]);

title('矩形脉冲周期信号'),xlabel('t'),ylabel('f(t)');

n=-30:30;

e=1;tao=2;zq=5;w=(2*pi)/zq;

xr=(e*tao/zq).*sinc(n.*tao./zq);

xi=zeros(61,1);

figure(2)

subplot(2,1,1),stem(n,xr,'.');grid;

xlabel('k'),ylabel('Real Part of X(k)');

subplot(2,1,2),stem(n,xi,'.');grid;

xlabel('k'),ylabel('Imaginary Part of X(k)');

n=-30:30;

e=1;tao=2;zq=5;w=(2*pi)/zq;

x=abs((e*tao/zq).*sinc(n.*tao./zq));

y=atan2(0,(e*tao/zq).*sinc(n.*tao./zq));

figure(3)

subplot(2,1,1),stem(n,x,'.');grid;

xlabel('k'),ylabel('Magnitude Part of X(k)');

subplot(2,1,2),stem(n,y,'.');grid;

xlabel('k'),ylabel('Phase Part of X(k)');

优质解答

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')%该信号频谱图

你再试一试看看,我运行了是正确的