| 查看: 448 | 回复: 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 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有263人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20