急求：含不可积积分的非线性拟合！！

lgq037

a=160.3;
b=4.15;
y=[0 15.15 43.94 63.63 65.15 103.02 101.51 113.63 131.05 144.68 172.73 181.82 193.18 196.97 203.03 209.09 218.18 224.24 225.75];
n=[1 0.93 0.84 0.81 0.80 0.74 0.73 0.71 0.69 0.68 0.66 0.65 0.64 0.64 0.64 0.63 0.62 0.62 0.61];


积分下限:1,上限:n。
想拟合参数:c和k
求高手帮忙解决,不胜感激!!!

liuyalong008

fit or nlinfit is not compatible with quad, so turning to primal numerical intergral formula maybe a soluble way

1. function f=nlinear_fit(beta,x)
   
2. a=160.3;
   
3. b=4.15;
   
4. f=zeros(1,length(x));
   
5. interval=zeros(length(x),100);
   
6. for j=1:length(x)
   
7.    interval(j,:)=linspace(1,x(j),100);
   
8.    h=interval(j,2)-interval(j,1);
   
9.     for i=1:100
   
10.       f(j)=f(j)+ interval(j,i).^(-beta(1)-1).*...
    
11.           exp((b-beta(1))./(beta(2)).*(1-interval(j,i).^beta(2)));  
    
12.     end
    
13.     f(j)=f(j)*(x(j)~=1);
    
14. end
    
15. f=-a.*f;   
    

复制代码

lgq037

非常感谢！！！！