| 查看: 416 | 回复: 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 |
回复此楼» 猜你喜欢
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有22人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
