写一个matlab洛伦兹函数拟合多个谱峰的代码并绘制曲线

以下是一个matlab洛伦兹函数拟合多个谱峰的示例代码，并绘制曲线：

% 生成示例数据
x = linspace(0,10,1000);
y = 1./(1+x.^2) + 0.05*randn(size(x));
num_peaks = 5;
true_params = zeros(num_peaks,3);
for i=1:num_peaks
    true_params(i,:) = [rand*10 randn rand/10];
    y = y + true_params(i,3) ./ ((x - true_params(i,1)).^2 + true_params(i,2)^2);
end

% 定义洛伦兹函数
lorentz = @(x,p) p(3) ./ ((x - p(1)).^2 + p(2)^2);

% 定义误差函数
err_fun = @(p) sum((y - arrayfun(@(x) sum(lorentz(x,p)), x)).^2);

% 初始猜测值
guesses = zeros(num_peaks,3);
for i=1:num_peaks
    [~,max_idx] = max(y);
    guesses(i,:) = [x(max_idx) 1 1];
    y(max_idx) = 0;
end

% 最小化误差函数
fit_params = fminsearch(err_fun, guesses);

% 绘制拟合曲线
figure;
plot(x,y,'k');
hold on;
for i=1:num_peaks
    plot(x,lorentz(x,fit_params(i,:)),'r');
end

运行以上代码，将会生成一个包含真实数据和拟合曲线的图像，其中每个红色线条代表一个谱峰的拟合曲线。