Question

考察序列x(n)=cos(0.48πn)+cos(0.52πn)分别取n=10和n=100两种情况下的DFT。Matlab程序问题？？

首先，确定x(n)的高密度频谱：
>>n=[0:1:99];x=cos(0.48*pi*n)+cos(0.52*pi*n);
>>n1=[0:1:9];y1=x(1:1:10);
>>subplot(2,1,1);stem(n1,y1);title(‘信号x(n),0<=n<=9’);
>>axis([0,10,-2.5,2.5]);text(10.2,-2.5,’n’)
>>Y1=dft(y1,10);magY1=abs(Y1(1:1:6));
>>k1=0:1:5;w1=2*pi/10*k1;
>>subplot(2,1,2);plot(w1/pi,magY1);title(‘DFT幅度的样本’);
>>xlabel(‘频率(单位:pi)’)
填90个零得到的密频谱：
>>n2=[0:1:99];y2=[x(1:1:10) zeros(1,90)];
>>subplot(2,1,1);stem(n2,y2);title(‘信号x(n),0<=n<=9+90zeros’);
>>Y2=dft(y2,100);magY2=abs(Y2(1:1:51));
>>k2=0:1:50;w2=2*pi/10*k2;
>>subplot(2,1,2);plot(w2/pi,magY2);title(‘DTFT幅度’);
>>xlabel(‘频率(单位:pi)’)
将x(n)的头100个样本来确定x(n)的高分辨率频谱：
>>subplot(2,1,1);stem(n,x);title(‘信号x(n),0<=n<=99’);
>>xlabel(’n’)
>>X=dft(x,100);magX=abs(X(1:1:51));
>>k1=0:1:50;w=2*pi/100*k;
>>subplot(2,1,2);plot(w/pi,magX);title(‘DTFT幅度’);
>>xlabel(‘频率(单位:pi)’)
麻烦高手帮忙修改一下程序！！现在只能运行出一个图，其它的都出错。

Answer 1

兄弟 你有没有写dft的函数