非周期信号的频谱课件.pptVIP

实验18 非周期信号的频谱 实验目的 学习Matlab的符号运算功能中计算傅里叶变换和反变换的方法。 掌握用Matlab绘制非周期信号频谱的数值方法和符号方法。 通过对非周期信号频谱的绘制。加深对傅里叶变换性质的理解。 实验原理与说明 傅里叶变换和反变换的符号运算 实验原理与说明 实验原理与说明 傅里叶变换的数值计算 实验原理与说明 用数值方法画出三角波的频谱。 计算示例1 计算示例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*Dirac(1,w)+i/w^2))-2*exp(-2*i*w)*(pi*Dirac(w)-i/w) F1=simple(F) F1 = -2*(cos(2*w)-cos(w))/w^2 计算示例3 计算示例4 实验内容 考虑非周期连续时间信号如门函数或单边指数，将其左移和右移后观察其频谱的变化。验证傅里叶变换的时移性质。 考虑非周期连续时间信号如门函数或三角波，将其用余弦调制。观察其频谱的变化。验证傅里叶变换的调制性质。 考虑非周期连续时间信号如门函数或单边指数，将其进行尺度变换后观察其频谱的变化。验证傅里叶变换的尺度变换性质。 实验内容 求下列信号的傅里叶变换。画出 和 的幅度谱和相位谱，并进行比较。（令A=2，a=3） 实验内容 实验内容 实验步骤与方法 由于经常在画如图18-2所示的时间波形、幅度频谱、相位频谱，可将exp18_2.m程序改编成函数。这样只要在画频谱图时调用该函数就可以了。 仿照例3的方法做实验1、2、3的内容。加深对傅里叶变换性质的理解。 仿照例4的方法，完成实验内容4的编程。上机调试程序，根据题目要求对幅度谱和相位谱加以比较。 仿照例1、2的符号计算方法，完成实验内容5、6的计算。并与理论计算进行比较。 实验报告要求 根据实验内容编写出的程序。以及绘出的各种波形图。对各种频谱图加以比较说明。 上机调试程序的方法。 根据实验归纳、总结出用Matlab绘制信号频谱图的方法。 心得体会及其他。 * * MATLAB的符号运算工具箱中，专门提供了傅里叶变换和反变换的函数。 正变换的调用格式为 F=fourier(f) 式中，f为时间函数的符号表达式，F为傅里叶变换式，也是符号表达式。 反变换的调用格式为 f=ifourier(F) 式中，F为傅里叶变换式的符号表达式，f为时间函数，是符号形式。 为了改善公式的可读性，MATLAB提供了pretty函数，调用格式为 Pretty(f) 式中，f为符号表达式。 如已知 ，求其频谱 。 syms t w 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 + i w f1=ifourier(F,t) f1 = exp(-2*t)*Heaviside(t) MATLAB提供了很多数值计算的工具，可以用来进行信号的频谱分析。quadl(即QUADL)是MATLAB中计算数值积分的函数。利用quadl函数可以计算非周期性连续信号的频谱。调