综合测控实验Simulink-2015重点解析.ppt

卷积的计算 由于滤波器为线性时不变因果系统，当输入是单位脉冲δ[n-k],则输出是h[n-k], 可见，系统的输出取决于输入和系统的脉冲响应h[n], 卷积事实上等于序列x[k]和和h[k]的反转位移序列相乘, 卷积的计算：逐点相乘后相加,这里要理解h[n-k]的含义, h[n-k]是h[k]关于纵轴的镜像； h[n-k]是h[-k]右移n个采样点所得的值。 离散系统的时域模型——差分方程 差分方程是递归定义的，函数f（t）在t=kT处第i阶差分的差分就是其第i+1阶差分，如：上述函数的二阶前向差分方程定义为： 离散系统中信号是以序列的形式出现的，因此只能用差分方程来分析。设f（t）为t的函数，在t=kT处取值为f（kT），此时f（t）在t=kT处的一阶差分方程定义为： 前向差分 后向差分 差分方程与卷积 差分方程与卷积均可以计算己知输入情况下的滤波器的输出，二者的结果对于同一个滤波器而言应该是相同的。 Z变换与差分方程求解 Z变换是对离散序列进行的一种数学变换。常用以求线性时不变差分方程的解。　 离散时间序列 x(n) 的Z变换: 逆Z变换：由Y(z)计算y[n] y[n] = Z ?1{Y ( z )} z-1延迟因子,不能把单步延迟项z-1和逆Z变换的符号Z ?1混淆。 ①线性性质： ②时延性质:延迟定理 用Z传输函数、脉冲响应、差分方程可以描述同一个数字系统。 数字滤波器设计 数字滤波器 经典滤波器 根据滤波特性来分：低通；高通；带通；带阻 (理想滤波器是不可实现的,只能在一定程度上去逼近) 从单位脉冲响应长度来分：IIR-DF；FIR-DF 现代滤波器 寻找在某种准则下的最优解 维纳滤波器，卡尔曼滤波器，自适应滤波器 从对信号处理的作用来分：选频；其他 1. ωp和ωs分别称为通带边界频率和阻带截止频率。 2. 从ωp到ωs称为过渡带，过渡带上的频响一般是单调下降的 3. 通带频率范围为0≤|ω|≤ωp，在通带中要求(1－δ1)|H(ejω)|≤1，阻带频率范围为ωs≤|ω|≤π，在阻带中要求|H(ejω)|≤δ2。 1 1+ ?p 1- ?p ?s ?p ?s ? 显然，?p 越小, 通带波纹越小，通带逼近误差就越小； ?s越大, 阻带波纹越小，阻带逼近误差就越小； ωp与ωs间距越小, 过渡带就越窄。所以低通滤波器的设计指标完全由通带边界频率ωp 、通带最大衰减?p 阻带边界频率ωs和阻带最小衰减?s确定。  1 1/A ? ?p ?s ? 归一化的数字低通滤波器幅度响应指标 FIR 滤波器可以设计为线性相位，并且总是稳定的。 在多数情况下，FIR滤波器的阶数 NFIR 显著大于具有等效幅度响应的IIR滤波器阶数NIIR 。NFIR/NIIR 通常为10的量级或更高. IIR 滤波器通常计算更简便。 在很多应用中，并不要求滤波器具有严格的线性相位，在这些情况下，通常会因计算简便而选择IIR滤波器。 求解得: 滤波器的特性由3db截止频率和阶数N确定. 附加条件1过渡比:通常用k表示,它反应过渡带的性能,过渡带越窄,k值趋近于1,对于低通滤波器k1. 附加条件2偏离参数:通常用k1表示, 设计要求 只要求简单的iir滤波器,有低通、高通等选项。 Fs=5120;dt=1.0/Fs;N=1024;T=dt*N; t=linspace(0,T,N); y=sin(2*3.14*100*t)+ 1/3*sin(9*2*3.14*100*t);plot(t,y); fp=280;fs=450; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标（低通滤波器的通、阻带边界频） [N,Wc]=buttord(wp,ws,rp,rs);%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc [b,a]=butter(N,Wc); %设计Butterworth低通滤波器 %[N,Wc]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp %[b,a]=ellip(N,rp,rs, Wc); %调用ellip计算椭圆带通DF系统函数系数向量b和a %[b,a]= =butter(N,Wc,‘high’); %设计Butterworth高通滤波器 figure(2); % 打开窗口2 freqz(b,a); %求数字低通滤波器的频率响应 %[h,f]=freqz(b,a); %求数字低通滤波器的频率响应 figure(3); % 打开窗口3 y1t=filter(b,a,y); %滤波器软件实现，应用差分方程 % conv()是卷积运算，