内容发布更新时间 : 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