Добавил:
korayakov
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:
‚«ЁпЁҐ ®ЄагЈ«ҐЁ©.
ЏаЁЎ«Ё¦Ґ®Ґ ўлзЁб«ҐЁҐ § 票© дгЄжЁ© Ї® д®а¬г«Ґ ’Ґ©«®а .
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
ЏаЁЎ«Ё¦Ґ®Ґ ўлзЁб«ҐЁҐ § 票© дгЄжЁ© Ї® д®а¬г«Ґ ’Ґ©«®а .
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