用MATLAB写PM调制解调有关仿真

在MATLAB中进行PM调制解调仿真，需要以下步骤：

1. 生成模拟信号：可以使用MATLAB中的函数生成一个正弦波或者其他信号作为调制信号。

2. 进行PM调制：使用MATLAB中的函数将调制信号进行PM调制。

3. 添加噪声：为了更加真实地模拟信道的影响，可以使用MATLAB中的函数向调制信号添加高斯白噪声或其他噪声。

4. 进行PM解调：使用MATLAB中的函数将接收到的PM调制信号进行解调，并还原原始调制信号。

5. 分析仿真结果：可以使用MATLAB中的函数分析仿真结果，例如计算误码率、绘制时域图和频谱图等。

下面是一个简单的PM调制解调仿真代码示例：

% 生成调制信号
fm = 10; % 调制信号频率
t = 0:0.001:1; % 时间序列
x = sin(2*pi*fm*t); % 调制信号

% PM调制
fc = 100; % 载波频率
kf = 2*pi*10; % 调制指数
y = pmmod(x,fc,kf,fs,'pem');

% 添加高斯白噪声
noise_power = 0.01;
y_noisy = awgn(y, noise_power, 'measured');

% PM解调
x_demod = pmdemod(y_noisy,fc,kf,fs,'pem');

% 绘制时域图和频谱图
figure;
subplot(2,1,1);
plot(t,x,'b',t,x_demod,'r');
legend('原始信号','解调信号');
xlabel('时间');
ylabel('幅值');
title('时域图');

subplot(2,1,2);
N = length(x_demod);
f = (-N/2:N/2-1)*(fs/N);
X = fftshift(fft(x_demod))/N;
plot(f,abs(X));
xlabel('频率');
ylabel('幅值');
title('频谱图');

这段代码生成一个频率为10Hz的正弦波作为调制信号，使用PM调制将其调制到100Hz的载波上，再添加高斯白噪声，最后进行PM解调并绘制时域图和频谱图。你可以根据需要修改参数和信号类型，进行更加复杂的仿真。