利用matlab做dsp实验,华工DSP实验一利用Matlab进行频谱分析

实验一利用Matlab进行频谱分析

一、对无限长时域离散信号进行频谱分析

1．代码：

function test1(N)

N = input('Type in the length = ');

k = 0:1:N-1;

subplot(2,1,1)

x = cos(pi*k/10)+sin(pi*k/6)+cos(2*pi*k/5);

X = fft(x);

stem(k-N/2,abs(fftshift(X)));

title('rectangular window')

xlabel('Frequency index k'); ylabel('Magnitue') subplot(2,1,2)

h=hamming(N)';

x = cos(pi*k/10)+sin(pi*k/6)+cos(2*pi*k/5).*h;

X = fft(x);

stem(k-N/2,abs(fftshift(X)));

title('hamming window')

xlabel('Frequency index k'); ylabel('Magnitude') 2．运行结果：

Type in the length = 16

Type in the length = 64