利用MATLAB结合双线性变换法设计一个IIR数字butterworth低通滤波器

麻烦请大大们请帮忙写一个MATLAB实现的设计方案。最近学校急着让交……

技术指标为：通带截止频率fp=1KHz,阻带截止频率fs=1.5KHz，通带衰减Rp=1.5dB,阻带衰减Rs=25dB，采样频率Fs=10kHz。
需要绘出滤波器的幅频特性曲线，判断设计是否符合要求。