m=[-4.1 -4 -3.7 -3.4 -2.7 -2 -1.6 0 1.6 2 2.9 3.5 3.9 4 4.2];
n=[10 9.1 6 4 2 0.8 0.5 0 0.5 0.8 2 4 6 8 10];
plot(m,n,-m,n)
f(-m)曲线在f(m)左侧,说明原曲线n=f(m)形心偏右。
yi=0:0.1:10;
xi=interp1(n(1:7),m(1:7),yi);
xi_=interp1(n(9:15),m(9:15),yi);
x_mean=mean([xi(6:end),xi_(6:end)]);
plot(m-x_mean,n,x_mean-m,n)