homework_one2151.asv

DSP信号处理

%题目2-1-5 演示程序：
%f1=1Hz，f2=3Hz，相对振幅delta=0.5，初相位phi=180，抽样间隔T=0.1秒，序列长length=10
clc;clf;
f1=1;f2=3; delta=0.5; phi=180; T=0.1; N=10;             %设置参数
n=0:N-1;                                              %设置序列长度
k=0:N-1;   

xn1=sin(2*pi*f1*T.*n);                                %产生序列xn1        
xn2=sin(2*pi*f2*T.*n+phi);                            %产生序列xn2   
xn=xn1+delta*xn2;                                     %产生复合函数序列xn

Xk1=fft(xn1,N);                                       %对序列xn1进行DFT变化
Xk2=fft(xn2,N);                                       %对序列xn2进行DFT变化
Xk=fft(xn,N);                                         %对复合序列xn进行DFT变化
 

%画图:比较Xk1+delta×Xk2与Xkd的实部、虚部、幅频特性及相频特性来讨论DFT的线性性质
subplot(2,4,1);stem(k,real(Xk1+delta*Xk2),'filled'); 
xlabel('k');ylabel('real（Xk1+delta*Xk2）');
title('Xk1+delta*Xk2 实部'); grid on;%Xk1+delta*Xk2 实部
subplot(2,4,5);stem(k,real(Xk),'filled'); 
xlabel('k');ylabel('real（Xk）');title('Xk实部'); grid on;%Xk实部
subplot(2,4,2);stem(k,imag(Xk1+delta*Xk2),'filled'); 
xlabel('k');ylabel('imag（Xk1+delta*Xk2）');
title('Xk1+delta*Xk2虚部'); grid on;%Xk1+delta*Xk2 虚部
subplot(2,4,6);stem(k,imag(Xk),'filled'); 
xlabel('k');ylabel('imag（Xk）');title('Xk虚部'); grid on;%Xk虚部
subplot(2,4,3);stem(k,abs(Xk1+delta*Xk2),'filled'); 
xlabel('k');ylabel('|Xk1+delta*Xk2|');
title('Xk1+delta*Xk2幅频特性'); grid on;%Xk1+delta*Xk2 幅频特性 
subplot(2,4,7);stem(k,abs(Xk),'filled'); 
xlabel('k');ylabel('|Xk|');title('Xk幅频特性'); grid on;%Xk幅频特性 
subplot(2,4,4);stem(k,angle(Xk1+delta*Xk2),'filled');
xlabel('k');ylabel('angle(Xk1+delta*Xk2)');
title('Xk1+delta*Xk2相频特性'); grid on;%Xk1+delta*Xk2 相频特性
subplot(2,4,8);stem(k,angle(Xk),'filled');
xlabel('k');ylabel('angle(Xk)');title('Xk相频特性'); grid on;%Xk相频特性