②y(k+2)-0.7y(k+1)+0.1y(k)=7f(k+2)-2f(k+1);f(k)= ε(t)
单位响应10edutilp5mA 0024681012141618205n (samples)冲激响应4x 10edutlip2mA002468101214161820零状Tim态e (sec)响应400200002468101214161820
三、实验心得
1.通过实验了解了连续时间系统的单位冲激响应，阶跃响应和零状态响应的意义和求解方法。 2.通过实验了解了离散时间系统的单位冲激响应，阶跃响应和零状态响应的求解方法。
3.了解了用MATLAB求解系统响应的方法。
4.知道了冲激响应的调用函数是impluse(),阶跃响应的调用函数是step(),零状态响应调用函数是Isim().
36
实验八、MATLAB在信号与系统的应用
----连续时间信号的频谱分析
一、实验目的
1.熟悉傅里叶变换的定义与性质。 2.熟悉常见信号的傅里叶变换。
3.了解傅里叶变换的MATLAB实现方法。
二、实验内容及实验结果
1.编成实现求下列信号的幅度频谱。
（1）求出f1(t)=ε(2t+1)-ε(2t-1)的频谱函数F1（jw）,请将它与上面门宽为2的门函数f(t)=ε(t+1)-ε(t-1)的频谱进行比较，观察两者的特点，说明两者的关系。 >> syms t w
>> Gt=sym('Heaviside(t+1)-Heaviside(t-1)'); >> Fw=simplify(fourier(Gt,t,w));
>> subplot(211),ezplot(Gt,[-2 2]);grid;
>> subplot(212),ezplot(Fw,[-10*pi 10*pi]);grid; >> axis([-10*pi 10*pi 0 2.2])
Heaviside(2 t+1)-Heaviside(2 t-1)10.50-2-1.5-1-0.500.5t2/w sin(1/2 w)11.5221.510.50-30-20-100w102030
37
Heaviside(t+1)-Heaviside(t-1)10.50-2-1.5-1-0.500.511.52t2/w sin(w)21.510.50-30-20-100102030w（2）三角脉冲f2(t)={1-|t| |t|<=1
{0 |t|>1 syms t w
>> Gt=sym('(1-abs(x))*(Heaviside(t+1)-Heaviside(t-1))'); >> Fw=fourier(Gt,t,w)；
>> subplot(211),ezplot(Gt,[-2 2])
>> subplot(212),ezplot(Fw,[-10*pi 10*pi])
(1-abs(x)) (Heaviside(t+1)-Heaviside(t-1)) = 021x0-1-2-2-1.5-1-0.500.511.52t2/w sin(w) (1-abs(x)) = 020x0-20-30-20-100102030w(3)单边指数信号f3(t)=exp(-t)ε(t) >> syms t w
>> Gt=sym('exp(-t)*Heaviside(t)') Gt =
exp(-t)*Heaviside(t)
>> Fw=simplify(fourier(Gt,t,w))
38
Fw =
1/(1+i*w)
(4)高斯信号f4(t)=exp(-t^2) >> syms t w
>> Gt=sym('exp(-1*t^2)') Gt = exp(-1*t^2)
>> Fw=fourier(Gt,t,w) Fw =
pi^(1/2)*exp(-1/4*w^2)
2.利用ifourier()函数求下列频谱函数的傅里叶反变换。 （1）F（jw）=-j(2w/(16+w^2)) (2)F(jw)=((jw)^2+5jw-8)/((jw)^2+6jw+5)
三、实验心得
1.通过实验了解了傅里叶变换的定义及性质。 2.通过实验了解了常见的傅里叶变换。 3.通过实验了解了傅里叶变换的方法。
4.通过实验知道了在MATLAB里实现傅里叶变换的函数为fourier(f),傅里叶反变换为ifourier(f)。
39
《MATLAB实践》实验报告
学院：电子信息学院专业：电子信息工程学号：
姓名：沈丽沙
0940302104
40