在MATLAB中实现下面例题的解调程序： 信源：m(t)=cos(4*pi*t)+ sin(4*pi*t)，载波中心频率fc=25Hz,调频器的压控震荡系数为6hz/v，载波平均功率1.5w，画出该调频信号的波形及其调频信号的震荡谱

在MATLAB中实现这个调制解调问题，你需要首先了解调频的基本步骤，包括调制、信号生成以及谱分析。这里我会给出一个简单的流程：

1. **创建信源信号**:
使用`cos`和`sine`函数创建给定的信源信号 `m(t) = cos(4*pi*t) + sin(4*pi*t)`。假设时间范围从0到1秒。

t = linspace(0, 1, 1000); % 创建时间向量
m = cos(4*pi*t) + sin(4*pi*t);

2. **载波和调制**:
- 定义载波频率 `fc = 25 Hz` 和压控震荡系数 `Kc = 6 Hz/V`。
- 调频信号可以通过乘法将信源 `m(t)` 与载波 `exp(j*2*pi*ftc*t)` 结合，其中 `ftc` 是载波频率。

fc = 25;
Kc = 6; % 调制系数 (Hz/V)
voltage_per_hz = sqrt(1.5 * 10^(-3)); % 平均功率转化为电压单位 (V)
freq_deviation = Kc * voltage_per_hz; % 载波频率偏移
fm_signal = m .* exp(j * 2 * pi * fc * t + freq_deviation * t);

3. **绘制波形**:
使用`plot`函数显示原始信源和调频后的信号。

subplot(2, 1, 1)
plot(t, m, 'b', t, fm_signal, 'r')
legend('Original Signal', 'FM Signal')
xlabel('Time (s)')
ylabel('Amplitude')

4. **谱分析**:
   使用`fft`函数计算幅度谱或功率谱密度，`abs(fft(fm_signal))` 或 `abs(fft(fm_signal)).^2` 可以得到调频信号的谱。

% 获取幅度谱
spectrogram = abs(fft(fm_signal));
% 或者获取功率谱密度
power_spectrum = abs(fft(fm_signal)).^2;

subplot(2, 1, 2)
imagesc(abs(spectrogram), 'XData', t, 'YData', [0 fc + freq_deviation])
xlabel('Time (s)')
ylabel('Frequency (Hz)')
colorbar
title('Spectrum of FM Signal')

这只是一个基本示例，实际应用中可能需要处理采样率、窗口函数和更复杂的处理。运行上述代码前请确保已经安装了MATLAB环境并配置好图形界面。