Approximate a System Frequency Response with a Filter in MATLAB
Convolve the input signal with (Real Part) harmonic signal $ \alpha {e}^{-2 \pi i {f}_{0} t + \phi} $.
Let's call this output as $ {y}_{1} $.
The parameters $ \alpha $ and $ \phi $ are from your reference system response.

This will give you the output of the reference system on your input signal only at your desired frequency.

Now convolve the input signal with (Real Part) harmonic signal $ 1 {e}^{-2 \pi i {f}_{0} t} $.
Let's call this output as $ {y}_{2} $.
This will extract this component of your input signal.

Calculate $ {y}_{3} = {y}_{1} - {y}_{2} $ and all you need is to add $ {y}_{3}  $ to your signal.

There are trickier and more efficient way to do so.
Yet I wanted to keep it simple and intuitive.