信号与系统分析第四章.ppt

依上述算法实现的相应C语言算法如下: floatrbg(float(*f)(),floata,floatb) { floath,sum,R［101］［101］,err; inti,j,k,n; h=(b-a); R［1］［1］=h*(f(a)+f(b))/2; n=100; for(i=1;i=n;i++) { sum=0; for(k=1;k=2*i-2;k++) sum=sum+f(a+(k-0.5)*h) h=h/2; R［i］［1］=0.5*(R［i-1］［1］+h*sum); for(j=2;j=i;j++) { err=(R［i］［j-1］-R［i-1］［j-1］)/(power(4,j-1)-1); R［i］［j］=R［i］［j-1］+err; err=(R［i］［j-1］-R［i-1］［j-1］)/(power(4,j-1)-1); R［i］［j］=R［i］［j-1］+err; if(abs(err)E)thenreturnR［i］［j］ } } return0 例4―21 计算抽样函数 的傅里叶变换。 解 式(4―132)抽样函数的傅里叶变换可写为 (4―132) 这样我们可利用龙贝格算法求出该函数的数值解,程序如下: floatt; floatD1(floatw) { returnsin(t)*cos(w*t)/t; } floatS1(floatw) { returnsin(t)*sin(w*t)/t; } main() { inti,a,b; a=LEFTVALUE; b=RIGHTVALUE; dt=STEP; t=STARTVALUE; for(i=1;in;i++) { D［i］=rbg(D1,a,b); S［i］=rbg(S1,a,b); t=t+dt; } } 这就是时域卷积定律，可简记为 （4-80） 2. 频域卷积性质 同时域卷积定律一样，我们也可以证明频域卷积 定律。在这里略去证明，只写出结论。 (4-81) (4-82) 例4―14 求图4.15所示梯形脉冲的傅里叶变换。 图4.15 梯形脉冲的傅里叶变换 解 梯形脉冲可看作是两个不同宽度的矩形脉冲 f1(t)与f2(t)的卷积，如图4.15所示。 f(t)=f1(t)*f2(t) 而矩形脉冲的傅里叶变换已在例4―3中求出,具体来说 图4.16 半波正弦脉冲 图4.17 三角形脉冲及其一、二街导的波形 表4―2 傅里叶变换的性质 4.4 线性非时变系统的频域分析 4.4.1 频域分析 ?在第二章线性非时变系统的时域分析中,我们已经指出线性非时变系统的零状态响应yf(t)是激励f(t)与冲击响应h(t)的卷积积分。即 (4―112) (4―113) 现在我们设h(t)、f(t)和yf(t)各自对应的傅里叶变换式为H(ω)、F(ω)和Yf(ω),　即 (4―114) (4―115) (4―117) (4―116) 图4.18 频域分析示意图 接下来我们讨论一下H(ω)。很显然在前面的讨论中,H(ω)已有了两个方面的含义:一是前面我们在式（4―114）中定义的H(ω)是与冲激响应h(t)对应的频谱函数;二是通过式（4―117）所表达的H(ω)的含义是零状态响应频谱函数Yf(ω)与激励函数F(ω)的比值,即 (4―118) 另外,考虑到式（4―115） 在看式（4-116）和（4-117） （4-119） （4-200） 例4―20如图4.19所示,试分析单位阶跃信号u(t)通过RC高通网络传输后的波形。 图 4.19 解 显然,当输入信号uS(t)为复指数信号e jωt时,如图有 则按H(ω)的定义有 对于单位阶跃信号u(t)而言,此时 最后一步考虑了冲激函数的取样性质。因此 4.4.2 无失真传输系统的频域分析 ? 无失真传输系统是指这样一个系统,它的输出信号与输入信号相比只有幅度的大小和出现时间的先后不同,而没有波形上的变化。系统的数学模型具有如下的形式: yf(t)=Kf(t+td) (4―121) 设输入信号f(t)的谱频函数为F(ω),输出信号yf(t)的谱频