我怎么能符合双指数曲线吗?

6视图(30天)

显示旧的评论

阿曼达Botelho阿马拉尔 2021年12月18日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1613675-how-can-i-fit-an-double-exponential-curve

评论道: 阿曼达Botelho阿马拉尔2021年12月20日

这就跟你问声好！我有一个问题拟合曲线。我有两个曲线,电流一电压,另,我需要Z = V /阻抗。我需要做的就是:找到适合的电压曲线。我发现拟合曲线和新新的阻抗。

我做了一个曲线的代码,我得到了一个不错的选择,然而,当策划新的阻抗不出来任何类似于“原始”(曲线没有调整)。

有人能帮助我吗?

这是我的代码和数据。

                          
                          负载(500 _constante_2u_v.mat);
                         
                          负载(500 _constante_2u_c.mat);
                         
                          负载(“t.MAt”);
                         
                          t = t (2:200000);
                         
                          % Impedancia
                         
                          Z_Trad_Primeira_500 = tensao_at_500. / corrente_at_500;
                         
                          图(1)
                         
                          情节(t * 10 ^ 3 Z_Trad_Primeira_500“r”。);在;
                         
                          %配合
                         
                          fcn1 = @ (b, b t) (1)。* exp (b (2)。* t) + b (3)。* exp (b (4)。* t);
                         
                          [f1, fval] = fminsearch (@ (b)规范(tensao_at_500 fcn1 (b、t)), (1;1;1;1);
                         
                          图(2)
                         
                          情节(t, tensao_at_500“p”)
                         
                          持有在
                         
                          情节(t, fcn1 (f1, t),“- - -”)
                         
                          持有从
                         
                          网格
                         
                          新Z %
                         
                          Znew = fcn1 (f1, t) / corrente_at_500;
                         
                          图(5)
                         
                          情节(t * 10 ^ 3 Z_Trad_Primeira_500“r”。);在;
                         
                          情节(t * 10 ^ 3 Znew“g——”);在;

16条评论
显示15个年长的评论隐藏15年长的评论

阿曼达Botelho阿马拉尔 2021年12月18日

编辑:阿曼达Botelho阿马拉尔 2021年12月18日

感谢你的回复!测试所有点,我找不到代码中的错误。我修改了方程,得到更好的近似。

fcn1 = @ (b, b t) (1)。* exp (b (2)。* t) + b (3)。* exp (b (4)。* t);

让我解释原因分析。正如我提到的,我想要分析阻抗值系统。我有一个默认的系统,我添加一些变量进行分析,所以对每一个改变我我将有一个新的曲线。

我现在的研究分为两个阶段:1)我必须检查曲线之间的差异,并检查每个参数的影响。(我使用dtw在这种情况下,如果你知道另一个分析我可以做我很感激)。2)试图让标准系统曲线的曲线。这是我在上面的代码分析。从方程的参数我可以检查每一个的变化根据参数的变化。

理解我在想什么吗?

马修•诺伊 2021年12月20日

你好,阿曼达

欢迎回来

据我了解的出版你参照,尝试做一个简单的两个参数指数适合一次电压信号,电流信号

这样做是在这里:

                                
                                负载(“tensao_at_500.MAT”);
                               
                                负载(“corrente_at_500.MAT”);
                               
                                负载(“t.MAt”);
                               
                                图(1)
                               
                                情节(t, tensao_at_500“b”t corrente_at_500“r”);
                               
                                标题(电压和电流的);
                               
                                包含(“时间(s)”);
                               
                                ylabel (“振幅”);
                               
                                %合适的电压信号
                               
                                %只使用衰减信号的一部分
                               
                                [val, indm] = max (tensao_at_500);
                               
                                tensao_fit = tensao_at_500 (indm:结束);
                               
                                t_fit = t (indm:结束);
                               
                                fcn1 = @ (b, b t) (1)。* exp (b (2)。* t);
                               
                                [f1, fval] = fminsearch (@ (b)规范(tensao_fit fcn1 (b, t_fit)), (val;1 e3;]);
                               
                                图(2)
                               
                                情节(t, tensao_at_500‘*’t_fit t_fit fcn1 (f1),“- - -”)
                               
                                标题(“电压”);
                               
                                传奇(“电压”,实验符合的);
                               
                                包含(“时间(s)”);
                               
                                ylabel (“振幅”);
                               
                                网格在
                               
                                %配合电流信号
                               
                                %只使用衰减信号的一部分
                               
                                [val, indm] = max (corrente_at_500);
                               
                                corrente_fit = corrente_at_500 (indm:结束);
                               
                                t_fit = t (indm:结束);
                               
                                fcn2 = @ (b, b t) (1)。* exp (b (2)。* t);
                               
                                [f2, fval] = fminsearch (@ (b)规范(corrente_fit fcn2 (b, t_fit)), (val;1 e3;]);
                               
                                图(3)
                               
                                情节(t, corrente_at_500‘*’t_fit t_fit fcn2 (f2),“- - -”)
                               
                                标题(“当前”);
                               
                                传奇(“当前”,实验符合的);
                               
                                包含(“时间(s)”);
                               
                                ylabel (“振幅”);
                               
                                网格在

阿曼达Botelho阿马拉尔 2021年12月20日

“只使用的腐烂的部分信号调整”

这正是我不认为!非常感谢!

登录置评。

在回答这个问题。

答案(0)

在回答这个问题。

类别

射频和混合信号混合信号Blockset 测量和Testbenches

找到更多的在测量和Testbenches在帮助中心和文件交换

s manbetx 845

MATLAB

释放

R2013a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!