0

我试图通过以下模型拟合以下正弦模型，这意味着确定性模型+白噪声：

function [x]=generate(N,m,A3)
f1=100;
f2=200;
T=1/f1;
for t=0:N*T/m:N*T
    WN(t)=rand(1,t)*2-1;
    x(t)=20*sin(2*pi*f1*t)+30*cos(2*pi*f2*t)+A3*WN(t);
end;

但是当我使用以下命令运行时：

generate(100,10,40)

我收到以下错误：

Subscript indices must either be real positive integers or logicals.

Error in generate (line 6)
    WN(t)=rand(1,t)*2-1;

那么有什么问题呢？据我所知，我可以为rand函数使用两个参数，对吗？这里有什么问题？

更新：这就是我为修复它所做的

function [ x]=generate(N,m,A3)
WN=[];
f1=100;
f2=200;
T=1/f1;
k=0;
for t=0:N*T/m:N*T
   % WN(k)=rand(1,1)*2-1;
   % x(k)=20*sin(2*pi*f1*t)+30*cos(2*pi*f2*t)+A3*WN(k);
    %k=k+1;
    WN=rand*2-1;
     disp(20*sin(2*pi*f1*t)+30*cos(2*pi*f2*t)+A3*WN);
end;

和结果

 generate(1000,10,40)
   25.4143

   -1.4678

   66.9518

   -9.6293

   51.9928

   55.3843

   59.4956

   -3.2451

   21.9826

   10.7896

   54.0055
4

3 回答 3

1

rand(1,t) 生成 1 by t 矩阵，要获得标量值，您可以使用 rand(1,1)

于 2013-03-21T10:52:23.223 回答
1

请注意， sin 和 cos 接受向量参数，因此不需要 for 循环，

T = 1 / fs;
t = (0:N-1)' * T; 
WN = rand(N, 1) * 2; 
WN = WN - mean(WN); 
x = 20 * sin(2*pi*f1*t) + 30 * cos(2*pi*f2*t) + A3 * WN;

此外，您的 T 值应该是 1 / 采样率，而不是 1 / f1（除非您在此频率下采样，在这种情况下，您不能期望看到正确的结果，因为您的频率都高于 fs / 2 的奈奎斯特频率，所以会有混叠）。此外，您应该减去整个向量 WN 的平均值，而不是盲目地假设平均值为 0.5（或 1，一旦乘以 2）。

于 2013-03-21T12:05:14.603 回答
1

我已解决如下：

function x = generate(N,m,A3)
f1 = 100;
f2 = 200;
T = 1./f1;
t = (0:(N*T/m):(N*T))';
wn = rand(length(t),1).*2 - 1;
x = 20.*sin(2.*pi.*f1.*t) + 30.*cos(2.*pi.*f2.*t) + A3.*wn;
end
于 2013-03-21T12:27:41.807 回答