谁能帮我讲解下这些MATLAB语言的意思.题目及语言如下已知模拟滤波器的系统函数为 ,分别用冲激响应不变法和双线性变换法将 转化为数字滤波器 ,并画出 和 的频率响应曲线.抽样频率分别为1

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/04 13:06:52
谁能帮我讲解下这些MATLAB语言的意思.题目及语言如下已知模拟滤波器的系统函数为 ,分别用冲激响应不变法和双线性变换法将 转化为数字滤波器 ,并画出 和 的频率响应曲线.抽样频率分别为1

谁能帮我讲解下这些MATLAB语言的意思.题目及语言如下已知模拟滤波器的系统函数为 ,分别用冲激响应不变法和双线性变换法将 转化为数字滤波器 ,并画出 和 的频率响应曲线.抽样频率分别为1
谁能帮我讲解下这些MATLAB语言的意思.题目及语言如下
已知模拟滤波器的系统函数为 ,分别用冲激响应不变法和双线性变换法将 转化为数字滤波器 ,并画出 和 的频率响应曲线.抽样频率分别为1000Hz和500Hz.
其程序如下,
clear;close all
b=1000;a=[1,1000];
w=[0:1000*2*pi];
[hf,w]=freqs(b,a,w);
subplot(2,3,1)
plot(w/2/pi,abs(hf));;title('模拟滤波器的幅频特性');grid on;
fs0=[1000,500];
for m=1:2
fs=fs0(m)
[d,c]=impinvar(b,a,fs)
[f,e]=bilinear(b,a,fs)
wd=[0:512]*pi/512;
hw1=freqz(d,c,wd);
hw2=freqz(f,e,wd);
subplot(2,3,2);plot(wd/pi,abs(hw1)/abs(hw1(1)));
hold on;grid on;title('脉冲响应不变法');
subplot(2,3,3);plot(wd/pi,abs(hw2)/abs(hw2(1)));hold on;title('双线性变化法');
end

谁能帮我讲解下这些MATLAB语言的意思.题目及语言如下已知模拟滤波器的系统函数为 ,分别用冲激响应不变法和双线性变换法将 转化为数字滤波器 ,并画出 和 的频率响应曲线.抽样频率分别为1
clear; %清除变量
close all %关闭所有窗口
b=1000;a=[1,1000]; %模拟滤波器分子分母的系数
w=[0:1000*2*pi]; %定义频带宽度
[hf,w]=freqs(b,a,w); %把频域转化到复频域
subplot(2,3,1) %分割窗口画图
plot(w/2/pi,abs(hf)); %画滤波器的幅频特性
title('模拟滤波器的幅频特性') %给这个图加个标题
grid on %打开网格
fs0=[1000,500]; %定义两个变化频率
for m=1:2
fs=fs0(m)
[d,c]=impinvar(b,a,fs) %脉冲响应不变法求滤波器系数
[f,e]=bilinear(b,a,fs) %双线性变化法求滤波器系数
wd=[0:512]*pi/512; %频率归一化
hw1=freqz(d,c,wd); %求脉冲响应不变法的数字滤波器的频率分量
hw2=freqz(f,e,wd); %求双线性变化法的数字滤波器的频率分量
subplot(2,3,2); %分割窗口画图
plot(wd/pi,abs(hw1)/abs(hw1(1))); %画滤波器的幅频特性
hold on %保持图形不清除
grid on %打开网格
title('脉冲响应不变法'); %给这个图加个标题
subplot(2,3,3)%分割窗口画图
plot(wd/pi,abs(hw2)/abs(hw2(1))) %画滤波器的幅频特性
hold on; %保持图形不清除
title('双线性变化法'); %给这个图加个标题
end