[工学]数字信号处理第一章.ppt

1.3 时域离散系统 模拟非因果系统是不能实现的，对于有些数字非因果系统，利用系统中的存储性能，可以近似实现，只是系统的输出有延时。 实际恢复： D/A转换器 D/AC方框图 解码：将数字信号转换成时域离散信号 零阶保持器：将前一个采样值进行保持，一直到下一个采样值来到，相当于进行常数内插。 作业：新书P29 2-（1），（3），（5） 3-（2） 4 5 6 7 11 12 13 %复指数序列 c = -(1/12)+(pi/6)*i; K = 2; n = 0:40; x = K*exp(c*n); subplot(2,1,1); stem(n,real(x)); grid on; xlabel(Time index n);ylabel(Amplitude); title(Real part); subplot(2,1,2); stem(n,imag(x)); grid on; xlabel(Time index n);ylabel(Amplitude); title(Imaginary part); % 实指数序列 n = 0:35; a = 1.2; K = 0.2; x = K*a.^n; stem(n,x); xlabel(Time index n);ylabel(Amplitude); % 正弦序列 n = 0:40; f = 0.1; phase = 0; A = 1.5; x = A*cos(2*pi*f*n - phase); clf; % Clear old graph stem(n,x); axis([0 40 -2 2]); grid on; title(Sinusoidal Sequence); xlabel(Time index n); ylabel(Amplitude); function [y,n] = seqadd(x1,n1,x2,n2) % 序列相加函数 % 实现y(n) = x1(n)+x2(n) % y = 在包含n1和n2的n点上求序列和, % x1 = 在位置向量n1上的第一序列 % x2 = 在位置向量n2上的第二序列(n2可与 n1不同) % y(n)的长度 n = min(min(n1),min(n2)) : max(max(n1),max(n2)); y1 = zeros(1,length(n)); y2 = y1; % 初始化 % 具有y的长度的x1 y1(find((n=min(n1)) (n=max(n1))))=x1; % 具有y的长度的x2 y2(find((n=min(n2)) (n=max(n2))))=x2; % 序列相加 y = y1+y2; function [y,n] = seqmult (x1,n1,x2,n2) % 序列相乘函数 % 实现y(n) = x1(n)+x2(n) % y = 在包含n1和n2的n点上求序列和, % x1 = 在位置向量n1上的第一序列 % x2 = 在位置向量n2上的第二序列(n2可与 n1不同) % y(n)的长度 n = min(min(n1),min(n2)) : max(max(n1),max(n2)); y1 = zeros(1,length(n)); y2 = y1; % 初始化 % 具有y的长度的x1 y1(find((n=min(n1)) (n=max(n1))))=x1; % 具有y的长度的x2 y2(find((n=min(n2)) (n=max(n2))))=x2; % 序列相加 y = y1 .* y2; function [y,ny] = seqshift(x,nx,n0) % 实现 y(n) = x(n-n0) % n0为平移样本数 ny = nx + n0; % 位置向量移位 y = x; % 序列的值不变 nx = 0:5; x = 0.5.^nx; n0 = 3; [y,ny] = seqshift(x,nx,n0); subplot(2,1,1); stem(nx,x); axis([0 10 0 1.2]); xlabel(nx); ylabel(x); subplot(2,1,2); stem(ny,y); axis([0 10 0 1.2]); xlabel(ny); ylabel(y);