关于BPSK、QPSK、MSK、QAM的调制与解调

%%  bpsk信号的产生
T = 1;%采样时间
f = 1e5;%采样率  fs=100000
t = 0:1/f:（T-1/f）; % 采样点
n = length（t）; % 采样点数=100000

binary_code = ‘1001100110‘;
len_code = length（binary_code）;  %得出二进制数的长度

% 一个码元对应的脉冲长度   采样点的个数和信号的长度一致
f_bpsk = 1000;     % 载波频率的频率为1000Hz，周期1ms
T_bpsk = 1/f_bpsk; % 载波频率的周期
A_bpsk = 10;       % 载波频率的幅值

t_mayuan = 0:1/f:（T_bpsk-1/f）;               % 一个码元所持续的时间内的采样点时刻
y_mayuan = A_bpsk*sin（2*pi*f_bpsk*t_mayuan）; % 一个码元的波形，用于表示0
t_bpsk = 0:1/f:（T_bpsk*len_code-1/f）;        %bpsk所持续的时间内的采样点时刻

y_bpsk = [];       %通过下面循环将其拼接在一起
for i=1:len_code
    tmp = str2num（binary_code（i））;  %str2num 表示 Convert string to number
    if tmp == 0
        y_bpsk = [y_bpsky_mayuan];
    end
   if tmp == 1
        y_bpsk = [y_bpsk-y_mayuan];
    end
end
% subplot（221）;plot（t_bpsky_bpsk‘green‘）;  %绘制对应的bpsk信号
% xlabel（‘t/s‘）;title（‘‘1001100110’的BPSK编码信号‘）;
% axis（[0 1e-2 -2*A_bpsk 2*A_bpsk]）;

%% bpsk信号叠加噪声
n=randn（size（y_bpsk））;  %产生随机噪声
% subplot（331）;plot（size（y_bpsk）randn）;%绘制噪声信号
y_bpsk=y_bpsk+n;%受燥后的信号
subplot（221）;plot（t_bpsky_bpsk‘green‘）;  %绘制对应的bpsk信号
xlabel（‘t/s‘）;title（‘‘1001100110’的BPSK编码信号（受rand）‘）;
axis（[0 1e-2 -2*A_bpsk 2*A_bpsk]）;

%% 对BPSK信号做FFT变换
Y=fft（y_bpsk）;%对信号做傅里叶变换

%%  幅度谱
%信号采集的频率 Fs = 1000 Hz，采集的数据点个数 N = 1000    信号长度N=1000和fs=1000采样频率要一一对应
L=length（t_bpsk）;   %disp（L）;  %L=1000
P1=2*abs（Y）/L;%幅度谱   要得到真实的振幅值的大小，只要将得到的变换后结果乘以2除以L即可。
P=P1（1:L）;%和下一行语句采样点的个数对应
% disp（P）;
f=（0:L-1）;
subplot（222）;plot（fP）;
xlabel（‘f/Hz‘）;title（‘BPSK信号的幅度谱‘）;

%% 相位图
% L=1000;

Pyy=[1:L];
for i=1:L
    Pyy（i）=phase（Y（i））; %计算相位
    Pyy（i）=Pyy（i）*180/pi; %换算为角度
end;
F=（[1:L]-1）*f_bpsk/L; %换算成实际的频率值
subplot（223）;plot（F（1:L）Pyy（1:L））;   %显示相位图
xlabel（‘f/Hz‘）;title（‘相位-频率曲线图‘）;

%%  载波频率
Fs = 2500;%信号采集的频率，Fs>2*f_bpsk采样定理的要求
T=1/Fs;
L = 2500;%采集的数据点个数    L和Fs必须相等  
t=（0:L-1）*T;
f_bpsk = 1000; % 载波频率1000Hz，周期1ms
y_mayuan = 10*sin（2*pi*f_bpsk*t）;
Y=fft（y_mayuan）;%对信号做傅里叶变换
P1=2*abs（Y）/L;%幅度谱
P2=P1（1:L/2）;%和下一行语句采样点的个数对应
f=（0:L/2-1）;
subplot（224）;plot（fP2）;
xlabel（‘f/Hz‘）;title（‘高频载波的幅度谱‘）;
fprintf（‘观察点1
‘）;

