作业帮 > 数学 > 作业

拟合_洒驾问题已知时间:tdata=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6

来源:学生作业帮 编辑:拍题作业网作业帮 分类:数学作业 时间:2024/05/02 05:50:23
拟合_洒驾问题
已知时间:tdata=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16];
酒精含量:vdata=[30 68 75 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4];
求式中的k1,k2,k3
tdata对就公式中的t,vdata对应公式中的v
造一个myfun.m
function F = myfun(k)
tdata=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16];
vdata=[30 68 75 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4];
D0=1;
F=D0.*k(1).*k(2)./(k(2)-k(1)).*(exp(-k(1).*tdata)./(k(3)-k(1))+exp(-k(2).*tdata)./(k(3)-k(2))-exp(-k(3).*tdata).*(1./(k(3)-k(1))+1./(k(3)-k(2))));
然后主程序内运行
k0=[.5 .6 .7]; % 这个你自己改成初始的参数估计
[k,resnorm] = lsqnonlin(@myfun,k0)