| 查看: 466 | 回复: 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 |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
物理学I论文润色/翻译怎么收费?
已经有139人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
200