QPSK的matlab

clear;
clc;

fc=1000;                      %载波频率
fs=10000;                     %采样频率
N=50;                         %信源个数
%********产生信息源**********
source=rand（150）>0.5;

for snr=1:10                 %信噪比
    %********串并转换************
    I1=source（1:2:N）;
    Q1=source（2:2:N）;

    %*****************采样**************
    sp=10; %一个数字波形采sp个点
    for i=1:length（I1）
        I2（（i-1）*sp+1:sp*i）=I1（i）;
        Q2（（i-1）*sp+1:sp*i）=Q1（i）;
    end

    %***************极性转化0->-11->1***********
    I3=2*I2-1;
    Q3=2*Q2-1;
    %*****************载波调制**********************
    t=0:1/fs:（length（I3）-1）/fs;
    s1=cos（2*pi*fc*t）;                          %载波形成
    s2=sin（2*pi*fc*t）;
    ss1=I3.*s1;                                 %I路调制
    ss2=Q3.*s2;                                 %Q路调制
    s=ss1+ss2;                                  %合成调制信号
    
    %***************计算频谱*******************
    nfft=2048;
    y=fft（snfft）;
    ps1=abs（y（1:nfft/2））.^2;
    p1=ps1/max（ps1）;
    f=fs*（0:（nfft/2-1））/nfft;


    %***************加高斯白噪声*************************
    s=awgn（ssnr）;
    %**********************************************

    %******************解调部分********************
    ss_d1=s.*s1;
    ss_d2=s.*s2;
    b1=fir1（102000/fs‘low‘）;        %低通滤波器设置
    sd1=filter（b11ss_d1）;        %滤波过程
    sd2=filter（b11ss_d2）;

    %********************抽样判决************
    sqpsk1=sd1（sp:sp:length（sd1））;
    for i=1:length（sqpsk1）
        if sqpsk1（i）>0
            sqpsk1（i）=1;