求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。
0 m; ]2 K' z afunction f=fh3(n)
, ^* ]. e; A) n( g' d* eg2=2^0.5-1;
8 h% c, x; q( `6 N# hmu0=4*pi*10^(-7);# l/ c4 g- h I* _
gaf=7.8*10^3;
Z* R% F+ f& ]5 N6 ^gacu=8.9*10^3;
) X0 k7 ~# I% G' [/ V4 ugat=4.5*10^3;
1 t9 _( R5 j6 ?; B$ k Mgaa=2.7*10^3;
- `8 {0 S. I/ @% Fp=2.7*10^4;2 |9 g4 B- n8 ~% x; c
q=3*10^(-5);
8 @3 P0 k+ F# L F! K% arou=0.7;
5 `* i* _2 Z, B! |kb=0.9;
. e# H" E2 K- c- V) ~%j0=8*10^6; @4 d* N, o: T% ~
j0=2*10^6;5 J7 [+ t% O6 H
kd=3;
% d0 e" t7 ]6 O- f V: u( ybm=2;
* F& Z4 W$ H6 }: `ht=2*10^(-3);) E7 u7 T' R! l
mu1=10^8/1.75;: R7 [6 @+ w# w, N
kt1=2*pi*gat*ht;
6 N# o# B: e3 mkt2=2*pi*gaa*ht;
$ D! e6 \9 @ bkr=1.02;" _, P1 J# o9 \9 Q- I% o
ku=(mu0/2/bm)^0.5;$ |/ n, U: R' u* Q! c, p% D0 F( M
kx=mu0^0.5*kb*kd^2*j0;
, k2 a) k6 C) [$ s m5 P- \9 _9 skp=(2*q/pi/p)^0.5;
. Z3 B7 j8 S- o( h! t& ^, k1 _, Ukg=mu0^0.5*kb*kd*j0*g2;
$ P8 I1 c# P* U, z# _$ P7 i$ mrx=ku*kx-kp*kg;0 F+ a3 n% G/ p4 _6 ^
ra=ku*kg;& h% ~) _7 `4 p) {- q2 ^
rc=kp*kx;
$ |. D& c6 B8 E }& [; s! Jkj=mu0*kb*kd*j0;: y: U! k3 z& I7 s/ s0 }6 S: G1 l
kj=kj/kx;/ V) C! S: L1 g. I4 o& a/ d
h0=3*ht;
6 w9 s+ N& u# ch1=2*(1+1/rou)*q/pi;5 T, C' W5 D$ m- N
h2=2*(kd-1);+ }0 |8 J1 Z& M
h3=1/bm*(mu0*p*q/(2*pi))^0.5;& p: {) r# k% u1 Z
h4=1/kj*(mu0*p*q/(2*pi))^0.5;0 B! @ O; }& w! h Z
hg=2*pi*kr*ht;
0 }( t4 S; _, p8 q. sw1=gaf/bm*(pi*p*q*mu0/2)^0.5;
3 Y* ~9 Y9 y' m6 sw1p=8*g2*kd*w1;
5 {' }7 F0 e/ s6 ?* V2 Nw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;+ }3 p1 @( H/ T* | \3 k( U: r
w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;% n/ I" u, y7 G+ t
w3pp=hg*h3;' n% m8 ?' |- q
w3=w3p+w3pp;
' r$ x- d( K- ~3 ~( \, Q: i" Kw4=2*kd*gacu/kj*(2*pi*p*q*mu0);
( y4 d3 [! Y0 F& z4 J/ Y" Y" R- U" Lw5=5*pi*ht*gat;7 o* S% q' {. R4 U2 c
w6=hg*h0*gat;5 t! A% O+ L% K+ [/ \
w7=hg*h1*gat;
( v4 C; E& _ s: F+ _0 @; Jw8=hg*h0*gat;
5 R1 S6 S5 b8 ?* _4 fw9=hg*h4*gat;3 X; [+ F5 ~/ x+ x! 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' Y) H1 q. o4 ^
|