信号处理教程_2012_24学时_第2章.ppt

【例2-14】* 求对应　　　　　　的时域信号　　。 解： 　　 　　从而有 【例2-15】* 已知信号　　的频谱为　　　　　　　　　， 　　　求　　。 解： 故 【例2-16】 已知信号　　的频谱为　　　　　　　　　　， 　　　求　　。 解： 其中： 即 所以 【例2-17】 已知信号　　的频谱为　　　　　　　　　　　 　　　　　为一常数，求　　。 解： 应用傅里叶变换频域卷积定理有 利用傅里叶变换性质和常见信号的傅里叶变换对 【例2-18】* 已知　　　　　　　　　　，求　　。 解： 根据傅里叶变换的时域卷积定理和时域微分特性有 　　即 　　故 2.6 Matlab 仿真 例 3-30 ：利用Matlab 呼出周期三角波函数的频谱。 t 0.5 1 f(t) 1 %三角形脉冲信号的傅里叶级数实现 N=10; n1=-N:-1;c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2;%计算n=-N到-1时的傅里叶级数 c0=0; %计算n=0时的傅里叶级数 n2=1:N;c2=-4*j*sin(n2*pi/2)/pi^2./n2.^2; %计算n=1到N时的傅里叶级数 cn=[c1 c0 c2];n=-N:N; subplot 211;stem(n,abs(cn));ylabel(Cn的幅度); subplot 212;stem(n,angle(cn));ylabel(Cn的相位); 例3-31：。 %周期脉冲函数的合成： Function [A_sym,B_sym]=CTFShchsym syms t n k y T=5;tao=0.2*T;a=0.5; if nargin4;Nf=6;end if nargin5;Nn=32;end % x=time_fun_s(t); h=1; x1=sym(Heaviside(t+0.5))*h; x=x1-sym(Heaviside(t-0.5))*h; A0=2*int(x,t,-a,T-a)/T; As=int(2*x*cos(2*pi*n*t/T)/T,t,-a,T-a); Bs=int(2*x*sin(2*pi*n*t/T)/T,t,-a,T-a); Nn=32;Nf=6; A_sym(1)=double(vpa(A0,Nn)); for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end if nargout==0 c=A_sym;disp(c) d=B_sym;disp(d) t=-8*a:0.01:T-a; f1=0.2/2+0.1871.*cos(2*pi*1*t/5)+0.*sin(2*pi*1*t/5);; f2=0.1514.*cos(2*pi*2*t/5)+0.*sin(2*pi*2*t/5);; f3=0.1009.*cos(2*pi*3*t/5)+0.*sin(2*pi*3*t/5);; f4=0.0468.*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5);; f5=-0.0312.*cos(2*pi*6*t/5)+0.*sin(2*pi*6*t/5);; f6=f1+f2; f7=f6+f3; f8=f7+f4+f5; subplot(2,2,1) plot(t,f1),hold on % y=time_fun_e(t) a=0.5;T=5;h=1;tao=0.2*T;t=-8*a:0.01:T-a; e1=1/2+1/2.*sign(t+tao/2); e2=1/2+1/2.*sign(t-tao/2); y=h.*(e1-e2); plot(t,y,r:) title(周期矩形波形的形成-基波) axis([-4,4.5,-0.1,1.1]) subplot(2,2,2) plot(t,f6),hold on % y=time_fun_e(t) a=0.5;T=5;h=1;tao=0.2*T;t=-8*a:0.01:T-a; e1=1/2+1/2.*sign(t+tao/2); e2=1/2+1/2.*sign(t-tao/2); y=h.*(e1-e2); plot(t,y,r:) title(周期矩形波形的形成-基波+2次谐波) axis([-4,4.5,-0.1,1.1]) subplot(2,2,3) plot(t,f7),hold on % y=time_fun_e(t) a=0.5;T=5;h=1;tao=0.2