function?t_charpt1_2
%?數值實驗1.2:誤差傳播與算法穩定性
%?輸入:遞推式選擇及遞推步數
%?輸出：各步遞推值及誤差結果，以及遞推值和誤差與遞推步數的關系圖
clc
promps={‘請選擇遞推關系式，若選E.1.7請輸入1，否則輸入2:‘};
result=inputdlg（promps‘charpt?1_2‘1{‘1‘}）;
Nb?=?str2num（char（result））;
if（（Nb~=1）&（Nb~=2））?errordlg（‘請選擇遞推關系式，若選E.1.7請輸入1，否則輸入2!‘）;return;end

result=inputdlg（{‘請輸入遞推步數n:‘}‘charpt?1_2‘1{‘10‘}）;
steps=str2num（char（result））;
if（steps<1）errordlg（‘遞推步數錯誤!‘）;return;end
result=zeros（1steps）;?err=result;
if（Nb==1）
????n=1;?result（n）=1/exp（1）;
????while（n????????result（n+1）=1-n*result（n）;
????????err（n+1）=abs（result（n+1）-func（n+1））;
????????n=n+1;
????end
elseif（Nb==2）
????n=steps;?
????err（n）=abs（result（n）-func（n））;
????while（n>1）
????????result（n-1）=（1-result（n））/n;
????????t=func（n-1）;?
????????err（n-1）=abs（result（n-1）-func（n-1））;
????????n=n-1;
????end
end
disp（[‘遞推值:‘num2str（result）]）;
disp（[‘誤差:‘num2str（err）]）;
plot（[1:steps]result‘-‘）;
grid?on
hold?on;
plot（[1:steps]err‘r--‘）;
xlabel（‘n‘）;?ylabel（‘En-?and?ERRn--‘）;
text（2err（2）‘\uparrow?err（n）‘）;
text（4result（4）‘\downarrow?En‘）;

%--------------------------------------------------------------------------
function?en=func（n）
%?計算En的精確值
if（n==1）
????en=1/exp（1）;
else
????en=1-n*func（n-1）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????270??2010-05-11?16:01??CH2_CarDoorDesign.m

?????文件???????1842??2010-05-10?17:14??CH2_Interpolation.m

?????文件?????438784??2010-05-27?17:12??數值實驗題.doc

?????文件???????2265??2010-05-11?10:06??CH3_CurveFitting.m

?????文件???????1385??2010-05-11?14:58??CH1_ErrorSpread.m

?????文件????????595??2010-05-11?11:32??CH1_illPosedProblem.m

-----------?---------??----------?-----??----

???????????????445141????????????????????6