《信号与系统》 实验三 连续信号的频域特性分析

参考

信号实验指导.pdf
信号实验三和实验四.pdf
《信号与系统》4.11.2系统函数的幅频特性和相频特性分析 MATLAB实现

题目一门函数的傅里叶变换

题目

求门]函数

$$g_4\left(t\right)$$
的傅立叶变换，并画出其频率特性曲线图。

code1 sf1函数的定义

function y=sf1(t,w)
y = (heaviside(t+2)-heaviside(t-2)).*exp(-1i*w*t);
end

注意放到同一文件夹比较好

code2

w = linspace(-6*pi,6*pi,512);
N=length(w);
F=zeros(1,N); %生成一个1行N列的全零矩阵
for k=1:N
    F(k)=quadl('sf1',-1, 2, [],[], w(k));
end
figure(1);
plot(w,real(F));
xlabel('\omega');
ylabel('F(j\omega)');

结果

题目二 频率特性函数幅频特性和相频特性

题目

已知频率特性函数为:
$$H\left(jw\right)=\frac{2{\left(jw\right)}^3+(jw)+4}{{(jw)}^4+3{(jw)}^3+2{(jw)}^2+5(jw)+2}$$
求其幅频特性和相频特性。

code

w = linspace(0,5,200)
b = [2 0 1 4];
a = [1 3 2 5 2];
H = freqs(b,a,w);
subplot(2,1,1);
plot(w,abs(H));
set(gca,'xtick',[0 1 2 3 4 5]);
set(gca,'ytick',[0 0.4 0.707 1]);
xlabel('\omega')
ylabel('|H(j\omega)|');
subplot(2,1,2);
plot(w,angle(H));
set(gca,'xtick',[0 1 2 3 4 5]);
xlabel('\omega')
ylabel('\phi(\omega)');

结果