s1 s2符合分布f(x)，怎么用matlab求y=|s1-s2|的分布函数

可以使用以下步骤来求解问题：

1. 首先，需要求出s1-s2的概率密度函数。假设f(x)为s1和s2的分布函数，则s1-s2的概率密度函数为：

g(x) = ∫∫f(s1)f(s2)δ(x-(s1-s2))ds1ds2

其中，δ(x)为Dirac delta函数，表示x=0时取值为无穷大，其余位置为0。在Matlab中，可以使用`integral2`函数来求解二重积分。

2. 然后，需要求解y=|s1-s2|的分布函数。由于|s1-s2|可以表示为max(s1,s2)-min(s1,s2)，因此可以考虑使用max和min函数来求解。

假设g(x)为s1-s2的概率密度函数，则y的概率密度函数为：

h(y) = 2∫g(x)F(y+x)dx

其中，F(x)为分布函数的累积分布函数。在Matlab中，可以使用`integral`函数来求解一重积分。

下面是一个具体的Matlab代码示例：

% 假设s1和s2符合正态分布
mu = 0;
sigma = 1;
f = @(x) normpdf(x, mu, sigma);

% 求解s1-s2的概率密度函数
g = @(x) integral2(@(s1, s2) f(s1).*f(s2).*dirac(x-(s1-s2)), -Inf, Inf, -Inf, Inf);

% 求解y的概率密度函数
h = @(y) 2*integral(@(x) g(x).*normcdf(y+x, mu, sigma), -Inf, Inf);

% 绘制概率密度函数图像
y = 0:0.1:5;
plot(y, h(y));
xlabel('y');
ylabel('h(y)');

在上述代码中，使用了`normpdf`函数求解正态分布函数的概率密度函数，使用了`normcdf`函数求解正态分布函数的累积分布函数。