返回顶部

收藏

数字低通滤波器(1阶)的实现

更多
function [ salida ] = pasabajos(nombre)
%Filtro pasabajos de 1er orden:

[X,Fs,bits] = wavread(nombre); %carga el wav
FX = fft(X);
Xmed=mean(X); %debe dar 0
Xvar=var(X);

pfiltro = [1 1]; %orden 1, promediador de 2 muestras: y[n]=(x[n]+x[n-1])/2

Y=1/2*filter(pfiltro,[1],X);
FY = fft(Y);

N=length(FX);
if ( mod(N,2) == 0 )
t=N/2;
else
t=floor(N/2)+1;
end
rectaFrec=linspace(0,pi*Fs/(2*pi),floor(N/2)+1);

FX = fftshift(FX);
FY = fftshift(FY);

subplot(2,1,1)
title(’Original’)
Espectro=abs(FX(t:N));
plot(rectaFrec,Espectro,’r'); %Grafico en funcion de la frecuencia en Hz

subplot(2,1,2)
title(’Pasabajos de 1er orden’)
Espectro=abs(FY(t:N));
plot(rectaFrec,Espectro,’b'); %Grafico en funcion de la frecuencia en Hz

salida = Y;
#该片段来自于http://outofmemory.cn

标签:shell,基础

收藏

0人收藏

支持

0

反对

0

»更多 您可能感兴趣的代码
  1. 2012-11-05 21:59:42java获得随机数代码 by 怪兽狂殴奥特曼
  2. 2014-03-09 17:30:52Memcached 多进程启停脚本 by 玖夏
  3. 2014-03-18 17:21:24无敌删除命令 by shixiaoliang
  4. 2014-04-21 21:41:17md5s Linux 主要源码 by Echo
  5. 2014-04-27 16:42:22yum 安装开发工具,开发库 by 屠夫
  6. 2014-05-17 10:14:31Ubuntu PPA Manager by 古殇
  7. 2014-05-18 09:54:22设置JDK环境变量 by 玖夏
  8. 2014-06-19 12:28:30nginx qps 监控 by xfk
  9. 2014-07-28 11:12:15windows batch 常用命令2 by 屠夫
  10. 2014-09-11 10:48:04Ubuntu 修改主机名/机器名 by 0晓风残月0
  11. 2018-11-05 16:27:56中文竖排 by 阮小七
相关聚客文章
  1. 博主 发表 2015-04-11 06:39:00 adb shell 基础: 看日志神器logcat
  2. 博主 发表 2015-04-11 06:10:00 adb 基础 ,重点在 adb shell ( adb basic )
  3. 博主 发表 2015-04-11 06:39:00 adb shell 基础: 看日志神器logcat
  4. 博主 发表 2015-04-11 06:10:00 adb 基础 ,重点在 adb shell ( adb basic )
  5. kingsamchen 发表 2015-06-27 09:16:20 使用cmder替换Windows自带的cmd
  6. 博主 发表 2015-07-14 02:57:17 AOS V0.2 发布,JavaEE 应用基础平台
  7. shang 发表 2015-07-31 09:27:05 linux通过shell脚本实现ftp自动更新
  8. 博主 发表 2015-08-19 00:00:00 fish 下使用 jenv
  9. DigDeeply 发表 2015-09-06 02:51:06 oh-my-zsh自动更改screen标题的问题解决
  10. 博主 发表 2015-07-15 03:24:44 常用Linux命令记录
  11. 果冻想 发表 2015-10-19 14:50:23 玩玩awk
  12. admin 发表 2015-11-03 13:55:43 写个添加svn帐户的脚本

发表评论