Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учёба / lab_mat-met

.pdf
Скачиваний:
8
Добавлен:
03.03.2016
Размер:
644 Кб
Скачать
☆

^bn_j_gp•cg• j•\gyggy lZdh` klZxlv g_e•g•cgbfb >ey pbo \biZ^d•\ lj_[Z g_e•g•cgm oZjZdl_jbklbdm Zijhdkbfh\m\Zlb qb •gl_jihex\Zlb gZ dh`ghfm djhp• ^bn_j_gpx\Zggy Hibk d_jmxqbo kb]gZe•\ lZ j•\gyggy ZijhdkbfZp•€ •gl_jiheyp•€ lj_[Z jhalZrh\m\Zlb \ nmgdp•€ j•\gygv a\¶yadm

Jha]eyg_fh ijbdeZ^

 

 

 

 

R

 

 

< _e_dljbqg•c ko_f• ^h RL-eZgpx]Z i•^dex-

 

 

 

 

 

 

 

qZ}lvky ^`_j_eh ihkl•cgh€ gZijm]b U 1g^mdlb\g•klv

 

 

 

 

 

 

 

 

i

 

 

 

dhlmrdb gZfZ]g_qm\Zggy } af•ggZ \_ebqbgZ ydZ aZ\^ZgZ

 

 

 

U

L

 

 

 

 

lZ[ebqghx nmgdp•}x Φ T

(iT ) lh[lh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L =

dΦ

= var

(12.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

iT

=

[iT1

, iT 2 , ..., iTn

] ;

 

N%

 

 

= [Φ

T1 , Φ T 2 , ..., Φ Tn ] .

 

Φ

T

 

 

AZ

11

aZdhghf

D•jo]hnZ

gZijm]Z

\

_e_dljbqg•c ko_f• \bagZqZ}lvky aZ nhjfmehx

iT

 

 

 

 

 

 

U =

iR +

L

di

.

(12.2)

 

 

 

 

 

 

 

dt

 

 

 

I•key i•^klZgh\db m hljbfZ}fh

U = iR +

dΦ

 

di

= iR +

dΦ

,

 

 

 

 

 

 

di dt

 

dt

Z[h

 

 

 

 

 

 

 

 

 

 

dΦ

= U −

iR – ^bn_j_gp•cg_ j•\gyggy m ghjfZevg•c nhjf• Dhrb

 

dt

 

 

 

 

 

 

 

 

 

 

i =

f (Φ )

– g_e•g•cg_ j•\gyggy a\yadm

Φ (0) = 0 – ihqZldh\• mfh\b

>ey \•^hfh€ lZ[ebqgh€ aZe_`ghkl• Φ T = f (iT ) fh`gZ jhajZom\Zlb Φ (i) ^ey [m^v-ydh]h agZq_ggy • aZ ^hihfh]hx ZijhdkbfZp•€ qb •gl_jiheyp•€ Ijb ZijhdkbfZp•€ djb\bo gZfZ]g_qm\Zggy a\bqZcgh \bdhjbklh\mxlv klmi•g_\• ihe•ghfb -]h Z[h -]h ihjy^dm h[h\¶yadh\h g_iZjgh]h lhfm sh nmgdp•y Φ (i) p_gljZevgh kbf_ljbqgZ

50

Nmgdp•y ydZ f•klblv m kh[• hibk g_e•g•cgbo ^bn_j_gp•cgbo j•\gygv lZ j•\gyggy a\yadm fh`_ fZlb gZklmigbc \b]ey^ m \biZ^dm

ZijhdkbfZp•€

•gl_jiheyp•€

function dy=frn_a(t,y)

function dy=frn_i(t,y)

global R U A

global R U it Ft

F=y(1);

F=y(1);

 

i=polyval(A,F);

i=interp1(Ft,it,F);

dy(1)=U-i*R;

dy(1)=U-i*R;

 

