My stochastic processes program

4 次查看(过去 30 天)
reem
reem 2011-5-30
Hello My teacher Matt
I have a program and When I asked my teacher if my solution is
true or not,he said to me that my solution is not correct
Please,could you help me
My program
N=400; %Define Number of samples((times))
f1=1; % Frequency of the cosinewave
FS=200; % Sampling Frequency
Q=0;
t=1:400;
n=0:N-1;
X1=0;
for i=1:400
for j=1:301
X=cos((2*pi*f1*(n/FS))+Q);
X1=X1+X;
Q= -2*pi/N;
end
Q=Q*2;
end
subplot(2,1,1); % Prepare the figure
plot(t,X1/(400*301)) % Plot mean value vs. t
grid;
title('mean value function of cosinewave ');
xlabel('Time index');
ylabel('Mean value');
Rxx=xcorr(X1); % Estimate its autocorrelation
subplot(2,1,2);
plot(Rxx); % Plot the autocorrelation
grid;
title('Autocorrelation function of the cosinewave');
xlabel('Time index');
ylabel('Autocorrelation');
My teacher at university said to me that I have mistake here and
he wrote this note but I did not understand what he mean!!!
X=cos((2*pi*f1*(n/FS))+Q); % Generate x(t)
Here X must be as X(i,j). Then it has 400 rows and 301 columns. Row number equals to the number of realizations
  1 个评论
Oleg Komarov
Oleg Komarov 2011-5-30
Can you please format the code properly: http://www.mathworks.com/matlabcentral/answers/7885-tutorial-how-to-format-your-question

请先登录,再进行评论。

回答(1 个)

Oleg Komarov
Oleg Komarov 2011-5-30
You basically calculate the same X all the time. Nothing changes inside the loops.
  4 个评论
reem
reem 2011-5-30
Yes,you are true
so I must apply that,isn't it???
Q1=Q1+Q
Oleg Komarov
Oleg Komarov 2011-5-30
Maybe, I don't know what is the stochastic process you have to reproduce.
Also, you probably have to change one of the input of X according to the indices i and j but the only vector you have in the expression is n (1 by 400) so you can use n(i), but this leaves the j index unused.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Vibration Analysis 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by