matlab滤音频 某一频率的噪声 利用带阻滤波器
%% 对一段音频滤去噪声%% 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滤音频 某一频率的噪声 利用带阻滤波器相关推荐
- matlab滤除某一频率的信号,设计一个matlab带通滤波器代码采样频率10Hz,滤除除0.9-1.1Hz之外的其他信号,信号输入为x1,输出为x2,不明白的问我可以追加条件,运行好使的追加分...
优质解答 % 用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器: %信号为0.5hz, 0.9hz, 1.1hz和1.5hz的正统信号叠加组成 %通带为[0.9,1.1] %频谱分辨率与信号实际长度 ...
- 利用matlab小波GUI工具箱对空挡噪声进行滤波
利用matlab小波GUI工具箱对空挡噪声进行滤波: 1.命令wavemenu调出matlab小波GUI工具箱, 2.File->load加载噪声音频文件,wavelet 1-D,选择小波函数, ...
- matlab与音频处理
MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 1 MATLAB 音频相关函数 1 MATLAB 处理音频信号的流程 2 音量标准化 2 声道分离合并与组合 3 数 ...
- MATLAB 与 音频处理 相关内容摘记
原文地址为: MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 MATLAB 与 音频处理 相关内容摘记 1 MATLAB 音频相关函数 1 MATLAB 处理音频信 ...
- 音频耦合电容容值大小如何确定 音频地环路造成的噪声分析
PART1:音频耦合电容容值大小如何确定 对于硬件工程师来说,或多或少都会接触到音频电路,其中的耦合电容是少不了的了. 提到这个耦合电容,大家都能知道其作用是隔直通交.然而我们发现不同的电路中,这个电 ...
- 使用Matlab提取ADC采样数据中的噪声
正常波形与噪声波形的区别 正常方波信号的特征:一个周期内的数据按时间单调变化,则仅有2种情况:① 数据在上升或下降边沿时,其两边的数据则是一边大一边小(如下图中C和D点):② 数据处于平稳时,则两边大 ...
- 利用matlab对音频做FFT进行频谱分析
利用matlab对音频的左右声道做频谱分析. syms y fs; syms in_put; in_put= 'C:\Users\getao\Desktop\立体声.wav' %in_put 为音频文 ...
- matlab 傅里叶变换 音频,傅里叶变换进行音频处理.pdf
傅里叶变换进行音频处理 2012年第 1期 福 建 电 脑 61 比较在时间缩放和俯仰平移算法中海明窗函数和 布菜克曼窗函数的效果 林志伟 ,达 力 ,王 浩 ,韩 炜 ,林 凡 (1.厦 门大学 计 ...
- MATLAB播放音频sound
MATLAB播放音频sound >> clear >> [y, Fs] = audioread('a_hisa.mp3'); %% y是音频数据m*n的矩阵,Fs是采样率.n是 ...
最新文章
- 机器学习知识点(十八)密度聚类DBSCAN算法Java实现
- 小明分享|sigmstar SSD201/SSD202 针对MIPI的LCD屏配置操作说明分享,更多Sigmastar相关技术文档及沟通讨论,可以给小明留言,小明邀请你进入我们的技术沟通群
- python if语句 for语句 while语句
- 信息科学 计算机 区别,电子信息科学技术和计算机科学技术有什么区别啊
- IC设计常用文件及格式介绍
- “杀死” APP 的留白设计!
- Atitit spring单元测试 注解 获取服务名 Spring文件单独放在一个文件夹,去掉dubbo配置,方便测试 里面包含的mybatis 找不到,只好设置成相对于class绝对路径可以了
- 计算机二级Java考试心得
- 基于Java+SpringBoot+vue+elementui社区疫情防控系统详细设计实现
- 银行家算法C++代码实现
- 线性代数之——消元法
- 运行项目页面显示“该网页无法正常运作”的解决办法 - 前端篇
- 电脑系统还原节点怎么创建
- AtCoder 001 A BBQ Easy
- 用“人话”解释不精确线搜索中的Armijo-Goldstein准则及Wolfe-Powell准则
- Win7使用技巧(持续更新...)
- 两个小时教你明白C++中类型转换关系
- Google Earth Engine(GEE)——Landsat 8TI/TOA/SR影像对比分析区别和去云即NDVI计算
- java vector的排序_java Vector排序及用法
- 树莓派编译安装FFmpeg(添加H.264硬件编解码器支持)