이전 질문 게시판은 새 글 쓰기를 막았습니다. [질문 게시판]을 이용바랍니다.
Date |
2009/05/27 21:39:15 |
Name |
시봉창 |
Subject |
매틀랩 코드 해석좀 부탁드립니다. |
Generate 100,000 observations of a random variable Y uniformly distributed between 0 and 1 using the Matlab function "rand( )". By comparing the estimated pdf and the real pdf, discuss if the observations were successfully generated.
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
N = 100000;
no_inv = 100;
st_v = 0; en_v = 1;
inv = (st_v+(en_v-st_v)/2/no_inv):(en_v-st_v)/no_inv:(en_v-(en_v-st_v)/2/no_inv);
rand('state',0);
y = rand(N,1);
h_y = hist(y, inv);
n_h_y = h_y/((en_v-st_v)*(N/no_inv));
bar(inv, n_h_y, 'y')
hold on, plot(inv, ones(no_inv), 'r')
------------------------------------------------------------
100,000개의 random variable을 생성하고 rand 함수를 이용하여 0과 1사이에 일정하게 분배(?)하여 얻은 pdf와 실제의 pdf를 비교하라. 뭐 이런 문제 같으데요..
st_v와 en_v는 0과 1사이에 분배하라고 했으니 그런 것 같구요.
그런데 inv에서 왜 행렬을 생성해주었는지 잘 모르겠습니다.
y=rand(N,1)은 100,000개의 randomvariable을 생성하는 과정이고.
h_y는 이렇게 생성된 randomvariable을 히스토그램의 bin에 넣는 과정인 것 같습니다.
n_h_y의 존재도 잘 모르겠구요.
|
통합규정 1.3 이용안내 인용
"Pgr은 '명문화된 삭제규정'이 반드시 필요하지 않은 분을 환영합니다.
법 없이도 사는 사람, 남에게 상처를 주지 않으면서 같이 이야기 나눌 수 있는 분이면 좋겠습니다."
|