%% 对一段音频滤去噪声%% Read Portion of Audio File
% Create a FLAC file from the example file |handel.mat|, and then read only
% the first 2 seconds.   %%
% Create a FLAC (|.flac|) file in the current folder.
load handel.matfilename = 'handel.flac';
audiowrite(filename,y,Fs);  %%
% Read only the first 2 seconds.
samples = [1,1*Fs];
clear y Fs
[y,Fs] = audioread(filename,samples);  %%
% Play the samples.
sound(y,Fs);
pause(1.5);%% 画y时域图 增加噪声得y2 FFT得Y2
t=(0:length(y)-1)/Fs;
t=t';
noise=0.3*sin(2*pi*1800*t); %1800hz噪声
y2=y+noise;
sound(y2,Fs);
pause(1.5);
figure;plot(t,y);
figure;plot(t,y2);
Y=fft(y);
Y2=fft(y2);
f=(0:length(Y2)-1)/length(Y2) * Fs;
figure;plot(f,abs(Y))
figure;plot(f,abs(Y2))
xlim([0,0.5*Fs])
ylim([0,350])%% 巴特沃斯带阻滤波器
Fs=8000; %采样频率
fp1=1500;fp2=2100;
fs1=1750;fs2=1850;
wp1=fp1/Fs*2*pi;   wp2=fp2/Fs*2*pi;   %通带截止频率
wp=[wp1,wp2];
ws1=fs1/Fs*2*pi;   ws2=fs2/Fs*2*pi;   %阻带截止频率
ws=[ws1,ws2];
Rp=1;As=30;
[n,wc]=buttord(wp/pi,ws/pi,Rp,As)
[b,a]=butter(n,wc,'stop') %求数字带阻滤波器系数
[H,w]=freqz(b,a);
dbH=20*log10(abs(H)/max(abs(H)));
figure;plot(w/2/pi*Fs,dbH,'r');
axis([0,4000,-40,2]);%% 滤波得y3
y3=filter(b,a,y2);
Y3=fft(y3);
figure;plot(f,abs(Y3));
xlim([0,0.5*Fs])
sound(y3,Fs);   

 

matlab滤音频 某一频率的噪声 利用带阻滤波器相关推荐

  1. matlab滤除某一频率的信号,设计一个matlab带通滤波器代码采样频率10Hz,滤除除0.9-1.1Hz之外的其他信号,信号输入为x1,输出为x2,不明白的问我可以追加条件,运行好使的追加分...

    优质解答 % 用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器: %信号为0.5hz, 0.9hz, 1.1hz和1.5hz的正统信号叠加组成 %通带为[0.9,1.1] %频谱分辨率与信号实际长度 ...

  2. 利用matlab小波GUI工具箱对空挡噪声进行滤波

    利用matlab小波GUI工具箱对空挡噪声进行滤波: 1.命令wavemenu调出matlab小波GUI工具箱, 2.File->load加载噪声音频文件,wavelet 1-D,选择小波函数, ...

  3. matlab与音频处理

    MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 1 MATLAB 音频相关函数 1 MATLAB 处理音频信号的流程 2 音量标准化 2 声道分离合并与组合 3 数 ...

  4. MATLAB 与 音频处理 相关内容摘记

    原文地址为: MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 1 MATLAB 音频相关函数 1 MATLAB 处理音频信 ...

  5. 音频耦合电容容值大小如何确定 音频地环路造成的噪声分析

    PART1:音频耦合电容容值大小如何确定 对于硬件工程师来说,或多或少都会接触到音频电路,其中的耦合电容是少不了的了. 提到这个耦合电容,大家都能知道其作用是隔直通交.然而我们发现不同的电路中,这个电 ...

  6. 使用Matlab提取ADC采样数据中的噪声

    正常波形与噪声波形的区别 正常方波信号的特征:一个周期内的数据按时间单调变化,则仅有2种情况:① 数据在上升或下降边沿时,其两边的数据则是一边大一边小(如下图中C和D点):② 数据处于平稳时,则两边大 ...

  7. 利用matlab对音频做FFT进行频谱分析

    利用matlab对音频的左右声道做频谱分析. syms y fs; syms in_put; in_put= 'C:\Users\getao\Desktop\立体声.wav' %in_put 为音频文 ...

  8. matlab 傅里叶变换 音频,傅里叶变换进行音频处理.pdf

    傅里叶变换进行音频处理 2012年第 1期 福 建 电 脑 61 比较在时间缩放和俯仰平移算法中海明窗函数和 布菜克曼窗函数的效果 林志伟 ,达 力 ,王 浩 ,韩 炜 ,林 凡 (1.厦 门大学 计 ...

  9. MATLAB播放音频sound

    MATLAB播放音频sound >> clear >> [y, Fs] = audioread('a_hisa.mp3'); %% y是音频数据m*n的矩阵,Fs是采样率.n是 ...

最新文章

  1. 机器学习知识点(十八)密度聚类DBSCAN算法Java实现
  2. 小明分享|sigmstar SSD201/SSD202 针对MIPI的LCD屏配置操作说明分享,更多Sigmastar相关技术文档及沟通讨论,可以给小明留言,小明邀请你进入我们的技术沟通群
  3. python if语句 for语句 while语句
  4. 信息科学 计算机 区别,电子信息科学技术和计算机科学技术有什么区别啊
  5. IC设计常用文件及格式介绍
  6. “杀死” APP 的留白设计!
  7. Atitit spring单元测试 注解 获取服务名 Spring文件单独放在一个文件夹,去掉dubbo配置,方便测试 里面包含的mybatis 找不到,只好设置成相对于class绝对路径可以了
  8. 计算机二级Java考试心得
  9. 基于Java+SpringBoot+vue+elementui社区疫情防控系统详细设计实现
  10. 银行家算法C++代码实现
  11. 线性代数之——消元法
  12. 运行项目页面显示“该网页无法正常运作”的解决办法 - 前端篇
  13. 电脑系统还原节点怎么创建
  14. AtCoder 001 A BBQ Easy
  15. 用“人话”解释不精确线搜索中的Armijo-Goldstein准则及Wolfe-Powell准则
  16. Win7使用技巧(持续更新...)
  17. 两个小时教你明白C++中类型转换关系
  18. Google Earth Engine(GEE)——Landsat 8TI/TOA/SR影像对比分析区别和去云即NDVI计算
  19. java vector的排序_java Vector排序及用法
  20. 树莓派编译安装FFmpeg(添加H.264硬件编解码器支持)

热门文章

  1. 【谈谈知识点】高代(Ⅰ)行列式矩阵
  2. Python 基础 — Python 简介
  3. 含参数的函数恒成立问题_Simplelife_新浪博客
  4. 使用MapBox自定义地图
  5. BLE 自适应跳频及算法实现
  6. 分步式数据库_创建真实数据科学项目的分步指南
  7. layui模板引擎的使用
  8. 基于asp.net737教学质量评价系统(评教)
  9. 沃尔玛申请了个物联网专利 能做到“补货”自动下单
  10. IPhone 手机屏幕 不显示了。 只有一个背景界面。其他都不显示。怎么回事。