紧急求助，利用Matlab对实验数据进行拟合求解参数。 - 计算模拟

做好每一件小事

21楼2013-04-16 19:37:17

tangwh

铁虫 (著名写手)

应助: 24 (小学生)
金币: 2830.8
红花: 3
帖子: 1055
在线: 106.4小时
虫号: 173851
注册: 2006-01-21
性别: GG
专业: 可再生与替代能源利用中的

没有Matlab编程基础，能不能帮我给编程贴出来呢

题外话

22楼2013-04-17 15:43:30

月浅浅

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2593.5
散金: 160
红花: 6
帖子: 1048
在线: 281.2小时
虫号: 934684
注册: 2009-12-29
性别: GG
专业: 光学

引用回帖:

9楼: Originally posted by 月只蓝 at 2013-04-15 08:55:43
function feixianxingnihe1
clear all;clc
format long
data1=;
tspan=data1(:,1); %r的数据，在此输入
xexp=data1(:,2); %E的数据，在此输入

k0=;
lb=-*1e5;
ub=*1e5;

%-------------------- ...

我参考您的程序，改了实验数据，但是一运行就报错：

Error using lsqnonlin (line 178)
LSQNONLIN requires the following inputs to be of data type double: 'LB'.

Error in myfeixianxingnihe (line 199)
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

不知道这个问题出在哪里啊？

23楼2015-07-15 16:22:55

月浅浅

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2593.5
散金: 160
红花: 6
帖子: 1048
在线: 281.2小时
虫号: 934684
注册: 2009-12-29
性别: GG
专业: 光学

引用回帖:

24楼2015-07-15 16:23:16

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

24楼: Originally posted by 月浅浅 at 2015-07-15 16:23:16
我参考您的程序，改了实验数据，但是一运行就报错：

Error using lsqnonlin (line 178)
LSQNONLIN requires the following inputs to be of data type double: 'LB'.

Error in myfeixianxingnihe (line 19 ...

请给出完整程序。

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

25楼2015-07-15 16:26:13

月浅浅

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2593.5
散金: 160
红花: 6
帖子: 1048
在线: 281.2小时
虫号: 934684
注册: 2009-12-29
性别: GG
专业: 光学

引用回帖:

25楼: Originally posted by 月只蓝 at 2015-07-15 16:26:13
请给出完整程序。...

这里回复不能贴图，我发表个主题，今天看了您发的两个拟合相关的帖子，里面可能有些东西没看太懂，非常希望您能给些意见

26楼2015-07-16 23:26:14

月浅浅

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2593.5
散金: 160
红花: 6
帖子: 1048
在线: 281.2小时
虫号: 934684
注册: 2009-12-29
性别: GG
专业: 光学

引用回帖:

25楼: Originally posted by 月只蓝 at 2015-07-15 16:26:13
请给出完整程序。...

请版主移步，给些意见

http://muchong.com/bbs/viewthread.php?tid=9175607&pid=1#pid1

27楼2015-07-16 23:48:52

月浅浅

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2593.5
散金: 160
红花: 6
帖子: 1048
在线: 281.2小时
虫号: 934684
注册: 2009-12-29
性别: GG
专业: 光学

引用回帖:

%-------------------------------------------------------------------------

function f = ObjFunc(k,tspan,xexp)
f=KineticsEqs(tspan,k)-xexp;

%------------------------------------------------------------------------
function xt = KineticsEqs(t,k)
xt=-(k(1)+k(2))+k(1)*exp(-t/k(3))+k(2)*exp(-t/k(4));

请问高手在这里function f = ObjFunc(k,tspan,xexp) 这个是起什么作用的？因为下面接着又有一个f=KineticsEqs(tspan,k)-xexp;

28楼2015-07-19 22:05:36