信号与系统上机实验1.docVIP

实验四 傅里叶变换(FT)及其性质 一、实验目的 1、学会运用Matlab求连续时间信号的傅里叶 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 二、实验原理及实例分析 （一）傅里叶变换的实现 例1：用Matlab符号运算求解法求单边指数信号的FT。 例2：用Matlab符号运算求解法求的IFT。 例3：用Matlab命令绘出例1中单边指数数信号的频谱图。 例4：用Matlab命令求图示三角脉冲的FT，并画出其幅度谱。 例5：用Matlab数值计算法求例3的三角脉冲幅度频谱图。 （二）FT的性质 1、尺度变换 例6：设矩形信号，利用Matlab命令绘出该信号及其频谱图。同时绘出的频谱图，并加以比较。 下面利用Matlab将常规矩形脉冲信号的频谱和其调制信号（课本例3-4信号）频谱进行比较。Matlab源程序如下： 傅里叶变换的其它性质可用类似的方法验证，希望大家课下练习。 三、实验内容 [注意：(1)写代码时j(i] 3、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图。 4、已知门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。 四、实验报告要求 实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。 五、开发环境 Matlab 7.0 六、实验内容： 1．（1） ft = sym(sin(2*pi*(t-1))/(pi*(t-1))); Fw = fourier(ft); subplot(2,2,1); ezplot(abs(Fw)); grid on; title(幅度谱); phase = atan(imag(Fw)/real(Fw)); subplot(2,1,2); ezplot(phase); grid on; title(相位谱); 2．（1） syms t; Fw = sym(10/(3+i*w) - 4/(5+i*w)); ft = ifourier(Fw,t); ezplot(ft,t); 3. f=sym((t-2)*(Heaviside(t+2)-Heaviside(t+1))+Heaviside(t+1)-Heaviside(t-1)+(-t+2)*(Heaviside(t-1)-Heaviside(t-2))); Fw = simplify(fourier(ft)); subplot(2,1,1); ezplot(abs(Fw)); grid on; title(幅度谱); phase = atan(imag(Fw)/real(Fw)); subplot(2,1,2); ezplot(phase); grid on; title(相位谱);