I understand that you want to derive an equivalent circuit model for the given impedance spectrum. From my understanding, you can try out the following workarounds to fit the data to an equivalent model.
The reference attached above discusses usage of two methods namely lsqcurvefit and fminunc to fit a nonlinear function to the provided data. The documentation for lsqcurvefit is attached below for your reference
In case you struggle to narrow down to an approach to fit the impedance data, it would be highly convenient if you could share your code so that I can go through it and address your issue in a more efficient manner.