clc;
clear all;
close all;
[X,Fs,bits] = wavread('C:\Users\xiaodan\Desktop\如果这就是爱情.wav');
X1 = X(:,1); %取出双通道中其中一个通道
%sound(X1,44100,16); %改变一下播放的频率或者比特数,会有不同的声音效果
i = length(X1);
t=(0:i-1)/Fs;
n = (1:i)';
v = 0.8*randn(i,1);
ar = [1,1/2];
v1 = filter(1,ar,v);
x = X1 + v1;
%sound(x,Fs,bits);
ma = [1,-0.8,0.4,-0.2];
v2 = filter(ma,1,v);
%sound(v2,Fs,bits);
L = 12;
mu = 0.000026;
dn = x ; % 期望响应
[ylms,W,elms] =filter_LMS(v2,x,L,mu);
%sound(ylms,Fs,bits);
%sound(elms,Fs,bits);
y = X1 - elms;
%sound(y,Fs,bits);
figure;
subplot(4,1,1);
plot(t,X1);grid;
ylabel('Amplitude');
xlabel('Time index(n)');
axis([0 20 -5 5]);
title('Original Signal(X1)');
wavwrite(X1,Fs,bits,'初始音频');
% 绘制滤波器输入信号
subplot(4,1,2);
plot(t,x);grid;
title('Desired input to the Adaptive Filter(X1+v1)');
wavwrite(x,Fs,bits,'含噪音频');
% 绘制自适应滤波器输出信号
%figure;
subplot(4,1,3);
plot(t,elms);grid;
title('Signal output from the Adaptive Filter(elms)');
wavwrite(elms,Fs,bits,'去噪音频');
subplot(4,1,4);
plot(t,y);grid;
title('Error(X1-elms)');
语音去噪LMS自适应滤波算法
3星 · 超过75%的资源 需积分: 32 69 浏览量
2014-05-11
18:49:36
上传
评论 21
收藏 1KB RAR 举报
小蛋0103
- 粉丝: 0
- 资源: 3
最新资源
- SRS-CentOS7-x86-64-5.0-r1.zip
- 全文检索使用ElasticSearch实现全文检索的详细说明和实践探索
- Python《基于TensorFlow的人脸识别(面部表情识别)卷积神经网络 》+源代码+文档说明
- Python《利用bert预训练的中文模型进行文本分类 数据集中文情感分析语料》+源代码+文档说明
- Keil MDK主题美化和代码美化
- 《单片机PID算法的恒温控制系统仿真与程序源码设计(DS18B20传感器) 》+源代码
- 基于matlab小波变换图像压缩系统代码15
- 基于matlab主成分分析图像压缩重建系统代码14
- C语言《基于AT89C52单片机搭建的PID直流电机控制程序,主要用于Proteus电路仿真》+源代码
- C++《基于规划和汽车动力学的自动驾驶项目仿真(基于ROS的自动驾驶项目仿真,使用DWA路径规划算法和双PID控制器)》+源代码
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
- 1
- 2
- 3
前往页