=heh\gZ ijh]jZfZ fh`_ fZlb gZklmigbc \b]ey^ ^ey \biZ^d•\

ZijhdkbfZp•€

•gl_jiheyp•€

global R U A

global R U it Ft

U=…

R=… y0=… tk=…

U=…

R=… y0=… tk=…

it=[…]

Ft=[…]

it=[…]

Ft=[…]

A=polyfit(it,Ft,5);

[t,y]=ode23(‘frn_i’,[0 tk],y0);

[t,y]=ode23(‘frn_a’,[0 tk],y0);

i=interp1(Ft,it,y);

i=polyval(A,y);

plot(t,y,t,i)

plot(t,y,t,i)

 

 

M k_j_^h\bs• Simulink ML ^ey •gl_jiheyp•€ \bdhjbklh\mxlv [ehd³Look-Up- Table´ a [•[e•hl_p• ³Functions&Tables´ 0/ ydbc a^•ckgx} dmkdh\h-e•g•cgm ZijhdkbfZp•x aZ\^Zgh€ lZ[ebqgh€ nmgdp•€ IZjZf_ljZfb pvh]h [ehdZ } \_dlhj Zj]mf_gl•\ lZ \_dlhj agZq_gv M a\¶yadm a pbf 6LPXOLQN-fh^_ev fh`_ fZlb gZklmi gbc \b]ey^

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

i

i

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Integrator

 

 

 

Look-Up

To

W orkspace

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Table1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G ain

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Look-Up

 

 

 

 

 

 

 

To W orkspace1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Table

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Clock

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

To W orkspace2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

51

AZ\^Zggy

