找一篇matlab用USRP进行调制解调的程序

以下是一个matlab用USRP进行调制解调的程序的示例：

调制部分：

% 设置参数
Fs = 1e6; % 采样率
fc = 2e6; % 载波频率
modulation_index = 1; % 调制指数
duration = 1; % 信号持续时间
waveform = 'fm'; % 调制类型

% 生成调制信号
t = linspace(0, duration, duration*Fs);
if strcmp(waveform, 'fm')
    m = sin(2*pi*200*t); % 调制信号
    modulated = fmmod(m, fc, Fs, modulation_index); % 调制
else
    m = randn(1, length(t)); % 调制信号
    modulated = ammod(m, fc, Fs, modulation_index); % 调制
end

% 连接到USRP
rx = comm.SDRuReceiver;
rx.Platform = 'B200mini';
rx.MasterClockRate = Fs*10;
rx.DecimationFactor = 10;
rx.CenterFrequency = fc;
rx.Gain = 20;
rx.SamplesPerFrame = length(modulated);
rx.OutputDataType = 'double';

% 接收数据
received = rx();

% 解调
if strcmp(waveform, 'fm')
    demodulated = fmdemod(received, fc, Fs, modulation_index); % 解调
else
    demodulated = amdemod(received, fc, Fs, modulation_index); % 解调
end

% 绘制结果
subplot(3, 1, 1);
plot(t, m);
title('调制信号');
subplot(3, 1, 2);
plot(t, modulated);
title('调制后信号');
subplot(3, 1, 3);
plot(t, demodulated);
title('解调后信号');

解调部分：

% 设置参数
Fs = 1e6; % 采样率
fc = 2e6; % 载波频率
modulation_index = 1; % 调制指数
duration = 1; % 信号持续时间
waveform = 'fm'; % 调制类型

% 连接到USRP
tx = comm.SDRuTransmitter;
tx.Platform = 'B200mini';
tx.MasterClockRate = Fs*10;
tx.InterpolationFactor = 10;
tx.CenterFrequency = fc;
tx.Gain = 20;
tx.SamplesPerFrame = duration*Fs;
tx.InputDataType = 'double';

% 生成调制信号
t = linspace(0, duration, duration*Fs);
if strcmp(waveform, 'fm')
    m = sin(2*pi*200*t); % 调制信号
    modulated = fmmod(m, fc, Fs, modulation_index); % 调制
else
    m = randn(1, length(t)); % 调制信号
    modulated = ammod(m, fc, Fs, modulation_index); % 调制
end

% 发送数据
tx(modulated);

% 接收数据
rx = comm.SDRuReceiver;
rx.Platform = 'B200mini';
rx.MasterClockRate = Fs*10;
rx.DecimationFactor = 10;
rx.CenterFrequency = fc;
rx.Gain = 20;
rx.SamplesPerFrame = length(modulated);
rx.OutputDataType = 'double';
received = rx();

% 解调
if strcmp(waveform, 'fm')
    demodulated = fmdemod(received, fc, Fs, modulation_index); % 解调
else
    demodulated = amdemod(received, fc, Fs, modulation_index); % 解调
end

% 绘制结果
subplot(3, 1, 1);
plot(t, m);
title('调制信号');
subplot(3, 1, 2);
plot(t, modulated);
title('调制后信号');
subplot(3, 1, 3);
plot(t, demodulated);
title('解调后信号');