绘制一组正弦波

Question

我有以下代码：

Fs = 1000;
T = 1/Fs;
L = 1000;
t = (0:L-1)*T;
k = 25:1:50;
m = 1:1:25;

其中k 和m 是对应的。我想绘制由以下原因产生的 25 个正弦波：

x = m*sin(2*pi*k*t);

我曾想过使用 for 循环来执行此操作，该循环每次从 m 和 k 中获取一个值，但我不确定如何继续。

Answer 1

下面是一个非常基本的绘图解决方案。您会注意到很难看到图中发生了什么，因此您可能需要考虑其他方式来呈现这些数据。

function q45532082
Fs = 1000;
T = 1/Fs;
L = 1000;
t = (0:L-1)*T;
k = 26:1:50;
m = 1:1:25;
%% Plotting
assert(numel(m) == numel(k)); % We make sure that the number of elements is the same.
figure(); hold on; % "hold" is needed if you want to see all curves at the same time.
for ind1 = 1:numel(m)
  plot(t,m(ind1)*sin(2*pi*k(ind1)*t));
end

这是结果：

请注意，您的代码中k和m中的元素数量不同，所以我不得不更改它。

Answer 2

使用plot 的功能，您还可以在没有循环的情况下绘制所有正弦波：

Fs = 1000;
T = 1/Fs;
L = 1000;
t = (0:L-1)*T;
k = 26:1:50;
m = 1:1:25;
x = m.*sin(2.*pi.*bsxfun(@times,t.',k)); %this results in an L*25 matrix, each column is data of one wave
% or, if you have version 2016b or newer:
% x = m.*sin(2.*pi.*t.'*k);
plot(t,x) % plot all sines at ones

作为@Dev-iL noted，我还必须更改k。

L = 1000 的结果太拥挤，所以我在这里用L = 50 绘制它：