数字信号处理(1) 代码
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
H=20*log10(abs(h)); %输出幅值转化为分贝单位
subplot(311); %把图像放到311位置
plot(t,x); %画出输入信号图像
[b,a]=butter(N,Wn,'s'); %巴特沃斯模拟滤波器设计
[b1,a1]=bilinear(b,a,0.5); %用双线性法，将模拟滤波器转换为数字滤波器
[h,omega]=freqz(b1,a1,256); %系统的绝对幅度响应（h包含了离散系统频响在 0~pi范围内256个频率等分点的值，omega则包含了范围内256个频率等分点）
ylabel('H (dB)'); %Y轴取值H
title('滤波器图像'); %图像题目
y=filter(b,a,x); %将信号进行滤波处理
title('滤波后的图像'); %图像题目
Rp=1; %幅度特性下降小于1DB
Rs=15; %幅度特性下降小于15DB
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); %巴特沃斯滤波器N估计
subplot(313); %把图像放到313位置
pl %网格布线
grid; %网格布线
title('输入信号'); %图像题目
subplot(312); %把图像放到312位置
x=sin(2*pi*20*t)+2*sin(2*pi*100*t)+5*sin(2*pi*200*t); %输入信号
Wp=tan(pi*30/Fs); %预畸变处理,Wp通带截至频率
Ws=tan(pi*60/Fs); %预畸变处理,Ws阻带截至频率
%抽样频率500，要求在频率小于30的通带内，幅度特性下降小于1DB，在频率大于60的阻带内，衰减大于15DB
Fs=500; %抽样频率
t=0:1/Fs:1; %采样周期
plot(omega,H); %画出滤波器示意图
grid; %网格布线
xlabel('\omega'); %X轴取值omega
subplot(311); %把图像放到311位置
plot(t,x); %画出输入信号图像
[b,a]=butter(N,Wn,'s'); %巴特沃斯模拟滤波器设计
[b1,a1]=bilinear(b,a,0.5); %用双线性法，将模拟滤波器转换为数字滤波器
[h,omega]=freqz(b1,a1,256); %系统的绝对幅度响应（h包含了离散系统频响在 0~pi范围内256个频率等分点的值，omega则包含了范围内256个频率等分点）
ylabel('H (dB)'); %Y轴取值H
title('滤波器图像'); %图像题目
y=filter(b,a,x); %将信号进行滤波处理
title('滤波后的图像'); %图像题目
Rp=1; %幅度特性下降小于1DB
Rs=15; %幅度特性下降小于15DB
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); %巴特沃斯滤波器N估计
subplot(313); %把图像放到313位置
pl %网格布线
grid; %网格布线
title('输入信号'); %图像题目
subplot(312); %把图像放到312位置
x=sin(2*pi*20*t)+2*sin(2*pi*100*t)+5*sin(2*pi*200*t); %输入信号
Wp=tan(pi*30/Fs); %预畸变处理,Wp通带截至频率
Ws=tan(pi*60/Fs); %预畸变处理,Ws阻带截至频率
%抽样频率500，要求在频率小于30的通带内，幅度特性下降小于1DB，在频率大于60的阻带内，衰减大于15DB
Fs=500; %抽样频率
t=0:1/Fs:1; %采样周期
plot(omega,H); %画出滤波器示意图
grid; %网格布线
xlabel('\omega'); %X轴取值omega