发布网友 发布时间:2022-05-01 14:06
共1个回答
热心网友 时间:2023-05-08 11:16
close all;%对数正态分布密度函数
) p(1)./x.*exp(-((log(x)-p(2))/p(3)).^2/2);
[maxy ind]=max(y);
%做非线性数据拟合
p=nlinfit(x,y,fun,[maxy*x(ind),log(x(ind)),1]);
%画拟合曲线
x1 = min(x):0.01:max(x);
yfit=fun(p,x1);
plot(x1,yfit,'k','linewidth',3);
[ymax i] = max(yfit); %找到极值点
xm = x1(i);
plot([xm xm],[0 ymax+10],'k','linewidth',3); %画竖线
xlim([min(x) max(x)]);
xlabel('BC浓度(ng/m^3)');
ylabel('频数');