哈工程信号处理实验三

信号处理实验三信号处理实验三 实验要求 研究抽样过程 分析产生混叠效应的原因 实现不同的重建方案 3 3 1 n 0 80 f0 300 fs 8000 x sin 2 pi f0 fs n pi 3 stem n x gtext f0 300 fs 8k plot n x gtext f0 300 fs 8k gtext Continue b n 0 80 fs 8000 subplot 221 x sin 2 pi 100 fs n pi 3 stem n x grid gtext f0 100 subplot 222 x sin 2 pi 225 fs n pi 3 stem n x grid gtext f0 225 subplot 223 stem n x grid gtext f0 350 subplot 224 x sin 2 pi 225 fs n pi 3 stem n x grid gtext f0 475 d n 0 80 fs 8000 subplot 221 x sin 2 pi 7525 fs n pi 3 stem n x grid gtext f0 7525 subplot 222 x sin 2 pi 7650 fs n pi 3 stem n x grid gtext f0 7650 subplot 223 x sin 2 pi 7775 fs n pi 3 stem n x grid gtext f0 7675 subplot 224 x sin 2 pi 7900 fs n pi 3 stem n x grid gtext f0 7900 e n 0 80 fs 8000 subplot 221 x sin 2 pi 32100 fs n pi 3 stem n x grid gtext f0 32100 subplot 222 x sin 2 pi 32225 fs n pi 3 stem n x grid gtext f0 32225 subplot 223 x sin 2 pi 32350 fs n pi 3 stem n x grid gtext f0 32350 subplot 224 x sin 2 pi 32475 fs n pi 3 stem n x grid gtext f0 32475 t 0 1 80000 1000 80000 x cos 2 pi 300 t pi 6 stem t x t 0 1 8000 1000 80000 x cos 2 pi 300 t pi 6 stem t x fmagplot x 1 80000 可以看到 当频率变化时 正弦信号并不是严格按信号频率变化的 而是呈周 期的 有时变密有时变疏 3 3 3 t 0 1 80000 1000 80000 x cos 2 pi 300 t pi 6 stem t x t 0 1 8000 1000 80000 x cos 2 pi 300 t pi 6 stem t x fmagplot x 1 80000 3 3 4 n 0 100 t n 80000 x cos 2 pi t 300 pi 6 stem t 10 x X W dtft x 900 plot W 10 pi abs X 3 3 5 fs 8000 fsim 80000 fcut 2 fs 2 fsim b a cheby2 9 60 fcut 调用雪比滤波器 c d freqz b a 1000 whole 1000 点的频率响应 d 501 1000 d 501 1000 2 pi d 的位置全部左移 2pi plot d 2 pi 1000 fsim abs c 3 3 6 fs 8000 fsim 80000 n 0 1000 l length n x cos 2 pi 300 n fsim pi 6 y x 1 80000 8000 l L length y 取长度避免后面的维度不相等 y2 zeros 1 length x 置零序列 for i 0 100 y2 i 10 1 y i 1 取 10 倍数中的数 相当于抽样置零 end t 0 0 001 1 plot t y2 fmagplot y2 0 0000125 3 3 7 fs 8000 fsim 80000 k fsim fs n 0 1000 l length n x cos 2 pi 2000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 2k fs 8000 fsim 80000 k fsim fs n 0 1000 l length n x cos 2 pi 6000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 6k x cos 2 pi 7000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 7k x cos 2 pi 9000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 9k x cos 2 pi 10000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 10k x cos 2 pi 15000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 15k x cos 2 pi 17000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 17k x cos 2 pi 18000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 18k x cos 2 pi 19000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for i 0 100 y2 i 10 1 y i 1 end t 0 0 001 1 y3 filter b a y2 subplot 211 fmagplot y3 0 0000125 subplot 212 plot t y3 gtext 19k x cos 2 pi 20000 n fsim pi 6 y x 1 k l L length y y2 zeros 1 length x for