信号系统实验2.docVIP

一、实验目的与要求： 熟悉连续非周期信号频谱特点及其分析方法； 掌握用MATLAB实现傅立叶变换的两种方法； 了解常用傅立叶变换性质的MATLAB实现方法； 二、实验内容、步骤 (一)实验内容 通过调用傅立叶变换函数fourier()，求出矩形脉冲信号的傅立叶幅度频谱，并画出信号的波形图和傅立叶幅度频谱。 设有一带开关的RC电路如图4-1所示，电容C=0.01F，电阻R=100(，开关未闭合之前，电容器上已充电到1V，试求当开关S闭合之后，电阻上电压的变化规律，并用MATLAB画出电压信号的波形及其幅度频谱。 图4－1 设一矩形脉冲，载波信号，试用傅立叶变换的数值解法实现调幅信号，并绘制出、及它们各自的频谱。 (二)实验步骤 （1）求出矩形脉冲信号的傅立叶幅度频谱 调用傅立叶变换函数fourier()，矩形脉冲信号。 syms t w; f=heaviside(t+1/2)-heaviside(t-1/2); // heaviside(t)阶跃信号生成函数 F=fourier(f); // 调用傅立叶变换函数 subplot(3,1,1); ezplot(f,[-5,5]); //绘图命令，绘制f(t)矩形脉冲信号 axis([-3,3,-2,2]); subplot(3,1,2); ezplot(F,[-20,20]); //绘图命令，绘制F(t)——矩形脉冲傅立叶变换函数 axis([-20,20,-1.1,1.1]); p=abs(F); // p=|F(t)| subplot(3,1,3); ezplot(p,[-20,20]); //绘图命令，绘制|F(t)| axis([-20,20,-1.1,1.1]); （2）RC电路输出电压信号的波形及幅度频谱 由图4—1分析RC电路，因为电容C=0.01F，电阻R=100(，开关未闭合之前，电容器上已充电到1V。 有i(t) =C*(dUC(t)/dt)，由回路定理，R* i(t)=-UC(t)， 所以有UC(t)+ RC*(dUC(t)/dt)=0，又RC=1， 所以UC(t)+ dUC(t)/dt=0，又系统零初始状态的输入UC(0)=1， 所以可调用dsolve函数如下表示出输出电压信号的波形。 syms t w y=dsolve(Dy+y=0,y(0)=1); //y= UC(t) figure(1) y1=-y; subplot(2,1,1); ezplot(y1); y1=y1*heaviside(t); //系统输出为因果信号，有始表示 Y=fourier(y1) //调用傅立叶变换函数 G=abs(Y); subplot(2,1,2); ezplot(G); （3）用傅立叶变换的数值解法绘制调幅信号频谱 矩形脉冲，调幅信号y(t) 的载波信号，。同实验内容（1），输入信号函数，进行波形间的乘积运算，然后调用有关函数（fourier()、ezplot()等）参与运算并绘制出波形图，程序代码略。 三、结论分析与实验体会 （一）结论分析 （1）矩形脉冲信号的傅立叶幅度频谱： 上图为矩形脉冲信号f(t)波形； 中图为矩形脉冲信号f(t)的傅立叶变换频谱图； 下图为矩形脉冲信号f(t)的傅立叶幅度频谱，p=|F(t)|。 （2）RC电路输出电压信号的波形及幅度频谱： 上图为输出电压信号的波形； 下图为其幅度频谱。 （3）用傅立叶变换的数值解法绘制调幅信号频谱： 右上图为矩形脉冲波形，左上图为其傅立叶变换频谱图； 右下图为调幅信号y(t)波形，左下图为其傅立叶变换频谱图。， 回答问题：结合频谱图，对傅立叶变换的频移特性进行说明。 由欧拉定理，=e^(j10πt)+e^(-j10πt)， 由傅里叶变换的频移特性 所以原矩形脉冲f(t)的傅立叶变换频谱向左向右分别平移10π个单位得到上图结果。 （二）实验体会 这次实验是对用MATLAB仿真软件对非周期信号进行进行频谱分析，实验中对于连续时间信号我们可以调用内置函数绘制出所要分析的各种信号的频谱，我们更直观的认识对于连续时间信号傅氏变换的定义，并对其信号的频谱有了更进一步的了解； 通过实验我们更熟悉了连续非周期信号频谱特点及其分析方法，通过波形图更好的了解傅里叶变化的有关性质，这次实验就观察到其频移特性的波形展现。 R C － + S