- 积分
- 23
- 注册时间
- 2010-11-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2012-2-18 22:29:23
|
显示全部楼层
来自 山东青岛
本帖最后由 liuyalong008 于 2012-2-18 22:30 编辑
fit or nlinfit is not compatible with quad, so turning to primal numerical intergral formula maybe a soluble way- function f=nlinear_fit(beta,x)
- a=160.3;
- b=4.15;
- f=zeros(1,length(x));
- interval=zeros(length(x),100);
- for j=1:length(x)
- interval(j,:)=linspace(1,x(j),100);
- h=interval(j,2)-interval(j,1);
- for i=1:100
- f(j)=f(j)+ interval(j,i).^(-beta(1)-1).*...
- exp((b-beta(1))./(beta(2)).*(1-interval(j,i).^beta(2)));
- end
- f(j)=f(j)*(x(j)~=1);
- end
- f=-a.*f;
复制代码
|
评分
-
1
查看全部评分
-
|