应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 急!!!麻烦大家帮我注释下这个matlab滤波器的程序,实在看不懂了...
11/1返回列表
查看: 468  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

pobm

新虫 (初入文坛)

[求助] 急!!!麻烦大家帮我注释下这个matlab滤波器的程序,实在看不懂了...

%**************     滤波器 *******************************%
%    Wp=0.2* pi; Ws=pi/2; Rs= 40; B=Ws- Wp;
%    beta= 0.5842*(Rs-21)^0.4 + 0.07886 * (Rs- 21);
%    N= ceil((Rs-8)/2.285/B);
%    Wc= (Wp+ Ws)/2/pi;
%    h=fir1(N,Wc,kaiser(N+1,beta));
%   
  
%xt信号的频域信号
  x2F=fft(x2,length(x2));
   %%%%%%%%滤波器设计
   fp=3000;fst=5550;fs=44100;
   wp=2*fp/fs;
   ws=2*fst/fs;
   deltaw=ws-wp;
   N0=ceil(6.6/deltaw);
   N=N0+mod(N0+1,2);
   n=N-1;
   wn=(ws+wp)/2;
   b=fir1(n,wn);
   
   [H,w]=freqz(b,1,length(x2F),'whole');
   mag=abs(H);
   n=0:length(mag)-1;
   fl=(0:length(x2F)-1)'*44100/length(x2F);
   figure(11)
   plot(fl,mag,'r-');
    title('滤波器幅频响应');
    xlabel('f'),ylabel('H');
    hold on;
    grid on
   
   
%%%通过滤波器之后的频域信号
for i=1:length(x2F)
    signalF(i)=x2F(i)*H(i);
end
fl=(0:length(x2F)-1)'*44100/length(x2F);
figure(12)
mag2=abs(signalF);
plot(fl,mag2,'r-');
title('信号滤波后幅频响应');
xlabel('f'),ylabel('x1’(f)');
hold on;
grid on

%%%%通过滤波器之后的时域信号
signalT=ifft(signalF,N1);
figure(13)
plot(t,signalT,'r-');
title('信号滤波后时域信号');
xlabel('f'),ylabel(' x1’(t)');
hold on;
grid on
回复此楼

» 猜你喜欢
1楼 2011-10-30 20:45:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 pobm 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