MATLAB与信号处理_3.pptVIP

实验三 信号变换及其MATLAB实现 一、实验目的 1.练习利用matlab 6.5中的命令或程序分析信号频谱。 2.进一步练习M-File的建立、保存、调用。 3.加深对几种傅里叶变换的理解。 二、实验原理与方法 连续时间信号傅里叶变换（FT）： 离散时间信号傅里叶变换（DTFT）： 离散傅里叶变换（DFT）： 快速傅里叶变换（FFT）： 三、实验内容及步骤 数学函数：1、 mod ；2、 length ；3、 sin 其他函数：fft 相关函数 复杂语句的理解 －?自己写出完成复杂功能的语句 要求 1.连续时间信号傅里叶变换分析： 例：分析信号x(t)=u(t)-u(t-1)的频谱。 r=0.005; t=-2:r:2; rt=[t0]-[t-10]; k=-10:0.005:10;w=pi*k; Rjw=rt*exp(-j *t‘ *w)*r; %（给出这之前语句的解释） R=abs(Rjw); ANGR=angle(Rjw); subplot(3,1,1);plot(t,rt); text(2.2,0,t); ylabel(r(t)=u(t)-u(t-1)) ; title(Analog Signal) ; subplot(3,1,2);plot(w/pi,R) ;text(10.5,0,w); ylabel(|Rt(jw)|) ; title(Continuous-time Fourier Transform-Magnitude Part) ; subplot(3,1,3);plot(w/pi,ANGR) ; text(10.5,-4,w); ylabel(ANGRt(jw)) ; title(Continuous-time Fourier Transform-Angle Part) ; 2.离散时间信号傅里叶变换分析： 例：分析信号x(n)=u(n)-u(n-n0)的频谱。n0=4或6或8。 n0=4; % 调整n0=6或8或其他值，分析并比较最后结果 N=256; n=0:N-1; k=-2:0.001:2;w=pi*k; %w=[-2pi,2pi] x=[nn0]; X = x * exp(-j*n*w); magX = abs(X); phaX = angle(X); subplot(2,2,1);plot(w/pi,magX); grid xlabel(Frequency in pi units); ylabel(|X(ejw)|) ; title(Discrete-time Fourier Transform---Magnitude) ; subplot(2,2,3);plot(w/pi,phaX/pi);grid on; xlabel(Frequency in pi units); ylabel(radians/pi) ; title(Discrete-time Fourier Transform---Angle Part) ; 3.离散傅里叶变换分析： 例：分析信号x(n)=u(n)-u(n-n0)的频谱。n0=4或6或8。 n0=4; % 调整n0=6或8或其他值,分析并比较最后结果 N=16; % 调整N=16或32或其他值,分析并比较最后结果 n=0:N-1; k=0:N-1; w=2*pi/N*k; x=[nn0]; X = x * exp(-j*n*w); magX = abs(X); phaX = angle(X); subplot(2,2,2);stem(k,magX, b); grid xlabel(Frequency in pi units); ylabel(|X(k)|) ; title(Discrete Fourier Transform---Magnitude) ; subplot(2,2,4);stem(k,phaX/pi, r);grid on; xlabel(Frequency in k=pi/N units); ylabel(radians/pi) ; title(Discrete Fourier Transform---Angle Part) ; function [Xk] = dft(xn,N) % Computes Discrete Fourier Transform % ----------------------------------- % [Xk] = dft(xn,N) % Xk = DFT coeff. array over 0 = k = N-1 % xn = N-point finite-duration sequence % N = Length of DFT % n = [0:1:N-1];