Mein MATLAB Forum

function[f_g]=faltungssumme(x_k,f_k,g_k,ges_k)

#Die Funktion bekommt später die zwei zu faltenden Zahlenfolgen übergeben
#Hier wird eine Funktion zur Generierung zunächst aufgerufen!
   nullvektor=zeros(1,length(x_k));
   schiebevektor=[g_k, nullvektor, nullvektor];
   x_start=x_k(1);
   x_end=x_k(length(x_k));
   x_diff=x_k(length(x_k))-x_k(1);
   x_schiebevektor=linspace(x_start-x_diff,x_end+x_diff,length(schiebevektor));
#Vorbelegung des Faltungsvektors zum Plotten
   gefaltet=zeros(1, length(schiebevektor));
#Nun kommt ein Schieberegister
   for(i=1:length(schiebevektor))
      schiebevektor=[0,schiebevektor];
      schiebevektor=schiebevektor(1:length(schiebevektor)-1);
      faltungsvektor=schiebevektor(length(x_k)+1:2*length(x_k));
#Hier wird nun Integriert!!!
      gefaltet(i)=sum(f_k.*faltungsvektor)*(x_diff/length(x_k));
      f_g=gefaltet(length(x_k)-((length(x_k)-1)*0.5):2*length(x_k)-((length(x_k)-1)*0.5)-1);
   endfor
   #plot(x_k,f_k,';f(x);',x_k,g_k,';g(x);',x_k,f_g,';f*g(x);',x_k,ges_k,';f*g gemessen(x);');
endfunction