找回密码
 注册
关于网站域名变更的通知
查看: 463|回复: 3
打印 上一主题 下一主题

lsqcurvefit 拟合

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-2-22 11:17 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
对于代码画出来的图如何使用lsqcurvefit 拟合?代码中c为X轴,FF为y轴
) H6 ]8 G6 H7 R' a6 i function FF0dian_TDdyTC
- @% h9 F) p  Q% xclc3 ]+ ^; C, k# E
TC=283:367;
% \) k: x. C- b% 0.025<vaDaRE<14  0.3<eps<0.55  0.29<dtad0<0.7  32<活化能<42  Co=UDa/eps9 j9 K7 l" p7 w7 V' p% ~: r
VADARE=14;) Y' Y) j8 A" v/ B" y/ M, H; ~
EPS=0.55;dtad0=0.29;UDA=10*EPS;
4 l* i4 c/ k- T6 o3 bRH=1;Y=36;
- n( B, {% f! FTR=300;TD=TC; $ m. f% N' F" y( F8 U, X  W
YD=TD/TR;
: n; J0 ]% W4 r0 b' s( C' i/ lfor i=1:85& e( [2 t* u, Y. E* n
start_t=0;end_t=2;. r1 u, I, t0 d* I
options=odeset('RelTol',1e-8,'AbsTol',[1e-9 1e-9]);& E- q9 m* m/ f8 w) i# \0 y
[T,y]=ode45(@SBR1,[start_t end_t],[0,277/300],options,TC(i));. H* d5 ?  m& b9 e$ W2 V2 ^
a=max(y(:,2));
5 W  x! E3 j  k- ^' V0 F  Ab=y(:,2)==a;- V/ X- C/ d" U& }$ ^; l) x
c=TC(i)./TR;
0 k6 E4 E" H0 @( e$ `& |7 A* IT(b);
4 }7 r0 t+ l7 H if T(b)<=1: W' g8 G+ L9 q; A/ J. k1 b
     YTA=c+1.05*dtad0/(UDA*(1+EPS*T(b))+EPS*RH);* b4 p3 }+ A5 a
     FF(i)=a-YTA;, H9 @3 t' W( f7 c
else
# l: y. n, B( Q: g9 ?2 p     YTA=c+1.05*dtad0/(UDA*(1+EPS)+EPS*RH);, n& n: |' N) q, x* o
     FF(i)=a-YTA;
) t& \" A% ?2 ^8 l# k' N; ~/ F3 [ end6 G/ F1 t( M1 f" B  x) j! R& C
end( ]6 M' T8 J5 H7 o  m- A3 h' S$ c
figure;
2 Z* u. ^5 n4 z! p  Pplot(TC/TR,FF,'.');hold on0 J# l8 S6 I7 |" E9 [7 O
end
2 U" Q7 l+ ]' s+ m- d
" V# z4 {) P* B* Y: Sfunction dy= SBR1(t,y,TC)
+ Z+ i/ G- h8 b: l8 K1 x" C& _dy=zeros(2,1);
$ j# c) M6 y. wVADARE=14;
+ o8 q" c5 O& C- z* W( o8 SEPS=0.55;dtad0=0.29;UDA=10*EPS;
1 q6 t, s* b, S) h4 o! PRH=1;Y=36;
/ A7 f! C7 p) u8 l! dTR=300;TD=TC; YD=TD/TR;3 Q* \$ G2 s6 R/ H
if t<=1
! K* L0 c2 j. `, _: E    YCE=((UDA*(1+EPS*t)*TC+RH*EPS*TD)/(UDA*(1+EPS*t)+RH*EPS))/TR;
; U/ V% X* r* p( D    dy(1) = VADARE*(t-y(1))*(1-y(1))*exp(Y*(1-1/y(2)));9 v7 U. t# b1 f- c$ Z1 O+ J
    dy(2) = ((dtad0*dy(1)-(UDA*(1+EPS*t)+RH*EPS)*(y(2)-YCE)))/(1+RH*EPS*t);7 I( g/ g' O$ ?7 p4 b3 _
end1 r* E' M4 K% a8 ^, c# y  R
if t>1
. I4 O8 O' R' j3 _   YCE=((UDA*(1+EPS)*TC+RH*EPS*TD)/(UDA*(1+EPS)+RH*EPS))/TR;
- t/ s; B% M( F2 n/ X8 ?   dy(1) = VADARE*(1-y(1))*(1-y(1))*exp(Y*(1-1/y(2)));1 T$ Y# y- L% `; d
   dy(2) = ((dtad0*dy(1)-(UDA*(1+EPS)+RH*EPS)*(y(2)-YCE)))/(1+RH*EPS);' |5 m- B; ~+ k6 Z- k" T
end  m, `; _6 l# q+ C
end
2 {) t1 v. _2 G& F  U% f, C
& q# V: @, A7 f7 F, R3 u$ k# N

该用户从未签到

2#
发表于 2021-2-22 13:02 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2021-2-22 13:04 | 只看该作者
请把待拟合的数据、哪些是待拟合参数等等都给出

该用户从未签到

4#
发表于 2021-2-22 13:15 | 只看该作者
所谓拟合,指的是已知形如 Output_Y = Function( Input_X, Parameters_C ) 这样的函数关系(其中Input_X为实验的若干控制变量,同时也是函数的输入;Parameters_C为待拟合参数;Output_Y为实验结果,同时也是函数的输出),给出Input_X与Output_Y 作为已知,求一组合适的Parameters_C数值使模型Function产生的误差最小。
3 M9 J9 L. N2 `1 P+ q+ ~* S, l! O& I; A$ X2 e$ x
最简单的例子,已知模型:电压U = 电阻R * 电流I,且不知道接在某简单电路上的电阻到底是多大。那么通过若干组实验,以电流I为输入(比如调整滑动变阻器来控制该电路中的电流数值),同时测量这个电阻两端的电压U作为输出,通过这若干组U与I的数值,可以估计出一个满足模型的电阻R值。控制输入I、测量输出U、以此二者来估计参数R。
1 Y. C  L' u' y3 n* E" A
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-6-16 03:06 , Processed in 0.093750 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表