信息处理 - 调频米(吨)m(t)在 MATLAB 中 - 吾爱随笔录

我在以下区间

m (t) = {\begin{cases} 1 & , 0 \leq t < \frac{t_{0}}{3} \\ - 2 & , \frac{t_{0}}{3} \leq t < \frac{2 t_{0}}{3} \\ 0 & , \frac{2 t_{0}}{3} \leq t < t_{0} \end{cases}

$m(t)=\begin{cases}1&, 0\leq t < \frac{t_0}{3}\\[2ex] -2&, \frac{t_0}{3}\leq t <\frac{2t_0}{3}\\[2ex] 0&, \frac{2t_0}{3}\leq t < t_0\end{cases}$

其中，采样频率。 $t_0 = 0.15\ \rm seconds$ $F_s=2\ \rm kHz$

我的 MATLAB 代码以构造 $m(t)$

t0= 0.15;
Fs= 2000;
t1= 0:1/Fs:t0/3- 1/Fs;
t2= t0/3:1/Fs:2*t0/3 - 1/Fs;
t3= 2*t0/3:1/Fs:3*t0/3-1/Fs;
m1= ones([1, length(t1)]);
m2= -2*ones([1, length(t2)]);
m3= zeros([1, length(t3)]);
m= [m1, m2, m3];
t= [t1, t2, t3];

我试图找到频率调制，并希望使用以下公式载波频率 Hz 和偏差常数。

x_{F M} (t) = \cos (2 π f_{c} t + 2 π k_{f} \int_{- \infty}^{t} m (x) d x)

$x_{FM}(t) = \cos\left(2 \pi f_c t + 2\pi k_f \int_{-\infty}^t m(x) dx\right)$

f_{c} = 200

$f_c=200$

k_{f} = 50

$k_f= 50$

当我们找到整数时，只允许使用 for 或 while 循环。为了整合，我做了这个for循环

result=0;
j=1;

for i= 0:1/Fs:t0- 1/Fs
   result= result + m(j);

    %% result2 will be used as we integrate it over 0 to 0.15
   result2(j) = result;
   j= j+1;
end

请告诉我，这种集成有问题吗？我绘制了这个，这对我来说似乎是正确的。

如果积分是真的，那么为什么在绘制这个信号时频率根本没有改变。

xFM= cos(2*pi*fc*t + 2*pi*kf*result)

完整代码在这里