如何找到曲线的拐点

56个观点(过去30天)

显示旧的评论

harith larabi 2022年10月4日14:22

0
链接

直接链接到这个问题

//www.ru-cchi.com/matlabcentral/answers/1817340-how-to-find-inflection-point-of-a-curve

评论道: 明星黾 2022年10月6日11:33

答:接受明星黾

你好,伙计们。

我在找一些信息。

也许你能帮我?

我需要找到曲线的拐点。

我知道我可以用二阶导数或正切但我不知道如何在matlab中应用或计算它?

在这里你会找到曲线

它代表了从弹性域到塑性域的过渡。

谢谢你！

问候

Harith

接受的答案

明星黾 2022年10月4日14:28

2
链接

直接链接到这个答案

//www.ru-cchi.com/matlabcentral/answers/1817340-how-to-find-inflection-point-of-a-curve#answer_1066475

编辑:明星黾 2022年10月4日20:40

首先,使用 sgolayfilt 函数来消除尽可能多的噪声，因为求导会放大它。使用的梯度函数来计算数值导数，例如:

framelen = 101;

Ys = sgolayfilt(y, 3, framelen);

Dydx = gradient(ys) ./ gradient(x);

选择合适的值 “framelen” 以得到想要的结果。

编辑- (2022年10月4日20:40)

的例子,

X = linspace(0.9,1.8, 250);%创建数据

y = exp (- (x - 1.3) ^ 2 * 50) / 8 + x + 0.005 * randn(大小(x));%创建数据

%图

%的阴谋(x, y)

%的网格

ydt =去趋势(y, 1);%趋势'y'以促进分析

framelen = 101;%选择合适的值

ydt = sgolayfilt(ydt, 3, framelen);%降噪信号

Dydx =梯度(ydts) ./梯度(x);%计算数值导数

[maxdydx, idxmax] = max(就要);插值指数下限

[mindydx, idxmin] = min(就要);插值指数上限

Idxrng = idxmax: idxmin;

Inflptx = interp1(dydx(idxrng)， x(idxrng)， 0)%找到拐点x值

inflptx = 1.3014

Inflpty = interp1(x, y, inflptx)找到拐点y值

inflpty = 1.4191

数字

