p_125_6_1.m

eda实例DSP上机程序

clc
clear
xn1=[1 0 0 0 0 0 0 0 0];
xn2=ones(1,9);
xn3=[1 2 3 4 5 4 3 2 1];
n=0:8;
xn4=sin(pi*n/4);
figure(1);
subplot(3,4,1),stem(n,xn1);xlabel('n'),ylabel('xn1');title('单位脉冲序列');
subplot(3,4,2),stem(n,xn2);xlabel('n'),ylabel('xn2');title('矩形序列');
subplot(3,4,3),stem(n,xn3);xlabel('n'),ylabel('xn3');title('三角序列');
subplot(3,4,4),stem(n,xn4);xlabel('n'),ylabel('xn4');title('正弦序列');
N=9;
Xk1=fft(xn1,N);
Xk2=fft(xn2,N);
Xk3=fft(xn3,N);
Xk4=fft(xn4,N);
subplot(3,4,5),stem(n,Xk1);xlabel('N'),ylabel('Xk1');title('单位脉冲序列DFT');
subplot(3,4,6),stem(n,Xk2);xlabel('N'),ylabel('Xk2');title('矩形序列DFT');
subplot(3,4,7),stem(n,Xk3);xlabel('N'),ylabel('Xk3');title('三角序列DFT');
subplot(3,4,8),stem(n,Xk4);xlabel('N'),ylabel('Xk4');title('正弦序列DFT');
xn11=ifft20(Xk1,N);
xn22=ifft20(Xk2,N);
xn33=ifft20(Xk3,N);
xn44=ifft20(Xk4,N);
subplot(3,4,9),stem(n,xn11);xlabel('N'),ylabel('xn11');title('单位脉冲序列');
subplot(3,4,10),stem(n,xn22);xlabel('N'),ylabel('xn22');title('矩形序列');
subplot(3,4,11),stem(n,xn33);xlabel('N'),ylabel('xn33');title('三角序列');
subplot(3,4,12),stem(n,xn44);xlabel('N'),ylabel('xn44');title('正弦序列');