用matlab写一个msk调制

可以使用以下代码实现MSK调制：

% 设置参数
fc = 1000; % 载波频率
fs = 10000; % 采样频率
T = 1/fs; % 采样时间间隔
N = 1000; % 采样点数
fc1 = 2000; % 频移1
fc2 = 4000; % 频移2
A = 1; % 幅度

% 生成调制信号
t = (0:N-1)*T; % 时间序列
s1 = A*cos(2*pi*fc*t); % 原始信号
s2 = A*cos(2*pi*(fc+fc1)*t + pi/2); % 频移1信号
s3 = A*cos(2*pi*(fc+fc2)*t + pi); % 频移2信号
s = [s1 s2 s3]; % 合并信号

% 绘制调制信号
figure;
plot(t, s1, 'b', t, s2, 'r', t, s3, 'g');
xlabel('时间');
ylabel('幅度');
title('MSK调制信号');

% 生成载波信号
t_c = (0:N*2-1)*T; % 时间序列
c1 = cos(2*pi*fc*t_c); % 载波1
c2 = sin(2*pi*fc*t_c); % 载波2
c = [c1 c2]; % 合并载波

% 进行调制
m = length(s); % 调制信号长度
n = length(c); % 载波长度
y = zeros(1, n); % 初始化调制信号
for i = 1:m
if s(i) > 0
y((i-1)*2+1:i*2) = c((i-1)*2+1:i*2); % 频移1
else
y((i-1)*2+1:i*2) = c((i-1)*2+2:-1:(i-1)*2+1); % 频移2
end
end

% 绘制调制后的信号
figure;
plot(t_c, y);
xlabel('时间');
ylabel('幅度');
title('MSK调制后的信号');