DB design log (10)

May 16

The last day. All personal work has been finished and we conbined the administrator and customer operation parts. Solve some small problem which were caused by the trivial database problem.

In the afternoon, we displayed our work to the tutor. Well, the whole developing process finished.

下面是一个简单的 MATLAB 程序,用于实现脉冲响应不变法设计 IIR 数字低通滤波器。 ```matlab % 清空工作区 clear all; % 设置通带和阻带参数 wp = 0.4*pi; % 通带边缘频率 ws = 0.35*pi; % 阻带边缘频率 Rp = 1; % 通带最大衰减 As = 10; % 阻带最小衰减 % 计算通带截止频率和阻带截止频率 wc = (wp + ws) / 2; % 通带截止频率 delta_w = wp - ws; % 过渡带宽度 delta_p = (10^(0.1*Rp)-1)/(10^(0.1*Rp)+1); % 通带最大波纹 delta_s = 10^(-0.1*As); % 阻带最小衰减 % 计算 IIR 滤波器的阶数和截止频率 N = ceil(0.5*log10((1/delta_s^2-1)/(1/delta_p^2-1))/log10(wc/(delta_w/2))); % 滤波器阶数 wc_1 = wc/(1+delta_p)^(-1/(2*N)); % 截止频率 % 计算 IIR 滤波器的极点和零点 k = 1:N; p = wc_1*cos(pi/2/N*(2*k-1))+1j*wc_1*sin(pi/2/N*(2*k-1)); % 极点 z = -ones(1, N); % 零点 % 计算数字低通滤波器的传递函数 [b, a] = zp2tf(z, p, 1); % 转换为分子和分母系数 Hd = freqz(b, a, 512); % 频率响应 % 绘制数字低通滤波器的幅频特性 w = linspace(0, pi, 512); mag = 20*log10(abs(Hd)); plot(w/pi, mag); grid on; xlabel('Normalized Frequency (\times\pi rad/sample)'); ylabel('Magnitude (dB)'); title('IIR Digital Lowpass Filter Design using Impulse Invariance Method'); axis([0 1 -60 10]); ``` 上述程序中,我们首先设置了所需的通带和阻带参数,然后计算出滤波器的阶数和截止频率。接下来,我们计算滤波器的极点和零点,并将其转换为分子和分母系数。最后,我们绘制了数字低通滤波器的幅频特性。 运行程序后,会得到以下的 IIR 数字低通滤波器的幅频特性图: ![IIR数字低通滤波器幅频特性图](https://img-blog.csdn.net/2018052216194825?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvYm9vazIwMTU=/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/q/80)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值