数字信号处理第3章导论.ppt

* 3.包含多重极点 所以 * 3.包含多重极点 例3.9-3 求 的拉氏逆变换。 解：F(s)展开式为 令 于是有 * 3.9.2 留数法 这是复变函数积分问题. 我们可以从 到 补足一条积分路径，构成一闭合围线积分， 沿该圆弧的积分应为零。由留数定理： 补足的这条路径CR是半径为?的圆弧， 若pi为一阶极点，则 若pi为一阶极点，则 若pi为k阶极点，则 * 3.9.2 留数法 F0(s)有一个单值极点 p1= –3，和一个二阶极点， p2= –1 例3.9-4 求 的拉氏逆变换。 解： F(s)不是真分式，首先将其分解为有理多项式与真分式之和。 记 它们的留数分别为 * 3.9.2 留数法 而 所以, F(s)的拉氏逆变换为 * 3.10 连续信号的频域与复频域的MATLAB分析 连续信号的频域和复频域表达式可以通过符号运算获得。其频谱的可视化可以用幅度谱和相位谱绘制。周期信号可以通过计算其傅里叶级数，画出它的幅度谱和相位谱；非周期性信号可以通过计算其傅里叶变换，画出它的幅度谱和相位谱。信号的复频域分析一般缺少可视化的直观表示，但可以用信号的拉普拉斯变换，绘制它的三维幅度曲面图和相位曲面图，来观察其复频域特征。 Stem () Plot() Ezplot() laplace() Ilaplace() * 3.10 连续信号的频域与复频域的MATLAB分析 周期性矩形脉冲信号如图3.10-1所示，画出它的幅度谱和相 位谱，以及前5次谐波叠加波形和前10次谐波叠加波形。 例3.10-1 解： * 3.10 连续信号的频域与复频域的MATLAB分析 解： n=-10:10; w1=0.4*pi; ％显示的谐波次数 n1=-10: -1;ft1=sin(0.2*pi*n1)./(pi*n1); ％计算负半轴的傅里叶级数 n2=1:10;ft2=sin(0.2*pi*n2)./(pi*n2); ％计算正半轴的傅里叶级数 ft=[ft1,0.2,ft2]; ％组合负半轴、零点和正半轴 的级数 fn = abs(ft);phase = angle(ft); ％计算幅度谱和相位谱 subplot(2,2,1);stem(n,fn);title(幅度谱); ％stem函数绘制离散序列 subplot(2,2,2);stem(n,phase);title(相位谱); syms t;s1=0.2;s2=0.2; ％直流分量 for k1=1:5 s1=s1+2*sin(k1*pi/5)*cos(w1*t*k1)/pi./k1; end for k2=1:10 s2=s2+2*sin(k2*pi/5)*cos(w1*t*k2)/pi./k2; end subplot(2,2,3);ezplot(s1);title(前5次谐波叠加); subplot(2,2,4);ezplot(s2);title(前10次谐波叠加); * 3.10 连续信号的频域与复频域的MATLAB分析 例3.10-2 用MATLAB分别绘制抽样信号 和矩形脉冲信号 的时域波形和频谱，并验证傅里叶变换的对偶性。 解： syms t;f1 = sin(t)/t; ％抽样函数f1(t)=Sa(t) f2 = pi*sym((Heaviside(t+1) -Heaviside(t-1))); ％计算门函数f2(t)= p G(t) F1=simple(fourier(f1));F2=simple(fourier(f2)); subplot(221);ezplot(f1,[ -10 10]); subplot(222);ezplot(F1,[ -2 2]);title(π[u(w+1) -u(w-1)]); subplot(223);ezplot(f2,[ -2 2]);title(π[u(t+1) -u(t-1)]); subplot(224);ezplot(F2,[ -10 10]); * 3.10 连续信号的频域与复频域的MATLAB分析 * * 3.7.3 典型信号的拉普拉斯变换 4．t的正幂信号 （ 是正整数） 所以 当 时 当 时 以此类推，得 * 3.8 拉普拉斯变换的基本性质 1．线性特性 若 ， ，则 例3.8-1 求