Up/Downsampling LMS

[Peter E. (eumelpeter)]

Hi, folgendes Blockschaltbild.. Dieses soll verdeutlichen wo Up and 
Downsampling stattfinden soll
Beim Downsampling werden ja zuerst alle Werte durch einen Tiefpass mit 
entsprechender Frequenz geschickt und anschließend nur jeder zweite Wert 
genommen.
Beim Upsampling soll an jeder zweiten Stelle eine 0 eingefügt werden und 
anschließend werden diese ja wieder durch einen Tiefpass gefiltert.
Die Theorie ist sehr klar und Up und Downsampling habe ich auch schon 
per Hand programmiert.Ist ja auch nicht schwer. Aber beim LMS hab ich 
Probleme.

Alles wird in Matlab programmiert.

Der LMS sieht folgendermaßen aus, ist auch korrekt implementiert.( Mit 
Anmerkungen wo das Up and Downsampling stattfinden soll, vielleicht 
liegt da ja auch bei mir ein Denkfehler das blockschaltbild liegt ja 
vor)

for i=1:length(x_ref);

   [Dieser Teil muss Downgesampelt werden]
    eingang=[x_ref(i,1);eingang(1:parameter-1)];
    y(i)=w_voll*eingang;

[Dieser Teil upgesampelt]
    e(i)=y(i)+d(i);

[Hier wieder ein Downsampling]
    w_voll=w_voll-(mu*e(i)*eingang)';
    
end

Nun zu meinem Code der zwar funktioniert, aber sehr sehr lange fürs 
Rechnen braucht . Über jegliche Hilfe wäre ich sehr dankbar! 
Wahrscheinlich ist die Implementierung auch falsch.

%Initialisierungen
w_para=zeros(1,parameter);
eingang=zeros(parameter,1);
y=zeros(1,length(x_ref));
e=zeros(1,length(x_ref));
 
v_parallel=filter(wtest,1,x_ref);

filt = fdesign.lowpass('N,Fc',12,10000,40000);
Hd = design(filt);

% Tiefpass gefiltertes Signal
x_ref_do=filter(Hd,x_ref);

mu=0.1;

for i=1:length(x_ref_do);

%Hier nur jeden zweiten Wert berechnen
if(mod(i,2)==0)
  eingang=[x_ref_do(i,1);eingang(1:parameter-1)];
    y(i)=w_para*eingang;

% Die Addition bei e muss in jedem Schritt gemacht werden oder? daher ist
% dieser Teil in jeder Schleife zu finden
    y_up(i)=filter(Hd,y(i));
    e(i)=y_up(i)+d(i)+v_parallel(i);
    e_do(i)=filter(Hd,e(i));

else

%Hier Upsampling, also mit 0 aufüllen jeden zweiten Wert
    y(i)=0; 

%Wert filtern für Upsampling
     y_up(i)=filter(Hd,y(i));

%addition der werte
    e(i)=y_up(i)+d(i)+v_parallel(i);

%e muss ja wieder downgesampelt werden
    e_do(i)=filter(Hd,e(i));

end

Danke für eure Hilfe und Bemühungen