| 查看: 441 | 回复: 0 | ||
[求助]
急!!!麻烦大家帮我注释下这个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 |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有133人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
求matlab求解一个方程组的问题
已经有6人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
- MATLAB 编程求解普通实验数据时间序列的最大李雅普诺夫指数！！！
已经有12人回复
- c语言函数（递归调用等） doc.
已经有10人回复
- William T. Freeman的基于例子超分辨率matlab程序跑不起来啊！
已经有6人回复
- MATLAB 画图小问题
已经有11人回复
- 求助：dushin程序输入文件无法读取！！！1
已经有10人回复
- Matlab人工神经网络工具箱的位置及使用
已经有5人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
