(6.5.1)--6利用MATLAB进行离散信号频域分析.pdfVIP

第 5 章 频域篇(二) 利用MATLAB进行离散信号频域分析 Signals and Systems 利用MATLAB进行离散信号频域分析 计算离散时间周期信号的频谱 计算离散时间周期信号的频谱  对于基波周期为N 的的信号，MATLAB 内部有两个非常高效的 计算函数：fft()和ifft()  若x是一个包含在单一周期 0  N  N 1 内的x[n ] 的N点向量， 那么x [n ] 的离散傅里叶级数就能用 a (1/ N ) fft(x) 计算出 a 0  k  N 1 a 这里N点向量 包含了 上的 k a  反之，若已知向量 是一个包含在单一周期 0  k  N 1 内的 a x N ifft(a) 傅里叶级数系数 ，则就能用 计算出时域序列 k 2 利用MATLAB进行离散信号频域分析 计算离散时间周期信号的频谱 计算离散时间周期信号的频谱  函数fft()和ifft()调用格式如下 a=fft(x，N) x=ifft(a，N)  x=ifft(a，N)_采样fft算法计算序列向量x 的N点离散傅里叶变换 N缺省时函数自打按x 的长度计算  x=ifft(a，N)_计算向量a的N点离散傅里叶反变换 3 利用MATLAB进行离散信号频域分析 计算离散时间周期信号的频谱 【例】假设是基波周期为32 的对称矩形脉冲序列。该信号在一个 1, n 0,1  周期内可以表示为 x[n ]  0, 其余n  命令行窗口 fx N=32; x=[1 1 zeros(1,30)]; a=(1/N)*fft(x); n=0:N-1; figure(1); subplot(3,2,3),stem(n,real(a),.); xlabel(k) ;ylabel( Re(ak)); subplot(3,2,4),stem(n,imag(a),.); xlabel(k) ;ylabel(Im(ak)); 4 利用MATLAB进行离散信号频域分析 计算离散时间周期信号的频谱 【例】假设是基波周期为32 的对称矩形脉冲序列。该信号在一个 1, n 0,1  周期内可以表示为 x[n ]  0, 其余n  命令行窗口 fx subplot(3,2,5),stem(n,abs(a),.); xlabel(k) ;ylabel(|ak|); title(幅频特性); subplot(3,2,6),stem(n,angle(a),. ); xlabel(k) ;ylabel(ak); title(相频特性); x1=N*ifft(a);