how to calculate the probability of an event with monte carlo simulation

조회 수: 11 (최근 30일)
I am trying to calculate the probability that the wind speed exp (Yt) will be below 15 kmh for more than 10 days in a 100 day period but when I run the code it tells me this "Unable to perform assignment because the left and right sides have a different number of elements", Can anyone help me to solve this problem? Thanks in advance.
function mycode()
n=1000;
x=zeros(n);
y=zeros(100);
y(1)=1.5 ;
y(2)=1.5 ;
for i= 1:n
for t= 3:100
y(t)=0.6*y(t-1)+0.4*y(t-2)+normrnd(0,0.1,[1,1]);
end
x(i)=sum(y<1.427); %1.427 is the wind speed
mean(x>10)
return
end
  댓글 수: 1
Rik
Rik 2020년 12월 11일
This time I edited your question for you. Next time, please use the tools explained on this page (and apply the standard indentation) to make your question more readable.

댓글을 달려면 로그인하십시오.

채택된 답변

Sibi
Sibi 2020년 12월 11일
mycode()
function mycode()
n=1000;
x=zeros(n,1);
y(1)=1.5 ;
y(2)=1.5 ;
for i= 1:n
for t= 3:100
y(t)=0.6*y(t-1)+0.4*y(t-2)+normrnd(0,0.1,[1,1]);
end
x(i)=sum(y<1.427);
end
mean(x>10)
return
end

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Energy Production에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by