| 查看: 431 | 回复: 0 | ||
ws89金虫 (著名写手)
悠然书生
|
[求助]
参数拟合程序，不知道是不是陷入了死循环，运行两天还在跑。求帮助
|
|
参数拟合程序，不知道是不是陷入了死循环，运行两天还在跑。师兄的程序，没有得到最后跟师兄一致的结果。求帮助。大神给看看是怎磨回事。 clc clear allT2=[-20 -10 -10 0 0 0 0 10 10 20]'; allf=1.0e+003 *[2.1304 0.9212 0.9212 0.4141 0.4141 0.4141 0.4141 0.193 0.193 0.093]'; allw2=[1 1 2 0.1 0.5 1 5 1 2 1]'; allG2=[23.46 8.062 14 1 2.97 4.18 13.483 2.13 3.128 1.08]'*10^6; allG1=[17 5.8 10 2.5 3.3 3.8 9.7 3.0 3.4 2.7]'*10^6; for i=1:10 f2(i,1)=allf(i); w2(i,1)=allw2(i); G1(i,1)=allG1(i); G2(i,1)=allG2(i); end x1=[f2 w2]; GG2=@(p,x)(p(1).*p(2).^2./p(4).^2+2.*p(1).*p(2)./p(4).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*p(5))+p(1).*x1(:,1).^(2.*p(5)).*x1(:,2).^(2.*p(5))+p(3).*p(2).^2./p(4)^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*cos... (pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*cos(pi/2.*p(5)).*cos(pi/2.*p(7))+p(3).*x1(:,1).^(2.*p(5)+p(7)).*x1(:,2).^(2.*p(5)+p(7)).*cos(pi/2.*p(7))+p(2).^2.... ./p(4).*cos(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*(p(5)+p(6)))); GG1=@(p,x)(p(3).*p(2).^2./p(4).^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*sin(pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*sin(pi/2.*p(5)).*sin(pi/2.*p(7))+p(3).*x1(:,1).^... (2.*p(5)+p(7)).*x1(:,2).^(2*p(5)+p(7)).*sin(pi/2.*p(7))+p(2).^2./p(4).*sin(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*sin(pi/2.*(p(5)+p(6)))); p0=[97.06669468363,272.111066749,60321.3626849,164.3298821730,0.517346991136036,0.20053385863749,0.26009044292186]; p=p0; GG22=p(1).*p(2).^2./p(4).^2+2.*p(1).*p(2)./p(4).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*p(5))+p(1).*x1(:,1).^(2.*p(5)).*x1(:,2).^(2.*p(5))+p(3).*p(2).^2./p(4)^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*cos... (pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*cos(pi/2.*p(5)).*cos(pi/2.*p(7))+p(3).*x1(:,1).^(2.*p(5)+p(7)).*x1(:,2).^(2.*p(5)+p(7)).*cos(pi/2.*p(7))+p(2).^2.... ./p(4).*cos(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*(p(5)+p(6))); %p0=p; GG11=p(3).*p(2).^2./p(4).^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*sin(pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*sin(pi/2.*p(5)).*sin(pi/2.*p(7))+p(3).*x1(:,1).^... (2.*p(5)+p(7)).*x1(:,2).^(2*p(5)+p(7)).*sin(pi/2.*p(7))+p(2).^2./p(4).*sin(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*sin(pi/2.*(p(5)+p(6))); a1=0.5;a2=1-a1; while(1) if (norm(GG11-allG1,inf)<1.0e6&&norm(GG22-allG2,inf)<1.0e6) break end options = optimset('MaxFunEvals',4000,'MaxIter',40000); p=fminsearch(@(p)(a1*(sum(((G1-GG1(p,x1))./G1).^2)+a2*sum(((G2-GG2(p,x1))./G2).^2))),p0); GG22=p(1).*p(2).^2./p(4).^2+2.*p(1).*p(2)./p(4).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*p(5))+p(1).*x1(:,1).^(2.*p(5)).*x1(:,2).^(2.*p(5))+p(3).*p(2).^2./p(4)^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*cos... (pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*cos(pi/2.*p(5)).*cos(pi/2.*p(7))+p(3).*x1(:,1).^(2.*p(5)+p(7)).*x1(:,2).^(2.*p(5)+p(7)).*cos(pi/2.*p(7))+p(2).^2.... ./p(4).*cos(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*cos(pi/2.*(p(5)+p(6))); GG11=p(3).*p(2).^2./p(4).^2.*x1(:,1).^p(7).*x1(:,2).^p(7).*sin(pi/2.*p(7))+2.*p(3).*p(2)./p(4).*x1(:,1).^(p(5)+p(7)).*x1(:,2).^(p(5)+p(7)).*sin(pi/2.*p(5)).*sin(pi/2.*p(7))+p(3).*x1(:,1).^... (2.*p(5)+p(7)).*x1(:,2).^(2*p(5)+p(7)).*sin(pi/2.*p(7))+p(2).^2./p(4).*sin(pi/2.*p(6))+p(2).*x1(:,1).^p(5).*x1(:,2).^p(5).*sin(pi/2.*(p(5)+p(6))); p0=p; end |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有15人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
20