| 查看: 378 | 回复: 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 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有201人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
求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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人