|
求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。
/ t4 p' i8 \5 f' q& Y1 ]6 cfunction f=fh3(n)
+ Y5 p% g2 b* ^g2=2^0.5-1;% l! \6 \. \) R* [& c
mu0=4*pi*10^(-7);3 u* \, I9 {, n) U- K
gaf=7.8*10^3;
- j, @& N7 ^% V1 egacu=8.9*10^3;- O6 b4 O' m5 r m
gat=4.5*10^3;/ K0 ?; |5 }" Z* ?6 @
gaa=2.7*10^3; O- n' W7 W" T# ]/ S5 W7 ?! F
p=2.7*10^4;) b1 M f( L# J
q=3*10^(-5);
* _% g3 `" ^ Q9 ~% srou=0.7;
% i q x0 y1 ]1 Z$ V: X/ Nkb=0.9;6 U! h" B+ ^0 N$ q. B0 k+ m
%j0=8*10^6;
: ] z/ n- t: K, @) q: K% Fj0=2*10^6;
8 }, L$ V& w5 i6 H. okd=3;
; n! u) s* [8 y6 ^8 X4 Dbm=2;7 m4 D1 ^+ V/ E/ v- S2 k; D- S K4 ^
ht=2*10^(-3);
0 s! v* O( f; Xmu1=10^8/1.75;
, m6 A5 e1 t+ @9 H: H0 qkt1=2*pi*gat*ht;
/ M& B+ Y e4 ]' Vkt2=2*pi*gaa*ht;; ]/ P1 @$ Y" j$ g% k; ~' d, a
kr=1.02;
9 F Y5 y8 Y, s1 z5 C$ n( sku=(mu0/2/bm)^0.5;
" t4 T7 g% O K: J! f4 Okx=mu0^0.5*kb*kd^2*j0;
+ N3 u; H" n5 [1 A/ T: X. Skp=(2*q/pi/p)^0.5;
" f: j& i( b8 R5 [kg=mu0^0.5*kb*kd*j0*g2;
( Q% T$ w& ^7 M9 _4 W* w) t( f R/ R& rrx=ku*kx-kp*kg;' ]' i1 G7 I4 G9 U0 ~
ra=ku*kg;
$ _: m3 \0 G1 S+ `. \: X. Nrc=kp*kx;
7 S' f, d0 q) Y9 R* r% V5 w) P3 lkj=mu0*kb*kd*j0;
& i# R" \; X5 vkj=kj/kx;& L" m9 S: M6 ^& [' f8 C" f
h0=3*ht;
3 H& u, P8 ~) n7 rh1=2*(1+1/rou)*q/pi;; E& h( ]( K% [: S7 ?
h2=2*(kd-1);
1 J' n6 W2 e; Y0 A. Bh3=1/bm*(mu0*p*q/(2*pi))^0.5;
~( N6 x8 I5 S6 \% \6 Ah4=1/kj*(mu0*p*q/(2*pi))^0.5;
) {) P, x; }$ E) x$ I# p2 ~8 E4 ohg=2*pi*kr*ht;
I# Y6 j4 s! R5 O6 o8 D6 f( p6 tw1=gaf/bm*(pi*p*q*mu0/2)^0.5;
# _& A) y2 z0 T. g2 n4 j7 Sw1p=8*g2*kd*w1;
7 r: z3 g. S9 S8 P8 Rw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;
' C1 H$ S3 Y/ j. j( d0 mw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;+ f$ f2 f$ r% z7 k
w3pp=hg*h3;+ U4 Y% u/ `* C& b6 o3 A+ d
w3=w3p+w3pp;3 D/ p4 ~+ m3 W
w4=2*kd*gacu/kj*(2*pi*p*q*mu0);
( W/ m' k. v6 }w5=5*pi*ht*gat;/ j; u$ n' ` \* a5 |6 t9 D9 }
w6=hg*h0*gat;6 P) g8 ]0 Z3 Z" o+ Z( r( ^" W
w7=hg*h1*gat;1 r# i7 {; s4 r, {
w8=hg*h0*gat;( @7 Q5 s& y7 ]9 C
w9=hg*h4*gat;+ a5 P, y/ h9 n+ R# o! ~6 z1 x
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);
6 t, }- s" O. R |