北京交通大学(数字信号处理研究生课程)chDFS概要.ppt

周期卷积的矩阵表示 例：N=4 Toeplitz 矩阵 定义: Toeplitz 矩阵的构成 1. 写出矩阵的第一列与第一行； 2. 从上至下逐行右移。 Matlab生成Toeplitz 矩阵 T = toeplitz(c,r) c矩阵的第一列 r矩阵的第一行 例： T=toeplitz([1 5 6 7],[1 2 3 4]) T = 1 2 3 4 5 1 2 3 6 5 1 2 7 6 5 1 利用Toeplitz 矩阵计算周期卷积 x=[1 2 3 4]; h=[5 6 7 8]; y=toeplitz(x,[x(1) fliplr(x(2:end))])*h y = 66 68 66 60 DFS的基本性质 4. 周期卷积定理 时域周期卷积定理： 频域周期卷积定理： 时域的周期卷积对应频域的乘积； 时域的乘积对应频域的周期卷积。 MATLAB计算DFS X=fft(x) x=ifft(X) X=fft([2 1 -3 1]) X = 1 5 -3 5 X=fft([0 1 0 -1]) X= 0 - 2.0i 0 2.0i h=fft(fft([1 2 3 4]))/4 h =1.0 4.0 - 0.0i 3.0 2.0+ 0.0i 信号Fourier级数的提出 利用 M点的滑动平均系统去噪 利用 M点的滑动平均系统去噪 % Signal Smoothing by Moving Average Filter N = 201; n = 1.0*rand(1,N)-0.5; k=0:N-1; s=2*k.*(0.9.^k)+2.0*cos(0.02*pi*k)+5.0; x=s+n; M =10; h= ones(M,1)/M; y = filter(h,[1],x); subplot(3,1,1); plot(k,x,b-, k,y,r-); title(x=2*k.*(0.9.^k)+2.0*cos(0.02*pi*k)+5.0+ns[k]); legend(x[k], y[k]); % Display input signal frequency spectrum w=linspace(-1.0*pi,pi,500); X=freqz(x,[1],w); subplot(3,1,2); plot(w/pi,abs(X)); grid on; title(Input signal frequency spectrum X(ejW)); % Display system frequency spectrum H=freqz(h,[1],w); subplot(3,1,3); plot(w/pi,abs(H)); grid on; title(System amplitude frequency spectrum H(ejW)); legend(M点滑动平均滤波(M=10)); 利用 M点的滑动平均系统去噪 江南Style_Psy时域波形与频谱 s 频谱 Hz 时域波形 据说此曲流行的一个原因是它与心脏的跳动韵律相似。 心脏跳动声音时频波形 s 频谱 Hz 时域波形 低频部分有点像！ 我和你_萨拉布莱曼时域波形与频谱 s 频谱 Hz 时域波形 她的声音空灵到了极致，仿佛带你漫步云端。 从频谱看和其他歌手的区别 基音很高，但更高频的分量很少！ 信号传输 为何引入信号的频域分析？ 信号的时域波形 加噪信号的时域波形 信号滤波 为何引入信号的频域分析？ 时间/秒 时间/秒 女生信号时域波形 男生信号时域波形 信号分析 为何引入信号的频域分析？ 电话拨号中的双音多频信号(DTMF) 信号识别 为何引入信号的频域分析？ 数字键 ’1’的波形 数字键 ’2’的波形 离散Fourier级数(DFS) 周期序列的DFS表示 周期序列DFS的推导 周期序列DFS的基本性质 周期序列的卷积 周期序列DFS 表示 周期为N的任意周期序列 其中： DFS系数的计算公式为 可用虚指数序列表示为 DFS系数 的周期为N。 周期序列DFS的推导 周期为N的单位脉冲序列定义为 周期序列DFS的推导 一般地，一个周期为N周期序列可表示为 如何用周期为4的单位脉冲序列表示下图序列? 周期序列DFS的推导 0, k不是N的整数倍 N, k是N的整数倍 周期为N的单位脉冲序列DFS 表示 周期序列DFS的推