24小时热门版块排行榜    

查看: 325  |  回复: 1

zhengmanzh

新虫 (小有名气)

[求助] 求高版本1stopt代做数据拟合

各位大神,跪求高版本做数据拟合。。。。。。感激涕零..............................................
parameters sigma0,u1,u2,sigma,u,p,fp_fmax;
parameterdomain=(0,1);
sharedmodel;
constant kc=20,gama=0.0712;
conststr pusai0=sigma0^2/pi*(arcsin(u1)+arcsin(u2));
conststr beita=2*pi*sigma0^2*1^2*(arcsin(u1)+arcsin(u2)-0.5*sin(2*arcsin(u1))-0.5*sin(2*arcsin(u2)))-4*pi*sigma0^2*1*(sqrt(1-u2^2)-sqrt(1-u1^2))^2/(arcsin(u1)+arcsin(u2));
conststr a=(2*pi*fp_fmax-(2*sigma0^2*1*(sqrt(1-u2^2)-sqrt(1-u1^2)))/pusai0)/sqrt(2*beita);
variable r,N_fmax,cdf;
function N_fmax=r*sqrt(2*beita)/pi^(3/2)/pusai0/1*
int((sqrt(1+gama/beita*(r*sigma/y)^2))/y/sqrt(2*pi)/sigma/y*exp(-1*(ln(y)-u)^2/2/sigma^2)*exp(-1*(p^2+r^2/y^2)/2/pusai0)*int(cosh(r*p*cos(sita)/y/pusai0)*(exp(-1*(a*p*sin(sita)/(sqrt(1+gama/beita*(r*sigma/y)^2)))^2)+sqrt(pi)*a*p*sin(sita)/(sqrt(1+gama/beita*(r*sigma/y)^2))*erf(a*p*sin(sita)/(sqrt(1+gama/beita*(r*sigma/y)^2)))),sita=0,pi/2),y=0,20);
cdf=int(int(x*exp(-1*(x^2+(p/sqrt(pusai0))^2)/2)*bessi0(p/sqrt(pusai0)*x),x=r/y/sqrt(pusai0),1000)*exp(-1*(ln(y)-u)^2/2/sigma^2)/sqrt(2*pi)/sigma/y,y=0,1000);
data;
0.0316    0.0700     0.9976
0.0355    0.0800     0.9947
0.0398    0.0920     0.9913
0.0447    0.1062     0.9873
0.0501    0.1207     0.9825
0.0562    0.1354     0.9768
0.0631    0.1580     0.9699
0.0708    0.1800     0.9615
0.0794    0.2000     0.9513
0.0891    0.2300     0.9389
0.1000    0.2450     0.9237
0.1122    0.2700     0.9052
0.1259    0.2920     0.8827
0.1413    0.3100     0.8555
0.1585    0.3278     0.8227
0.1778    0.3300     0.7835
0.1995    0.3240     0.7373
0.2239    0.3230     0.6835
0.2512    0.3130     0.6220
0.2818    0.2860     0.5551
0.3162    0.2500     0.4821
0.3548    0.2210     0.4050
0.3981    0.1920     0.3267
0.4467    0.1530     0.2512
0.5012    0.1140     0.1822
0.5623    0.0800     0.1233
0.6310    0.0600     0.0769
0.7079    0.0400     0.0436
0.7943    0.0230     0.0220
0.8913    0.0115     0.0097
1.0000    0.0040     0.0039
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengmanzh

新虫 (小有名气)

内容已删除
2楼2014-06-07 22:41:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhengmanzh 的主题更新
信息提示
请填处理意见