概述
2.1 基础知识回顾
对线性离散时间系统,若 y 1 [ n ] y_1[n] y1[n]和 y 2 [ n ] y_2[n] y2[n]分别是输入序列 x 1 [ n ] x_1[n] x1[n]和 x 2 [ n ] x_2[n] x2[n]的响应,则输入 x [ n ] = α x 1 [ n ] + β x 2 [ n ] x[n]=αx_1[n]+βx_2[n] x[n]=αx1[n]+βx2[n]的输出响应为 y [ n ] = α y 1 [ n ] + β y 2 [ n ] y[n]=αy_1[n]+βy_2[n] y[n]=αy1[n]+βy2[n]。上式的叠加性质对任意常量α和β以及任意输入 x 1 [ n ] x_1[n] x1[n]和 x 2 [ n ] x_2[n] x2[n]都成立。若存在一组非零的α和β值,或者一组非零的输入序列 x 1 [ n ] x_1[n] x1[n]和 x 2 [ n ] x_2[n] x2[n],上式不成立,则为非线性系统。
对于离散时不变系统,若 y 1 [ n ] y_1[n] y1[n]是 x 1 [ n ] x_1[n] x1[n]的响应,则输入 x [ n ] = x 1 [ n − n 0 ] x[n]=x_1[n-n_0] x[n]=x1[n−n0]的输出响应为 y [ n ] = y 1 [ n − n 0 ] y[n]=y_1[n-n_0] y[n]=y1[n−n0],式中 n 0 n_0 n0是任意整数。上面的输入输出关系对任意输入序列及其响应的输出成立。若对至少一个输入序列及其相应的输出序列不成立,则为时变系统。
若对于任意有界输入序列 x [ n ] x[n] x[n],其输出 y [ n ] y[n] y[n]也是一个有界序列,则该离散时间系统是有界输入有界输出(BIBO)稳定的,即若 ∣ x [ n ] ∣ < B x |x[n]|<B_x ∣x[n]∣<Bx(对所有n值),则相应的输出 y [ n ] y[n] y[n]也有界,即 ∣ y [ n ] ∣ < B y |y[n]|<B_y ∣y[n]∣<By(对所有n值),这里 B x B_x Bx和 B y B_y By都是有限常数。
数字滤波器对单位样本序列 ∣ y [ n ] ∣ |y[n]| ∣y[n]∣的响应称为单位样本响应,简称为冲激响应,用 h [ n ] {h[n]} h[n]表示。相应地,离散时间系统对单位阶跃序列 μ [ n ] {μ[n]} μ[n]的响应用 s [ n ] {s[n]} s[n]表示,称为单位阶跃响应,简称为阶跃响应。
输入信号 x [ n ] x[n] x[n]的冲激响应 h [ n ] h[n] h[n]可用来描述线性时不变离散系统的响应 y [ n ] y[n] y[n],表示为 y [ n ] = ∑ k = − ∞ ∞ h [ k ] x [ n − k ] y[n]=sum_{k=-∞}^∞h[k]x[n-k] y[n]=∑k=−∞∞h[k]x[n−k],通过简单的变量代换可表示成 y [ n ] = ∑ k = − ∞ ∞ h [ n − k ] x [ k ] y[n]=sum_{k=-∞}^∞h[n-k]x[k] y[n]=∑k=−∞∞h[n−k]x[k],上面两式中的和称为序列 x [ n ] x[n] x[n]和 h [ n ] h[n] h[n]的卷积和,表示为 y [ n ] = h [ n ] ∗ x [ n ] y[n]=h[n]*x[n] y[n]=h[n]∗x[n]。
2.2 滑动平均滤波器
滑动平均滤波器: y [ n ] = 1 M ∑ k = 0 M − 1 x [ n − k ] y[n]=frac{1}{M}sum_{k=0}^{M-1}x[n-k] y[n]=M1∑k=0M−1x[n−k]
% 用滑动平均滤波器去除噪声
R = 50;
d = rand(R,1)-0.5;
m = 0:1:R-1;
s = 2*m.*(0.9.^m);
x = s + d';
plot(m,d,'r-',m,s,'b--',m,x,'g:')
xlabel('Time index n'); ylabel('Amplitude')
legend('d[n]','s[n]','x[n]');
pause
M = input('Number of input samples = ');
b = ones(M,1)/M;
y = filter(b,1,x);
plot(m,s,'r-',m,y,'b--')
legend('s[n]','y[n]');
xlabel ('Time index n');ylabel('Amplitude');
最后
以上就是彩色月饼为你收集整理的【MATLAB数字信号处理】离散时间系统的时域分析(一)的全部内容,希望文章能够帮你解决【MATLAB数字信号处理】离散时间系统的时域分析(一)所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复