"Вычисления, визуализация и программирование в среде MATLAB 5.x" - читать интересную книгу автора (Мартынов Н. Н., Иванов А. П.)

0.7391
' ╧ЁютхЁшь, эрёъюы№ъю сышчъю ъ эєы■ чэрўхэшх ЇєэъЎшш MyFunctionl т Єюўъх
0.7391 (ёь. Ёшё. 4.3).
*> MATLAB Command Window и╚° ┼▌
File Edit__Windo* hjelp
DS1' " ╔ " S
╗ MyFimctionl(0.7391) ж aiis =
-2.4881e-005
-2J________________________________ _______jJs_ ╨шёєэюъ 4.3
┼ёыш эрщфхэю рсёюы■Єэю Єюўэюх чэрўхэшх ъюЁэ , Єю чэрўхэшх ЇєэъЎшш т ¤Єющ Єюўъх Ёртэю эєы■. ╥ръшь юсЁрчюь, тхышўшэр ЇєэъЎшш т яЁшсышцхээю эрщфхээюь эєых ъюётхээю їрЁръЄхЁшчєхЄ яюуЁх°эюёЄ№ Ёхчєы№ЄрЄр. ╫Єюс√ єяЁрты Є№ яюуЁх°эюёЄ№■, эєцэю юёє∙хёЄты Є№ т√чют ЇєэъЎшш fzero ё ЄЁхь  рЁуєьхэЄрьш:
fzero ( name, xO, tol )
уфх ярЁрьхЄЁ tol чрфрхЄ ЄЁхсєхьє■ тхышўшэє яюуЁх°эюёЄш (ю°шсъш). ╧ютЄюЁшь яЁхф√фє∙шх т√ўшёыхэш , яюЄЁхсютрт сюы№°хщ ЄюўэюёЄш ЁрёўхЄют (Єю хёЄ№ ьхэ№°хщ яюуЁх°эюёЄш):
format long
ї = fzero('MyFunctionl',pi/2, le-8)
x =
0.73908513263090 MyFunctionl(x) ans = 9.778322596076.805e-010
юЄъєфр тшфэю, ўЄю фхщёЄтшЄхы№эю фюёЄшуэєЄр сюы№°р  ЄюўэюёЄ№ эрїюцфхэш  эєы  ЇєэъЎшш.
┼∙х Ёрч яюфўхЁъэхь, ўЄю ЇєэъЎш  fzero эрїюфшЄ эєыш Єюы№ъю тх∙хёЄтхээю-чэрўэ√ї ЇєэъЎшщ юфэющ тх∙хёЄтхээющ яхЁхьхээющ. ╬фэръю ўрёЄю с√трхЄ эхюсїюфшью эрщЄш ъюьяыхъёэ√х ъюЁэш тх∙хёЄтхээючэрўэ√ї ЇєэъЎшщ, юёюсхээю т ёыєўрх ьэюуюўыхэют.
MATLAB 5.x. ┬√ўшёыхэш , тшчєрышчрЎш , яЁюуЁрььшЁютрэшх
135
─ы  ¤Єющ Ўхыш т ёшёЄхьх MATLAB ёє∙хёЄтєхЄ ёяхЎшры№эр  ЇєэъЎш  roots, ъюЄюЁющ т ърўхёЄтх рЁуєьхэЄр яхЁхфрхЄё  ьрёёшт ъю¤ЇЇшЎшхэЄют ьэюуюўыхэр. ═ряЁшьхЁ, фы  ьэюуюўыхэр ї└4- ╟їы3 + 3їы2- ╟ї + 2, шьх■∙хую фтр тх∙хёЄтхээ√ї (1 ш 2) ш фтр ъюьяыхъёэ√ї ъюЁэ  (i ш -i), эєцэю ёэрўрыр ёЇюЁьшЁютрЄ№ ьрёёшт хую ъю¤ЇЇшЎшхэЄют:
Coef = [ 1, -3, 3, -3, 2 ] яюёых ўхую т√чтрЄ№ ЇєэъЎш■ roots:
у = roots( Coef )
у =
2.00000000000000
0.00000000000000 + l.OOOOOOOOOOOOOOi
0.00000000000000 - l.OOOOOOOOOOOOOOi
1.00000000000000
┬ чрфрўх ю эрїюцфхэшш эєыхщ ЇєэъЎшш ёыюцэ√ь ьюьхэЄюь  ты хЄё  эрїюцфхэшх эрўры№эюую яЁшсышцхэш  ъ эєы■ ЇєэъЎшш, р Єръцх ряЁшюЁэр  юЎхэър шї ъюышўхёЄтр. ╧ю¤Єюьє трцэю ярЁрыыхы№эю ё яЁшьхэхэшхь ЇєэъЎшщ Єшяр roots шыш fzero тшчєрышчшЁютрЄ№ яютхфхэшх шёъюь√ї ЇєэъЎшщ эр Єюь шыш шэюь юЄЁхчъх чэрўхэшщ рЁуєьхэЄр. ╠ръёшьры№эюх ёюфхщёЄтшх т ¤Єюь ьюцхЄ юърчрЄ№ ЇєэъЎш  fplot (name, [x0, xl]), ёЄЁю ∙р  уЁрЇшъ ЇєэъЎшш ё шьхэхь name эр юЄЁхчъх юЄ ї╬ фю xl.
═ряЁшьхЁ, фы  Ёрэхх ЁрёёьюЄЁхээющ ЇєэъЎшш MyFunctionl ё яюью∙№■ т√чютр ЇєэъЎшш
fplot( 'MyFunctionl', [0,pi/2] )-
єсхцфрхьё , ўЄю эр юЄЁхчъх юЄ 0 фю pi/2 фхщёЄтшЄхы№эю ёє∙хёЄтєхЄ хфшэёЄтхээ√щ эєы№ ¤Єющ ЇєэъЎшш (ёь. Ёшё. 4.4).
* Figure No. 1
Eile Edit Window Jdelp 1
╨шёєэюъ 4.4
╘єэъЎш  fplot ЄЁхсєхЄ юЄ яюы№чютрЄхы  х∙х ьхэ№°х яюфуюЄютшЄхы№эющ ЁрсюЄ√, ўхь фрцх шчєўхээр  эрьш Ёрэхх ЇєэъЎш  plot (ёь. уы. 2). .
136
├ыртр 4. ╨х°хэшх Єшяют√ї чрфрў рыухсЁ√ ш рэрышчр
╧юшёъ ьшэшьєьр ЇєэъЎшш
┬ ёшёЄхьх MATLAB шьх■Єё  ёяхЎшры№э√х ЇєэъЎшш фы  яюшёър ьшэшьєьют чрфрээ√ї ЇєэъЎшщ. ╧Ёш ¤Єюь тючьюцхэ яюшёъ ьшэшьєьр ъръ фы  ЇєэъЎшш юфэющ тх∙хёЄтхээющ яхЁхьхээющ, Єръ ш фы  ЇєэъЎшщ ьэюушї яхЁхьхээ√ї.
─ы  ЇєэъЎшщ юфэющ яхЁхьхээющ шї ьшэшьєь√ Ёрч√ёъштрхЄ ЇєэъЎш  f min:
fmin( name, xO, xl )
╟фхё№ name яЁхфёЄрты хЄ шь  ЇєэъЎшш, є ъюЄюЁющ эрїюф Єё  ьшэшьєь√, р ї╬ ш xl чрфр■Є юЄЁхчюъ яюшёър. ╚ыы■ёЄЁшЁютрЄ№ ЁрсюЄє ¤Єющ ЇєэъЎшш сєфхь эр яЁшьхЁх ЇєэъЎшш hump (яхЁхтюфшЄё  ъръ луюЁс╗), ёяхЎшры№эю яюёЄрты ■∙хщё  ё ёшёЄхьющ MATLAB т фхьюэёЄЁрЎшюээ√ї Ўхы ї. ▌Єр ЇєэъЎш  чрфрхЄё  ЇюЁьєыющ
є = 1./((ї - 0.3).╦2 + 0.01)+1./((ї - 0.9).ы2 + 0.04) - с; р хх уЁрЇшъ ыхуъю яюыєўрхЄё  ё яюью∙№■ т√чютр ЇєэъЎшш
fplot( 'humps', [0,3] )
╚ч Ёшёєэър 4.5 тшфэю, ўЄю ыюъры№э√щ ьшэшьєь ¤Єющ ЇєэъЎшш ёє∙хёЄтєхЄ эр юЄЁхчъх юЄ 0.5 фю 1.0. ╧юяЁюсєхь эрщЄш Єюўъє ьшэшьєьр ёыхфє■∙шь т√чютюь ЇєэъЎшш fmin:
ї = fmin( 'humps', 0.5, 1.0 ) ї =
0.63701067459059
╨шёєэюъ 4.5
MATLAB 5.x. ┬√ўшёыхэш , тшчєрышчршш , яЁюуЁрььшЁютрэшх 137