首页 > 其他分享 >采样定理

采样定理

时间:2023-04-30 09:11:32浏览次数:37  
标签:采样 infty frac 定理 theta pi omega mathrm

信号\(x(t)\)的频谱为 \(X(\mathrm{j}\omega)\)。
对信号使用周期单位冲激串采样得到采样信号 \(x_p(t)\):

\[x_p(t) = x(t)p(t) \]

其中,\(p(t)\)为采样函数, 是周期为 \(T\)的周期单位冲激串, 并且 \(p(t)\) 的基波频率 \(\omega_s=\frac{2\pi}{T}\)。 采样函数如下:

\[p(t) =\sum_{n=-\infty}^{+\infty}\delta(t-nT) \]

周期为\(T\)的单位冲激串函数\(p(t)\) 的傅里叶变换为强度为 \(\frac{2\pi}{T}\) 周期为 \(\frac{2\pi}{T}\)的周期冲激串:

\[\begin{aligned} P(\mathrm{j}\omega)&=\frac{2\pi}{T}\sum_{k=-\infty}^{+\infty}\delta(\omega-\frac{2\pi k}{T}) \\ &=\frac{2\pi}{T}\sum_{k=-\infty}^{+\infty}\delta(\omega-k\omega_s) \end{aligned} \\ 其中,\omega_s=\frac{2\pi}{T} \]

根据冲激函数的采样性质,\(x_p(t)\) 为一个冲激串,每个冲激的强度等于 \(x(t)\) 在以 \(T\) 为间隔处的样本值。采样过程如下图:
img
采样得到信号为:

\[x_p(t) = \sum_{n=-\infty}^{+\infty}x(nT)\delta(t - nT) \]

根据傅里叶变换的相乘性质:

\[r(t)=s(t)p(t)\longleftrightarrow R(\mathrm{j}\omega)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} S(\mathrm{j}\theta)P(\mathrm{j}(\omega-\theta)) \mathrm{d}\theta \]

可得:

\[x_p(t)=x(t)p(t)\longleftrightarrow X_p(\mathrm{j}\omega)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} X(\mathrm{j}\theta)P(\mathrm{j}(\omega-\theta)) \mathrm{d}\theta \]

因为冲激函数的采样性质, \(X(\mathrm{j}\omega) P(\mathrm{j}(\omega-\theta))\) 的积分只在 \(\omega=\theta\) 的位置不为0,因此根据上式有:

\[\begin{aligned} X_p(\mathrm{j}\omega)&=\frac{1}{2\pi}\int_{-\infty}^{+\infty} X(\mathrm{j}\theta)P(\mathrm{j}(\omega-\theta)) \mathrm{d}\theta \\ &=\frac{1}{2\pi}\int_{-\infty}^{+\infty} X(\mathrm{j}\omega)P(\mathrm{j}(\omega-\theta)) \mathrm{d}\theta \end{aligned} \]

已知周期单位冲激串的傅里叶变换为:

\[P(\mathrm{j}\omega)=\frac{2\pi}{T}\sum_{k=-\infty}^{+\infty}\delta(\omega-k\omega_s) \]

因此, \(X_p(\mathrm{j}\omega)\) 就是 \(X(\mathrm{j}\omega)\) 与上式中每个分量冲激的卷积结果之和。
而信号与一个冲激的卷积结果为信号在冲激位置的位移,即 \(x(t)*\delta(t-T)=x(t-T)\), 那么有:

\[X(\mathrm{j}\omega)*\delta(\omega-k\omega_s)=X(\mathrm{j}(\omega-k\omega_s))\]

结合上面三式,有:

