Question

我有带噪声的正弦信号，其中每次振荡的噪声点数在我的代码中应为每次振荡的5个点相同），我想每次改变振荡数：2、3、4….15 （在我的代码中更改了可变的“随机”）

在每个振荡次数下，我都将提取幅度作为频率的函数<不幸的是，由于很少振荡，FFT不能正常工作，所以信号中的点太少，因此我必须将信号（带噪声的正弦）拟合为正弦波，以便将新信号的频率与正弦波的频率进行比较

请参阅代码，如何进行拟合，以便提取信号的频率？

  %my code
  random=40;
  f=5; % the frequency of the sine wave also the number of points per 
     oscillation

  %the number of oscillation is random/f

  t = (1:random)';
  X = ones(random,2);


  y_1= sin((2*pi)/f*t);
  X(:,2) = y_1;
  y=y_1+randn(random,1);
  y = y(:);
  beta = X\y;
  yhat = beta(1)+beta(2)*sin((2*pi)/f*t);
 figure
 plot(t,y,'.b','markersize',12);
 hold on
 plot(t,yhat,'r','linewidth',2);

Answer 1

这是一个常见问题。请尝试下面的链接或Stackoverflow上的其他链接。就您而言，您有很多异常值，因此您可以使用RANSAC之类的方法将其排除。

https://www.mathworks.com/matlabcentral/answers/121579-curve-fitting-to-a-sinusoidal-function

https://www.mathworks.com/matlabcentral/answers/195371-sine-curve-fitting-for-the-given-data

在不使用FFT的情况下查找信号的频率

1 个答案: