再次麻烦了
我的问题是: f=-v/2*sqrt(pi)*quad(t^(z^2-1/2)*(1+2*z^2*log(t))*exp(-v^2*z^2*t^(2*z^2)),t,0,1)=0;这是我的等式方程。其中,v 是自己给定的。v=10:0.5:32.5;要求解z的最小值。试用FMINBND函数来解决。并返回上述方程 f的最小值。我想z也要自定义一个区间来解决此问题。实验数据如下:
f
1.6082e-3
3.5517e-4
7.8242e-5
1.7214e-5
3.7843e-6
8.3157e-7
1.8268e-7
4.0125e-8
V
10
12.5
15
17.5
20
22.5
25
27.5
Z*
0.2854
0.2561
0.2341
0.2169
0.2030
0.1914
0.1816
0.1731
第二个问题:min v=fv(w,z)
Subject g(v,z)=0
其中:g(v,z)=w-B=fw(v,z)-B;B为常数。erfc(y)=1-
w= fw(v,z)=1/2quad(erfc(v*z*exp(-2*x*z^2))*exp(-x),x,0,inf)
diff(w,z)=-v/sqrt(pi)*quad(exp(-v^2*z^2*exp(-4x*z^2)-x-2*x*z^2)*(1-4*x*z^2),x,0,1/4*z^2)- v/sqrt(pi);quad(exp(-v^2*z^2*exp(-48x*z^2)-x-2*z^2)*exp(1-4*x*z^2),x, 1/4*z^2,inf);
我们的目的是要求解:
g(v,z)=0
diff(w,z)=-v/sqrt(pi)*quad(exp(-v^2*z^2*exp(-4x*z^2)-x-2*x*z^2)*(1-4*x*z^2),x,0,1/4*z^2)- v/sqrt(pi);quad(exp(-v^2*z^2*exp(-48x*z^2)-x-2*z^2)*exp(1-4*x*z^2),x, 1/4*z^2,inf)=0;
这个方程组。此问题:B=[1e-3,1e-4,1e-5,1e-6,1e-7,1e-8,1e-9,1e-10];
最后MATLAB实验数据:v*,z*,为最优解
w
1e-3
1e-4
1e-5
1e-6
1e-7
1e-8
1e-9
1e-10
V*
10.7871
14.5947
18.3964
22.1957
25.9939
29.7915
33.5879
37.3787
Z*
0.27529
0.2373
0.2116
0.1927
0.1781
0.1664
0.1567
0.1485
我想对这个过程进行MATLAB编程