求大神解释下段MATLAB程序并注释一下,急求!
答案:1 悬赏:40 手机版
解决时间 2021-03-23 11:28
- 提问者网友:相思似海深
- 2021-03-23 01:36
求大神解释下段MATLAB程序并注释一下,急求!
最佳答案
- 五星知识达人网友:鸽屿
- 2021-03-23 02:12
clear all;
close all;
fs=2500; %采样率
fp=99,fc=105; %通带边界频率和截止频率
wp=2*pi*fp/fs; %频率做归一化处理
ws=2*pi*fc/fs;
Bt=ws-wp; %过渡带宽
N0=ceil(6.6*pi/Bt); %确定滤波器长度
N=N0+mod(N0+1,2); %保证为奇数
wc=(wp+ws)/2/pi; %过渡带中心频率
hn=fir1(N-1,wc,hamming(N)); %汉明窗滤波器
freqz(hn,1,512); %绘制频率响应图
%%构造测试信号
f1=30;f2=45;f3=80;f4=93;f5=99;f6=130; %各频率分量(Hz)
t=0:0.01:2; %时长
x=0.5*cos(2*pi*f1*t)+0.3*cos(2*pi*f2*t)+0.2*cos(2*pi*f3*t)+0.7*cos(2*f4*pi*t)+0.8*cos(2*pi*f5*t)+0.6*cos(2*pi*f6*t); %测试信号
figure(5);
subplot(211)
plot(t,x);
y=filter(hn,1,x); %时域滤波
title('滤波前信号');
subplot(212)
plot(t,y);
title('滤波后信号');
close all;
fs=2500; %采样率
fp=99,fc=105; %通带边界频率和截止频率
wp=2*pi*fp/fs; %频率做归一化处理
ws=2*pi*fc/fs;
Bt=ws-wp; %过渡带宽
N0=ceil(6.6*pi/Bt); %确定滤波器长度
N=N0+mod(N0+1,2); %保证为奇数
wc=(wp+ws)/2/pi; %过渡带中心频率
hn=fir1(N-1,wc,hamming(N)); %汉明窗滤波器
freqz(hn,1,512); %绘制频率响应图
%%构造测试信号
f1=30;f2=45;f3=80;f4=93;f5=99;f6=130; %各频率分量(Hz)
t=0:0.01:2; %时长
x=0.5*cos(2*pi*f1*t)+0.3*cos(2*pi*f2*t)+0.2*cos(2*pi*f3*t)+0.7*cos(2*f4*pi*t)+0.8*cos(2*pi*f5*t)+0.6*cos(2*pi*f6*t); %测试信号
figure(5);
subplot(211)
plot(t,x);
y=filter(hn,1,x); %时域滤波
title('滤波前信号');
subplot(212)
plot(t,y);
title('滤波后信号');
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