用MATLAB写出1、随机变量X服从参数为试验次数20，概率为0.25的二项分布。 （1）生成X的概率分布律； （2）产生18个随机数（3行6列）； （3）已知 ，求 （4）求 （5）画出X的分布律和分布函数图形。

好的，下面是MATLAB的代码实现：

1. 生成X的概率分布律

n = 20; % 试验次数
p = 0.25; % 概率
x = 0:n; % 可能的取值
pmf = binopdf(x, n, p); % 生成概率分布律
disp('X的概率分布律：');
disp([x; pmf]);

2. 产生18个随机数

r = binornd(n, p, 3, 6); % 产生18个随机数
disp('随机数矩阵：');
disp(r);

3. 已知 X=10，求 P(X=10)

x0 = 10; % 已知 X=10
px0 = binopdf(x0, n, p); % 计算 P(X=10)
disp(['P(X=' num2str(x0) ')=' num2str(px0)]);

4. 求 P(X≥15)

x1 = 15; % X的取值大于等于15
px1 = 1 - binocdf(x1-1, n, p); % 计算 P(X≥15)
disp(['P(X≥' num2str(x1) ')=' num2str(px1)]);

5. 画出X的分布律和分布函数图形

figure;
subplot(2,1,1);
stem(x, pmf, 'filled');
title('X的概率分布律');
xlabel('X');
ylabel('P(X)');
grid on;

subplot(2,1,2);
cdf = binocdf(x, n, p); % 生成分布函数
plot(x, cdf);
title('X的分布函数');
xlabel('X');
ylabel('F(X)');
grid on;

注意：在第5步中，为了更好的展示分布函数，我使用了生成分布函数的函数 `binocdf`，而不是手动计算。