% #########################################################
% 
% Faltung von Verteilungsfunktionen - Versuch
% 
% #########################################################

clear all; close all;
x = -16:.1:16;                   % Wertebereich der x-Achse
pMerkmal = 'Vorderkantendicke d_{VK} [%]'

% #########################################################
% Funktion Fertigungsstreuung
% #########################################################
mittelNew=1.2;
stabwNew=1.7;
typeNew='Normal'
inputAnew = mittelNew
inputBnew = stabwNew
inputCnew = 0

fktFertigungsstreuung = pdf(typeNew,x,inputAnew,inputBnew,inputCnew)


% #########################################################
% Funktion Messwerte "Used Parts"
% #########################################################
mittelUsed=0;
stabwUsed=1;
typeUsed='Normal'
inputAused = mittelUsed
inputBused = stabwUsed
inputCused = 0

fktUsed = pdf(typeUsed,x,inputAused,inputBused,inputCused)


% #########################################################
% Funktion : Verschleiß
% #########################################################
%[fktVerschleiss,res] = deconv(fktFertigungsstreuung,fktUsed);
fktVerschleiss = conv(fktFertigungsstreuung,fktUsed,'same')


% #########################################################
% Diagrammerstellung
% #########################################################

plot(x,fktFertigungsstreuung,'-k',x,fktUsed,'--k',x,fktVerschleiss,'-.k')
hleg1 = legend('Verteilung Neu','Verteilung Gebraucht','Verteilung Verschleiß');
xlabel(pMerkmal,'fontsize',12,'fontweight','b');
ylabel('Häufigkeit [%]','fontsize',12,'fontweight','b');
hold all