信号处理工具箱.pptVIP

利用在MATLAB设计IIR数字滤波器可分以下几步来实现 (1)按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标； (2)根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc； (3)利用最小阶数N产生模拟低通滤波原型； (4)利用截止频率Wc把模拟低通滤波器原型转换成模拟低通、高通、带通或带阻滤波器； (5)利用冲激响应不变法或双线性不变法把模拟滤波器转换成数字滤波器。 例7.24 设计一个在通带内的最大衰减为3 dB，在阻带内的最小衰减为40 dB的4阶低通模拟椭圆滤波器原型。 MATLAB程序如下： n=4; rp=3; rs=40; [z, p, k]=ellipap(n, rp, rs); [b, a]=zp2tf(z, p, k); w=logspace(-1, 1); freqs(b, a) 第7章 信号处理箱 7.1 信号及其表示 7.2 信号的基本运算 7.3 线性时不变系统 7.4 傅里叶(Fourier)变换 7.5 IIR数字滤波器的设计方法 7.6 FIR数字滤波器设计 7.7 信号处理的图形用户界面工具 7.1 信号及其表示 7.1.1工具箱中的信号产生函数 函数名 功能 函数名 功能 sawtooth 产生锯齿波或三角波信号 pulstran 产生冲激串 square 产生方波信号 rectpule 产生非周期的方波信号 sinc 产生sinc函数波形 tripuls 产生非周期的三角波信号 chirp 产生调频余弦信号 diric 产生Dirichlet或周期sinc函数 gauspuls 产生高斯正弦脉冲信号 gmonopuls 产生高斯单脉冲信号 vco 电压控制振荡器 7.1.2离散时间信号的表示 在MATLAB中，离散时间信号x(n)的表示：需用一个向量x表示序列幅值，用另一个等长的定位时间变量n，才能完整地表示一个序列。 [例7-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]) 运行结果如图7.10所示。 图 7.10 离散时间信号图形 7.1.3几种常用离散时间信号的表示 1．单位脉冲序列 直接实现：x=zeros(1,N); x(1,n0)=1; 2．单位阶跃序列 直接实现：n=[ns:nf]; x=[(n-n0)=0]; 3．实指数序列 直接实现：n=[ns:nf]; x=a.^n; 4．复指数序列 直接实现：n=[ns:nf]; x=exp((sigema+jw)*n); 5．正(余)弦序列 直接实现：n=[ns:nf]; x=cos(w*n+sita); 7.2 信号的基本运算 7.2.1信号的相加与相乘 y(n)=x1(n)+x2(n) y(n)=x1(n)×x2(n) MATLAB实现：y=x1+x2; y=x1.*x2 7.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) 7.2.3 序列翻褶与序列累加运算 序列翻褶：y(n)=x(-n)。MATLAB可实现: y=fliplr(x) 序列累加的数学描述为： MATLAB实现：y=cumsum(x) 7.2.4 两序列的卷积运算 两序列卷积运算： MATLAB实现：y=conv(x1,x2)。序列x1(n)和x2(n)必须长度有限。 7.2.5 两序列的相关运算 两序列相关运算： 。MATLAB实现：y=xcorr(x1,x2)。 7.2.6 信号的能量和功率 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; 7.3.1．二次分式模型 连续系统： 离散系统： 7.3 线性时不变系统 在MATLAB中，用sos、ss、tf、zp分别表示二次分式模型、状态空间模型、传递函数模型和零－极点增益模型。其中sos