| 查看: 882 | 回复: 3 | ||
云凌霜1993金虫 (小有名气)
|
[求助]
matlab 曲线拟合的问题，拜托了
|
|
求助大家，通过实验数据得出了一系列C和F，如何通过matable拟合出指定方程C= 1-exp(-k/F)里头的K值呢？从别人那里要来的代码计算总出现错误，真心看不懂。 Undefined function or method 'eval' for input arguments of type 'double'. Error in ==> main>zuixiao at 40 y(i)=eval(C1)*100;% (corresponding output) Error in ==> fmincon at 540 initVals.f = feval(funfcn{3},X,varargin{:}); Error in ==> main at 11 k1 = fmincon(@zuixiao,[0],[],[],[],[],0,[],[]); Caused by: Failure in initial user-supplied objective function evaluation. FMINCON cannot continue. ??? Undefined function or method 'eval' for input arguments of type 'double'. 代码贴在这里了，拜托大家了！！！ function main() clear all close all warning off syms k F C data=[0.0153 99.91 0.0244 99.72 0.0312 99.31 0.0463 97.82];% (figure) C= 1-exp(-k/F) ;% (defined formula) k1 = fmincon(@zuixiao,[0],[],[],[],[],0,[],[]); C=subs(C,k,k1) ;% current value disp(['The corresponding expressions for：' ]) disp(vpa(C,4)) x=0.01:0.001:0.05; for i=1:length(x) Temp=subs(C,F,x(i)) ; y(i)=eval(Temp)*100; end figure('color',[1 1 1]) scatter(data(:,1),data(:,2),'.') hold on plot(x,y) xlabel('Flux(mL/cm^2s)') ylabel('Conversion %') legend('Experimental','Fitted curve') box on function error=zuixiao(x) syms k F C data=[0.0153 99.91 0.0244 99.72 0.0312 99.31 0.0463 97.82];% (figure) C=1-exp(-k/F);% (defined formula) C=subs(C,k,x);% (current value) m=length(data(:,1)); error=0;% (initial error) for i=1:m C1=subs(C,F,data(i,1)); y(i)=eval(C1)*100;% (corresponding output) error =error+(y(i)-data(i,2))^2;% (square error) end |
回复此楼» 猜你喜欢
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有18人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有163人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有17人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于matlab微分、及非线性拟合问题
已经有5人回复
|
为什么要拟合呢？看公式k=-Fln(1-C)。直接代数据进去算完后再统计平均不可以吗？matlab也很好实现 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
云凌霜19932楼2017-06-09 10:46:45
云凌霜1993
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1620.2
- 红花: 1
- 帖子: 86
- 在线: 95.7小时
- 虫号: 2341170
- 注册: 2013-03-12
- 性别: MM
- 专业: 电化学
3楼2017-06-09 15:22:09
zhangjiangbo
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3293.3
- 散金: 50
- 红花: 5
- 帖子: 740
- 在线: 70.9小时
- 虫号: 309683
- 注册: 2006-12-23
- 专业: 电气科学与工程
4楼2017-06-10 22:32:01
10