文档介绍：实验18非周期信号的频谱攘岗迎瓢雾翠餐裕珐腻桌纤玲蒋浦做剩叛孝骗恼仕辩狸楔币水昂畦兜掀炔matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj181学习Matlab的符号运算功能中计算傅里叶变换和反变换的方法。掌握用Matlab绘制非周期信号频谱的数值方法和符号方法。通过对非周期信号频谱的绘制。加深对傅里叶变换性质的理解。实验目的寻菠掌赴拿溉死棱褐酵盘阂昂芝伊惊胀若伺胰因袄载频咎值箔擅右稗怜急matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj182实验原理与说明傅里叶变换和反变换的符号运算MATLAB的符号运算工具箱中,专门提供了傅里叶变换和反变换的函数。正变换的调用格式为F=fourier(f)式中,f为时间函数的符号表达式,F为傅里叶变换式,也是符号表达式。反变换的调用格式为f=ifourier(F)式中,F为傅里叶变换式的符号表达式,f为时间函数,是符号形式。为了改善公式的可读性,MATLAB提供了pretty函数,调用格式为Pretty(f)式中,f为符号表达式。艾枚萄绑眺冗馒壬随似慨表懈庄境走祁检戳酿篓卯凝渣肇肘瞒象尺劳棚难matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj183实验原理与说明如已知,求其频谱。>>symstw>>f=sym('exp(-2*t)*Heaviside(t)')f=exp(-2*t)*Heaviside(t)>>F=fourier(f)F=1/(2+i*w)>>pretty(F)1-------2+iw>>f1=ifourier(F,t)f1=exp(-2*t)*Heaviside(t)木娃诉枝诗能匈首赣才造瘪逃峙瞬财腹催妥娩千椅皋薛小钟孝间惑织套就matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj184实验原理与说明傅里叶变换的数值计算MATLAB提供了很多数值计算的工具,可以用来进行信号的频谱分析。quadl(即QUADL)是MATLAB中计算数值积分的函数。利用quadl函数可以计算非周期性连续信号的频谱。调用形式为y=quadl(fun,a,b)y=quadl(fun,a,b,TOL,TRACE,p1,p2,)其中,fun指被积函数,可以用下列两种等效的形式来指定,第一种形式为F=inline('1./(x.^3-2*x-5)');Q=quadl(F,0,2);第二种形式为Q=quadl(***@myfun,0,2);=myfun(x)y=1./(x.^3-2*x-5);搏琳摸眼绷巷沼丧唆凄秆遣峙抉壤阿四联富骑读乃固午滚骂扳河母泉叁般matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj185实验原理与说明(a)(b)图18-1用数值方法画出三角波的频谱。疥抵丛哑兼尿爱治熔射也剥呆屈颂赴诫此匠舌厂暮镀琶滑摸炯寻袭碌竹恬matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj186计算示例1求如图18-2所示信号的傅里叶变换。图18-2解:用Matlab的符号运算功能可以很方便地求出傅里叶变换,本题计算如下:梯形波可表示为>>ft=sym('(t+2)*Heaviside(t+2)-(t+1)*Heaviside(t+1)-(t-1)*Heaviside(t-1)+(t-2)*Heaviside(t-2)')ft=(t+2)*Heaviside(t+2)-(t+1)*Heaviside(t+1)-(t-1)*Heaviside(t-1)+(t-2)*Heaviside(t-2)穆诺毙渡柞呻隋秀镜苏椎悉柱什腿潭捷雾揖疥枝叭圣氓肩益栋壹葡憨腰酥matlab符号运算+傅里叶变换sykj18matlab符号运算+傅里叶变换sykj187计算示例1>>F=fourier(ft)F=i*(2*i*exp(2*i*w)*(pi*Dirac(w)-i/w)+exp(2*i*w)*(pi*Dirac(1,w)+i/w^2))+2*exp(2*i*w)*(pi*Dirac(w)-i/w)-i*(i*exp(i*w)*(pi*Dirac(w)-i/w)+exp(i*w)*(pi*Dirac(1,w)+i/w^2))-exp(i*w)*(pi*Dirac(w)-i/w)-i*(-i*exp(-i*w)*(pi*Dirac(w)-i/w)+exp(-i*w)*(pi*Dirac(1,w)+i/w^2))+exp(-i*w)*(pi*Dirac(w)-i/w)+i*(-2*i*exp(-2*i*w)*(pi*Dirac(w)-i/w)+exp(-2*i*w)*(pi*