1、溫州大學物理與電子信息工程學院mat i ab仿真及其應用實驗報告課程名稱：matlab仿真及其應用班級：07通信姓名：周偉學號：07110003137實驗地點：5b309日期：2009-11-14實驗名稱：實驗三matlab基本編程操作實驗目的：l熟悉matlab的腳本編寫2、熟悉matlab的函數編寫3、書寫matlab的變量使用4、熟悉matlab的程序控制結構實驗內容：1. 分別選用if或switch結構實現函數表示-1 x<-a/(x) = < -a<x<aa1 x>-ay2xh2. 根據+ x + + . + +近似計算指數，當與指數函數的誤差小于2!

2、n0.01時停止，分別用for和while結構實現。3. 迭代計算兀屮=丄，給出可能的收斂值，并給出不同收斂對應的處置£+2范圍。4在第3題的代碼中增加try和catch控制塊，以避免出現£=-2的情況5、從鍵盤輸入數值，迭代計算耳+】=一實驗結果及分析：1、分別選用if或switch結構實現函數表示-1 x<-ax/(%) = < -a<x<aa1 x>-aedit d«huc r承 top iindow h«lp以 katlabm-futqpgclrl*0ei cureclo” ccmtnd windwyaritblei

3、mport d«tayorktote asctrl*sgw編寫程序3 editor- d:>atlabvorkfunctionl.beile editlex t cell tols dehug resk top window helpd 3 1function y=functionl (x a)2 _if x<=-a3 -y=t;4 _elseif -a<x&&x<al5 -y=x/a;6 _else7 -y=i;8 _endfunction y二function l(x，a)if x<=-aelseif -a<x&&

4、;x<ay=i;end保存后在命令窗口鍵入y=functionl(l,5)運行結果:couand tindot>> jr=fimctionl (lj 5)y =0. 2000»switch 實現：function y=functionl(x,a)switch xcase x<=-ay 二-1;case -a<x&&x<ay=x/a;otherwisey=i;end保存后在命令窗口鍵入y=function 1(1,5)運行結果：» y=functionl( 1,5)y =0.2000r22、根據= l + x + +近似計算

5、指數，當與指數函數的誤差小于 2!n0.01時停止，分別用for和while結構實現。用for結構實現：用for結構：function y=function3(x) y= 1 ;t= 1 ;a=x;m=y+x;for i=l:infif (m-y)>=0.00001 y=y+a/functionl(t); a=a*x;t=t+1; m=y+a/function 1 (t);elsebreak;endendendfunction y=functionl(n)for i=l:l:n y=y*i;endend運行結果：» y=f2(3) y 二20. 0855 >> y=

6、f2 (6) y 二403.4288用while結構實現：function y=functionl(x) y=l;i=l;a=l;while abs(a*(x/(i+1)-1 )>0.0000l|i<10 a=l;for i=l :1 :ia=a*x; a=a/i;endy 二 y+a；i=i+l;end富 editor - untitled*eile edit wet cell tols debug resktop window helpd & s 抽swa 老1 屯尼123456789101112function y=functionl(x)y=l;i=l;a=l；wh

7、ile abs(a*(x/(i+l)-l)>0. 00001|i<10 a=l；for i=l:l:ia=a*x;a=a/i;endy=y+a; i=i+l;end運行結果:>> y=functionl(2)7. 3891»>> y=functionl(4)54.59813、迭代計算£+=一，給出可能的收斂值，并給出不同收斂對應的處置 "+2范圍。function y=function3(x)x0=x/2;xl=3/(x0+2);while abs(xo-x l)>=0.00001xo=xl;xl=3/(x0+2);end

8、4、在第3題的代碼中增加try和catch控制塊，以避免出現£=-2的情況function y=function4(x)x0=x/2;xl=3/(x0+2);while abs(xo-x 1 )>=0.00001xo=x1;xl=3/(x0+2);tryx0=-2;catchdisp('eirof);endend35、從鍵盤輸入數值，迭代計算耳二亠兀+2function function5()x=input('input a number:*);x0=x/2;xl=3/(x0+2);while abs(xo-x 1 )>=0.00001xo=x1;xl=3/(x0+2);xlend運行結果» function5()input