MATLAB课程PPT第八章.pptVIP

8.7.1 模拟滤波器原型设计 例1：设计模拟原型 Butterworth低通滤波器： Ωp=0.2 π Rp=7dB Ωs=0.3 π As=16dB 例2：设计ChebyshevⅠ型模拟原型低通滤波器： Ωp=0.2 π Rp=1dB Ωs=0.3 π As=16dB 8.7.2 模拟到数字滤波器变换 1、冲激不变法： 函数impinvar：冲激响应不变法实现模拟到数字的滤波器变换。 [bz,az]=impinvar(b,a,Fs) [bz,az]=impinvar(b,a) 扩展函数imp_invr：冲激响应不变法模拟到数字滤波器变换 [b,a]=imp_invr(c,d,T) 2、双线性变换: MATLAB中用函数bilinear实现双线性变换 * * 8.1.1 基本信号表示 单位阶跃序列 采用扩展函数 stepseq.m可生成单位阶跃序列 例: 产生序列x(n)并绘制离散图 0≤n≤20 x(n)=u(n-4)-3u(n-20) 3 实指数序列 x(n)=an n=0:N-1; x=a.^n; 8.1.1 基本信号表示 4 复指数序列  x(n)=e(?+jwo)n 5 正余弦序列 x(n)=cos(ω0*n+Θ) 6 随机序列  由函数 rand(1,N) 和 randn(1,N) 生成 例：产生序列x(n), 并绘制其离散数据图。 X(n)=cos(0.4π n)+0.2w(n) 0≤n≤50 7 周期序列  x(n)=x(n+N) 例：产生序列 =[x1 x1 x1 x1]并绘制出离散图 .x1=[1 2 3] 8.1.2 序列操作 信号加 应用扩展函数 sigadd.m 实现信号加 信号乘 应用扩展函数 sigmult.m 实现 改变比例 y=a*x 移位 采用扩展函数 sigshift.m 实现 5 折叠 采用扩展函数 sigfold.m 实现 8.1.2 序列操作 取样和 采用命令 y=sum(x(n1:n2)) 实现 7　取样积 采用命令 y=prod(x(n1:n2)) 实现 信号能量　 Ex=sum(x.*conj(x)) Ex=sum(abs(x).^2) 信号功率 Px=sum(abs(x).^2)/N 10 奇偶综合 扩展函数 evenodd.m可将任一给定的序列分解成偶对称部分xe(n)和奇对部分xo(n) 8.1.2 序列操作 例：已知序列x(n)=u(n)-u(n-10)，要求将它分解成奇偶序列 8.1.3 　卷积 求卷积可直接采用MATLAB中的函数 conv，其调用格式是：y=conv(x,h) 扩展函数 con_m 可求出带下标的序列卷积． 例1：设线性时不变(LTI)系统的冲激响应为：h(n)=(0.9)nu(n),输入序列为x(n)=u(n)-u(n-10)，求系统的输出序列y(n)。（绘制出x、h及y的离散数据图） 例2：设离散系统可由下列差分方程表示： y(n)-y(n-1)+0.9y(n-2)=x(n) (1)计算n=-20,…100上的单位冲激响应； (2)计算n=-20,…100上的单位阶跃响应。 §8.2 序列的傅里叶变换 例1：设x(n)=2n，-5≤n≤5，求相应的X(ejω) 若序列x(n)绝对可和，则x(n)的傅里叶变换为： 解：根据定义 记矢量x=[x(0) x(1)…x(n)]，取 即将[0, π]均匀分成M+1点，记矩阵 矢量 ，则上式可写成X=Wx，若记矢量 n=[-5,-4,…5],k=[0,1,2,…M]，则有 §8.2 序列的傅里叶变换 例2：实序列 求出 X(ejω)的实部和虚部，同时求出奇偶部分相应的Xe(ejω)和Xo(ejω). 即