matlab表达式随频率变化,信号与系统实验(MATLAB版)实验11连续LTI系统的频域分析.ppt...

一、实验目的　　(1) 掌握连续时间信号傅里叶变换和傅里叶逆变换的实现方法，以及傅里叶变换的时移特性、傅里叶变换的频移特性的实现方法; 　　(2) 了解傅里叶变换的特点及其应用; 　　(3) 掌握函数fourier和函数ifourier的调用格式及作用;　　(4) 掌握傅里叶变换的数值计算方法，以及绘制信号频谱图的方法。 二、实验原理　　1. 系统的频率特性　　连续LTI系统的频率特性又称频率响应特性，是指系统在正弦信号激励下稳态响应随激励信号频率的变化而变化的情况，又称系统函数H(ω)。对于一个零状态的线性系统，如图11.1所示。 　　其中，|H(ω)|随ω变化的规律称为系统的幅频特性; φ(ω)随ω变化的规律称为系统的相频特性。　 频率特性不仅可用函数表达式表示，还可用随频率f(或ω)变化的曲线来描述，如图11.2所示。　　当频率特性曲线采用对数坐标描述时，又称为波特图，如图11.3所示。  三、 涉及的MATLAB函数　　1. fourier函数　　功能：实现信号f(t)的傅里叶变换。　　调用格式：　　F=fourier(f)：是符号函数f的傅里叶变换，默认返回函数F是关于ω的函数。　　F=fourier(f，v)：是符号函数f的傅里叶变换，返回函数F是关于v的函数。　　F=fourier(f，u，v)：是关于u的函数f的傅里叶变换，返回函数F是关于v的函数 　　2. ifourier函数　　功能：实现信号F(jω)的傅里叶逆变换。　　调用格式：　　f=ifourier(F)：是函数F的傅里叶逆变换，默认的独立变量为ω，默认返回是关于x的函数。　　f=ifourier(F，u)：返回函数f是u的函数，而不是默认的x的函数。　　f=ifourier(F，v，u)：是对关于v的函数F进行傅里叶逆变换，返回关于u的函数f。 四、 实验内容与方法　　1. 验证性实验(参考程序)　　1) 编程实现信号的傅里叶变换和傅里叶逆变换　　(1) 傅里叶变换。　　已知连续时间信号f(t)=e-2|t|，通过程序完成信号f(t)的傅里叶变换 MATLAB程序： syms t;  f=fourier(exp(-2*abs(t)));  ezplot(f); 信号[WTBX]f(t)的傅里叶变换如图11.4所示。 　　(3) 傅里叶变换数值计算。　　已知门函数f(t)=g2(t)=ε(t+1)-ε(t-1)，试采用数值计算方法确定信号的傅里叶变换F(jω)。MATLAB程序： R=0.02; t=-2：R：2; f=stepfun(t，-1)-stepfun(t，1); W1=2*pi*5; %频率宽度N=500; %采样数为N k=0：N; W=k*W1/N; %W为频率正半轴的采样点F=f*exp(-j*t′*W)*R; %求F(jw)F=real(F); W=［-fliplr(W)，W(2：501)］; F=［fliplr(F)，F(2：501)］; subplot(2，1，1); plot(t，f); xlabel(′t′); ylabel(′f(t)′); axis(［-2，2，-0.5，2］); title(′f(t)=u(t+1)-u(t-1)′); subplot(2，1，2); plot(W，F); xlabel(′w′); ylabel(′F(w)′); title(′f(t)的傅里叶变换′); 信号的傅里叶变换如图11.6所示。 (4) 连续函数的傅里叶变换。MATLAB程序：clf;dt=2*pi/8; w=linspace(-2*pi，2*pi，2000)/dt; k=-2：2; f=ones(1，5); F=f*exp(-j*k′*w); f1=abs(F); plot(w，f1); grid; 连续函数的傅里叶变换如图11.7所示。 (5) 连续周期信号的傅里叶级数。MATLAB程序：clf; %计算连续周期信号的傅里叶级数 N=8; n1=- N：-1; %计算N为负数时的傅里叶级数c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2; c0=0; %计算N为零时的傅里叶级数 n2=1：N; %计算为N正数时的傅里叶级数c2=-4*