matlab程序设计 课件 第9章.pptVIP

9.3 MATLAB在信号处理中的应用 9.3.1 离散傅里叶变换 如果离散序列{x(n),n= -∞,…,0,1,…,∞}是绝对可加的，即 则其离散时间傅里叶变换为 （9.22） 对应的离散时间傅里叶反变换为 （9.23） 9.3.1 离散傅里叶变换 有限长离散序列{x(n),n=0,1,…,N-1}的傅里叶变换称为离散傅里叶变换（DFT），其正逆变换形式为 (9.24) (9.25) 其中， (9.26) 9.3.1 离散傅里叶变换 对于求取离散信号连续频谱的MATLAB实现，可以从离散时间傅里叶变换的定义（9.22）式出发。设置一系列较密的频率wi，求出一系列X(wi)=x*exp(-j*n*wi)，其中，x和n是等长的行向量。于是，x与exp(-j*n*wi)的乘积就是逐项相乘后的连加演算，得出一个标量X(wi)，即该频点上的频率响应。把不同的wi写成一个行向量w，代入上式，即可分别得出相应wi处的响应，得出X(w)=x*exp(-j*n*w)。 【例9.30】 求离散信号的连续频谱。 取一个周期的正弦信号 ，完成下列计算： （1）作N点采样，求它的连续频谱； （2）对该信号进行K个周期延拓，求其连续频谱。 编写M脚本文件exam9_30如下： 9.3.1 离散傅里叶变换 % exam9 _30 clear disp(取一个周期的正弦信号x(t)=sin(2*pi*t)，对其进行N点采样) N=input(N=); % 从键盘输入采样点数 t=[1:N]/N; xt=sin(2*pi*t); % 对一个周期的正弦信号进行N点采样 dt=2*pi/N; % 步长 w=linspace(0,2*pi,1000)/dt; % 设置频率 Xw=xt*exp(-j*[1:N]*w)*dt; % 根据DFT定义求得频率响应X(w) % 周期延拓 disp(重复K个周期的N点信号离散傅里叶变换) K=input(K=); % 输入延拓周期 xt1=reshape(xt*ones(1,K),1,K*N); % 延拓后的时域信号 Xw1=xt1*exp(-j*[1:length(xt1)]*w)*dt; % 延拓后的频率信号 % 显示图形 subplot(2,1,1),plot(w,abs(Xw),k),grid title(一个周期的N点正弦信号x(t)=sin(2*pi*t)的频率响应), xlabel(w),ylabel(X(w)) subplot(2,1,2),plot(w,abs(Xw1),k),grid title(周期延拓后的频率响应), xlabel(w),ylabel(X(w)) exam9_30的运行结果和图9.35 9.3.1 离散傅里叶变换 9.3.1 离散傅里叶变换 DFT实际上是对有限长序列频谱的离散化，使时域有限长序列与频域有限长序列相对应，从而可在频率域对信号进行处理。可以把(9.24)式写成矩阵形式