设计一个椭圆带通滤波器，它的指标如下： %椭圆低通滤波器的设计2113[b，a]=ellip(5，0.1，40，5*2/Fs)；关键是ellip这个函数的5261使用，第一个参数表示4102滤波器的阶数，第二个参数表示通带1653波纹度，单位dB，第三个参数表示阻带衰减，单位dB，第四个参数表示归一化的低通滤波器截止频率，计算公式如下，假设截止频率为F，则有：第四个参数=F/(Fs/2)=F*2/Fs由上面的描述可知，该滤波器的截止频率为5Hz，所以会把15和30Hz滤掉，不过我觉得截止频率设为10Hz会更好一点，你说呢。剩下的比较简单，有不懂的欢迎追问。
这个matlab程序怎么编？（设计低通滤波器） 1.1 实验目的21131．了解数字信号处理系统的一般构成5261；2．掌握奈奎4102斯特抽样定理。1.2 实验1653仪器1．YBLD智能综合信号源测试仪 1台2．双踪示波器 1台3．MCOM－TG305数字信号处理与现代通信技术实验箱 1台4．PC机（装有MATLAB、MCOM－TG305配套实验软件）1台1.3 实验原理一个典型的DSP系统除了数字信号处理部分外，还包括A/D和D/A两部分。这是因为自然界的信号，如声音、图像等大多是模拟信号，因此需要将其数字化后进行数字信号处理，模拟信号的数字化即称为A/D转换。数字信号处理后的数据可能需还原为模拟信号，这就需要进行D/A转换。一个仅包括A/D和D/A两部分的简化数字信号处理系统功能如图1所示。A/D转换包括三个紧密相关的过程，即抽样、量化和编码。A/D转换中需解决的以下几个重要问题：抽样后输出信号中还有没有原始信号的信息？如果有能不能把它取出来？抽样频率应该如何选择？奈奎斯特抽样定理（即低通信号的均匀抽样定理）告诉我们，一个频带限制在0至fx以内的低通信号x(t)，如果以fs≥2fx的抽样速率进行均匀抽样，则x(t)可以由抽样后的信号xs(t)完全地确定，即xs(t)包含有x(t)的成分，可以通过适当的低通滤波器不失真地恢复出x(t)。最小。
高分求：用双线性变换法设计原型低通为椭圆型的数字IIR带通滤波器 p1=400；p2=500；通带边界频率s1=350；s2=550；阻带截止频率Ap=1；通带最大衰减As=40；阻带最小衰减Ft=2000；抽样频率T=2；wp1=2*pi*p1/Ft；wp2=2*pi*p2/Ft；ws1=2*pi*s1/Ft；ws2=2*pi*s2/Ft；Wp1=(2/T)*tan(wp1/2)；Wp2=(2/T)*tan(wp2/2)；Ws1=(2/T)*tan(ws1/2)；Ws2=(2/T)*tan(ws1/2)W0=Wp1*Wp2；w0=sqrt(W0)；BW=Wp2-Wp1；带通滤波器的通带宽度lp=1；归一化处理ls=Ws1*BW/(W0-Ws1^2)；[N，Wn]=ellipord(lp，ls，Ap，As，'s')；[B，A]=ellip(N，1，40，Wn，'s')；[BT，AT]=lp2bp(B，A，w0，BW)；[num，den]=bilinear(BT，AT，0.5)；[z，p，k]=tf2zp(num，den)；figure(1)；zplane(z，p)；title('零极点')[h，w]=freqz(num，den，512)；figure(2)plot(w/pi，20*log10(abs(h)))；axis([0 1-100 1])；title('频谱特性曲线')gridn=0：800；k=n/8000；通过滤波器3f1=2*pi*450；f2=2*pi*6000；x=sin(f1*k)+sin(f2*k)；y=filter(num，den，x)；x1=sin(f1*k)figure(3)plot(x1)；x1图形输出axis([0，100*pi，-5，5])；title('x1(t)')；x2=sin(f2*k)；figure(4)plot(x2)；x2图形输出axis([0，100*pi，-5，5])；title('x2')；figure(5)plot(x)；axis([0，100*pi，-5，5])；title('输入信号')；figure(6)plot(y)；axis([0，100*pi，-5，5])；title('输出信号')。
