Something About Homework 2.doc

Something About Homework 2 信院03607班 冯春芳 28号 031878 CONTENT 1.Preface 2.目录 Page 1 3. DFT各种现象的图形 4.混叠现象参数设置的疑问 5.程序运行小问题 模拟原始信号求N点DFT的幅度谱和相位谱： 有混叠现象的信号求N点DFT的幅度谱和相位谱： 有频谱泄漏现象的信号求N点DFT的时域图和频域图： 加长矩形长度改善频漏现象的信号求N点DFT的时域图和频域图： Kaiser函数改善频漏现象的信号求N点DFT的时域图和频域图： 栅栏现象和信号末尾补零改善的信号求N点DFT图： 观察混叠现象时不同参数的图： N=64; n=0:63; t=0.05*n; q=n*2*pi/N; x=cos(3*pi/20*t);%生成信号序列 y=fft(x,N); subplot(3,1,1); plot(t,x); title(source signal);%绘制源曲线信号 subplot(3,1,2); plot(q,abs(y)); %幅度谱 title(magnitude); subplot(3,1,3); plot(q,angle(y)); %相位谱 title(phase); 以上是源程序，当t取以下值时图发生的变化并不是我起先想得那样，有点疑惑， 在FFT算法中，一个信号用间隔Ts进行抽样所得到的N个样本表示，最后得到的结果是频域区间[0,fc]内Xd(f)的长度为N的多个样本。采样频率fs=1/Ts，频率抽样间隔df=fs/N为傅立叶变换的频率分辨率。 对连续信号采样，要求连续信号是带限的，采样频率要足够高。Fs应满足Nyquist采样定理才不产生混叠。 这里t就是采样间隔，那么采样频率fs=1/t2fh=2*0.075=0.15Hz,那么t6.7s,理想情况下幅度谱应是在f=0.075Hz有一个冲激，但根据 信号带宽大致估计：带宽正比于信号持续时间的倒数，为可靠起见，可取信号持续时间倒数的10倍 那么带宽又是另一个值，t的取法就没那么简单了 t=0.01*n时求N点DFT的幅度谱和相位谱： t=0.1*n时求N点DFT的幅度谱和相位谱： t=0.5*n时求N点DFT的幅度谱和相位谱： t=1*n时求N点DFT的幅度谱和相位谱： t=7*n时求N点DFT的幅度谱和相位谱： t=15*n时求N点DFT的幅度谱和相位谱： t=50*n时求N点DFT的幅度谱和相位谱： 随着t的增大，频漏现象越严重，但特别奇怪的是当t=7*n时两个冲激突然很靠近一位更大时会更靠近，但当t=15*n时却又分开了，且冲激更陡了，限于本人有限的知识对其百思不得其解。 另一个问题就是在运行最后的总文件my.m文件时出现下列的错误，上课时问老师，说获得的m文件和gui文件与m文件里的fuction名（function varargout = my(varargin)）必须一致，但是怎么生成m文件后没有gui文件？如上，虽然改fuction名，但仍没有解决问题，不过这并不影响结果。还有就是在m文件中回调时直接把程序粘上是不是比写函数名运行速度要慢？关于以上问题打扰老师敬请赐教！ “??? Error using == struct2handle Error while evaluating axes CreateFcn.”