用matlab仿真瑞利信道下qpsk的误码率，不能用matlab自带的瑞利函数，瑞利信道h = sqrt(0.5)*( randn(1,1)+ sqrt(-1)*randn(1,1) )，matlab是2016a版本

在MATLAB 2016a版本中，你可以通过以下步骤模拟QPSK在瑞利信道下的误码率：

首先，你需要创建一个伪随机数据序列作为QPSK信号，然后添加瑞利噪声到这个序列上。以下是基本步骤：

1. **生成QPSK信号**:

   % 设置QPSK参数
   M = 4;           % QPSK的点数 (4 for QPSK)
   data = randi([0 M-1], N, 1); % 创建一个包含0到M-1整数的数据序列

   % 将二进制数据转换为QPSK符号
   qam_symbols = pskmod(data, M, 'qpsk');

2. **创建瑞利信道模型**:
照你的描述，我们可以直接生成标准正态分布的加复共轭噪声，并乘以sqrt(0.5)来得到瑞利信道影响：

   snr_db = 10;      % 设定信噪比，例如10dB
   noise_power = 10^(-snr_db/10);  % 计算噪声功率

   h = sqrt(0.5) * [randn(1, N) + j*randn(1, N)];  % 瑞利噪声
   noisy_qam = qam_symbols .* h;                 % 应用噪声

3. **解调并计算误码率**:
使用QPSK解调器对噪声信号进行解调，通常可以使用`demodulate`函数。然后计算误码率，比如使用伯努利试验或累积错误计数：

   decoded_data = demodulate(noisy_qam, 'qpsk');  % 解调符号
   ber = sum(decoded_data ~= data) / N;       % 比较原数据和解调后的数据，计算误码率

4. **重复多次并求平均**:
如果想获得更准确的结果，可以增加迭代次数，每次解调后累加误码并除以总次数，最后取平均值。

ber_avg = 0;
for i = 1:num_trials
    % ... (以上步骤)
    ber_avg = ber_avg + ber;
end
ber_avg = ber_avg / num_trials;  % 平均误码率