\[X_p(\mathrm{j}\omega)=\frac{1}{T}\sum_{k=-\infty}^{+\infty}X(\mathrm{j}(\omega-k\omega_s) \]

结果是采样得到的信号 \(x_p(t)\) 的频谱是原信号 \((t)\) 的频谱以 \(\omega_s\)为周期重复的函数,即频谱被周期性地复制在间隔为 \(\omega_s\) 的原频谱两侧,同时幅度变为原来的 \(\frac{1}{T}\)。
img
从上图可见,两个重复频谱的间隔为 \(\omega_s-2\omega_M\), 其中 \(\omega_M\) 为原信号频谱的截止频率(带限),当 \(\omega_s>2\omega_M\) 时,复制的频谱没有重叠,此时将采样得到的信号 \(x_p(t)\) 经过一个通带 \(\omega_c\)在 \(\omega_M\)和\(\omega_s-\omega_M\)范围内的低通滤波器,并对幅度加以 \(T\) 倍的增益,那么,输出的信号\(x_r(t)\)将完全还原为 \(x(t)\),即 \(X_r(\mathrm{j}\omega)=X(\mathrm{j}\omega)\), 这就是采样定理
img

标签:采样,infty,frac,定理,theta,pi,omega,mathrm
From: https://www.cnblogs.com/SPDH/p/17364902.html

相关文章

  • 中国剩余定理(CRT)学习笔记
    约定\(A\perpB\)表示\(\gcd(A,B)=1\)。\(A\midB\)表示\(B\equiv0\pmod{A}(A\neq0)\)。引入考虑以下这道题:有物不知其數,三三數之剩二,五五數之剩三,七七數之剩二。問物幾何?——《孫子算經》也就是说,求出下列关于\(x\)方程组的最小整数解:\[\begin{cases}x\equi......
  • 【学习笔记】拓展中国剩余定理
    若干方程组:\(\begin{cases}x\equivc_1\quad(\modp_1)\\x\equivc_2\quad(\modp_2)\\···\\x\equivc_m\quad(\modp_m)\end{cases}\)求x但不保证p互质。采用两两方程合并的形式。\(\begin{cases}x\equivc_1\quad(\modp_1)\\x\equivc_2\quad(\modp_2)\......
  • P1495 【模板】中国剩余定理(CRT)/ 曹冲养猪
      #include<iostream>#include<algorithm>#include<cstring>#include<cmath>usingnamespacestd;#defineintlonglongintn,a[20],M[20],Mi[20];intgcd(inta,intb,int&x,int&y){if(b==0){x=1;y=0......
  • 中国剩余定理(CRT)(待完善)
    中国剩余定理(CRT)求同余方程组\(\left\{\begin{aligned}x\equiva_1(\modm_1)\\x\equiva_2(\modm_2)\\\cdots\\x\equiva_n(\modm_n)\end{aligned}\right.\)的解,满足\(m_1,m_2,\cdots,m_n\)两两互质。\[设M=\prod_{i=1}^nm_i,Mi=\fracM{m_i},t_i是线性同余方程M_it......
  • Lucas定理——定义、证明、实现、运用
    目录什么是Lucas定理证明Lucas定理Lucas定理求解组合数的C++实现什么是Lucas定理这是一个有助于分解组合数来求解的定理,适合模数小,数字大的问题。有质数\(p\),对于\(n,m\),如果\(n=k_1p+b_1,m=k_2p+b_2\),有\[C_n^m\equivC_{k_1}^{k_2}C_{b_1}^{b_2}\pmodp\]由此可以分解成......
  • 【IT老齐010】CAP定理
    【IT老齐010】CAP定理分布式架构的基本理论。指的是在一个分布式系统中,一致性(Consistency)、可用性(Availability)、分区容错性(Partitiontolerance)。C:更新操作成功后,所有节点在同一时间的数据完全一致。(复习:事务的一致性:事务前后的数据完整性保持一致)A:用户访问数据时,系统能......
  • hdu 5446 长春区域赛网络赛1010 Unknown Treasure(lucas定理+中国剩余定理+移位乘法)
    题目链接:hdu5446题目大意:求出Cmn%M,M=p1⋅p2⋯pk题目分析:首先对于每个质数pi我们,我们可以利用Lucas定理求出Cmn%pi的值,Lucas定理如下:Cmn%p=Cm/pn/p⋅Cm%pn%p%p然后我们可以利用中国剩余定理求取最后答案:M=∏i=1kpi,Mi=M/piCmn%M=∑i=1kCmn%pi⋅Mi⋅inv[Mi]因为做乘法......
  • A*B Problem 485 (数学题+九余数定理)
    A*BProblem1000ms |          内存限制:655352设计一个程序求出A*B,然后将其结果每一位相加得到C,如果C的位数大于等于2,继续将C的各位数相加,直到结果是个一位数k。例如:6*8=48;4+8=12;1+2=3;输出3即可。第一行输入一个数N(0<N<=1000000),表示N组测试数据。......
  • Eddy's digital Roots 1163 (数学+九余数定理)
    Eddy'sdigitalRootsTimeLimit:2000/1000MS(Java/Others)   MemoryLimit:65536/32768K(Java/Others)TotalSubmission(s):5278   AcceptedSubmission(s):2952ProblemDescriptionThedigitalrootofapositiveintegerisfoundbysumming......
  • 中国剩余定理
        #include<bits/stdc++.h>usingnamespacestd;typedeflonglongLL;LLexgcd(LLa,LLb,LL&x,LL&y){if(!b){x=1,y=0;returna;}LLd=exgcd(b,a%b,y,x);y-=a/b*x;re......