MATLAB曲线相减

问题

现有两条曲线($\vec{x}1$,$\vec{y}1$),($\vec{x}2$,$\vec{y}2$)，两条曲线的坐标不一定对齐，点数不一定相同，即 $\vec{x}1\ne \vec{x}2$，现要求两条曲线的差值。

例

要计算下面两个曲线的差值

x1 = 0:0.5:2*pi;
y1 = sin(x1);
x2 = 0:0.45:2.5*pi;
y2 = cos(x2);

plot(x1,y1,':o');
hold on
plot(x2,y2,':*');

解决方法

在网上找了一圈，没找到相关的函数，自己写了一个函数

% 求y1-y2的差值
function [x,y] = CurveSubtraction(x1,y1,x2,y2)
fprintf('the length of x1 and x2 is：[%d] and [%d]\n',length(x1),length(x2));

% 排序
[x1,index] = sort(x1);
y1 = y1(index);
[x2,index] = sort(x2);
y2 = y2(index);

% 删除重复元素
[x1,ia] = unique(x1);
y1 = y1(ia);
[x2,ia] = unique(x2);
y2 = y2(ia);

min_x = max(min(x1),min(x2)); %求出横坐标的公共部分
max_x = min(max(x1),max(x2));

x1_cut_index = find(x1>=min_x&x1<=max_x); %裁剪到起点和终点后的坐标索引
x2_cut_index = find(x2>=min_x&x2<=max_x);
len_x1 = length(x1_cut_index); %裁剪之后的向量长度
len_x2 = length(x2_cut_index);
fprintf('the range of common subsequence:[%d-%d]. x1 has [%d] elements，x2 has [%d] elements\n',min_x,max_x,len_x1,len_x2)
x1_cut = x1(x1_cut_index); % 裁剪到相同起点和终点后的横纵坐标
y1_cut = y1(x1_cut_index);
x2_cut = x2(x2_cut_index);
y2_cut = y2(x2_cut_index);

if len_x1>len_x2 %如果y1点数多余y2点数，把y2插值到和y1一样多，然后相减
    inter_y2 = interp1(x2_cut,y2_cut,x1_cut);
    x = x1_cut;
    y = y1_cut - inter_y2;
    fprintf('Interpolating y2');
else % 否则，把y1插值到和y2一样多，然后相减
    inter_y1 = interp1(x1_cut,y1_cut,x2_cut);
    x = x2_cut;
    y = inter_y1 - y2_cut;
    fprintf('Interpolating y2');
end
end

测试

x1 = 0:0.5:2*pi;
y1 = sin(x1);
x2 = 0:0.45:2.5*pi;
y2 = cos(x2);

plot(x1,y1,':o');
hold on
plot(x2,y2,':*');
%% 
[x,y] = CurveSubtraction(x1,y1,x2,y2);
hold on
plot(x,y,'-^');
line([0 8],[0 0]);
legend('y1','y2','y1-y1')