如何构造存储交点坐标的矩阵。(数据数组和定义曲线)

1查看(最近30天)
Brenna狼
Brenna狼 在2021年1月28日
编辑: 马特·J 在2021年1月28日
我有一个包含大量数据的文本文件,我已经将其加载并命名为vect。我需要知道阻尼正弦数据与数据平均值相交的所有坐标。
所有intersctions将发生在Y =平均水平,但我需要的时间每个路口。(试图确定oscilation的周期)
我现在有代码选择替代的采样速率和曲线数据,并平均,但在尝试使用P =相交(VECT,zavg)不输出响应。
raw_data =负载(500 c.txt);
frequency_of_choice = 15;
X = raw_data(:,1);
X1 = X(frequency_of_choice:frequency_of_choice:端,:);
结束
y = raw_data (:, 2);
Y1 = Y(frequency_of_choice:frequency_of_choice:端,:);
结束
x0 = raw_data (1,1);
y0 = raw_data(1、2);
x = (x0, x1);
y = (y0, y1);
VECT = [X,Y];
avg (x) =意味着
L =长度(x);
XZ =零([L 1]);
zavg = XZ +平均;
XMIN = 2.785;
XMAX = 6.48;
P =相交(VECT,zavg)
数字
情节(Y,X)
xlim ([xmin xmax])
持有
情节(zavg“r”
持有离开
1评论
马特·J
马特·J 在2021年1月28日
所有intersctions将发生在Y =平均水平,但我需要的时间每个路口。(试图确定oscilation的周期)
请注意,这是获得振荡周期的一个非常近似的方式。我建议你这样做只是为了开发一个初始参数猜测适当的曲线拟合,例如,使用

登录评论。

答案(1)

马特·J
马特·J 在2021年1月28日
编辑:马特·J 在2021年1月28日
x = linspace(0 8 *π,1000);%例x, y数据
y = exp (-0.1 * x)。* sin (x) + 3;
xc = x (diff(符号(y-mean (y))) ~ = 0);%交点
情节(x, y xc 0 * xc + (y),“罗”);
6个评论
马特·J
马特·J 在2021年1月28日
似乎还好我,
raw_data =负载(websave (500 c.txt“//www.tatmou.com/matlabcentral/answers/uploaded_files/502108/500c.txt”));
Y = raw_data(:,1);
x = raw_data (:, 2);
xc = x (diff(符号(y-mean (y))) ~ = 0);%交点
情节(x, y xc 0 * xc + (y),“罗”);
xlim ([4.5, 5.8])

登录评论。

下载188bet金宝搏


释放

R2020b

社区寻宝

在MATLAB中心找到宝藏,并发现社区如何可以帮助你!

开始狩猎!