matlab数字化调制代码

以下是Matlab中的数字化调制代码示例：

1. 二进制相移键控（BPSK）调制

% 设置参数
N = 1000; % 信号长度
fc = 100; % 载波频率
Rb = 10; % 比特率
Tb = 1/Rb; % 每个比特的持续时间
fs = 10*Rb; % 采样频率
Ts = 1/fs; % 采样时间间隔
t = (0:N-1)*Ts; % 时间轴

% 生成随机比特序列
data = randi([0,1],1,N);

% BPSK调制
s = sqrt(2/Tb)*cos(2*pi*fc*t + pi*(data-0.5));

% 绘制时域波形
subplot(2,1,1)
plot(t,s)
xlabel('时间')
ylabel('幅度')
title('BPSK信号时域波形')

% 绘制频谱图
S = fftshift(fft(s))/N;
f = (-N/2:N/2-1)*(fs/N);
subplot(2,1,2)
plot(f,abs(S))
xlabel('频率')
ylabel('幅度')
title('BPSK信号频谱图')

2. 四相移键控（QPSK）调制

% 设置参数
N = 1000; % 信号长度
fc = 100; % 载波频率
Rb = 10; % 比特率
Tb = 1/Rb; % 每个比特的持续时间
fs = 10*Rb; % 采样频率
Ts = 1/fs; % 采样时间间隔
t = (0:N-1)*Ts; % 时间轴

% 生成随机比特序列
data = randi([0,1],1,2*N);

% 将二进制比特序列转换为QPSK符号序列
symbols = 2*data(1:2:end)-1 + 1i*(2*data(2:2:end)-1);

% QPSK调制
s = sqrt(2/Tb)*real(symbols).*cos(2*pi*fc*t) - sqrt(2/Tb)*imag(symbols).*sin(2*pi*fc*t);

% 绘制时域波形
subplot(2,1,1)
plot(t,s)
xlabel('时间')
ylabel('幅度')
title('QPSK信号时域波形')

% 绘制频谱图
S = fftshift(fft(s))/N;
f = (-N/2:N/2-1)*(fs/N);
subplot(2,1,2)
plot(f,abs(S))
xlabel('频率')
ylabel('幅度')
title('QPSK信号频谱图')

以上是二进制相移键控（BPSK）调制和四相移键控（QPSK）调制的Matlab示例代码，可以根据自己的实际应用进行修改。