对于初学者,您可以像这样对方程进行矢量化:
t = 0:0.001:10;
input = cos(t);
xam = (Ac + input).*cos(2*pi*fc*(1:length(t)));
我怀疑您每次看到完全相同的输出曲线的原因是您的fc
值是整数。如果fc
是一个整数,cos(2*pi*fc*t)
对于 的所有整数值计算为 1,t
调制后输入信号保持不变。当您使用时,1/fc
您会得到一个非整数值并且输出被调制。
我相信你想要做的是以下几点:
t = 0:0.001:10; % Time in seconds
f = 1; % Frequency of input signal in rad/s
fc = 3; % Frequency of carrier wave in rad/s
input = cos(2*pi*f*t); % Input signal
xam = (Ac + input).*cos(2*pi*fc*t); % Amplitude modulated signal
注释显示了每个参数的含义。如果您只想将采样率传递给fs
您的函数,则解决方案如下所示:
N = 1001;
fs = 1000;
t = (0:1:N)*(1/fs);
f = 1;
fc = 3;
input = cos(2*pi*f*t);
你的函数看起来像这样:
function xam = modulate(input, fs, fc, Ac)
t = (0:1:length(input))*(1/fs);
xam = (Ac + input).*cos(2*pi*fc*t);