首页 > 其他分享 >01. 数字信号处理概述

01. 数字信号处理概述

时间:2023-11-30 16:13:13浏览次数:37  
标签:infty 01 frac limits 信号处理 概述 pi Omega sum

符号规定

\(\Omega\) 模拟角频率
\(\omega\) 数字角频率
\(T\) 采样周期
\(\Omega_{s}\) 采样角频率
\(\Omega_N\) 连续信号最大角频率

连续信号的采样

对于一个连续信号\(x_{c}(t)\),其经过冲击采样之后获得之后变成一个冲击串函数\(x_{s}(t)\),即:

\[x_{s}(t) = \sum\limits_{n=-\infty}^{\infty}x_{c}(t)\delta(t-nT) = \left \{ \begin{matrix} x_{c}(nT) &t=nT \\ 0 &其他 \end{matrix} \right. \]

显然,对于\(x_{s}(t)\)任然是一个连续信号,我们需要经过量化之后获得我们所需要的离散信号\(x[n]\),有\(x[n] = x_{c}{nT}\)
流程如图:
image

频域分析

对\(x_{s}(t)\)进行连续傅里叶变换有:

\[X_{s}(j\Omega) = \frac{1}{2\pi}X_{c}(j\Omega)*S(j\Omega) \]

其中

\[S(j\Omega) = \frac{2\pi}{T}\sum\limits_{k=-\infty}^{\infty}\delta(\Omega-k\Omega_{s}) \]

所以有:

\[X_{s}(j\Omega) = \frac{1}{T}\sum\limits_{k=-\infty}^{\infty}X_{c}(j(\Omega-k\Omega_{s})) \]

显然,采样后的频域是对原本连续信号频域的周期延拓,这里面要求有:

\[\Omega_{s} > 2\Omega_{N} \]

这里\(\Omega_{N}\)称之为奈奎斯特频率,\(2\Omega_{N}\)称之为奈奎斯特率

在对\(x[n]\)进行DTFT,有:

\[X(e^{j\omega}) = \frac{1}{T} \sum\limits X[n]e^{-j\omega n} \]

对比上述式子,有:

\[X_{s}({j\Omega}) = X(e^{j\omega})|_{\omega= \Omega T} \]

所以上面冲击串函数的在数字频域上可以表示:

\[X_{s}(j\omega ) = \frac{1}{T}\sum\limits_{k=-\infty} ^{\infty}X_{c}(j(\frac{\omega}{T}-2\pi \frac{k}{T})) \]

即对采样之后在数字角频域进行了归一化

例1

假设\(T=1/600\),模拟连续信号\(x_{c}(t) = cos400\pi t\)

可以计算出:\(\Omega_{s}=1200\pi, \Omega_{N} = 800\pi\)
由时域采样定理知:可以用一组离散信号表示原连续信号
此时\(X_{s}(j\Omega)\)图像如下:
image

在进行归一化有:
image

连续信号的复原

对于数字序列\(x[n]\),我们可以构造一个连续的冲击串函数\(x_{s}{t}\),有:

\[x_{s}(t) = \sum\limits_{n=-\infty}^{\infty}x[n]\delta(t-nT) \]

对于这个函数,我们已知其在频域上等于原连续信号的周期延拓,所以我们希望使用一个低通滤波器对其进行滤波再添加一个T倍的增益,该门函数的截止频率一般取\(\Omega_N\)的一般。


对于截止频率为\(\Omega_{H}\),频率增益为A的门函数,其时域形式为:

\[h(t) = \frac{1}{\pi}A\Omega_{H}sa(t\Omega_{H}) \]


所以有:

\[h_{r}(t) = sa(\frac{t\Omega_{N}}{2}) = \frac{\sin(\pi t/T)}{\pi t/T} \]

所以在频域上将\(X_{s}(j\Omega)与H(j\Omega)\)相乘,在时域上表现为卷积,有:

\[x_{r}(t)= x_{s}(t)*h_{r}(t) = \sum\limits_{n=-\infty}^{\infty}x[n]h_{r}(t-nT) = \sum\limits_{n=-\infty}^{\infty}x[n]\frac{sin(\pi (\frac{t}{T}-n)}{\pi(\frac{t}{T}-n)} \]

则可以抽象成一个D/C转换器
image

完整的步骤

时域连续信号在进过采样之后变成数字序列,而DSP处理的对象即为数字序列,经过输入到系统中后输入一个响应,将输出响应再进行复原成时域连续信号。

后面我们将围绕中间系统部分详细展开。

标签:infty,01,frac,limits,信号处理,概述,pi,Omega,sum
From: https://www.cnblogs.com/C-qian/p/17867624.html

相关文章

  • vulhub中spring的CVE-2018-1273漏洞复现
    ​ 1.影响版本SpringDataCommons1.13-1.13.10(IngallsSR10)SpringDataREST2.6-2.6.10(IngallsSR10)SpringDataCommons2.0to2.0.5(KaySR5)SpringDataREST3.0-3.0.5(KaySR5)2.漏洞原理SpringData是一个用于简化数据库访问,并支持云服务的开源框......
  • P1084 [NOIP2012 提高组] 疫情控制
    首先军队可以原地不动,时间越多越容易合法,先套上二分。在不回到根的情况下,军队深度肯定越小越好。所以军队能往上移就移,如果能回到根就暂时在根对应的儿子那里驻扎。这个过程用树上倍增优化。做完这一步后,我们找出需要军队驻扎的根的儿子(向下不经过军队就能到达叶子),现在就是要让......
  • P1081 [NOIP2012 提高组] 开车旅行
    题目有点长,一步一步来。预处理出每座城市两人分别会选择的下一座城市用set即可实现。倍增优化DP令\(f_{i,j}\)表示从城市\(j\)出发,行驶\(2^i\)天会到达的城市。令\(ga_{i,j}\)表示从城市\(j\)出发,行驶\(2^i\)天,小A行驶的路程。\(gb_{i,j}\)同理。答案枚......
  • HPM5301EVKlite windows下上手使用
    HPM5301是先辑半导体出的一款300Mhz的riscV芯片,售价7.99,开发板售价39.99,性价比超高资源获取evklite基本不需要原理图,下载SDK和ENV即可SDK连接:https://github.com/hpmicro/hpm_sdkENV连接:https://github.com/hpmicro/sdk_env占坑待续。。。......
  • P2495 [SDOI2011] 消耗战
    题意给定一棵有边权的无根树。\(q\)次询问,每次询问\(k\)个点。求断边使得根节点\(1\)与\(k\)个点不连通的最小边权。Sol虚树。\(n^2\)dp是trivial的。考虑优化。注意到其中很多点都是无用的。考虑保留有效点。不难发现,有效点集为询问点两两\(lca\)的集合......
  • [SUCTF2018]followme
    [SUCTF2018]followme看起来像是在做渗透测试的时候抓的包,导出所有HTTP对象大部分是类似这种的文件login%3f_=6975b9a9f7a359d322e06c0e28db112blogin.php%3freferer=http%3A%2F%2F192.168.128.145%2Fbuild%2Fadmin%2F<!DOCTYPEHTML><htmlclass="no-js"><head>......
  • 初中英语优秀范文100篇-012 My Experience of Being a Volunteer - 我的一次志愿者经
    PDF格式公众号回复关键字:SHCZFW012记忆树1Lastyear,Ipaidavisittothehomefortheagedwithmyclassmatesasvolunteers.翻译去年,我和我的同学作为志愿者去老年人之家探望了老人们。简化记忆探望老人句子结构这个句子可以分为四个主要部分:1状语短语:“La......
  • 聪明办法学python—task05&选学01
    条件语句if语句:1.elseif——>elif2.后面是冒号3.其余与c语言相同循环结构while循环while条件:bodystatmentsfor-in循环for变量in字符串/范围/集合:statements结束循环breakcontinuereturn与c语言相同 ......
  • #Y1010. 平衡括号计数
    平衡字符串的性质将左括号看出1,右括号看成-1(1)最后加起来是0(2)中间一直要是正数所以代码如下usingnamespacestd;intn;intdx[2]={1,-1};intans=0;voiddfs(intcur,intx){ if(x<0)return; if(cur==n){ if(x==0)ans++; return; } for(inti=0;i<2;i++){ dfs......
  • 聪明办法学Python Task05 & 选学01
    聪明办法学Python学习笔记Chapter4条件Conditionalsif语句ifstatementConditionalsMakeDecisionsif语句流程一个例子deff(x):print("A",end="")ifx==0:print("B",end="")print("C",end=&......