/Matrix 차원을 사용하는 중에 오류가 발생했습니다.

Question

익명의 기능을 사용하여 데이터 용 로렌츠 배포판을 만들려고합니다. Gaussian 분포는 아무런 문제없이 생성되지만 Lorentzian 분포는이 오류를 만난다.

오류 :

Error using/

Matrix dimensions must agree.

Error in @(p,x)1/(pi*p(2)*(1+((x-p(1))/(p(2))).^2))

코드 : (1) 1 차원의 벡터이므로 행렬 분할이 유효하지

% Variables necessary for histograms 
num_bins = 100; 
low_range = linspace(min(low_data), max(low_data), num_bins); 
low_counts = histc(low_data, low_range); 
low_err = sqrt(low_counts); 


% Creating and testing fits 
gauss_fun_form = @(p, x) 1/(p(2)*sqrt(2*pi))*exp(-(x-p(1)).^2/(2*p(2).^2)); 
loren_fun_form = @(p, x) 1/(pi*p(2)*(1+((x-p(1))/(p(2))).^2)); 

low_avg = mean(low_data); 
low_std = std(low_data);  

low_gauss_param = [low_avg low_std]; 
low_gauss = gauss_fun_form(low_gauss_param, low_range); 
low_gauss_scale = max(low_counts)/max(low_gauss); 

low_loren_param = [low_avg 2*log(2)*low_std]; 
low_loren = loren_fun_form(low_loren_param, low_range); 
low_loren_scale = max(low_counts)/max(low_loren); 

Answer 1

(pi*p(2)*(1+((x-p(1))/(p(2))).^2)) generate vector of 100 dimensions 

동안 ./ 적용되어야한다. 다음과 같이 loren_fun_form의 정의를 변경하십시오.

loren_fun_form = @(p, x) 1./(pi*p(2)*(1+((x-p(1))/(p(2))).^2));