时域采样频谱分析.docVIP

数字信号处理研讨 关键词： 采样定理运用 不同采样频率下的信号幅频特性 matlab分析信号幅频特性 FFT变换运用求幅频特性 1、根据离散时间信号与模拟信号的时域、频域关系，求序列的频谱。 设： 求出它的傅里叶变换，画出幅频特性曲线； 分别以采样频率和对它进行采样得到序列，画出的幅频特性曲线和的幅频特性曲线，根据所得结果分析同一个模拟信号以不同采样频率采样后所得序列频谱的变化情况。 解答 采样定理知识 采样定理是DSP技术的一个里程碑，也称为香农采样定理（Shannon sampling theorem），或奈奎斯特采样定理（Nyquist sampling theorem）。采样定理证明了采样频率必须为信号频率的两倍以上，采样数据才可以正确的重构出原始的模拟信号。 奈奎斯特频率（Nyquist frequency） 采样频率的一半 模拟信号中高于奈奎斯特频率以上的信息是无法包含在数字信号中。但是这些频率信息会混淆到数字信号所能表示的信息中。 任何高于奈奎斯特频率的频率都会对应到0到奈奎斯特频率之间 具体知识见： /view/b1c8f42a4b73f242336c5fbd.html?from=relatedhasrec=1 1）画出原信号幅频特性 Matlab代码： clc,clear Fs = 5000; % Sampling frequency T = 1/Fs; % Sample time L = 1000;% Length of signal t = (0:L-1)*T;% Time vector % Sum of a 50 Hz sinusoid and a 120 Hz sinusoid x=sinc(1000*t); subplot(211) plot(Fs*t(1:50),x(1:50)) title(x=sin(1000*pi*t)./(1000*pi*t)) xlabel(time) NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(x,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2+1); %Plotsingle-sidedamplitude spectrum. k=abs(Y(1:NFFT/2+1)); subplot(212) plot(6*f,k) title(幅频特性曲线) xlabel(频率 (Hz)) ylabel(|X(jw)|) 2）做离散信号分析并画出频谱 Matlab代码： clear all; T1=1/2000; n=0:64/(1000*T1); xn1=sinc(1000*n*T1); Xk1=fft(xn1,1024); subplot(221);stem(n,xn1,.); grid on xlabel(n);ylabel(x(n));title(2000Hz采样 x(n)); k=0:1023;wk=2000*pi*2*k/1024; subplot(222);plot(wk,abs(Xk1));grid on; xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|);title(FT[x(n)]); T2=1/5000; n=0:64/(1000*T2); xn2=sinc(1000*n*T2); Xk2=fft(xn2,1024); subplot(223);stem(n,xn2,.); grid on xlabel(n);ylabel(x(n));title(5000Hz采样 x(n)); k=0:1023;wk=2000*pi*2*k/1024; subplot(224);plot(wk,abs(Xk2));grid on; xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|);title(FT[x(n)]); 2、设： 求出它的傅里叶变换，画出幅频特性曲线； 分别以采样频率和对它进行采样得到序列，画出的幅频特性曲线和的幅频特性曲线，根据所得结果阐述数字角频率与模拟频率的关系。 参考第一题的解答过程 结论 同一模拟信号下的抽样连续时间序列和离散时间序列的幅频特性是一样的，只决定傅里叶变换的频谱幅度和频谱周期延拓的间隔。 数字角频率和模拟角频率的关系如下： ，分析得知同一信号模拟角频率不变，抽样频率决定了数字频率的变化。