资源描述：
《【精品】连续时间信号频域分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验三连续时间信号的频域分析实验人：Miwan1实验目的（1）掌握连续时间信号傅里叶变换的实现方法；（2）了解傅里叶变换的特点及其应用；（3）掌握傅里叶变换的数值计算方法以及绘制信号频谱图的方法；2实验原理（1）傅甲叶级数的三介形式（2）傅甲叶级数的指数形式（3）非周期信号的傅里叶变换3涉及的Matlab函数（1）fourier函数：实现信号的傅里叶变换。调用格式：F=fourier（D是符号函数f的傅里叶变换，默认返回函数F是关丁•w的函数。说明：该函数假定两个序列都从零开始。（2）quad8函数：用來计算非周期信号的频谱。调用格式：y=quad8Cf
2、unc^,a,b）y=quad8（tfunc,,a,b,TOL,TRACE,pl,p2,...）其中：func是一个字符串，表示被积函数的.m文件名；a,b分别表示定积分的下限和上限;TOL表示指定允许的相对或绝对积分误差，非零的TRACE表示以被积瓯数的点绘图形式来跟踪该quad8函数牛成的返回值，如果TOL和TRACE均赋以空矩阵，则两者均口动使用默认值；'pl,p2,'表示被积函数所需要的多个额外输入参数。（3）quadl函数：计算非周期信号的频谱。调川格式：y=quad1（Tunc',a,b）y=quadl（＜func,,a,b,TOL,TRAC
3、E,p1,p2,...）其中func是一个字符串，表示被积函数4实验内容与方法（1）周期信号的傅里叶级数的Matlab实现①利用Matlab画出周期三角波信号的频谱。源程序如下：N=10;nl=-N:-1;cl=-4*j*sin(nl*pi/2)/piA2./nl.A2;c0=0;n2=l:N;c2=-4*j*sin(n2*pi/2)/piA2./nl.A2;cn=[clcOc2];n=-N:N;subplot211;stem(nzabs(cn));ylabel('cnfl勺幅度');subplot212;stem(n,angle(cn));ylabe丄
4、(1cn0勺扌U位1);xlabel(1omega/omega_0')程序运行结果如下图1所示图1①周期矩形脉冲幅度为1,脉冲宽度为tao=l,周期T=5*tao,试求出该信号的三角形式的傅里叶系数，并绘出各次谐波叠加的傅里叶综合波形图。源程序如下：%周期矩形脉冲函数的合成function[A_symzB_sym]=CTFSchsymsymstnkx;T=5;tao=0.2*T;a=0.5;ifnargin<4;Nf=6;endifnargin<5;Nn=32;endx=time_fun_x(t);AO=2*int(x,t,-a,T-a)/T;As=
5、int(2*x*cos(2*pi*n*t/T)/Tzt,-a,T-a);Bs=int(2*x*sin(2*pi*n*t/T)/Tzt,-a,T-a);A_sym(1)=doub丄e(vpa(AO,Nn));fork=l:NfA_sym(k+1)=double(vpa(subs(AsAnzk),Nn));B_sym(k+1)=double(vpa(subs(Bs,nzk),Nn));endifnargout==0c=A_sym;disp(c)d=B_sym;disp(d)t=-8*a:0.01:T-a;fl=0.2/2+0.1871.*cos(2*pi*l
6、*t/5)+0.*sin(2*pi*l*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壮/5);f4=0.0468.*cos(2*pi*4壮/5)+0・*sin(2*pi*4*t/5);f5=-0.0312.*cos(2*pi*6壮/5)+0.*sin(2*pi*6壮/5);f6=fl+f2;f7=f6+f3;f8=f7+f4+f5;subplot221plot(tzfl)fholdony=time_fun_e(t);p
7、lot(tzy,*r:')title(•周期矩形波的形成一一基波-)axis([-4,4.5,-0.1,1.1])subplot222plot(t,f6)Aholdony=time_fun_e(t);plot(tzy,'r:')title(•周期矩形波的形成——基波+二次谐波*)axis([-4,4.5,-0.1,1.1])subplot223plot(tAf7)zholdony=time_fun_e(t);plot(t,yz'r:')t让h(•基波+二次谐波+二次谐波1)axis([-4,4.5,-0.1,1.11)subplot224plot(t,f
8、8),holdony=time_fun_e(t);plot(t,yz1r:1)t
