那应该如何修改呢，你能给我修改完后的代码吗

当 `r` 的大小为 `50`，而 `t` 的大小为 `1000` 时，可以猜测是在瑞利衰落信道模拟时出现了问题。可以将代码中的瑞利衰落信道模拟部分修改为如下代码：

% 瑞利衰落信道模拟
fc = 900e6; % 载波频率
v = 50; % 移动速度，单位：km/h
fd = v / 3.6 * fc / 3e8; % 多普勒频移
h = sqrt(0.5) * (randn(size(t)) + 1i*randn(size(t)));
r = r .* exp(1i*2*pi*fd*t) .* h;

这里使用 `t` 的大小来生成瑞利衰落信道，从而保证 `r` 和 `t` 的大小相同。另外，使用 `randn` 函数生成的高斯随机数向量 `h` 的大小也应该与 `t` 的大小相同，从而保证 `r` 和 `h` 的大小相同。修改后的完整代码如下：

% 随机离散信号生成
x = randi([0 1], 1, 1000);

% 调制
fc = 10e6; % 载波频率
t = (0:length(x)-1) / fc;
s = 2*x-1; % BPSK调制
y = s .* cos(2*pi*fc*t);

% 采样
fs = 20*fc; % 采样频率
ts = 1/fs; % 采样时间间隔
r = y(1:fs/fc:end);

% 瑞利衰落信道模拟
fc = 900e6; % 载波频率
v = 50; % 移动速度，单位：km/h
fd = v / 3.6 * fc / 3e8; % 多普勒频移
h = sqrt(0.5) * (randn(size(t)) + 1i*randn(size(t)));
r = r .* exp(1i*2*pi*fd*t) .* h;

% 接收端采样
fs = 10*fc; % 采样频率
ts = 1/fs; % 采样时间间隔
r = r(1:fs/fc:end);

% 解调
s_hat = real(r) > 0;

% 比特误码率计算
ber = sum(x ~= s_hat) / length(x);
fprintf('比特误码率：%.4f\n', ber);

希望这次修改能够解决问题。