matlab提取振动视频的位移时程曲线的代码

Matlab提取振动视频的位移时程曲线可以使用光流法实现。以下是一个简单的Matlab代码示例：

% 读取视频文件
vid = VideoReader('video.mp4');

% 设置参数
winSize = 15;
maxLevel = 2;
criteria = struct('type', 'Count+EPS', 'maxCount', 10, 'epsilon', 0.03);

% 初始化
prev_frame = readFrame(vid);
prev_gray = rgb2gray(prev_frame);
prev_pts = detectMinEigenFeatures(prev_gray, 'MinQuality', 0.01);
prev_pts = prev_pts.Location;

% 创建空数组存储位移量和时间
displacements = [];
times = [];

% 循环读取视频帧
while hasFrame(vid)
    frame = readFrame(vid);

    % 转换为灰度图像
    gray = rgb2gray(frame);

    % 通过光流法获取位移量
    [next_pts, status] = ...
        opticalFlowLK(prev_gray, gray, prev_pts, 'Window', winSize, 'MaxLevel', maxLevel, 'Criteria', criteria);

    % 计算位移量
    displacement = mean(next_pts - prev_pts);

    % 将位移量添加到数组中
    displacements = [displacements; displacement];

    % 计算时间
    time = vid.CurrentTime;
    times = [times; time];

    % 更新下一帧
    prev_gray = gray;
    prev_pts = next_pts;
end

% 绘制时程曲线图
plot(times, displacements(:, 1), 'r-', 'LineWidth', 2);
hold on;
plot(times, displacements(:, 2), 'b-', 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Displacement (pixels)');
legend('x', 'y');
grid on;

其中，`opticalFlowLK`是光流法函数，用于计算两帧图像之间的位移量；`detectMinEigenFeatures`函数用于检测图像中的角点。代码中的`winSize`、`maxLevel`和`criteria`参数可以根据实际情况进行调整。