Regression Fit a SinWave to a dataset?

2 次查看(过去 30 天)
Ok..I'm almost there. Yesterday I received some help but still have questions (I'm an ultra beginner).
I love baseball and have a data set with the amount of runs scored by the Pittsburgh Pirates for their last 501 games. The data set is sampled by 1. I completed a spectral analysis on that data set and a 3.5678 game cycle is identified in the data set. I would like to know exactly how to: 1.) generate the sin/cos wave with a period of 3.5678 games 2.) fit the newly created sin/cos wave to the Score data set.
Any code or other insight will be MOST appreciated.
  1 个评论
Daniel Shub
Daniel Shub 2012-5-1
http://www.mathworks.com/matlabcentral/answers/6200-tutorial-how-to-ask-a-question-on-answers-and-get-a-fast-answer
http://www.mathworks.com/matlabcentral/answers/8626-how-do-i-get-help-on-homework-questions-on-matlab-answers

请先登录,再进行评论。

采纳的回答

Wayne King
Wayne King 2012-5-1
Something is probably off in your spectral analysis in identifying the frequency. You have 501 points, which means the spacing between frequencies is 1/501. The closest you can actually get to a period of 3.5678 is 0.2794 cycles/game. So let me use that
y = score(:); %make sure score is a column vector
X = ones(501,3);
n = 0:500;
X(:,2) = cos(2*pi*0.2794*n)';
X(:,3) = sin(2*pi*0.2794*n)';
betahat = X\y;
Fit is
yhat = betahat(1)+betahat(2)*cos(2*pi*0.2794*n)+betahat(3)*sin(2*pi*0.2794*n);
If you have the Statistics Toolbox
[b,bint,r,rint,stats] = regress(y,X);
returns useful information in stats for example:
the R-square statistic, the F statistic and p value for the full model, and an estimate of the error variance.
If you have the Statistics Toolbox and R2012a, you can use
mdl = LinearModel.fit(X,y)
which outputs a table. You can get the coefficients with:
mdl.Coefficients.Estimate
  2 个评论
Clifford Shelton
Clifford Shelton 2012-5-1
how did you get the .2794 number?
Is there a way to calculate that?
Sorry if my questions seem retarded...I'm so new to this!
Clifford Shelton
Clifford Shelton 2012-5-1
Also..for some reason..I cannot complete
%betahat = X/y;
without getting the error
"??? Error using ==> mldivide
Matrix dimensions must agree."
what am I doing wrong?

请先登录,再进行评论。

更多回答(2 个)

Wayne King
Wayne King 2012-5-1
Because when you do a spectral analysis, the frequency spacing in your case is 1/501. So you only get frequency estimates at 0, 1/501, 2/501,3/501, 4/501. So I was assuming that your spectral analysis really revealed a peak at the closest Nyquist frequency (a multiple of 1/501).
If you want, you can use 0.2803 which is 1/3.5678 for the frequency. That should not change the results much I would think.
  2 个评论
Dr. Seis
Dr. Seis 2012-5-1
He gave you the hint... it is a multiple of 1/501
140/501 = 0.2794...
He chose this because ABS(140/501 - 1/3.5678) < ABS(141/501 - 1/3.5678), which means 140/501 is closer to the frequency you are after.
Clifford Shelton
Clifford Shelton 2012-5-1
Honestly, I'm very new to this...so that "hint" didn't register with me.
I'm still having trouble with the error:
??? Error using ==> mldivide
Matrix dimensions must agree.
After looking at the values of my variables it seems that the X variable has a value of <501x3 double> while the y variable has a value of <501x1 double>. Could the two different matrix values be the problem?
When I copy and past your code it works. But when I try to use the data from my dataset i get the error.
Not sure why.

请先登录,再进行评论。


Wayne King
Wayne King 2012-5-1
You are not using what I gave you. You are using
%betahat = X/y;
Look at the example I gave you:
y = score(:); %make sure score is a column vector
X = ones(501,3);
n = 0:500;
X(:,2) = cos(2*pi*0.2794*n)';
X(:,3) = sin(2*pi*0.2794*n)';
betahat = X\y;
It's very important which slash you use.
  2 个评论
Clifford Shelton
Clifford Shelton 2012-5-1
Ok. I don't fully understand how you knew that the number needed to caclulatea cycle of 3.5678 games was .2794.
So my question would be what is the easiest way for me to caclulate what I would need to use for cycles identified by my spectral analysis as being:
31.19552 games
5.87064 games
8.00192 games
4.64491 games
Do I just use the same 1/501 calculation? That doesnt' seem tob e right. Please help! Thanks!!

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Strategy & Logic 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by