|
求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。 2 y) h0 r. `+ g% n5 L. N# g
function f=fh3(n). y! w! w- y6 o
g2=2^0.5-1;, Q& H j/ u- i9 H' ?1 ~: h1 t
mu0=4*pi*10^(-7);
3 n5 o; {7 U5 w* Igaf=7.8*10^3;
7 D1 p) O7 C# y: n5 W- j1 o/ q6 Jgacu=8.9*10^3;
3 {# P. S9 w. x4 K# Xgat=4.5*10^3;
) [, Q& D( N' F, b: qgaa=2.7*10^3;
9 B# i) `6 |7 L+ `+ s' bp=2.7*10^4; Q& R- k& q/ @) q) H4 C
q=3*10^(-5);) R8 [" m" H. r; B u
rou=0.7;
4 e3 c9 L( B( ]- @kb=0.9;
* h, Y2 n4 z; ]( t& N% F r%j0=8*10^6;
9 S M! _: @0 l0 C% Hj0=2*10^6;7 e# h* C! y! F x
kd=3;
" `5 |5 X! s8 \& q+ Abm=2;
; |7 b' U" G2 v- Fht=2*10^(-3);
2 Z2 a6 }: z8 N; c6 ^. Pmu1=10^8/1.75;
% `2 n3 Y7 O$ d# hkt1=2*pi*gat*ht;
4 O, k7 t8 Z& P* Akt2=2*pi*gaa*ht;
6 o2 N5 R- n& w* l9 Pkr=1.02;7 t. o, I" h& d) @7 e+ o
ku=(mu0/2/bm)^0.5;8 C5 V: v u) b. D
kx=mu0^0.5*kb*kd^2*j0;
' N g Z4 Y) Q- q, z6 f8 Ekp=(2*q/pi/p)^0.5;5 K1 V7 n; g7 Z( _0 ~7 H% P2 d
kg=mu0^0.5*kb*kd*j0*g2;
. F' X* q1 u% N0 urx=ku*kx-kp*kg;
. B1 B- \+ w* _2 Pra=ku*kg;3 L) `& t( y8 F @" X
rc=kp*kx;) \( R; Y+ y( f* j: d( \' Z' a
kj=mu0*kb*kd*j0;7 R1 z5 I9 O. M4 W9 {0 g. d
kj=kj/kx;
) d2 q1 f8 S6 Z4 `h0=3*ht;
; ^$ P' n8 y4 W- zh1=2*(1+1/rou)*q/pi;
4 E: ?% k2 p6 r1 ^8 ?5 r# X) w% Xh2=2*(kd-1);
: `1 W/ i( @! ?! o& v( ?h3=1/bm*(mu0*p*q/(2*pi))^0.5;
4 {8 k4 K# {! fh4=1/kj*(mu0*p*q/(2*pi))^0.5;. T( g5 T3 D# z; u ~
hg=2*pi*kr*ht;
+ E& p! ^2 _6 R, q7 Kw1=gaf/bm*(pi*p*q*mu0/2)^0.5;2 @: V/ b/ n0 }* c! C; x! ]
w1p=8*g2*kd*w1;8 K% N. p. F; ^: {# T( Q
w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;
/ e& d( K% B6 | ~2 ?8 Y+ aw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;: k* R2 G `7 p" W& f; K5 @
w3pp=hg*h3;
1 C$ e0 q1 z* ]- Q: U' Hw3=w3p+w3pp;
; o# L4 ~, ? P8 i [9 G) Zw4=2*kd*gacu/kj*(2*pi*p*q*mu0);* l) T& S- R5 f9 h- T
w5=5*pi*ht*gat;
4 [2 V1 i3 o- ]w6=hg*h0*gat;
1 R/ G5 i8 ], D1 e8 q; |6 aw7=hg*h1*gat;( {8 G3 @- x5 @, g: `- r
w8=hg*h0*gat;2 S Y$ B: @" W, t4 @: `% `
w9=hg*h4*gat;/ H; p0 m( v; s" w1 z3 S7 s
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);
' Z* F% N+ h9 s5 d |