2014-02-27 148 views
0

我正在尝试使用匿名函数为我的数据创建一个洛伦兹分布。高斯分布创建时没有问题,但洛伦兹分布遇到此错误。错误使用/矩阵尺寸必须同意

错误:

Error using/

Matrix dimensions must agree.

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

代码:

% 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); 

回答

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

而1为1个维的向量,所以分割矩阵是无效的,./应适用。改变loren_fun_form的定义如下应该工作

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