Arcct-CSDN博客

原创 Digital Signal Processing Using Matlab读书笔记第三章（七）（P63-P67）

定理3 取样原则一个频带宽度为F0的带限信号xa(t)，如过它的取样频率Fs大于频带宽度F0的两倍，则可以通过它的取样值x(n)重建。否则将会导致x(n)信号的错误识别，一个模拟带宽信号的取样频率为2F0的话，该频率称为奈奎斯特频率。应该注意的是，在xa(t)取样后，x(n)代表的最高模拟频率是Fs/2（或...

2021-05-10 11:49:10 184

原创 Digital Signal Processing Using Matlab读书笔记第三章（六）（P59-P63）

例3.16>> b=[0.0181,0.0543,0.0543,0.0181];>> a=[1.0000,-1.7600,1.1829,-0.2781];>> m=0:length(b)-1;>> l=0:length(a)-1;>> K=500;k=0:1:K;>> w=pi*k/K;>> num=b*exp(-j*m'*w);>> den=a*exp(-j*l'*w);>> H

2021-04-24 14:29:52 109

原创 Digital Signal Processing Using Matlab读书笔记第三章（五）（P55-P59）

例3.13>> w=[0:1:500]*pi/500;>> H=exp(j*w)./(exp(j*w)-0.9*ones(1,501));>> magH=abs(H);angH=angle(H);>> subplot(2,1,1);plot(w/pi,magH);>> grid>> subplot(2,1,2);plot(w/pi,angH/pi);grid例3.14因为输入不是绝对可加的，DTFT不能计算全

2021-04-21 17:23:41 118

原创 Digital Signal Processing Using Matlab读书笔记第三章（四）（P51-P55）

例3.10>> n=-5:10;x=rand(1,length(n))+j*rand(1,length(n));>> k=-100:100;>> w=(pi/100)*k;>> X=x*(exp(-j*pi/100)).^(n'*k);>> y=conj(x);>> Y=y*(exp(-j*pi/100)).^(n'*k);>> Y_check=conj(fliplr(X));>> error=

2021-04-21 16:13:16 156

原创 Digital Signal Processing Using Matlab读书笔记第三章（三）（P47-P51）

DTFT的属性在前面的章节中，我们为了展示图像讨论了两个重要的属性，现在我们讨论剩下的有用的属性。应用如果序列x(n)是实数且是偶对称，那么X也是实数偶对称。复数序列的一个类似的性质将会在问题3.7中讨论。在下面几个例子中，我们将会使用有限区间序列证明其中一些性质例3.7...

2021-04-13 09:41:30 146

原创 Digital Signal Processing Using Matlab读书笔记第三章（二）（P43-P47）

如果x(n)是有限区间，那么matlab可以数字上地用来计算任何频率下的X，方法就是直接应用(3.1)。除此之外，我们通过等间隔分割[0:pi]的方法来估计X，然后(3.1)可以运用在矩阵向量乘积运算，为了理解这个，让我们假设序列x(n)在n1<=n<=nN有N个取样信号(不一定在[0,N-1])，我们以下面这一方式估计X式子中表示在[0,pi]之间(M+1)个等间隔频率，3.1可以写成下式当x(n)和X各自表示成列向量x和X时，有W是一个（M+1)*N的矩阵除此之外，如果我们.

2021-04-12 10:03:18 233

原创 Digital Signal Processing Using Matlab读书笔记第三章（一）（P40-P42）

离散时间傅里叶分析我们已经看到一个线性时不变系统是如何通过它的单位冲激响应表示的，单位冲激响应使得我们可以通过卷积计算任意的输入x(n)对应的响应该卷积的表示形式是基于能通过按比例以及时延的单位信号的线性结合，每一个基础信号提供了一个新的信号形式，每一形式在不同的系统下既有优点，也有缺点，但是，当系统是线性以及时不变的，只有一个表示形式是最有用的，这就是离散时间傅里叶变换。离散时间傅里叶变换（DTFT）如果x(n)是绝对可和的，即那么它的DTFT如下X的反向离散时间傅里叶变换(IDTFT)

2021-04-10 16:47:40 142 1

原创 Digital Signal Processing Using Matlab读书笔记第二章（六）（P32-P35）

例2.5

2021-04-06 12:33:28 245

原创 Digital Signal Processing Using Matlab读书笔记第二章（五）（P27-P32）

回顾序列相关性如果我们将卷积操作（2.11）和两序列相关性的定义（2.6）进行比较，我们可以发现二者有紧密的相关性如果序列是有限的，这些相关性的计算可以用conv函数实现例2.8这个例子我们会阐述互相关序列的应用...

2021-04-03 21:48:29 175

原创 Digital Signal Processing Using Matlab读书笔记第二章（四）（P19-P27）

几何级数一个单边指数信号，形如α是一个任意的常数，被称为几何级数，级数的和在数字信号处理中有许多应用序列的相关性

2021-04-03 19:59:53 228 2

原创 Digital Signal Processing Using Matlab读书笔记第二章（三）（P15-P19）

例2.2a题n=-2:10;[x11,n11]=sigshift(x,n,5);x=[1:7,6:-1:1];[x11,n11]=sigshift(x,n,5);[x12,n12]=sigshift(x,n,-4);[x1,n1]=sigadd(2*x11,n11,-3*x12,n12);stem(n1,x1);b题>> n=-2:10;x=[1:7,6:-1:1];>> [x11,n11]=sigfold(x,n);>> [x12,n12]=

2021-03-28 22:03:08 181

原创 Digital Signal Processing Using Matlab读书笔记第二章（二）（P10-P15）

对序列的操作1.信号相加在matlab中可以用"+“符号但是X1(n)和X2(n)的长度必须相等，如果序列长度不相等，或者相同长度序列的取样位置不同，我们不能直接用”+"符号，我们必须首先扩充X1和X2，以便于它们有相同的位置向量n，这需要格外注意matlab的索引操作，特别的，逻辑操作符“&”，关系操作符“<=”和“==”以及查找函数（find函数作用是查询非零元素的位置）都需要X1和X2的长度相等，下面的函数称为sigadd函数，解释了这些操作function [y,n]=siga

2021-03-27 21:15:17 249

原创 Digital Signal Processing Using Matlab读书笔记第二章（一）（P7-P10）

第二章-离散时间信号与系统一个模拟信号可以用x(t)表示，这个变量t代表任何物理量，但我们通常假设它为秒来代表时间，一个离散信号用x(n)来表示，该变量n是整数值，代表时间上的离散量，，所以也被称为离散时间信号，该离散时间信号是一个数字序列，可以通过下列形式表示向上的箭头示意n=0时的值。在matlab我们可以通过由适当的值构成的行适量表示一个有限区间的序列。但是，这样一个矢量没有任何信息关于该序列的位置，所以一个x(n)的正确表示需要两个矢量，一个用于x一个用于n，举个例子，一个序列x(n)={2

2021-03-26 09:36:11 390

原创定点数与浮点数与单精度数

定点数定点数分为定点整数与定点小数，运算中不考虑小数点的问题。N=S*2^p，S是N的全部有效数字，称为N的尾数，p是N的阶码，表示小数点的位置。浮点数B=（+/-）S*2^（+/-）pp是阶码（常用补码），p前的（+/-）是阶符，S是尾数（常用原码），S前的（+/-）是尾数符号这个1表示的是S在源码状态下的规格化表示法例子解释规格化S在源码状态下的规格化表示法设字长8位，阶码2位，位数4位110.1=0.11012^3,Pf=0，Sf=0，规格化表示为01101101110.1

2021-03-09 22:43:12 599

Arcct的博客