homework_one2112.m

DSP信号处理

%题目2-1-1  演示程序：a=0.9，length=10
clc;clf;

a=0.9;N=10;                 %设置参数
n=0:N-1;                    %设置序列长度范围
xn=a.^n;                    %产生序列xn

Xk=fft(xn,N);               %对序列xn进行DFT变化
k=0:N-1;    

%画图
figure(1);
stem(n,xn,'filled');xlabel('n');ylabel('x(n)'); 
title('实指数序列x(n)：a=0.9, length=10 ');grid on;%序列xn（时域）

figure(2);
subplot(2,2,1);stem(k,real(Xk),'filled');
xlabel('k');ylabel('real(Xk)');title('频域（DFT）实部'); grid on; %实部

subplot(2,2,2);stem(k,imag(Xk),'filled'); 
xlabel('k');ylabel('imag(Xk)');title('频域（DFT）虚部'); grid on; %虚部

subplot(2,2,3);stem(k,abs(Xk),'filled'); 
xlabel('k');ylabel('|Xk|');title('频域（DFT）幅频特性'); grid on;%幅频特性 

subplot(2,2,4);stem(k,angle(Xk),'filled');
xlabel('k');ylabel('angle(Xk)');title('频域（DFT）相频特性'); grid on;%相频特性