如何找到最适合的曲线

4次观看(过去30天)
Stashu Kozlowski.
Stashu Kozlowski. 2021年3月20日
评论: 星节演唱 2021年3月20日
我有以下一组数据。
从它的外观,它沿着日志正常分布的线路遵循一些东西。但我难以拟合一个符合物的日志规范。我怎么能适应它?我也已附加了数据文件以进行参考。
更新:我已经意识到数据未归一化,因此您无法符合正确的日志规范曲线。所以我首先需要找到数据下的区域。

答案(2)

星节演唱
星节演唱 2021年3月20日
使用完全不同的方法(系统识别工具箱):
D = readmatrix (“RLC1.xls”);
t = D(:,1);
S = D(:,2);
Ts = =平均值(差异(t));
tsd = std(diff(t));
FS = 0.25;
[SR,TR] =重组(S,T,FS);
fn = fs / 2;
数字
情节(TR,SR)
网格
L =元素个数(tr);
FTS = FFT(s)/ l;
fv = linspace(0,1,修复(L / 2)+1)* Fn;
4 = 1:元素个数(艘);
数字
情节(阵线,图像放大(FTs (Iv)))
网格
ylim ([1] * 0.0005)
u = [1;0 (l - 1, 1)];
TDD = IDDATA(SR,U,1 / FS);
sys =特遣部队(tdd、3、2)
数字
比较(tdd系统)
产生此传输功能数据:
sys =
输入“U1”输出“Y1”:
0.007132 S ^ 2 + 0.0004604 s + 1.881e-06
-------------------------------------
S ^3 + 0.01617 S ^2 + 9.041 -05 S + 7.589e-08
连续时间确定传递函数。
参数化:
数字杆子:3零数量:2
数字自由系数:6
使用“tfdata”“getpvec”“getcov”为了参数和他们的不确定因素。
地位:
估计的使用TFEST在时域数据“TDD”上。
适合估计数据:95.47
FPE:4.55E-06,MSE:4.304E-06
还有这个 相比 阴谋:
4评论
星节演唱
星节演唱 2021年3月20日
如果你想适合它 “lnf” 函数,您可能需要添加第三个“缩放”参数,以便它适合曲线:
lnf = @(p,x)p(3)。* exp( - (log(x)-p(1))。^ 2./(2*p(2).^2))./(x。* P(2)* SQRT(2 * PI));
p = fminsearch(@(p)规范(S-min(s) - lnf(p,t)),rand(3,1));
一个参数集 [5,0.5,75] 产生一个合理的适合,所以用那些代替 兰特 调用作为初始参数估计。

登录评论。


沃尔特罗伯森
沃尔特罗伯森 2021年3月20日
2度的两个多项式的比率相当合适,具有近似系数
P = [0.00213762295171613 56.4867705538076 45.003715200478]
Q = [1 -51.094359290793 32989.9993982203]
Yfit = polyval(p,x)./ polyval(q,x);
这是使用CFTOOL,'rat22'合身。

标签

s manbetx 845


释放

R2020B.

社区宝藏狩猎

找到Matlab Central中的宝藏,并发现社区如何帮助您!

开始狩猎!