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;