| 查看: 3278 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助用matlab 多元非线性拟合2 已有4人参与
|
|||
|
2个问题 第一个问题， 2个自变量X, Y, 一个因变量Z 数据如下，希望用matlab 拟合一个X, Y, Z三维的图，并给出由数据 ，拟合出的方程 最好把编写的程序也发上来，我以后可以 照瓢画葫芦 模仿着做。 X Y Z 200 50 1.044 500 250 1.12 1000 500 1.25 第二个问题， 3个自变量 a ， b， c 一个因变量 f（a,b,c） 不知道， 这种情况不知道是 几维的图了？？ 希望能 给出由数据 拟合出的方程，请给出 编程的具体步骤 f(a,b,c) a b c 1.04 200 200 200 1.23 500 500 500 1.6 1000 1000 1000 |
回复此楼» 猜你喜欢
华师大读博
已经有5人回复
-
国自然面上和省基金B类撒花
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有5人回复
-
0854求调剂
已经有28人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
一志愿A区211，22408 321求调剂
已经有8人回复
-
295分求调剂
已经有6人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有13人回复
-
085600材料与化工调剂
已经有6人回复
|
超级菜鸟求助，能不能帮我把这道题转化为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
2楼2013-12-14 15:58:30
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68132.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2013-12-14 17:16:50
4楼2013-12-15 08:45:07
20