|
求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。
* ` w+ ^' Z3 r: I( C4 G: Cfunction f=fh3(n)
" d9 o5 r, V# o, D4 ? g& u9 ^ |g2=2^0.5-1; Y5 I1 w6 a8 o+ J- x" ~
mu0=4*pi*10^(-7);
# s# e2 O6 z. }/ _' [+ n7 Q1 lgaf=7.8*10^3;
# S* Z/ Z& l6 ?+ R5 K' Wgacu=8.9*10^3;
( M- G6 n+ W3 z( o J+ A2 bgat=4.5*10^3;3 @" s& L( A' S# ~6 f; t+ a, v0 r
gaa=2.7*10^3;
4 m/ k3 ~# @$ P! {* Op=2.7*10^4;2 A+ O) i9 K: s8 i, |
q=3*10^(-5);
2 L0 w$ X) c3 Erou=0.7;3 E) W' P" b) [9 O* m; w' v" C
kb=0.9;) `5 t6 p9 q' `- j6 K3 N
%j0=8*10^6;
$ _- I" l$ D6 H5 m: a2 N; d2 Q6 J' Cj0=2*10^6;
4 Y- W- d$ D' \4 S# ^kd=3;
: T# b6 s( I3 G& k) ] ebm=2;
, n( |0 r3 w3 G0 ~) \. ^- Dht=2*10^(-3);
" f2 u% `) u( A1 _mu1=10^8/1.75;( g, {4 G/ E" X, r$ v
kt1=2*pi*gat*ht;3 y+ F5 b" S2 Z6 `* q+ W
kt2=2*pi*gaa*ht;% p$ P4 H0 F, u' B9 G+ d" ?/ d W3 t
kr=1.02;
& U e& \2 w- M0 k" m* ~ cku=(mu0/2/bm)^0.5;
" W: j$ z/ N% {+ W! [( \: Okx=mu0^0.5*kb*kd^2*j0;# b+ T3 o. k$ H; J6 U# t: C6 n& ]
kp=(2*q/pi/p)^0.5;. ? t4 I% J& ~6 P; A
kg=mu0^0.5*kb*kd*j0*g2;
- ?9 v; E8 O5 |9 hrx=ku*kx-kp*kg;/ g! y# h8 K+ R9 A) |/ Q
ra=ku*kg;3 ~) r# r) `4 e( c9 M
rc=kp*kx;
+ }5 Y+ }1 `7 V9 v5 D% j* Zkj=mu0*kb*kd*j0;1 c2 x( w }, o
kj=kj/kx;
' | w% ?9 b0 Yh0=3*ht;
7 @: p+ X6 b3 F& r8 b) jh1=2*(1+1/rou)*q/pi;
- [3 K# R' H. jh2=2*(kd-1);5 Z( {* D" g: Z0 g+ p' ]
h3=1/bm*(mu0*p*q/(2*pi))^0.5;& |; Q" Q6 @' e& ~5 Z2 V: w9 W
h4=1/kj*(mu0*p*q/(2*pi))^0.5;7 W$ E7 J$ }3 o
hg=2*pi*kr*ht;
. W' ]2 B1 x% N# o0 p7 l, B" dw1=gaf/bm*(pi*p*q*mu0/2)^0.5;
c) y9 i& W. A7 Iw1p=8*g2*kd*w1;; q$ S# o) Q9 _, Y. F6 d! z
w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;
( r8 Y! m- [; i+ G/ Ew3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;8 ?0 N, Q$ d9 x4 t8 F1 z# f H( u# ^
w3pp=hg*h3;9 V' z: S" g. U& @' _3 F% |
w3=w3p+w3pp;8 A7 u! ?0 o6 C7 x- u5 n6 d) Z
w4=2*kd*gacu/kj*(2*pi*p*q*mu0);0 q0 Y# y9 B; y6 u! X {" i, p
w5=5*pi*ht*gat;
( d% [3 Q% J, T( B2 d* [: Zw6=hg*h0*gat;" t ?- u) H& K/ L5 u+ b0 J
w7=hg*h1*gat;& H: V' N. E" S: f5 u) H
w8=hg*h0*gat;
5 a2 a& e( c% Z$ iw9=hg*h4*gat;
) Q) f8 S: s+ w2 \" uf=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);
0 b( u( P+ Z! K7 \3 X7 r: { |