习题一 绘制典型信号及其频谱图电子工程学院202班单边指数信号的理论表达式为 figure(4);调整,将a 分别等于1、5、10等值,观察时域波形和频域波形。
由于波形 较多,现不失代表性地将a=1和a=5时的各个波形图列表如下进行对比, 其 他a 值的情况类似可推知。
单边指数信号信号 名称 单边 时间函数f t 频谱函数F ■指数 脉冲Ee% t a对提供的MATLAB 程序作了一些说明性的补充, MATLAB 程序为%单边指数信号 clc; close all ; clear all ;E=1;a=1; %调整a 的值,观察不同a 的值对信号波形和频谱的影响 t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); figure(1); plot(t,f);xlabel( 't' );ylabel( 'f(t)' );title( '信号时域图像’);figure(2);plot(w,abs(F));xlabel( '\omega'特性'); figure (3);plot(w,20*log10(abs(F)));xlabel( );ylabel( '|F(\omega)|'); ti tle('幅频'\omega' );ylabel( '|F(\omega)| indB' );title(幅频特性/dB'); plot(w,a ngle(F)*57.29577951);xlabel( )/(°)' );title( '相频特性’);'\omega' );ylabel( '\phi(\omega50 90 80.60 40 22a O/b0 J200升|I IIg C.5-/dB 1\ / \ 1 \ f 1 f11 1>IJ、f1 ■l! «. ■i11 -10 -20- -10ID203加特性出a IB0 -2fJ..j:645-15■24-W■22分析:由上表中a=1和a=5的单边指数信号的波形图和频谱图的对比可以发现, 当a 值增大时,信号的时域波形减小得很快,而其幅频特性的尖峰变宽,相 频特性的曲线趋向平缓。
矩形脉冲信号的理论表达式为MATLAB 程序为:% 矩形脉冲信号 clc;相 频 特 性矩形脉冲信号信号 名称 矩形 脉冲时间函数f t频谱函数F ■2EE Sa2 二■ Sin2close all ; clearall ;E=1; %矩形脉冲幅度 width=2;%对应了时域表达式中的 taot=-4:0.01:4; w=-5:0.01:5; f=E*rectpuls(t,width);F=E*width*si nc(w.*width/2); figure(1); 't' );ylabel( 'f(t)' );title( '信号时域图像’); figure(2);plot(w,abs(F));xlabel( 频特性'); '\omega' );ylabel('|F(\omega)|');title(figure (3);plot(w,20*log10(abs(F)));xlabel( '\omega');ylabel('|F(\omega)| indB' );title( '幅频特性 /dB');figure(4);plot(w,a ngle(F));xlabel( '\omega' );ylabel( '\phi(\omega)' );title(相频特性’);调整,将 分别等于1、4等值,观察时域波形和频域波形。
由于波形较多, 现不失代表性地将a=1和a=4时的各个波形图列表如下进行对比,其他 值的情况类似可推知。
70.S0 2%MATLAB 中的矩形脉冲函数, width 即是tao , t 为时间plot(t,f);xlabel(C幅频特性/dB相频特性幅频特叫HBY V-1J0.Lin-1 ■-5 --3> j016-3 -2-10123分析:由以上的图标对比可知,(1) 解释“幅值特性/dB ”中许多向下跳变的尖峰这是由于求取分贝数要用lg 函数,IgO 为负无穷,所以出现 了图像中的很多向下跳变的尖峰。
实际上,矩形脉冲信号一般不 看以分贝为单位的幅频特性曲线。
升余弦信号的理论表达式为: 名称MATLAB 程序为:% 升余弦信号 clc; close w=-5:0.01:5;f1=E*rectpuls(t,width); %MATLAB 中的矩形脉冲函数, width 即是 tao , t 为时间f=0.5*(1+cos(2*pi.*t/width)).*f1; %用矩形脉冲函数乘以因子得到升余弦函数F=E*width*si nc(w.*width/2)*0.5./(1-(w*width*0.5/pi).A 2);figure(1);三、 升余弦脉冲信号信号 时间函数f t 频谱函数F ■升余 弦脉习+ cosE Saall clear allE=1; width=2;t=-4:0.01:4; %对应了时域表达式中的 taoplot(t,f);xlabel( figure (2);'t' );ylabel( 'f(t)' );title('信号时域图像’); plot(w,abs(F));xlabel( 频特性');figure (3);'\omega' );ylabel( '|F(\omega)|');title('幅plot(w,20*log10(abs(F)));xlabel( '幅频特性/dB');'\omega' );ylabel( '|F(\omega)| indB' );title( f igure(4); plot(w,angle(F));xlabel( 相频特'\omega');ylabel( '\phi(\omega)' );title(调整,将分别等于1、4等值,观察时域波形和频域波形。
由于波形较多,现不失代表性地将=1和=4时的各个波形图列表如下进行对比,其他值的情况类似可推知。
分析:(1) 首先解释 时,幅值谱中出现的极大值的原因如下所示,生余弦脉冲的时域频域表达式如下所示。
由生余弦函数的傅 立叶变换表达式可知,当分母等于0时,幅值就会变为无穷。
图中的极大值 即是值接近极点,使得幅值跳变到了很大的值。
升余(2) 解释“幅值特性/dB ”中许多向下跳变的尖峰这是由于求取分贝数要用lg 函数,IgO 为负无穷,所以出现了图像中的很多向下跳变的尖峰。
实际上,升余弦信号一般不看以分贝为单位的幅频特弦脉匚丿用律,吒丿E Sa 2丿21 - v2 n>2相 特 性频性曲线。
MATLAB 程序为:%三角脉冲信号 clc; close调整,将分别等于2、4等值,观察时域波形和频域波形。
由于波形较多, 现不失代表性地将=2和=4时的各个波形图列表如下进行对比,其他 值的情况类似可推知。
四、 三角脉冲信号 三角脉冲信号的理论表达式为:信号 名称 时间函数f t三角 脉冲8E .—2— sin I - o T I 4 丿all clear allE=1;width=4; t=-4:0.01:4;%对应了时域表达式中的 tao w=-5:0.01:5; f=E*tripuls(t,width); %MATLAB 中的三角脉冲函数,F=0.5*E*width*(s in c(w.*width/4).A 2);figure(1); width 即是tao , t 为时间plot(t,f);xlabel( figure (2);'t' );ylabel('f(t)');title('信号时域图像’);plot(w,abs(F));xlabel( 特性'); figure (3);plot(w,20*log10(abs(F)));xlabel( '幅频特性/dB'); '\omega');ylabel( '|F(\omega)|');title('\omega' );ylabel( '|F(\omega)| in'幅频dB' );title( figure(4);plot(w,a ngle(F)*57.29577951);xlabel( )/ (°)' );title('相频特性’);'\omega' );ylabel( '\phi(\omega廷肖茫叫分析:(1)首先对比和4时的结果,可以明显看到三角脉冲宽度变宽之后其频域的幅频特性曲线反而变窄了,这与理论公式的结果相一致。
(2)解释“幅值特性/dB”中许多向下跳变的尖峰这是由于求取分贝数要用lg函数,IgO为负无穷,所以出现了图像中的很多向下跳变的尖峰。
实际上,升余弦信号一般不看以分贝为单位的幅频特性曲线。
(3)相频特性曲线中,可以明显看到其相频特性曲线的角度一直为0°这是因为三角脉冲的傅立叶变换表达式一直为实数,这与公式也是相符合的,是三角脉冲的特点。
11 /11。