4 Star 11 Fork 2

叶夜笙歌 / FastICA

加入 Gitee
与超过 1200万 开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)
免费加入
该仓库未声明开源许可证文件(LICENSE),使用请关注具体项目描述及其代码上游依赖。
克隆/下载
test.m 2.95 KB
一键复制 编辑 原始数据 按行查看 历史
叶夜笙歌 提交于 2014-03-30 17:35 . FastICA盲源分离算法DEMO
close all;
clear all;
clc;
% x1 = wavread('C:\Users\sound\Desktop\sound\wav_data11\channel1.wav');
% x2 = wavread('C:\Users\sound\Desktop\sound\wav_data1\channel2.wav');
% x3 = wavread('C:\Users\sound\Desktop\sound\wav_data1\channel3.wav');
% x4 = wavread('C:\Users\sound\Desktop\sound\wav_data11\channel4.wav');
% x5 = wavread('C:\Users\sound\Desktop\sound\wav_data1\channel5.wav');
% x6 = wavread('C:\Users\sound\Desktop\sound\wav_data1\channel6.wav');
% x7 = wavread('C:\Users\sound\Desktop\sound\wav_data1\channel7.wav');
% x8 = wavread('C:\Users\sound\Desktop\sound\wav_data11\channel8.wav');
% x4 = wavread('C:\Users\sound\Desktop\wav_data1\channel4.wav');
x1 = load('C:\Users\sound\Desktop\sigArray1.txt');
x2 = load('C:\Users\sound\Desktop\sigArray2.txt');
x3 = load('C:\Users\sound\Desktop\sigArray3.txt');
x4 = load('C:\Users\sound\Desktop\sigArray4.txt');
x5 = load('C:\Users\sound\Desktop\sigArray5.txt');
x6 = load('C:\Users\sound\Desktop\sigArray6.txt');
% N=length(x1);
% x(1,1:N-8)=x1(9:N);
% x(2,1:N-9)=x2(10:N);
% x(3,1:N-10)=x3(11:N);
% x(1,:)=x1;
% x(2,:)=x2;
% x(3,:)=x3;
% x1 = wavread('001000010mix1.wav');
% x2 = wavread('001000010mix2.wav');
% x1 = wavread('C:\Users\sound\Desktop\\111111111mix1.wav');
% x2 = wavread('C:\Users\sound\Desktop\\111111111mix2.wav');
% x3 = wavread('C:\Users\sound\Desktop\\111111111mix3.wav');
% x4 = wavread('C:\Users\sound\Desktop\\111111111mix4.wav');
% x5 = wavread('C:\Users\sound\Desktop\\111111111mix5.wav');
% x6 = wavread('C:\Users\sound\Desktop\\111111111mix6.wav');
% x7 = wavread('C:\Users\sound\Desktop\\111111111mix7.wav');
% x8 = wavread('C:\Users\sound\Desktop\\111111111mix8.wav');
% x9 = wavread('C:\Users\sound\Desktop\\111111111mix9.wav');
% x3 = wavread('source3.wav');
% x4 = wavread('source8.wav');
% xx = x3 +x4;
x(1,:)=x1;
x(2,:)=x2;
% x(3,:)=x3(40000:50000);
% x(4,:)=x4(40000:50000);
% x(5,:)=x5(40000:50000);
% x(6,:)=x6(40000:50000);
% x(7,:)=x7(40000:50000);
x(3,:)=x3;
x(4,:)=x4;
x(5,:)=x5;
x(6,:)=x6;
% x(9,:)=x9;
% x(5,:)=x5;
% x(6,:)=x6;
% x(7,:)=x7;
% x(8,:)=x8;
icasig = fastica(x);
% icasig1 = fastica(icasig);
figure;
subplot(5,2,1);
plot(x1);
subplot(5,2,2);
plot(icasig(1,:));
subplot(5,2,3);
plot(icasig(2,:));
subplot(5,2,4);
plot(icasig(3,:));
% subplot(5,2,5);
% plot(icasig(4,:));
% subplot(5,2,6);
% plot(icasig(5,:));
% subplot(5,2,7);
% plot(icasig(6,:));
% subplot(5,2,8);
% plot(icasig(7,:));
% subplot(5,2,9);
% plot(icasig(8,:));
% subplot(5,2,10);
% plot(icasig(9,:));
% figure;
% subplot(3,1,1);
% plot(x3);
% subplot(3,1,2);
% plot(x4);
% subplot(3,1,3);
% plot(xx);
% sound(x1,8000);
sound(icasig(1,:),8000);
sound(icasig(2,:),8000);
sound(icasig(3,:),8000);
sound(icasig(4,:),8000);
sound(icasig(5,:),8000);
sound(icasig(6,:),8000);
% sound(icasig(3,:),8000);
% sound(icasig(4,:),8000);
% sound(icasig(5,:),8000);
% sound(icasig(6,:),8000);
% sound(icasig(7,:),8000);
% sound(icasig(8,:),8000);
% sound(icasig(9,:),8000);
Matlab
1
https://gitee.com/a_bad_geek/FastICA.git
git@gitee.com:a_bad_geek/FastICA.git
a_bad_geek
FastICA
FastICA
master

搜索帮助