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

lsqcurvefit 拟合

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
对于代码画出来的图如何使用lsqcurvefit 拟合?代码中c为X轴,FF为y轴% w5 J" _% J) s4 g+ O' V
function FF0dian_TDdyTC
  ?. D" S4 d! {+ z* G; j5 O! V8 bclc
7 L- U5 F+ Z% ]% |7 y, M! WTC=283:367;; B/ g  X- m- y$ u# Y( n. |; s
% 0.025<vaDaRE<14  0.3<eps<0.55  0.29<dtad0<0.7  32<活化能<42  Co=UDa/eps( |* q$ \2 ]% _: `2 C. T
VADARE=14;
% k1 B/ ]" G5 `EPS=0.55;dtad0=0.29;UDA=10*EPS;2 d, ?. Q; p: {6 V/ E; u& x
RH=1;Y=36;. K& b7 m" I( [4 A
TR=300;TD=TC;
- K4 ]: ~. G4 H" vYD=TD/TR;" L+ T+ J, L6 M$ j: F7 z! {
for i=1:85; w" h& c4 N# h
start_t=0;end_t=2;6 N7 H! F' Y7 \: h
options=odeset('RelTol',1e-8,'AbsTol',[1e-9 1e-9]);2 l2 w, j( E( w
[T,y]=ode45(@SBR1,[start_t end_t],[0,277/300],options,TC(i));+ G/ s* z& O) f2 F& h
a=max(y(:,2));
/ T: D6 N; i* zb=y(:,2)==a;
& C0 |" l9 R+ D7 n2 Jc=TC(i)./TR;
  h1 T# R' T. ]9 R% f- h) e+ I% IT(b);; d8 K' Y5 F$ Y1 }6 v
if T(b)<=1' N! `* }) [/ {: |% W5 t! K
     YTA=c+1.05*dtad0/(UDA*(1+EPS*T(b))+EPS*RH);- Q, y) t3 E. P# o  ^+ u
     FF(i)=a-YTA;3 b$ f* R1 k$ R
else0 o4 ?& P- q$ w- o. [5 X- Y
     YTA=c+1.05*dtad0/(UDA*(1+EPS)+EPS*RH);
& x8 A0 z" K1 d# p& h; d5 u8 b     FF(i)=a-YTA;
$ V9 b% s. V0 A" [2 l end
2 H6 s! f6 H# _+ Zend
0 a! }2 f3 L1 P( p2 t% I5 w/ c( hfigure;5 k' M6 h" a) r, k% @( f) Y7 n
plot(TC/TR,FF,'.');hold on! G4 }' D) M0 T, r
end+ e7 ^* {' r7 I  ^0 A- X- }

, u0 Q# ^4 {7 c2 N) Q$ L8 }function dy= SBR1(t,y,TC)
* O3 V& K4 d4 v  Jdy=zeros(2,1);3 U8 O2 ^, b5 k9 r; W
VADARE=14;
) \4 f0 ]5 v# p- q! QEPS=0.55;dtad0=0.29;UDA=10*EPS;
/ K6 n" T2 v1 u9 Y- }* wRH=1;Y=36;
# p& u/ |( [0 ^/ U1 @' YTR=300;TD=TC; YD=TD/TR;- L, T$ q) E) E3 D5 \
if t<=10 C/ a2 U4 m& Z* X. e& g/ ~- |
    YCE=((UDA*(1+EPS*t)*TC+RH*EPS*TD)/(UDA*(1+EPS*t)+RH*EPS))/TR;
1 n: C6 k9 t: [5 O/ \$ l# ^    dy(1) = VADARE*(t-y(1))*(1-y(1))*exp(Y*(1-1/y(2)));
1 S6 g* ?) _7 m' h1 z# [: J4 [    dy(2) = ((dtad0*dy(1)-(UDA*(1+EPS*t)+RH*EPS)*(y(2)-YCE)))/(1+RH*EPS*t);
. H, _# w: E7 O+ z/ f2 t1 Dend
9 ~" p# c1 n$ _9 Vif t>1
5 M; s' |5 w" s   YCE=((UDA*(1+EPS)*TC+RH*EPS*TD)/(UDA*(1+EPS)+RH*EPS))/TR;
" Z7 ?) N& Y7 i   dy(1) = VADARE*(1-y(1))*(1-y(1))*exp(Y*(1-1/y(2)));$ r2 P/ f3 Q8 Q4 N
   dy(2) = ((dtad0*dy(1)-(UDA*(1+EPS)+RH*EPS)*(y(2)-YCE)))/(1+RH*EPS);2 R6 k# B1 h6 S4 s% F+ K7 A" F1 C
end- C3 C0 s# i7 f% z, P; R5 v$ o
end0 M* I* {% J( y: I4 f; f$ f

2 L5 N& p0 a0 O3 \3 ~) Q$ b+ A: r$ }& \3 G" T

该用户从未签到

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产生的误差最小。
7 K8 N5 Q+ g$ ]9 m4 G, S
% x( S" ]& j4 D$ h( D6 {  `/ t6 O最简单的例子,已知模型:电压U = 电阻R * 电流I,且不知道接在某简单电路上的电阻到底是多大。那么通过若干组实验,以电流I为输入(比如调整滑动变阻器来控制该电路中的电流数值),同时测量这个电阻两端的电压U作为输出,通过这若干组U与I的数值,可以估计出一个满足模型的电阻R值。控制输入I、测量输出U、以此二者来估计参数R。/ C/ }6 B6 e, X' d
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-8-5 03:03 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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