基于matlab时域采样和频域采样验证毕业设计.doc

基于matlab时域采样和频域采样验证毕业设计（论文） 提示频域采样用以下方法容易变程序实现 ① 直接调用MATLAB函数fft计算就得到在的32点频率域采样 ② 抽取的偶数点即可得到在的16点频率域采样即 当然也可以按照频域采样理论先将信号x n 以16为周期进行周期延拓取其主值区16点再对其进行16点DFT FFT 得到的就是在的16点频率域采样 四思考题 如果序列x n 的长度为M希望得到其频谱在上的N点等间隔采样当N M时 如何用一次最少点数的DFT得到该频谱采样 五实验报告及要求 1 运行程序打印要求显示的图形 2 分析比较实验结果简述由实验得到的主要结论 　 3 简要回答思考题 4 附上程序清单和有关曲线 六程序清单和信号波形 1时域采样理论的验证 程序清单 时域采样理论验证程序 Tp 641000 观察时间Tp 64微秒 产生M长采样序列x n Fs 1000T 1Fs Fs 1000T 1Fs M TpFsn 0M-1 f nFsM A 444128alph pi50205omega pi50205 xn Aexp -alphnT sin omeganT Xk Tfft xnM M点FFT[xnt ] subplot 311 plot fabs Xk xlabel fHz ylabel x1 jf title x1 n 的幅度特性 Fs 300Hz Tp 641000 观察时间Tp 64微秒 产生M长采样序列x n Fs 1000T 1Fs Fs 300T 1Fs M TpFsn 0M-1 f nFsM A 444128alph pi50205omega pi50205 xn Aexp -alphnT sin omeganT Xk Tfft xnM M点FFT[xnt ] subplot 311 plot fabs Xk xlabel fHz ylabel x1 jf title x1 n 的幅度特性 Fs 200Hz Tp 641000 Fs 200T 1Fs M TpFsn 0M-1 A 444128alph pi50205omega pi50205 xnt Aexp -alphnT sin omeganT Xk Tfft xntM yn xa nT subplot 325 tstem xntyn box on title a Fs 1000Hz k 0M-1fk kTp subplot 326 plot fkabs Xk title a TFT[xa nT ]Fs 1000Hz xlabel f Hz ylabel 幅度 axis [0Fs012 abs Xk ] 信号波形 频域采样理论的验证 程序清单 M 27N 32n 0M 产生M长三角波序列x n xa 0floor M2 xb ceil M2 -1-10 xn [xaxb] Xk fft xn1024 1024点FFT[x n ] 用于近似序列x n 的TF X32k fft xn32 32点FFT[x n ] x32n ifft X32k 32点IFFT[X32 k ]得到x32 n X16k X32k 12N 隔点抽取X32k得到X16 K x16n ifft X16kN2 16点IFFT[X16 k ]得到x16 n subplot 322 stem nxn box on title b 三角波序列x n xlabel n ylabel x n axis [032020] k 01023wk 2k1024 subplot 321 plot wkabs Xk title a FT[x n ] xlabel \omega\pi ylabel X ej\omega axis [010200] k 0N2-1 subplot 323 stem kabs X16k box on title c 16点频域采样 xlabel k ylabel X_1_6 k axis [080200] n1 0N2-1 subplot 324 stem n1x16n box on title d 16点IDFT[X_1_6 k ] xlabel n ylabel x_1_6 n axis [032020] k 0N-1 subplot