请给出一个完整matlab程序，再尽量不调用程序内函数时使其满足一下条件：给定一个倾斜界面的叠加剖面； 根据射线偏移理论对该叠后剖面进行偏移处理自激自收得到的反射信息对应的反射点可能来自以vt/2为半径，以自激自收点为圆心的圆弧上的任意一点。得到来自界面上两个点的反射波旅行 t1 和 t2 及波速v时，分别以两个点为圆心，对应的时间和速度乘积为半径画圆，做两个圆的公切线可以得到反射点的位置和反射段位置。其中，画出六个道集画圆结果； 对比分析偏移结果和叠加剖面。

由于题目中提到了许多细节，这里给出的程序仅供参考，实际应用时需要根据具体情况进行修改和调整。

% 假设已经读入了叠加剖面和射线信息，分别为stack和rays
% stack为一个n x m的矩阵，表示n条道集上的m个采样点
% rays为一个n x 2的矩阵，表示n条道集上的自激自收点和反射点的位置

% 设置参数
vt = 2000;  % 界面上的波速
n = size(stack, 1);  % 道集数
m = size(stack, 2);  % 采样点数

% 初始化反射段和反射点位置
refl_segs = zeros(n, 4);  % 每条道集上的反射段位置，用左右两个端点表示
refl_pts = zeros(n, 2);  % 每条道集上的反射点位置

% 对每条道集进行偏移处理
for i = 1:n
    % 获取自激自收点和反射点位置
    sx = rays(i, 1);
    gx = rays(i, 2);
    
    % 获取对应的旅行时间和波速
    t1 = abs(sx - refl_segs(i, 1)) / vt;
    t2 = abs(gx - refl_segs(i, 2)) / vt;
    v1 = refl_segs(i, 3);
    v2 = refl_segs(i, 4);
    
    % 计算圆心和半径
    cx1 = (t1 * v1^2 - t2 * v2^2) / (2 * (v1^2 - v2^2)) + sx;
    cx2 = (t2 * v2^2 - t1 * v1^2) / (2 * (v2^2 - v1^2)) + gx;
    cy1 = sqrt((vt/2)^2 - (cx1 - sx)^2);
    cy2 = sqrt((vt/2)^2 - (cx2 - gx)^2);
    r1 = v1 * t1;
    r2 = v2 * t2;
    
    % 计算公切线斜率
    k = (cy2 - cy1) / (cx2 - cx1);
    
    % 计算反射点位置
    refl_x = (r2^2 - r1^2 - cx2^2 + cx1^2 - cy2^2 + cy1^2) / (2 * (cx1 - cx2));
    refl_y = sqrt(r1^2 - (refl_x - cx1)^2);
    refl_pts(i, :) = [refl_x, refl_y];
    
    % 计算反射段位置
    x1 = refl_x - refl_y / k;
    y1 = 0;
    x2 = refl_x + (m - refl_y) / k;
    y2 = m;
    refl_segs(i, :) = [x1, y1, x2, y2];
    
    % 画圆和公切线
    figure;
    plot(1:m, stack(i, :));
    hold on;
    axis([1 m -1000 1000]);
    circle(cx1, cy1, r1);
    circle(cx2, cy2, r2);
    line([sx, gx], [0, m], 'Color', 'r');
    line([refl_segs(i, 1), refl_segs(i, 3)], [refl_segs(i, 2), refl_segs(i, 4)], 'Color', 'g');
    line([refl_x - refl_y / k, refl_x + (m - refl_y) / k], [0, m], 'Color', 'b');
end

% 绘制反射点和反射段位置
figure;
for i = 1:n
    subplot(n, 2, 2*i-1);
    plot(1:m, stack(i, :));
    hold on;
    axis([1 m -1000 1000]);
    plot(rays(i, :), [0, m], 'Color', 'r');
    plot(refl_pts(i, 1), refl_pts(i, 2), 'x', 'Color', 'g');
    title(sprintf('Stack and ray path for trace %d', i));
    
    subplot(n, 2, 2*i);
    plot(1:m, stack(i, :));
    hold on;
    axis([1 m -1000 1000]);
    line([refl_segs(i, 1), refl_segs(i, 3)], [refl_segs(i, 2), refl_segs(i, 4)], 'Color', 'g');
    plot(refl_pts(i, 1), refl_pts(i, 2), 'x', 'Color', 'r');
    title(sprintf('Stack and reflection segment for trace %d', i));
end

% 绘制偏移结果和叠加剖面
figure;
for i = 1:n
    subplot(n, 2, 2*i-1);
    plot(1:m, stack(i, :));
    hold on;
    axis([1 m -1000 1000]);
    plot(refl_pts(i, 1), refl_pts(i, 2), 'x', 'Color', 'g');
    title(sprintf('Stack and reflection point for trace %d', i));
    
    subplot(n, 2, 2*i);
    plot(1:m, stack(i, :));
    hold on;
    axis([1 m -1000 1000]);
    line([refl_segs(i, 1), refl_segs(i, 3)], [refl_segs(i, 2), refl_segs(i, 4)], 'Color', 'r');
    title(sprintf('Stack and migrated trace for trace %d', i));
end

% 画圆函数
function circle(x, y, r)
    ang = 0:0.01:2*pi;
    xp = r*cos(ang);
    yp = r*sin(ang);
    plot(x+xp, y+yp);
end