随笔分类 - MATLAB
摘要：代码： 运行结果： 第1、2小题的图这里不放了。 相位延迟phase-delay为0.01时对应的a 的值0.9802 此时1阶全通系统的留数、极点为 系统零极点图 该系统部分分式展开后，求逆z变换得脉冲响应 由下图知，两种方法得到的系统脉冲响应h的幅度谱、相位谱、群延迟大致类似（ω=π时不同）。
阅读全文
摘要：代码： 运行结果： D=1，单个滤波器 这里取D=4，单个重复4次，系统函数部分分式展开， 第2、3小题不会。
阅读全文
摘要：代码： 运行结果： 增益系数K 单个iir高通滤波器，零极点图 L=6阶梳状，系统函数部分分式展开 零极点的模和幅角 直接形式分子、分母系数直接求幅度谱、相位谱和群延迟 依据梳状滤波器系统函数部分分式展开，求逆z变换得到脉冲响应序列，由其求幅度谱、相位谱和群延迟，如下 直接形式求法和逆z变换求法，幅
阅读全文
摘要：代码： 运行结果： 先计算单个IIR低通， 零极点 L=6阶梳状低通，系统函数部分分式展开如下 梳状低通滤波器零极点图 幅度谱、相位谱、群延迟 可以看出，在0到2π范围内，单个低通重复出现了6次，原来的谱压缩到六分之一。
阅读全文
摘要：今天是六一儿童节，陪伴不了家人，心里思念着他们，看着地里金黄的麦子，远处的山，高高的天 代码： 运行结果： 陷波滤波器，ω0=0.5π，引入极点r=0.7 系统函数部分分式展开 系统零极点如下图 幅度谱、相位谱、群延迟 零点位于原点位置，相当于去掉零点，如下 去掉零点后，陷波滤波器的幅度谱、相位谱和
阅读全文
摘要：代码： 运行结果： 系统函数部分分式展开，查表求逆z变换就可得到h(n) 零极点的模和幅角 将脉冲序列当成输入得到h_chk(n)，系统函数求逆z变换得到h(n)， 二者幅度谱、相位谱、群延迟对比如下，可见，幅度谱一样，相位谱和群延迟有所不同。
阅读全文
摘要：代码： 运行结果： 这里的G为增益系数，使得幅度谱在共振频率处最大，等于1 系统函数部分分式展开 系统函数的零极点图 h_chk是由将脉冲序列当成系统输入而得到的，h是由系统函数部分分时展开后查表求逆z变换得到的， 二者幅度谱一致，但是相位谱和群延迟稍有不同。 r=0.9和0.99的结果这里就不放了
阅读全文
摘要：代码： 运行结果： 系统函数部分分式展开， 零极点的模和幅角： 用脉冲序列当输入得到脉冲响应序列h_chk(n)，系统函数H(z)取逆z变换得h(n)，二者如下图 h_chk(n)的幅度谱、相位谱、群延迟 h(n)的幅度谱、相位谱、群延迟 r=0.9、0.99的图这里就不放了。
阅读全文
摘要：代码： 运行结果： 运用Parks-McClellen算法，滤波器长度M=27，比较合适 脉冲响应、幅度谱和误差函数 振幅谱 下面是P7.9的结果，可看出，P-M法得到的长度（M=27)比窗函数（M=43）得到的小得多。 第7章终于弄完了，内心小欢喜一下，放一张田地的照片放松放松，明天开始第8章！
阅读全文
摘要：代码： 运行结果： 先用firpmord求出N、f、m、weights，这些参数作为firpm函数参数，再求出脉冲响应序列 求出脉冲响应序列的幅度谱、误差响应函数 振幅谱，按照题目要求画了网格线 依前面经验，Parks-McClellan方法要比窗函数法、频率采样法设计的滤波器长度小，节约硬件资源，
阅读全文
摘要：代码： 运行结果： 窗函数法，使用了矩形窗、三角窗、Hann窗、Hamming窗、Blackman窗、Kaiser窗， 1、Rectangular窗 2、Bartlett三角窗 3、Hann、Hamming窗、Blackman窗的图这里不放了，直接放Kaiser窗的结果 4、频率采样方法 5、PM法
阅读全文
摘要：代码： 运行结果： 30阶时，稍有欠缺 增大滤波器长度，M=35时， 幅度谱、误差函数如下 振幅谱如下
阅读全文
摘要：代码： 运行结果： 用P-M法设计出来滤波器长度M=23 幅度谱如下 振幅谱如下： 第2个输入，及对应的输出 第2个输入输出的谱（DTFT）对比： 从图中看，输入0.5π分量幅度有所增大，几乎达到原来的2倍。
阅读全文
摘要：代码： 运行结果： 设计一个50阶（即长度M=51）的线性相位FIR，通带宽度不超过0.02π，阻带衰减达到30dB， 最后要把输入中的高斯噪声过滤掉。 As=33dB，满足设计要求。 用P-M方法设计的脉冲响应幅度谱， 振幅谱 输入信号 滤波前后，输入数出对比 输入输出的谱： 右下图可见，随即噪声
阅读全文
摘要：参照Example7.27，因为0.1π=2πf1 f1=0.05，0.9π=2πf2 f2=0.45 所以0.1π≤ω≤0.9π，0.05≤|H|≤0.45 代码： 运行结果： 根据线性相位FIR性质，differentiator为第3类线性相位FIR，下图为脉冲响应、幅度谱和振幅谱。 脉冲响应和
阅读全文
摘要：代码： 运行结果： 滤波器长度M=47，阻带衰减满足设计指标。 幅度谱和相位谱 振幅谱，把阻带、通带放大，数数极值点的个数。 下图，9个极值点 下图，8个极值点 下图，9个极值点 总共有9+8+9=26个极值点，M=47，L=(M-1)/2=23，0到π上，最多L+3=26个极值点。 输入输出序列
阅读全文
摘要：代码： 运行结果： 阶数N=22时，阻带衰减43dB，不满足要求。 所以增大N，当N=24时，As=52dB，滤波器长度M=25
阅读全文
摘要：又是一年五一节，朋友圈都是晒名山大川的，晒脑袋的，我这没钱的待在家里上网转转吧 频率采样法设计带通滤波器，过渡带中有一个样点 代码： 运行结果： 频率采样法，得到的脉冲响应、振幅响应和幅度响应： 幅度响应、相位响应和群延迟 采用fir2函数求脉冲响应序列 matlab自带的fir2函数比频率采样法得
阅读全文
浙公网安备 33010602011771号