应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求助用matlab 多元非线性拟合2
222/3返回列表上一页123下一页
查看: 3079  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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倍及其以上，才能拟合出来，否则不行。
拟合的原理是用最小二乘法解多元方程。
赞一下(1人) 回复此楼
请大家去我的版面看看，正在收集开源期刊，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 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