信号与系统实验二 下载本文

内容发布更新时间 : 2024/12/28 12:24:18星期一 下面是文章的全部内容请认真阅读。

图片如图2所示。

图5

(3)高斯信号f3?e?t

答:程序:

程序如下所示。

syms t w;

f3=sym('exp(-t^2)'); F3=fourier(f3); FFP3=abs(F3);

ezplot(FFP3,[-10*pi 10*pi]); ylabel('|F(j\\omega)|'); grid on;

axis([-15 15 0 1.8]);

图片:

图片如图2所示。

2

图6

题目三:利用ifourier()函数求下列频谱函数的傅式反变换

2w(jw)2?5jw?8(1)F(jw)??j(2)F(jw)?

16?w2(jw)2?6jw?5

答:程序:

程序如下所示。

syms t w;

Fw1=sym('-j*2*w/(16+w^2)'); ft1=ifourier(Fw1,w,t)

syms t w;

Fw2=sym('(j*w)^2+5*j*w-8/((j*w)^2+6*j*w+5)'); ft2=ifourier(Fw2,w,t)

运行结果:

ft1 =

(j*(pi*heaviside(-t)*exp(4*t)*i+ (pi*exp(-4*t)*dirac(t)*i)/4- (pi*exp(4*t)*dirac(t)*i)/4-pi*exp(-4*t)*heaviside(t)*i))/pi

ft2 =

-(2*pi*j^2*dirac(t, 2) - (2*(pi*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*i - pi*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*i + pi*sign(3*imag(1/j) - 2*imag((1/j^2)^(1/2)))*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*i - pi*sign(3*imag(1/j) + 2*imag((1/j^2)^(1/2)))*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*i -

pi*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*heaviside(t)*2*i + pi*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*heaviside(t)*2*i))/(j^2*(1/j^2)^(1/2)) + pi*j*dirac(t, 1)*10*i)/(2*pi)

题目四:设H(jw)?1,试用MATLAB画出该系统的幅频20.08(jw)?0.4jw?1特性|H(jw)|和相频特性?(w)。

答:程序:

程序如下所示。

w=0:0.001:20; b=[1];

a=[0.08,0.4,1]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H)); grid on;

xlabel('\\omega(rad/s)'); ylabel('|H(j\\omega)|'); title('H(jw)的幅频特性'); subplot(2,1,2); plot(w,angle(H)); grid on;

xlabel('\\omega(rad/s)'); ylabel('\\phi(\\omega)'); title('H(jw)的相频特性');

图片:

图片如图2所示。

图7