第4章MAT六B在信号处理中的应用研究报告.pptVIP

第 4 章 MATLAB在信号处理中的应用 第4章 MATLAB在信号处理中的应用 4.1 信号及其表示 4.2 信号的基本运算 4.3 信号的能量和功率 4.4 线性时不变系统 4.5 线性时不变系统的响应 4.6 线性时不变系统的频率响应 4.7 傅里叶(Fourier)变换 4.8 IIR数字滤波器的设计方法 4.9 FIR数字滤波器设计 4.1 信号及其表示 4.1.1连续时间信号的表示 连续时间信号：时间变化连续。如y=x(t) 离散时间信号(序列)：时间离散，如x(nT)=x(t)|t=nT. 4.1.2工具箱中的信号产生函数 函数名 功能 函数名 功能 sawtooth 产生锯齿波或三角波信号 pulstran 产生冲激串 square 产生方波信号 rectpule 产生非周期的方波信号 sinc 产生sinc函数波形 tripuls 产生非周期的三角波信号 chirp 产生调频余弦信号 diric 产生Dirichlet或周期sinc函数 gauspuls 产生高斯正弦脉冲信号 gmonopuls 产生高斯单脉冲信号 vco 电压控制振荡器 4.1.3离散时间信号的表示 在MATLAB中，离散时间信号x(n)的表示：需用一个向量x表示序列幅值，用另一个等长的定位时间变量n，才能完整地表示一个序列。 [例4-10] 绘制离散时间信号的棒状图。其中x(-1)=-1, x(0)=1, x(1)=2, x(2)=1, x(3)=0, x(4)=-1。MATLAB源程序为： n=-3:5; %定位时间变量 x=[0,0,-1,1,2,1,-1,0,0]; stem(n,x); grid; % 绘制棒状图 line([-3,5],[0,0]); %画x轴线 xlabel(n); ylabel(x[n]) 运行结果如图4.11所示。 图 4.11 离散时间信号图形 4.1.4几种常用离散时间信号的表示 1．单位脉冲序列 直接实现：x=zeros(1,N); x(1,n0)=1; 2．单位阶跃序列 直接实现：n=[ns:nf]; x=[(n-n0)=0]; 4.2 信号的基本运算 4.2.1信号的相加与相乘 y(n)=x1(n)+x2(n) y(n)=x1(n)×x2(n) MATLAB实现：y=x1+x2; y=x1.*x2 4.2.2序列移位与周期延拓运算 序列移位：y(n)=x(n-m)。MATLAB实现：y=x; ny=nx-m 序列周期延拓：y(n)=x((n))M，MATLAB实现：ny=nxs:nxf；y=x(mod(ny,M)+1) 4.2.3 序列翻褶与序列累加运算 序列翻褶：y(n)=x(-n)。MATLAB可实现: y=fliplr(x) 序列累加的数学描述为： MATLAB实现：y=cumsum(x) 4.2.4 两序列的卷积运算 两序列卷积运算： MATLAB实现：y=conv(x1,x2)。序列x1(n)和x2(n)必须长度有限。 4.2.5 两序列的相关运算 两序列相关运算： 。MATLAB实现：y=xcorr(x1,x2)。 4.3 信号的能量和功率 1.信号能量 数字定义： MATLAB实现: E=sum(x.*conj(x)); 或 E=sum(abs(x).^2); 数字定义： 2. 信号功率 MATLAB实现: P=sum(x.*conj(x))/N; 或 E=sum(abs(x).^2)/N; 4.4 线性时不变系统 4.4.1 系统的描述 1．常系数线性微分/差分方程 2．系统传递函数 3．零－极点增益模型 连续系统： 连续系统： 离散系统： 离散系统： 4．极点留数模型 离散系统： 连续系统： 5．二次分式模型 连续系统： 离散系统： 6．状态空间模型 连续系统： 离散系统： 4.4.2 系统模型的转换函数 在MATLAB中，用sos、ss、tf、zp分别表示二次分式模型、状态空间模型、传递函数模型和零－极点增益模型。其中sos表示二次分式，g为比例系数，sos为L×6的矩阵，即 （4－15） 1．ss2tf函数 格式：[num, den]=ss2tf(A,B,C,D,iu) 功能：将指定输入量iu的线性系统(A,B,C,D)转换为传递函数模型[num,den]。 2．zp2tf函数 格式：[num,den]=zp2tf(z,p,k) 功能：将给定系统的零－极点增益模型转换为传递函数模型，z、p、k分别为零点列向量、极点列向量和增益系数。