資源描述：
《實(shí)驗(yàn)九 數(shù)值微積分與方程數(shù)值求解.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、實(shí)驗(yàn)九數(shù)值微積分與方程數(shù)值求解一、實(shí)驗(yàn)?zāi)康?、掌握求數(shù)值導(dǎo)數(shù)和數(shù)值積分的方法2、掌握代數(shù)方程數(shù)值求解的方法3、掌握常微分方程數(shù)值求解的方法二、實(shí)驗(yàn)內(nèi)容1、求函數(shù)在指定點(diǎn)的數(shù)值導(dǎo)數(shù)。symsxy=[x,x.^2,x.^3;1,2*x,3*x.^3;0,2,6*x];dy=diff(y)subs(dy,x,1)subs(dy,x,2)subs(dy,x,3)dy=[1,2*x,3*x^2][0,2,9*x^2][0,0,6]ans=123029006ans=14120236006ans=16270281006>>2、用數(shù)值方法求定積分。（1）的近似值g=inline('sqrt(c
2、os(t.^2)+4*sin((2*t).^2)+1)');I=quadl(g,0,2*pi)I=6.7992+3.1526i（2）g=inline('log(1+x)./(1+x.^2)');I=quadl(g,0,1)I=0.27223、分別用3種不同的數(shù)值方法解線性方程組。第一種（求逆矩陣方法）：A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]';x=inv(A)*bx=0.6667-1.00001.50000.0000第二種（左除運(yùn)算符）：A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9
3、,0,2];b=[-4,13,1,11]';x=Abx=0.6667-1.00001.5000-0.0000第三種（LU分解）：A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]';[L,U]=lu(A);x=U(Lb)x=0.6667-1.00001.5000-0.00004、求非齊次線性方程組的通解。A=[2731;3522;9417]A=273135229417>>b=[642]'b=642>>null(A,'r')ans=0.0909-0.8182-0.45450.09091.0000001.0000>>
4、x0=Abx0=-0.18180.9091005、求代數(shù)方程的數(shù)值解。（1）附近的根functionf=fz(x)f=3*x-sin(x*pi/180)-exp(x);fzero('fz',-1.5)ans=0.6287（2）在給定的初值，求方程組的數(shù)值解functionF=myfun(X)x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y^2+log(z)-7;F(2)=3*x+2^y-3*z+1;F(3)=x+y+z-5;X=fsolve('myfun',[1,1,1],optimset('Display','off'))X=0.21362.43592.
5、35056、求函數(shù)在指定區(qū)間的極值。（1）在（0，1）內(nèi)的最小值。functionf=f(x)f=(x.^3+cos(x)+x.*log(x))./exp(x)[x,fmin]=fminbnd('f',0,1)x=0.5223fmin=0.3974（2）在[0，0]附近的最小值點(diǎn)和最小值。functiong=g(u)x1=u(1);x2=u(2);g=2*x1.^3+4*x1.*x2-10*x1.*x2+x2.^2;[u,gmin]=fminsearch('g',[0,0])u=3.00009.0000gmin=-27.00007、求微分方程的數(shù)值解，并繪制解的曲線。funct
6、iondy=vdp5(x,y)dy=zeros(2,1)dy(1)=y(2)dy(2)=5*y(2)/x-y(1)/x[X,Y]=ode23('vdp5',[0,100],[0,0]);plot(X,Y(:,1),'-')8、求微分方程組的數(shù)值解，并繪制解的曲線。functiondy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);[T,Y]=ode45('rigid',[012],[011]);>>plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(
7、:,3),'+')
