fastica去噪matlab_帮帮忙，这个地震去噪的fastICA程序我怎么改不对啊

clear

clc

f=30;

fs=1000;

number=100;

t=-number/2+1:number/2;

x=1:100;

a=(1-2*(pi*f*t/fs).^2).*exp(-(pi*f*t/fs).^2);

subplot(221);plot(x,a);

title('雷克子波');

y1=zeros(1,61);m1=(rand(1)-0.5).*2 ;                                                                                                                         ;

y2=zeros(1,90);m2=(rand(1)-0.5).*2;

y3=zeros(1,13);m3=(rand(1)-0.5).*2;

y4=zeros(1,24);m4=(rand(1)-0.5).*2;

y5=zeros(1,54);m5=(rand(1)-0.5).*2;

y6=zeros(1,52);m6=(rand(1)-0.5).*2;

y=[y1,m1,y2,m2,y3,m3,y4,m4,y5,m5,y6,m6];

subplot(222);

plot(y);

title('反射系数序列y(t)');

s1=conv(y',a);

s2=s1((ceil(length(a)/2):size(s1)-floor(length(a)/2)));

s=s2&