情节(x, y,“DisplayName的”，“数据”）

持有在

情节(inflptx inflpty,“+ r”，“MarkerSize”10“DisplayName的”，“拐点”）

持有从

网格

包含(“X”）

ylabel (“Y”）

传奇(“位置”，“最佳”）

．

8的评论
显示隐藏 7年长的评论

harith larabi 2022年10月5日13:40

B0Inflection.xlsx

你好,谢谢你。

在变量管理器x有数据的和为dydx。

你会发现附加的。xlxs，我运行的代码和abce是一样的

                                   X = linspace(0.9,1.8, 250);%创建数据
                                  
                                   y = exp (- (x - 1.3) ^ 2 * 50) / 8 + x + 0.005 * randn(大小(x));%创建数据
                                  
                                   %图
                                  
                                   %的阴谋(x, y)
                                  
                                   %的网格
                                  
                                   ydt =去趋势(y, 1);%趋势'y'以促进分析
                                  
                                   framelen = 101;%选择合适的值
                                  
                                   ydt = sgolayfilt(ydt, 3, framelen);%降噪信号
                                  
                                   Dydx =梯度(ydts) ./梯度(x);%计算数值导数
                                  
                                   [maxdydx, idxmax] = max(就要);插值指数下限
                                  
                                   [mindydx, idxmin] = min(就要);插值指数上限
                                  
                                   Idxrng = idxmax: idxmin;
                                  
                                   Inflptx = interp1(dydx(idxrng)， x(idxrng)， 0)%找到拐点x值
                                  
                                   Inflpty = interp1(x, y, inflptx)找到拐点y值
                                  
                                   数字
                                  
                                   情节(x, y,“DisplayName的”，“数据”）
                                  
                                   持有在
                                  
                                   情节(inflptx inflpty,“+ r”，“MarkerSize”10“DisplayName的”，“拐点”）
                                  
                                   持有从
                                  
                                   网格
                                  
                                   包含(“X”）
                                  
                                   ylabel (“Y”）
                                  
                                   传奇(“位置”，“最佳”）

谢谢你,Harith

明星黾 2022年10月5日14:33

有重复的 “x” 需要首先删除的值。

结果表明，对信号进行去噪和对数值导数-分别去噪后，会出现两个拐点

T1 = readtable (“//www.ru-cchi.com/matlabcentral/answers/uploaded_files/1146155/B0Inflection.xlsx”）

T1 =296×7表

Var1 Var2 Var3 Var4 Var5 Var6 Var7 ＿＿＿＿ _____ _____ _________ ＿＿＿＿ ＿＿＿＿ ____________0.2 66.25 -0.42 0.92014南南{' X =行D '} 0.21 66.75 -0.43 0.92708南南{' Y =行'}0.21 66.75 -0.44 0.92708南南{0×0字符}0.23 67 -0.44 - 0.93056南南{0×0字符}0.23 67.25 -0.46 0.93403南南{0×0字符}0.25 67.25 -0.47 0.93403南南{0×0字符}0.26 67.5 -0.47 0.9375南南{0×0字符}0.25 67.5 -0.49 0.9375南南{0×0字符}0.27 67.75 -0.48 0.94097南南{0×0字符}0.27 67.75 -0.49 0.94097南南{0×0字符}0.28 68.25 -0.51 0.94792南南{0×0字符}0.28 - 68.25-0.5 0.94792 NaN NaN {0×0 char} 0.3 68.5 -0.51 0.95139 NaN NaN {0×0 char} 0.31 68.75 -0.52 0.95486 NaN NaN {0×0 char} 0.31 68.75 -0.52 0.95486 NaN NaN {0×0 char} 0.33 69 -0.54 0.95833 NaN NaN {0×0 char}

[x,第九]=独特(T1 {: 4});

第九T1 y = {1};

数字

情节(x, y)

网格

包含(“X”）

ylabel (“Y”）

标题(“原始数据”）

ydt =去趋势(y, 1);%趋势'y'以促进分析

数字

情节(x, ydt)

网格

包含(“X”）

ylabel (“Y”）

标题(“去趋势数据”）

framelen = 101;%选择合适的值

ydt = sgolayfilt(ydt, 3, framelen);%降噪信号

Dydx =梯度(ydts) ./梯度(x);%计算数值导数

Dydxs = sgolayfilt(dydx, 3, framelen);去噪数值导数

数字

情节(x, ydts,“DisplayName的”，“平滑和去趋势数据”）

持有在

情节(x,就要,“DisplayName的”，平滑数值导数的）

持有从

网格

包含(“X”）

ylabel (“Y”）

传奇(“位置”，“最佳”）

zxi =找到(diff(签署(就要)));

为k = 1:元素个数(zxi)

idxrng = max (1, zxi (k) 1): min(元素个数(x) zxi (k) + 1);

Inflptx (k) = interp1(dydxs(idxrng)， x(idxrng)， 0);

Inflpty (k) = interp1(x(idxrng)， y(idxrng)， inflptx(k));

结束

InflectionPoints = [inflptx;inflpty]

InflectionPoints = 2×2

1.0117 1.3536 0.5737 4.0925

数字

情节(x, y,“DisplayName的”，“数据”）

持有在

情节(inflptx inflpty,“+ r”，“MarkerSize”10“DisplayName的”，“拐点”）

持有从

网格

包含(“X”）

ylabel (“Y”）

传奇(“位置”，“最佳”）

．

明星黾 2022年10月6日11:33

我的荣幸!

如果我的回答能帮你解决问题，请接受它!

．

登录发表评论。

类别

MATLAB 语言基础知识数据类型表

2022世界杯八强谁会赢？

MATLAB

世界杯预选赛小组名单社区寻宝

在MATLAB中央找到宝藏，并发现社区可以如何帮助你!世界杯预选赛小组名单

开始狩猎!