Convolution of functions using data?

Hi everyone.

I would like to find the convolution between two datasets and graph the result. I tried with several tools on matlab but I didn't obtain what I'm searching.

I have two spectra that overlap and I would like to find the convolution between them. The spectra are in vector form, defined by two vectors x and y for each spectrum.

I already used polyfit and interp1 to get the function for my dataset but they can't fit them optimally.

for the first spectrum: x='l544dim' and y='a544dim_n'

for the second spectrum: x='l_laser' and y='a_laser_norm'

Thanks in advance!