clear;clc;close all
%3.5.1 载入语音信号
N=10^4;
% N=8000;
load chirp; s1=y(1:N); %应为列向量
load gong; s2=y(1:N);
%chirp与gong是matlab自带的语音文件，前者是敲钟声，后者是嬉笑声。
%3.5.2 对输入信号中心化
s1=s1/std(s1);s2=s2/std(s2); %std()函数是用于标准差
%===================================
% k=1:100;
% s1=sin(2*pi*k/32);
% s2=square(2*pi*k/8);
%==========================
%3.5.3 组合信号、混合信号、定义解混矩阵
s=[s1,s2];
figure(1)
subplot(211)
plot(s1)
subplot(212)
plot(s2)
M=2;
maxiter=100;
eta=10^(-3);
A=randn(M,M);
x = s*A;
figure(2)
subplot(211)
plot(x(:,1))
subplot(212)
plot(x(:,2))
