24小时热门版块排行榜    

查看: 3079  |  回复: 21

wyliu_as

铁杆木虫 (著名写手)

引用回帖:
8楼: Originally posted by maguohua at 2013-12-15 10:07:34
那请仁兄帮我 把二组数据都拟合一下,并给出方程,有金币奉上...

稍等,我正在下载origin,我家里电脑没有该软件,搞定后马上给你
请大家去我的版面看看,正在收集开源期刊,http://emuch.net/bbs/myf.php?ffid=282
11楼2013-12-15 10:25:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maguohua

新虫 (小有名气)

引用回帖:
11楼: Originally posted by wyliu_as at 2013-12-15 10:25:50
稍等,我正在下载origin,我家里电脑没有该软件,搞定后马上给你...

OK 等你的好消息
12楼2013-12-15 13:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超人不怕

银虫 (小有名气)

引用回帖:
3楼: Originally posted by 月只蓝 at 2013-12-14 17:16:50
问题一

function threeD_
clear all;clc
data=;
x=data(:,1);y=data(:,2);z=data(:,3);
=griddata(x,y,z,linspace(min(x),max(x),50)',linspace(min(y),max(y),50),'v4');%插值
figure,surf(X,Y,Z);shad ...

超级菜鸟求助,能不能帮我把这道题转化为lsqcurvefit来解?急!
function KineticsEst1_Diff  
clear all
clc

% 动力学数据
t = [0  20  40  60  120  180  300];
CAm = [10  8  6  5  3  2  1];

% 用最小二乘样条拟合法计算微分dCA/dt--使用不经过实验点的B样条插值函数
knots = 3;
K = 3;                  % 三次B样条
sp = spap2(knots,K,t,CAm);
pp = fnder(sp);         % 计算B样条函数的导函数
dCAdt = fnval(pp,t)    % 计算t处的导函数值
rAm = dCAdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
CAi = fnval(sp,ti);
plot(t,CAm,'ro',ti,CAi,'b-')
xlabel('t')
ylabel('C_A')
legend('实验值','B样条拟合')

% 非线性拟合
beta0 = [0.0053 1.39];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
     lsqnonlin(@OptObjFunc,beta0,[],[],[],rAm,CAm);         
ci = nlparci(beta,residual,jacobian);

% 参数辨识结果
fprintf('Estimated Parameters:\n')
fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))
fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2))
fprintf('  The sum of the squares is: %.1e\n\n',sum(residual.^2))

% 绘制反应速率拟合曲线
figure
plot(t,rAm,'ro',t,Rate(CAm,beta),'b*')
xlabel('t')
ylabel('dC_Adt')
legend('Experiment','Kinetic Model')

% ------------------------------------------------------------------
function f = OptObjFunc(beta,rAm,CAm)
rAc = Rate(CAm,beta);
f = rAc - rAm;

% ------------------------------------------------------------------
function rA = Rate(CA,beta)
rA = -beta(1)*CA.^beta(2);   % -rA = -dCA/dt = k*CA^n, 其中k=beta(1), n=beta(2)
大菜鸟
13楼2013-12-15 16:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxiang29

银虫 (正式写手)

引用回帖:
6楼: Originally posted by wyliu_as at 2013-12-15 09:03:04
我不会用mat,不过origin有很强大的拟合功能,可以最多做9个自变量的拟合,可以给出拟合方程,还可以自己制定方程进行拟合。不知道楼主会不会用origin,很简单的,希望有用。
三元的可以画出你所需的三维立体图像, ...

这么厉害啊,我现在有x1,x2,x3三个自变量,一个因变量Y,有15组数据,但一直不会拟合方程以及分析,请问您是怎么用origin搞出来的呀???
smile to the world
14楼2013-12-16 22:13:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maguohua

新虫 (小有名气)

引用回帖:
14楼: Originally posted by chenxiang29 at 2013-12-16 22:13:10
这么厉害啊,我现在有x1,x2,x3三个自变量,一个因变量Y,有15组数据,但一直不会拟合方程以及分析,请问您是怎么用origin搞出来的呀???...

好像3个变量的 origin 就不行了
15楼2013-12-19 16:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maguohua

新虫 (小有名气)

引用回帖:
7楼: Originally posted by 月只蓝 at 2013-12-15 09:11:29
f(a,b,c)=p1+p2*ln(a/3+b/3+c/3)+p3/(a/3+b/3+c/3));
式中,
p1=              -532.920658380386
p2=              2376.8235508824
p3=              665.287941422274
拟合结果评价系数:
Correlation Coef. (R): 1
R-Square: ...

还得请教您  
f(a,b,c)=p1+p2*ln(a/3+b/3+c/3)+p3/(a/3+b/3+c/3));
式中,
p1=              -532.920658380386
p2=              2376.8235508824
p3=              665.287941422274

我拥有了a,b,c的值,用这个方程进行计算
用matlab 可以直接算吗?
还是要用excel  算!?  
用matlab 算是不是 还要编程??
16楼2013-12-19 16:28:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
16楼: Originally posted by maguohua at 2013-12-19 16:28:04
还得请教您  
f(a,b,c)=p1+p2*ln(a/3+b/3+c/3)+p3/(a/3+b/3+c/3));
式中,
p1=              -532.920658380386
p2=              2376.8235508824
p3=              665.287941422274

我拥有了a,b,c的值, ...

方程和参数都有 用MATLAB和excel都能算。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
17楼2013-12-19 16:52:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

扯吧,就这么几个数据点,还不控制变量,居然想得到拟合方程?
18楼2013-12-19 20:28:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyliu_as

铁杆木虫 (著名写手)

引用回帖:
11楼: Originally posted by wyliu_as at 2013-12-15 10:25:50
稍等,我正在下载origin,我家里电脑没有该软件,搞定后马上给你...

我经过尝试,发现用于拟合的数据个数必须是自变量的2倍及其以上,才能拟合出来,否则不行。
拟合的原理是用最小二乘法解多元方程。
请大家去我的版面看看,正在收集开源期刊,http://emuch.net/bbs/myf.php?ffid=282
19楼2013-12-20 08:14:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyliu_as

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
14楼: Originally posted by chenxiang29 at 2013-12-16 22:13:10
这么厉害啊,我现在有x1,x2,x3三个自变量,一个因变量Y,有15组数据,但一直不会拟合方程以及分析,请问您是怎么用origin搞出来的呀???...

选用regression功能,将因变量放在第一列,自变量放在后面,全选中自变量,点多重回归,就运算出来出来了
请大家去我的版面看看,正在收集开源期刊,http://emuch.net/bbs/myf.php?ffid=282
20楼2013-12-20 20:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maguohua 的主题更新
信息提示
请填处理意见