用MATLAB帮我写代码，读取给定振动视频，利用基于直线检测的位移提取算法对给定视频进行处理，提取结构多个测点的动态变形曲线,以无振动时提取位移的均方根值作为评价算法的指标或者多种算法的提取结果相互比较

好的，下面是MATLAB代码，实现基于直线检测的位移提取算法并提取结构多个测点的动态变形曲线：

% 读取给定振动视频
video = VideoReader('vibration_video.mp4');

% 设定参数
window_size = 21; % 窗口大小
sigma = 3; % 高斯核标准差
threshold = 0.2; % 非极大值抑制阈值
min_length = 30; % 直线最小长度

% 初始化变量
frame_num = video.NumFrames;
displacement = zeros(frame_num, num_points); % 多个测点的位移
xy_points = cell(num_points, 1); % 多个测点的坐标

% 循环处理每一帧
for k = 1:frame_num
    % 读取视频帧
    frame = read(video, k);
    % 灰度化
    gray_frame = rgb2gray(frame);
    % 高斯滤波
    smooth_frame = imgaussfilt(gray_frame, sigma);
    % 边缘检测
    edge_frame = edge(smooth_frame, 'canny');
    % 直线检测
    [H, theta, rho] = hough(edge_frame);
    peaks = houghpeaks(H, 10, 'threshold', ceil(0.3*max(H(:))));
    lines = houghlines(edge_frame, theta, rho, peaks, 'FillGap', window_size, 'MinLength', min_length);
    
    % 计算位移
    for i = 1:num_points
        % 获取当前测点坐标
        x = xy_points{i}(1);
        y = xy_points{i}(2);
        % 循环处理每条直线
        for j = 1:length(lines)
            % 获取直线端点坐标
            x1 = lines(j).point1(1);
            y1 = lines(j).point1(2);
            x2 = lines(j).point2(1);
            y2 = lines(j).point2(2);
            % 判断直线是否经过当前测点
            if abs((y2-y1)*x - (x2-x1)*y + x2*y1 - y2*x1) / sqrt((y2-y1)^2 + (x2-x1)^2) < threshold
                % 计算位移
                displacement(k, i) = displacement(k, i) + sqrt((x-x1)^2 + (y-y1)^2) - sqrt((x-x2)^2 + (y-y2)^2);
            end
        end
    end
end

% 计算均方根值
rms = sqrt(mean(displacement.^2, 1));

在上面的代码中，我们首先读取给定的振动视频，然后设定了一些参数，例如窗口大小、高斯核标准差、非极大值抑制阈值和直线最小长度等。接下来，我们循环处理每一帧，将其灰度化并进行高斯滤波和边缘检测。然后，我们使用霍夫变换检测直线，并计算每个测点的位移。最后，我们计算多个测点的位移的均方根值作为评价算法的指标。

如果需要对多种算法的提取结果相互比较，可以在上面代码的基础上再添加其他算法的实现，并将它们的结果与上面的基于直线检测的位移提取算法的结果进行比较。