%% 求3dB带宽 
%找到频谱的最大值
disp（P）;
max=0;  %将max的值初始化
N=length（P）;  %  求出P的长度为1000
% fprintf（‘观察点2
‘）;
%P（11）  检测P（1，i）的值
for i=1:N
    if P（1i）>=max
        max=P（1i）;
    else 
        max = max;
    end
end
fprintf（‘幅频的最大值为 %g
‘max）;
%求出最大值对应的频率
global fh;
for i=1:N
    if P（1i）==max
        fh=i;
    end
end
fprintf（‘最大值对应的频率为 %gHz
‘fh）;

%求出下降3dB对应的幅度大小
global m;
for i=1:N
    if P（1i）>=0.7*max &  P（1i）<=0.8*max
        m=P（1i）;
    end
end
fprintf（‘幅频的0.707*最大值为 %g
‘m）;
%求出下降3dB对应的幅度所对应的频率
global fl;
for i=1:N
    if P（1i）==m
        fl=i;
    end
end
fprintf（‘下降3dB对应的幅度所对应的频率为 %gHz
‘fl）;

%求出3dB带宽
global B;
B=fh-fl;
fprintf（‘3dB带宽为 %gHz
‘B）;

%% 码元速率的计算
Rb=B/2;
fprintf（‘码元速率为 %g波特
‘Rb）;

%% 信号的功率     计算公式为 pow=sum（xn.^2）/（length（x）/Fs），xn为信号序列，Fs为采样率
% disp（y_bpsk）;  %显示y_bpsk序列的各点的值
%disp（leng

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-05-25 10:26  PSK调制
     目录           0  2020-05-26 18:07  PSK调制BPSK
     文件         551  2020-05-25 10:05  PSK调制BPSKDTDT_bianhuan.m
     文件        3761  2020-05-25 10:19  PSK调制BPSKpsk_rand.m
     文件        5181  2020-05-25 09:42  PSK调制BPSKpsk_tiaozhi_jietiao.m
     文件        3930  2020-05-25 10:16  PSK调制BPSKcompute_Rb.m
     文件        3129  2020-05-25 10:22  PSK调制BPSKcompute_Rb_way2.m
     文件        2994  2020-05-26 17:17  PSK调制BPSKerjingzhi_tiaozhi.m
     文件         861  2020-05-25 09:28  PSK调制BPSKfft_FFT.m
     文件         644  2020-05-14 09:10  PSK调制BPSKgenerate_bpsk.m
     文件        2955  2020-05-15 09:39  PSK调制BPSKsample_5.m
     文件        1028  2020-05-25 10:00  PSK调制BPSK est_2.m
     文件        3697  2020-05-18 13:27  PSK调制BPSK源代码1.txt
     文件        3126  2020-05-18 08:33  PSK调制BPSK码元速率的估计.txt
     目录           0  2020-05-27 16:12  PSK调制MQAM
     文件        1466  2020-05-27 16:10  PSK调制MQAMMQAM_error.m
     文件        3359  2020-05-27 15:56  PSK调制MQAMQAM16.m
     文件        1775  2020-05-27 16:00  PSK调制MQAMQAM_16_tiaozhi_jietiao.m
     文件         230  2020-05-27 16:04  PSK调制MQAMsample_1.m
     目录           0  2020-06-09 10:10  PSK调制MSK
     文件       14181  2020-06-07 07:20  PSK调制MSK1.docx
     文件        3490  2020-06-10 22:33  PSK调制MSKMSK_MSK_1.m
     文件        2458  2020-06-09 10:06  PSK调制MSKMSK_MSK_2.m
     文件        2548  2020-06-01 12:49  PSK调制MSKMSK_signal.m
     文件        2930  2020-06-09 09:43  PSK调制MSKMSK_signal_1.m
     文件        3156  2020-06-05 11:03  PSK调制MSKMSK_signal_2.m
     文件        9391  2020-06-02 16:13  PSK调制MSKUntitled2.m
     文件        2861  2020-05-30 15:10  PSK调制MSKmsk.m
     文件        1679  2020-05-31 13:23  PSK调制MSKmsk_pingpu.m
     文件        2434  2020-06-01 10:20  PSK调制MSKmsk_pingpu_1.m
     文件         111  2020-06-01 19:47  PSK调制MSK est.m
............此处省略17个文件信息