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

Лабы / Лабы 2 / LAB22

.TXT
Скачиваний:
19
Добавлен:
17.04.2013
Размер:
1.66 Кб
Скачать
‚«Ёп­ЁҐ ®ЄагЈ«Ґ­Ё©.
ЏаЁЎ«Ё¦Ґ­­®Ґ ўлзЁб«Ґ­ЁҐ §­ зҐ­Ё© дг­ЄжЁ© Ї® д®а¬г«Ґ ’Ґ©«®а .


1. ‚л§ў вм Є®­бв ­вг eps (¬ иЁ­­л© ­г«м). ‘®§¤ вм Їа®жҐ¤гаг
(㤮Ў­ҐҐ ў ўЁ¤Ґ m-д ©« ):
epsilon=1;
k=0;
while (1+epsilon)>1,
epsilon=epsilon/2;
k=k+1;
end
disp('epsilon='), epsilon
disp('зЁб«® ЁвҐа жЁ© а ў­®'),k

2. €§гзЁвм б Ї®¬®ймо Є®¬ ­¤л help Є®¬ ­¤л ®ЄагЈ«Ґ­Ё© fix, round, ceil, floor
‡ ¤ вм зҐвлаҐе¬Ґа­л© ўҐЄв®а x=[2.36 -2.35 -3.78 4.53] Ё ®ЄагЈ«Ёвм ҐЈ®
гЄ § ­­л¬Ё Є®¬ ­¤ ¬Ё. Џа® ­ «Ё§Ёаг©вҐ १г«мв в.

3. €§ Є в «®Ј  d:/matlab бЄ®ЇЁа®ў вм ў d:/vm1/temp д ©« l22.m
Џа®ўҐаЁвм ҐЈ® а Ў®вг. Љ Є®ў б¬л᫠१г«мв в®ў? (‚лзЁб«Ґ­ЁҐ cos Ё tg ¤«п
Ў®«миЁе §­ зҐ­Ё©  аЈг¬Ґ­в  Ї® д®а¬г«Ґ ’Ґ©«®а )

4. ‚ Є в «®ЈҐ d:/vm1/temp ЇҐаҐЁ¬Ґ­г©вҐ д ©« l22.m, ­ §ў ў ҐЈ® Є ЄЁ¬-«ЁЎ®
¤агЈЁ¬ г­ЁЄ «м­л¬ Ё¬Ґ­Ґ¬ (ў ¤ «м­Ґ©иҐ¬ ‚л Ўг¤ҐвҐ аҐ¤ ЄвЁа®ў вм ᮤҐа¦ ­ЁҐ
Ё ЁбЇ®«м§®ў вм ¤ ­­л© д ©« ­  б«Ґ¤го饬 § ­пвЁЁ, Ї®н⮬㠯® ®Є®­з ­ЁЁ а Ў®вл
­Ґ § Ўг¤м⥠ЇҐаҐЇЁб вм ҐЈ® ­  «Ёз­л© бҐвҐў®© ¤ЁбЄ). €§¬Ґ­ЁвҐ ЇҐаҐЁ¬Ґ­®ў ­­л©
д ©« в Є, зв®Ўл Ї®«гзЁвм Їа®Ја ¬¬г, ўлзЁб«пойго cos(x) Ї® д®а¬г«Ґ ’Ґ©«®а .
‚ Є зҐб⢥ ®Ўа §ж  ­Ё¦Ґ ЇаЁў®¤Ёвбп Їа®Ја ¬¬ , ўлзЁб«пой п §­ зҐ­ЁҐ sin(x)

n=input('n='); % Є®«-ў® з«Ґ­®ў ў д®а¬г«Ґ ’Ґ©«®а 
x=input('x='); %  аЈг¬Ґ­в
a=x; % ⥪гйЁ© з«Ґ­ ап¤ 
y=0; % ⥪гй п ­ Є®Ї«Ґ­­ п б㬬  ап¤ 
c=-x*x; % ўбЇ®¬. ЇҐаҐ¬Ґ­­ п
for k=1:n,
y=y+a; % k- п з бвЁз­ п б㬬  ап¤  ’Ґ©«®а 
a=a*c/((2*k+1)*2*k); % (k+1)-© (б«Ґ¤гойЁ©!) з«Ґ­ ап¤  ’Ґ©«®а 
end
disp('б㬬  n з«Ґ­®ў ап¤  ’Ґ©«®а '),y
err=abs(sin(x)-y);
disp('ЁбвЁ­­ п ®иЁЎЄ '), err
r=abs(a);
disp(' ­ «ЁвЁзҐбЄ п ®жҐ­Є  ®иЁЎЄЁ'),r
Соседние файлы в папке Лабы 2