38. 已知系统的零极点分别如下：
z1?2.2,z2??1?j,z3??1?j,z4?1.4p1?3.7?j2,p2?3.7?j2,p3??2.1?j,p4??2.1?j
求系统的系统函数H(z)。
39. 设计椭圆IIR数字低通滤波器，其性能指标为：通带截止频率为1000Hz，阻带截止频率为1250Hz，通带波纹为0.4dB，最小阻带衰减为45dB，抽样频率为5000Hz。绘制所设计的滤波器增益响应。
40. 编写总体均值滤波器程序。原始未受干扰的序列为：s[n]=3[n(0.8)n], 加性噪声信号d[n]为随机序列，幅度0.6，受干扰的序列为：x[n]= s[n]+ d[n]，绘制噪声序列和60次检测结果的总体平均的序列。
d = rand(R,1)-0.5; b = ones(M,1)/M;
b = round(a); % Generate impulse noise y = medfilt1(x,3); % Median filtering y = filter(b,1,x); axis([0 50 -1 8]); c = conv(a, b);
input('Denominator coefficients = '); h = freqz(num, den, w); real(h) imag(h) abs(h) angle(h)
[z,p,k] = buttap(4)
[num,den] = butter(N,Wn,'s');
[num,den] = cheby1(N,Rp,2*pi*Fp,'s'); [num,den] = ellip(N,Rp,Rs,2*pi*Fp,'s'); U = fft(u,M); zplane(num,den)
[num, den] = zp2tf(z, p, k);z,p 是列向量 [num, den] = residuez(r,p,k); [r,p,k] = residuez(num,den);
[y,t] = impz(num,den,L);
disp('Coefficients of the power series expansion');