主题:处理一个等式问题
我的问题是: 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编程。
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编程。