JhajZom\Zlb i_j_o•^g• ijhp_kb m ko_fZo gZ^Zgbo gZ jbk IZjZf_ljb ko_f m lZ[ebp• oZjZdl_jbklbdb g_e•g•cgbo _e_f_gl•\ L i R gZ^Zg• m lZ[e m klZ[•e•ljhgZ oZjZdl_jbklbdZ •^_ZevgZ jbk Ih[m^m\Zlb ]jZn•db kljmf•\ m ]•edZo _e_dljbqgh€ ko_fb m nmgdp•€ qZkm

 

R

 

 

R /2

R /3

 

R

 

 

R

 

R

 

L

 

 

 

 

 

 

 

U

 

U

R

L

U

C

 

U

C

 

 

 

 

 

 

1

 

 

2

 

 

3

 

 

4

 

R V D 1 V D 2

 

R

 

 

L9'

 

 

 

 

 

 

 

V D 1

 

 

8kl

89'

 

 

 

 

 

L

 

 

L

 

 

 

U

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V D 2

 

 

 

 

 

 

 

5

 

 

 

6

 

 

7

 

 

 

 

 

5,6) U = 0.1Em sin125t

 

 

 

 

 

 

 

 

 

 

 

Jbkmghd

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LZ[ebpy

U=Em,B

‹

‹ R Hf

‹

‹ C fdN

‹ ‹ ko

Ukl,

 

 

\Zj

ko

 

\Zj

ko

 

\Zj

 

<

200

1

 

1

21

7

3

150

13

5

12

220

2

 

1

23

8

3

170

14

5

18

240

3

 

1

25

9

3

190

15

5

24

250

4

 

2

22

10

4

100

16

6

16

270

5

 

2

24

11

4

120

17

6

20

300

6

 

2

26

12

4

140

18

6

22

 

 

 

 

 

 

 

 

 

 

LZ[ebpy

N <[

0

0.181

0.33

0.451 0.551 0.632

0.699

0.753

0.798

0.835

0.865

i :

0

2

4

6

8

10

12

14

16

18

20

U <

0

66.4

118

158

190

214

233

248

259

268

275

i :

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

 

 

 

 

 

 

52

 

 

 

 

 

EZ[hjZlhjgZ jh[hlZ ‹

J1R?GGY LJ:GKP?>?GLGBO L: :E=?;J:1QGBO J1<GYGV

P•ev jh[hlb gZ\qblbky jha\¶yam\Zlb ljZgkp_^_glg• lZ Ze]_[jZ•qg• j•\gyggy

L_hj_lbqg• \•^hfhkl• Qbk_evgbc jha\¶yahd j•\gyggy

F(x)=0

(13.1)

ih^Ley}lvky gZ ^\Z _lZib: \•^^•e_ggy dhj_g•\ lZ mlhqg_ggy €o ihqZldh\bo gZ[eb`_gv •l_jZp•cgbfb f_lh^Zfb

GZc[•evr jhaih\kx^`_gbfb f_lh^Zfb mlhqg_ggy dhj_g•\ } f_lh^b [•k_dp•c ohj^ ^hlbqgbo lZ f_lh^ ijhklbo •l_jZp•c

M iZd_l• ML j•r_ggy g_e•g•cgbo j•\gygv \bdhgm} nmgdp•y

x = fzero ( Fun, x0, Options, p1, p2…)

Fun – kljhdh\Z af•ggZ •f¶y m-nZceZ m ydhfm f•klblvky hibk nmgdp•€ Z[h Matlab-\bjZa ^ey €€ jhajZomgdm

x0 – ijb[ebag_ agZq_ggy dhj_gy fh`_ [mlb kdZeyj Z[h \_dlhj jhaf•jhf M ^jm]hfm \biZ^dm j•r_ggy rmdZ}lvky gZ aZ\^Zghfm •gl_j\Ze•

Options – \_dlhj hip•c yd• d_jmxlv ijhp_khf h[qbke_ggy lZ \b\h^m j_amevlZl•\

p1, p2– \bdhjbklh\mxlv ydsh nmgdp•y fZ} [•evr h^gh]h Zj]mf_glZ F(x, p1, p2…);

Options, p1, p2 … – g_h[h\'yadh\• iZjZf_ljb

>ey j•r_ggy Ze]_[jZ•qgbo j•\gygv \ iZd_l• ML } nmgdp•y roots.

X = roots (P) – h[qbkex} \_dlhj dhj_g•\ X ihe•ghfZ a dh_n•p•}glZfb P.

F_lh^ [•k_dp•c

F_lh^ [•k_dp•c Z[h f_lh^ iheh\bggh]h ^LeHggy kdeZ^Z}lvky a ihke•^h\gh]h ih^•e_ggy \•^j•adZ ydbc f•klblv dhj_gv gZ\i•e

x=(a+b)/2,

53

^_ Z • b – e•\Z lZ ijZ\Z f_`• dhj_gy lh[lh

 

b>Z,

(13.2)

F(a)*F(b)<0.

(13.3)

>ey dh`gh]h gZklmigh]h ^•e_ggy h[bjZ}lvky lZ iheh\bgZ \•^j•adm gZ d•gpyo ydh€ nmgdp•y fZ} ijhlbe_`gbc agZd Ijb pvhfm •gl_j\Ze •kgm\Zggy dhj_gy a\m`m}lvky aZ jZomghd af•gb h^g•}€ a ch]h f_`: e•\h€ Z o) Z[h ijZ\h€ b o).

1l_jZp•cgbc ijhp_k aZd•gqm}lvky ydsh \bdhgm}lvky mfh\Z:

 

E − D ≤ ε,

(13.4)

^_ ε - aZ\^ZgZ lhqg•klv h[qbke_ggy dhj_gy ε<<1 1gdheb

\bfZ]Zxlv sh[

h^ghqZkgh \bdhgm\ZeZky mfh\Z

 

 

) [

 

≤ ε.

(13.5)

 

 

F_lh^ [•k_dp•c – ijhklbc lZ gZ^•cgbc kihk•[ ihrmdm dhj_g•\ j•\gyggy F(x)=0 <•g a[•]Z}lvky ^ey [m^v-ydbo [_ai_j_j\gbo nmgdp•c F(x), m lhfm qbke• g_^bn_j_gpx}fbo R\•^dbklv a[•`ghkl• g_\_ebdZ AZjZ^b ^hky]g_ggy lhqghkl• ε

g_h[o•^gh \bljZlblb

 

 

 

E −

D

 

1 = ORJ

 

 

 

 

(13.6)

ε

 

 

 

 

 

 

•l_jZp•c P_ ihagZqZ} sh aZ^ey ihemq_ggy dh`gbo \•jgbo ^_kylbqgbo agZd•\ g_h[o•^gh ajh[blb [ebavdh •l_jZp•c

Ydsh gZ \•^j•adm [a,b] agZoh^blvky d•evdZ dhj_g•\ lh ijhp_k a[•]Z}lvky ^h h^gh]h a gbo F_lh^ g_ijb^Zlgbc ^ey ihrmdm djZlgbo dhj_g•\ iZjgh]h ihjy^dm

F_lh^ ohj^

F_lh^ ohj^ Z[h f_lh^ ijhihjp•cgbo qZklbg f•klblvky m ihke•^h\ghfm ih^•e_gg• \•^j•adZ >a,b@ ydbc fZ} dhj_gv gZ qZklbgb ijhihjp•cg• agZq_ggyf

nmgdp•€ gZ d•gpyo \•^j•adm

 

 

) D

=

D −

[

,

(13.7)

 

 

E −

 

 

) E

[

 

 

a\•^db

54

[ =

D) E −

E) D

.

(13.8)

) E −

 

 

) D

 

=_hf_ljbqgh p_ _d\•\Ze_glgh aZf•g• ]jZn•dZ nmgdp•€ F(x ohj^hx ydZ ijhc^_ q_j_a lhqdb a,F(a lZ b,F(b)).

>ey aZd•gq_ggy •l_jZp•cgh]h ijhp_km aZf•klv mfh\b \bdhjbklh\mxlv

mfh\m

 

[L+ − [L ≤ ε,

(13.9)

^_ xi+1, xi – \•^ih\•^gh hklZgg•c h[qbke_gbc • ihi_j_^g•c ^h gvh]h gZ[eb`_ggy dhj_gy M j_rl• p_c f_lh^ ihoh^blv gZ f_lh^ [•k_dp•c Ze_ aZ[_ai_qm} [•evr r\b^dm a[•`g•klv

F_lh^ ^hlbqgbo

F_lh^ ^hlbqgbo Z[h f_lh^ GvxlhgZ ihey]Z} m ihke•^h\g•c ZijhdkbfZp•€ nmgdp•€ F(x ^hlbqgbfb ^h djb\h€ m lhqp• ihi_j_^gvh]h gZ[eb`_ggy ([L ) [L ) yd• i_j_lbgZxlv hkv Z[kpbk m lhqp• gZklmigh]h gZ[eb`_ggy o• 1 sh \bagZqZ}lvky aZ nhjfmehx

[

 

= [ −

)([L )

.

(13.10)

+

) ([ )

L

L

 

 

 

 

 

L

 

 

Ihke•^h\g•klv a[•]Z}lvky ^h ^•ckgh]h agZq_ggy dhj_gy j•\gyggy F(x ydsh gZqZevg_ gZ[eb`_ggy dhj_gy e_`blv m •gl_j\Ze•>a,b] (F(a)*F(b gZ ydhfm iho•^g• ) ([) lZ ) ([) ljbfZxlv \eZkgbc agZd lZ \bdhgm}lvky mfh\Z

) [ ) [ > .

(13.11)

1l_jZp•€ ijbibgyxlv ijb \bdhgZgg• mfh\ lZ Z[h

F_lh^ GvxlhgZ _n_dlb\gbc ydsh \•^hf_ ]Zjg_ ihqZldh\_ gZ[eb`_ggy ^ey dhj_gy lZ ih[ebam dhj_gy ]jZn•d nmgdp•€ fZ} \_ebdm djmlbagm. M [eZ]hijb}fgbo \biZ^dZo qbkeh \•jgbo ^_kylbqgbo agZd•\ m q_j]h\hfm gZ[eb`_gg• ih^\hx}lvky, lh[lh ijhp_k a[•]Z}lvky ^m`_ r\b^dh

55

G_^he•dhf f_lh^Z ^hlbqgbo } g_h[o•^g•klv jhajZoh\m\Zlb m dh`g•c lhqp• g_ l•evdb agZq_ggy nmgdp•€, Ze_ • agZq_ggy iho•^gbo

F_lh^ ijhklbo •l_jZp•c

F_lh^ ijhklbo •l_jZp•c ihey]Z} m aZf•g• \bo•^gh]h j•\gyggy F(x)=0 _d\•\Ze_glgbf chfm j•\gyggyf

[ = ϕ [ (13.12) lZ h[q•kex\Zggy ihke•^h\g•kl•

[L + = ϕ [L ,

(13.13)

(i=1 ydZ a[•]Z}lvky ijb L → ∞ ^h lhqgh]h j•r_ggy. 1l_jZp•x ijbibgyxlv ijb \bdhgZgg• mfh\b

>hklZlgvhx lZ g_h[o•^ghx mfh\hx a[•`ghkl• f_lh^Z }

 

ϕ [

 

< .

(13.14)

 

 

R\b^d•klv a[•`ghkl• ajhklZ} a• af_gr_ggyf ϕ [ .

AZ\^Zggy

H[qbkeblb i_jrbc ^h^Zlgbc dhj_gv ljZgkp_^_glgh]h j•\gyggy • mk• ^•ckg• dhj_g• Ze]_[jZ€qgh]h j•\gyggy a lZ[ebp• a lhqg•klx -3 , 10-5 • -7 .

F_lh^bqg• j_dhf_g^Zp•€

GZqZevg• gZ[eb`_ggy dhj_g•\ o0 \bagZqblb \•amZevgh ih[m^m\Z\rb ]jZn•d nmgdp•€ F(x).

I_jr g•` \bdhjbklh\m\Zlb f_lh^b ^hlbqgbo Z[h ijhklbo •l_jZp•c i_j_\•jblb €o a[•`g•klv

>ey dhgljhex ijZ\bevghkl• j•r_gv \b\_klb gZ _djZg ihke•^h\g•klv gZ[eb`_gv dhj_g•\ lZ agZq_ggy nmgdp•€ m pbo lhqdZo

Hp•gblb r\b^d•klv a[•]m j•agbo f_lh^•\

56

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LZ[ebpy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

‹

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(x)

 

F_lh^ j•r_ggy

 

\Zj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

3

4

 

1

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

4,2x3-31,92x2+74,3x-51,87

;bk_dp•c

 

 

VLQ [

+

 

 

 

−

 

[ +

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

FRV([ + ) −

 

 

[ =

 

 

 

3,6x3-172,8x2+5,184x-237,32

Ohj^

 

3

 

sin3 (x) −

 

x +

1 = 0

 

 

5,8x3-47,56x2+121,2x-97,02

;bk_dp•c

 

4

 

DUFWJ([) −

[ +

 

=

 

 

 

6,1x3-90,28x2+388,2x-506,2

Ohj^

 

5

([ +

) FRV( [)

 

 

−

 

=

 

 

 

3,6x3-39,96x2+12,17x+426,4

;bk_dp•c

 

6

0.25 + arctg(0,58x +

0,1) −

x2 = 0

 

2.7x3-37,26x2+16,71x-202,7

Ohj^

 

7

 

OQ([) −

 

 

 

 

 

 

 

 

=

 

 

1,3x3-5,98x2-1,09x+13,76

Ijhklbo •l_jZp•c

 

 

 

[ +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

1

 

 

 

 

 

 

 

−

 

x2 =

0

 

4,5x3-26,1x2+176,6x-112,4

Ijhklbo •l_jZp•c

 

 

 

cos(1,05 −

x)

 

 

 

 

 

 

9

 

OQ([) −

 

[

+

 

=

 

 

5,1x3-62,22x2+142,7x+109,2

>hlbqgbo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

OQ([) −

 

 

 

 

 

=

 

 

 

 

 

1,6x3-3,04x2-29,18x+8,98

>hlbqgbo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VLQ( [)

[

 

 

 

 

 

 

 

 

 

-2,3x3+0,23x2+17,05x+13,48

 

 

11

 

 

− [ =

 

 

;bk_dp•c

 

12

 

 

−

 

−

−

=

 

 

1,6x3-14,24x2+38,13x-29,02

Ohj^

 

13

 

FRV( [) +

=

 

 

 

5,3x3-36,04x2+12,25x+28,05

Ohj^

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

+

 

−

([ −

 

 

) =

 

 

-2,6x3+4,68x2+14,38x+3,822

;bk_dp•c

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

−

−

−

[ =

 

 

 

-1,5x3-14,25x2-37,98x-22,03

>hlbqgbo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

57

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ijh^h\`_ggy lZ[ebp•

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

16

 

[ FRV( [) +

=

 

 

 

3,4x3-46,58x2+127,3x-60,34

 

>hlbqgbo

 

 

 

 

 

 

 

 

 

 

 

 

 

17

 

[ −

[− =

 

 

 

 

2,8x3-25,76x2+6,18x+107,4

 

;bk_dp•c

 

18

 

VLQ([) −

 

[ =

 

 

 

 

-1,4x3-10,78x2-22,54x-11,85

 

Ohj^

 

19

DUFWJ

([

)

+

 

[

=

 

 

 

3,1x3-62,6x2+414,7x-898,9

 

Ohj^

 

 

 

 

−

 

 

 

 

 

 

 

20

 

([ −

 

) −

[

=

 

 

 

1,6x3-12,48x2+25,04x-8,12

 

;bk_dp•c

 

 

 

 

 

 

 

 

 

 

21

 

[ − VLQ([) =

 

 

5,4x3-54x2+140,6x-73,8

 

>hlbqgbo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

H [ −

[ −

=

 

 

 

 

2,7x3-17,6x2-45,4x+123

 

;bk_dp•c

 

 

 

 

 

 

 

 

 

 

 

 

23

FRV([ +

) −

[ =

 

 

 

-1,8x3-5,58x2+1,5x+119

 

;bk_dp•c

 

 

 

 

 

 

 

 

 

24

DUFWJ([) − [ =

 

 

-2,5x3+8,25x2+61,9x-117

 

Ijhklbo •l_jZp•c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

 

H−

[

 

+

[ −

=

 

 

 

 

1,7x3-16,43x2+21,4x-12

 

;bk_dp•c

 

26

 

H

−

H + +

=

 

 

 

3,7x3-17,9x2-55,1x+105

 

Ohj^

 

27

 

VLQ( [) −

[ =

 

 

-3,1x3+9,55x2+71x-111

 

>hlbqgbo

 

28

DUFWJ

([)

−

 

(

 

[)

 

 

(

[)

=

 

3,4x3-46,58x2+127,3x-60,34

 

Ohj^

 

 

 

FRV

 

VLQ

 

 

 

 

 

29

 

H−

 

 

−

[ =

 

 

 

 

 

5,3x3-36,04x2+12,25x+28,05

 

;bk_dp•c

 

30

 

FRV ([) −

[ =

 

 

 

4,2x3-31,92x2+74,3x-51,87

 

Ijhklbo •l_jZp•c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58

EZ[hjZlhjgZ jh[hlZ ‹

JHA<¶YAM<:GGY KBKL?F G?E,G1CGBO J,<GYGV

P•ev jh[hlb gZ\qblbky jha\¶yam\Zlb kbkl_fb g_eLg•cgbo gbfb f_lh^Zfb

L_hj_lbqgL\L^hfhklL

Kbkl_fZ n jL\gygv a n g_\L^hfbfb fZ} \b]ey^

I

([ [ [Q ) =

 

I

 

([ [ [Q ) =

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

Q

) =

 

 

[

[

 

I

 

([

 

 

:[h m \_dlhjghfm \b]ey^•

F(X) = 0, X =[x1, x2, …xn].

jL\gygv Ll_jZpLc-

(14.1)

>ey jLr_ggy lZdbo g_eLg•cgbo kbkl_f \bdhjbklh\mxlvLl_jZpLcgLf_lh^b M iZd_l• ML j•r_ggy kbkl_f g_e•g•cgbo j•\gygv \bdhgm} nmgdp•y

X = fsolve ( Fun, X0, Options, p1, p2…)

Fun – kljhdh\Z af•ggZ sh f•klblv •f'y m-nZceZ m ydhfm aZibkZgZ kbkl_fZ g_e•g•cgbo j•\gygv

X0 – \_dlhj gZqZevgbo gZ[eb`_gv dhj_g•\

Options – \_dlhj hip•c yd• d_jmxlv ijhp_khf h[qbke_ggy lZ \b\h^m j_amevlZl•\ p1, p2,.. - \bdhjbklh\mxlv ydsh Fun(X, p1,p2...) f•klylv [•evr h^gh]h Zj]mf_glZ Options, p1, p2 - g_h[h\'yadh\• iZjZf_ljb

[X, Fval] = fsolve (Fun, X0,...) - ih\_jlZ} dj•f dhj_g•\ agZq_ggy nmgdp•€ m lhqdZo O.

Jha]eyg_fh ^_yd• •l_jZp•cg• f_lh^b j•r_ggy kbkl_f e•g•cgbo j•\gygv

F_lh^ a\bqZcgboLl_jZpLc

>ey aZklhkm\Zggy pvh]h f_lh^Z lj_[Z ihqZldh\m kbkl_fm jL\gygv i_j_l\hjblb ^h \b]ey^m

59

Соседние файлы в папке Учёба