homework_one212.m

DSP信号处理

%题目2-1-2  演示程序：a=0.5, b=0.8, length=10
clc;clf;
a=0.5;b=0.8;p=a+j*b;N=10;                 %设置参数
n=0:N-1;                                  %设置序列长度范围
xn=p.^n;                                  %产生序列xn

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

%画图
figure(1);
subplot(2,1,1);stem(n,real(xn),'filled');xlabel('n');ylabel('real（xn)'); 
title('复指数序列（时域）实部x(n)');grid on;        %复指数序列实部xn（时域）
subplot(2,1,2);stem(n,imag(xn),'filled');xlabel('n');ylabel('imag（xn）'); 
title('复指数序列（时域）虚部x(n)');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;%相频特性