MATLAB举例.ppt

第二章例子 例2.2.2画出以下各序列在给定区间的波形 (a) x(n)=1.5*δ(n+1) - δ(n-3) -4≤n≤5 (b) x(n)=n[u(n) – u(n-8)] - 10exp(-0.3(n-10))[u(n-10) – u(n-16)] 0≤n≤20 (c) x(n)=cos(0.07πn) + 0.2w(n) 0≤n≤30其中w(n)为具有零均值的单位方差的高斯随机序列 (d)求(b)中函数 x(n)的能量E(n) 程序hc222，需要调用impseq.M与stepseq.M 第二章例子 例2.3.4和2.3.5 设x = [3,-3,7,0,-1,5,2]; nx = [-4:2]; h = [2,3,0,-5,2,1]; nh = [-1:4]; 求其卷积y(n)。 程序hc234,hc235，调用convwthn.M 其中hc234为手工计算，画出计算过程 其中hc235为利用MATLAB程序中的卷积函数计算卷积 注意：卷积结果的长度为N1＋N2－1 第三章例子 例3.2.5 求序列x(n)=[1,3,5,3,1]的DTFT，画出它在ω=-8~8范围内的频率特性。并把x(n)的位置零点左右移动，讨论时移对DTFT的影响。 计算公式为： 程序hc325 说明： 1、序列的DTFT是连续函数； 2、序列的DTFT是周期函数 ； 3、实序列的DTFT具有对称性。 4、信号在时间轴上的平移不影响其DTFT的幅频特性，只影响它的相频特性。 5、时域对称的序列。它具有相位随频率线性变化的特点。对称中心所处的位置决定了相频特性的斜率的大小。 第三章例子 例3.5.1 求x=[2,-1,1,1]的离散傅立叶变换。并画到图上，与其DTFT比较。 其DTFT已经在程序hc324中求出，本例程序hc351 例3.5.2：重做例3.5.1，频谱取64点。（尝试改变补零的点数） 程序为hc352 第四章例子 例4.3.2 考虑长度为11的矩形窗函数(关于Y轴对称，需循环移位)序列，计算其频谱。 解：假如选N=20作为重复周期，则要在序列后面补9个零。使用FFT时，我们必须按N=20的周期延拓序列中取从n=0到19的主值部分，因此FFT的输入为 x=[ones(1,6),zeros(1,N-11),ones(1,5)] 在编写程序时，要准备给出不同的N进行比较。 所以程序hc432的写法应能适用于不同的N。 第四章例子 例4.5.4 计算定义在全部t上的xa(t)=cos5t的频谱，它的理论频谱是： 它包含了权重为π的位于Ω=±5上的两个脉冲函数。 解：在计算机计算中，正余弦函数必须截断为有限长度L。信号＝cos5t的带宽限制于5。纯理论地看，只要采样周期小于π/5=0.63秒，就不会发生频率混叠。然而如果把cos5t截断为长L的信号，则它的频谱就不再是有限带宽了，所以必须采用更小的采样周期，任选T＝0.1，并选N＝50，得到L＝TN＝5。按此来截断信号。 * DFT的Goertzel算法 (9.4.9)中的递推关系对n=0，1，……，N重复N+1次，每次只需要计算一次实数乘和两次实数加，而带有复数运算的方程(9.4.10)仅在n＝N时刻计算一次。所以，对实数序列x(n)，由于对称性，用这种算法求出X(k)和X(N-k)的值只需要N次实数乘法和一次复数乘法运算。 实际上，不必计算复数值X(k)，只要求出幅度平方值|X(k)|2 ，故计算的最后一步还可以简化，从而完全避开复数运算。 * Goertzel算法子程序 按(9.4.7)式编成子程序gfft。它根据输入序列x，和规定的DFT样本序号k，计算DFT样本X。 function X=gfft(x,k) % 用Goertzel算法计算序号为k的DFT样本 N=length(x);x1=[x,0]; % 递推要N+1次， d1 = 2*cos(2*pi*k/N); %滤波中间项系数 v = filter(1,[1,-d1,1],x1);% 用滤波实现卷积 W = exp(-i*2*pi*k/N); % 为下一步求W X = v(N+1) - W*v(N); % 求出第k个DFT * 在双音拨号中的应用 DTMF信号参数选择的考虑因素： (1).为了抗干扰，要检测8个频率及其倍频； (2).为了能分辨清这16个频率，要求采集数据的长度NFs/D，已知Fs=8000，而根据最小相邻基频间隔（73Hz）大于2D的要求，D37，得到N应在200以上。而为了提高检