《科學(xué)與工程計(jì)算方法及其應(yīng)用》課程上機(jī)實(shí)驗(yàn)報(bào)告姓名：班級(jí)：自動(dòng)化2班學(xué)號(hào)：20110703217日期：2013/4/23指導(dǎo)老師：劉小燕本次實(shí)驗(yàn)題號(hào)：實(shí)驗(yàn)1程序流程圖牛頓迭代法程序代碼及注釋求和function[sum]=fun(k)s=1;sum=0;forj=1:kfori=1:js=s*i;endsum=sum+s;s=1;end分段函數(shù)functionfx=fun(x)ifx>1fx=xA2+1;elseifx>0&&x<=1fx=2*x;elsefx=xA3;end二分法clearclcsymsRy=inline(exp(-0.005*R)*cos(sqrt(2000-0.01*RA2*0.05))-0.01);a=250;b=650;err=0.000001;while(b-a)>errs=(a+b)/2;if(y(s)*y(a)<0)b=s;elsea=s;endenddisp('k=');disp(k);4.NewtonclearclcsymsRy=inline(exp(-0.005*R)*cos(sqrt(2000-0.01*RA2*0.05))-0.01);dy=inline(diff(sym(y)));k=0;b=650;x1=x-y(x)/dy(x);err=0.000001;disp(‘k=，)disp(k)disp(‘x=，)disp(x)k=k+1;whileabs(abs(x1)-abs(x2))>errdisp(‘k=，)disp(k)disp(‘x=，)disp(x)x=x1;x1=x-y(x)/dy(x);k=k+1;enddisp(‘k=，)disp(k)disp(‘x=，)disp(x)運(yùn)算結(jié)果求階乘求函數(shù)值牛頓迭代法/MATLAB文件?編輯（目調(diào)試受）艇（②窗口姓幫助（旦） DI<X>昭唱G口|梅百「§1當(dāng)前且錄0。I口：'叩亂濟(jì)血呻招""I...|ft7|快捷療式S如何束加便］最近更新福留通版亶相甌福留通版亶相甌▼I堆閡基本二I

宕禰匕i I類型ans 605.'o"o' doubledf <1x1 inline> inlinef <1x1 inline〉 inlinex <1x1 syrn> syrn田ra??>>sjnnsx>>f=inline(exp(-0.005*R)*cos(sqrt(2000-0.01*R'2*0.05))-0.01)???Undefinedfwictionorvariable'R'.>>f=inline(exp(-0.005*x)*cos(sqrt(2000-0.01*x'2*0.05))-0.01)f=Inlinefunction:f(x)=exp(-1./200.*x).*cos(1./100.*(20000000-5.*x.'2).*(l./2))-l./100>>df=inline(diff(sjnn(f)));?niudun(f,df,200,0.000001)Workspace|CommandHistoiy ?xWorkspace|CommandHistoiy ?x—~J ~—clc—diffCfJ)—clc—diff(sym(f))—df=inline(diff(sjnn(f)))—df=inline(；diff ')—df=inline(diff(sjrm(7f?)))—clc—symsk—f=inline(exp(-0.005*R)*cos(sqrt(2000-0.01*R*2*0.05))-0.01)—f=inline(exp(-0.005*x)*cos(sqrt(2000-0.01*x*2*0.05))-0.01)—df=inline(diff(syrn(f))):-^niudun(f,df,200,0.000001)Li i 凸事開(kāi)始冬)I605.8841:ans=605.8841??mows?1.