homework_one213.asv

DSP信号处理

%题目2-1-3  演示程序：f=1Hz，初始相位delta=0，抽样间隔T=0.1秒，序列长length=10
clc;clf;
f=0.5; delta=0; T=0.1; N=10;              %设置参数
t=0:pi/30:N*T;                            %设置步长
n=0:N-1;                                  %设置序列长度
x(t)=sin(2*f*t+delta);                    %产生函数x(t)
xn=sin(2*f*T+delta);                      %对函数x(t)进行抽样
Xk=fft(xn,N);                             %对序列xn进行DFT变化
k=0:N-1;    

%画图
figure(1);
subplot(2,1,1);stem(t,x(t));xlabel('t');ylabel('函数x(t)'); 
title('时域上：抽样前的函数x(t)');grid on;  %画函数x(t)
subplot(2,1,2);stem(n,xn),'filled');xlabel('n');ylabel('xn'); 
title('时域上：抽样后的序列x(n)');grid on;  %画抽样序列x(n)

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;%相频特性