|
求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。
: [. B1 }# @. X6 bfunction f=fh3(n) j$ t+ x- E! i2 m* f3 i: U
g2=2^0.5-1;8 p# e; `. q: Z
mu0=4*pi*10^(-7);% {" \2 E. S" Q1 v
gaf=7.8*10^3;
: {4 ^$ _4 u! Y5 ?, C- _gacu=8.9*10^3;7 S1 i6 z C8 Z1 ?9 [: e
gat=4.5*10^3;7 P3 N6 p# J" c$ v a
gaa=2.7*10^3;
9 Y# A* B# c& k- H" ^$ U- Pp=2.7*10^4;, G n( \0 [" o6 K) Q+ J: n( Y
q=3*10^(-5);. I4 H, l* E: @& d
rou=0.7;
. G1 N x- B1 k. C2 k" f# \: ekb=0.9;
0 O; c+ o2 T/ |% M) N z%j0=8*10^6;
/ P' r$ ~" ]2 d/ `7 ^& K" W( G$ }j0=2*10^6;
4 C' w0 t1 j4 Q! _kd=3;( b, W0 Y' D: P
bm=2;" \" k4 f, z7 v' ?0 R6 c9 t# b2 ~, ~
ht=2*10^(-3);. r5 d+ a# y! B- ^. J |4 G
mu1=10^8/1.75;: ]& ?; E7 m! L5 P- E% X
kt1=2*pi*gat*ht;
& E6 s5 g8 o* Y/ }kt2=2*pi*gaa*ht;0 Q6 d9 k, R8 p9 a
kr=1.02;
1 W/ m& |/ Q4 L% a$ d0 bku=(mu0/2/bm)^0.5;( c9 f/ h5 c' S/ g# b6 D' x( N
kx=mu0^0.5*kb*kd^2*j0;6 J. Z, D0 G, v* g. R/ {- {
kp=(2*q/pi/p)^0.5;, E( D& }9 }- D* }# V
kg=mu0^0.5*kb*kd*j0*g2;( ~8 A" E5 a; B8 N0 G
rx=ku*kx-kp*kg;
1 H, \: E( D% @) Q3 K7 }+ nra=ku*kg;7 c: z' O7 H: I" W
rc=kp*kx;) L5 f( ], K) y9 A2 ?9 E+ k
kj=mu0*kb*kd*j0;
0 x+ ^' ?- g& y- W1 Skj=kj/kx;
' j# ^) b( s/ o) ]% p% c; e/ ]3 T; ^h0=3*ht;
% W+ Y% t! m% m; d- E3 W' vh1=2*(1+1/rou)*q/pi;
$ i: L1 }5 m2 Ch2=2*(kd-1);
# A) E; P1 f4 {* }h3=1/bm*(mu0*p*q/(2*pi))^0.5;- w' w+ f' Y; ^
h4=1/kj*(mu0*p*q/(2*pi))^0.5;1 x1 P+ [9 E: |
hg=2*pi*kr*ht;
- [. ^8 O8 m% n) Uw1=gaf/bm*(pi*p*q*mu0/2)^0.5;
! E" g1 W; a5 h+ P# E! G0 `w1p=8*g2*kd*w1;
9 Y) p2 }0 G: k7 \( q) _& ~w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;' |0 d4 j p5 q: g. L
w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;3 M3 M( k9 Z4 T0 w
w3pp=hg*h3;
- L% Y% W+ t6 e+ V) R4 I, e* qw3=w3p+w3pp;) D6 y( Y! {8 y$ c
w4=2*kd*gacu/kj*(2*pi*p*q*mu0);2 F4 J7 n& Y5 s
w5=5*pi*ht*gat;2 K `0 L2 V/ @5 v2 f% r
w6=hg*h0*gat;
1 `! g/ O1 \* hw7=hg*h1*gat;$ A/ J: A, w9 [' Y1 F
w8=hg*h0*gat;& b' e* k: V9 }" }4 L0 l& r; X, e0 [
w9=hg*h4*gat;4 f5 y! z1 R! w7 `/ _* V
f=w1p*r*x^0.5+w2*x^1.5/(g2*r-kd*x)+w3*r*x^(-0.5)+w4*x^0.5+w5*r^2+w6*r+w7/r+w8*r*x+w9*r*x^(-0.5)/(g2*r-kd*x);
+ x" k- Q) Y8 i4 Z h' @8 M6 C# _5 R |