数字信号处理——离散信号与系统分析

一、离散信号的产生。

数字信号处理，一般是建立在下图的架构中，由FPGA/DSP为核心，以ADC和DAC为桥梁，来处理外部的模拟信号。

其中FPGA/DSP处理的都是数字信号，而数字信号属于离散信号。

数字信号：时间、幅度都为离散的。

离散信号：时间离散，幅度可以是离散、也可以是连续的。

二、基本的离散信号。

2.1、单位脉冲序列。

clear;close all;clc;
N = 15;
n = -N : N;
x = zeros(1,2*N+1);
x(N+1) = 1;
stem(n,x);

2.2、单位阶跃序列。

clear;close all;clc;
N = 15;
n = -N : N;
x = zeros(1,2*N+1);
x(N+1:2*N+1) = 1;
stem(n,x);

2.3、矩形序列。

clear;close all;clc;
N = 15;
M = 5;
n = -N : N;
x = zeros(1,2*N+1);
x(N+1:N+1+M) = 1;
stem(n,x);

2.4、复指数序列。

clear;close all;clc;
N = 15;
n = -N : N;
x = exp((-0.1+0.5i)*pi*n/8);
xReal = real(x);
xImag = imag(x);
xAbs = abs(x);
xAngle = angle(x);
subplot(221)
stem(n,xReal);
subplot(222)
stem(n,xImag);
subplot(223)
stem(n,xAbs);
subplot(224)
stem(n,xAngle);

将x改为x = exp(0.5i*pi*n/8);时

2.5、正弦序列。

clear;close all;clc;
N = 15;
n = -N : N;
x = sin(0.5*pi*n/8);
stem(n,x);

三、序列的基本运算。

3.1 序列翻转。

clear;close all;clc;
N = 15;
n = -N : N;
x = zeros(1,length(n));
for i=1:5
	x(i+N) = i;
end
subplot(211)
stem(n,x);
subplot(212)
stem(n,flip(x));

3.2 序列平移。

clear;close all;clc;
N = 15;
n = -N : N;
x = zeros(1,length(n));
for i=1:5
	x(i+N) = i;
end
subplot(211)
stem(n,x);
subplot(212)
stem(n,circshift(x,3));

3.3 序列内插。

clear;close all;clc;
N = 15;
n = -N : 1 : N;
m = -N : 0.5 : N;
x = zeros(1,length(n));
for i=1:3
	x(i+N) = i;
end
subplot(211)
stem(n,x);
subplot(212)
y = interp1(n,x,m,'nearst'); % linear,spline
stem(m,y);

3.4 序列抽取。

clear;close all;clc;
N = 15;
M = 3;
n = -N : N;
m = -N : M : N;
x = zeros(1,length(n));
for i=1:10
	x(i+N) = i;
end
subplot(211)
stem(n,x);
subplot(212)
y = downsample(x,M);
stem(m,y);

3.5 序列卷积。

clear;close all;clc;
N = 15;
n = -N : N;
m = -2*N : 2*N;
x = zeros(1,length(n));
y = zeros(1,length(n));
for i=1:3
	x(i+N) = i;
	y(i+N) = 1;
end
subplot(311)
stem(n,x);
subplot(312)
stem(n,y);
subplot(313)
stem(m,conv(x,y));
axis([-15,15,0,10]);

3.6 序列自相关。

clear;close all;clc;
N = 15;
n = -N : N;
x = zeros(1,length(n));
for i=1:3
	x(i+N) = i;
end
subplot(211)
stem(n,x);
subplot(212)
[a,b] = xcorr(x,'unbiased');
stem(b,a)
axis([-15,15,0,1]);

3.7 序列互相关。

互相关可以从一个混合信号中，判断其是否包含某种信号。如：从一段音乐中，判断其是否包含1kHz的正弦信号。

clear;close all;clc;
N = 15;
n = -N : N;
m = -2*N : 2*N;
x = zeros(1,length(n));
y = zeros(1,length(n));
for i=1:3
	x(i+N) = i;
	y(i+N) = 1;
end
subplot(311)
stem(n,x);
subplot(312)
stem(n,y);
[a,b] = xcorr(x,y,'unbiased');
subplot(313)
stem(b,a)
axis([-15,15,0,1]);

四、离散系统的分类。

4.1 线性与非线性系统。

系统同时具有均匀特性与叠加特性，才是线性系统。

4.2 时变与非时变系统。

非时变系统：对于固定的输入，其输出也是固定的，不与输入的作用时刻有关。

4.3 因果与非因果系统。

有输入然后才有输出的，是因果系统。输入不能超前于输出，有因才有果。

4.4 稳定与不稳定系统。

对于任意有界输入，其输出也是有界的，是稳定系统。

五、工程应用。

输入信号x，作用于系统H，得到输出信号y。那么，

已知x和H，求y。是分析。用于分析系统的特性。

已知y和H，求x。是控制。用于闭环控制系统，可以参考《自动控制原理》。

已知x和y，求H。是设计。用于设计一种满足需求的系统，如：数字滤波器。