matlab作业题答案.docVIP

作业题 1 求函数在指定点的导数值 syms a x y; a=[x x^2 x^3;1 2*x 3*x^2;0 2 6*x]; y=det(a); dy=diff(y); x=[1 2 3]; eval(dy) ans = 6 24 54 2 符号法求下列函数的导数或积分 y=x10+10x+logx10,求y’ y=sym(x^10+10^x+log(10)/log(x)) Dy=diff(y) y = log(10)/log(x) + 10^x + x^10 Dy = 10^x*log(10) + 10*x^9 - log(10)/(x*log(x)^2) y=ln(1+x), 求y’’∣x=1 syms x y; y=log(1+x); D2y=diff(y,2); x=1; eval(D2y) ans = -0.2500 y=ex/cosx,求y’ syms x y; y=exp(x)/cos(x); Dy=diff(y) Dy = exp(x)/cos(x) + (exp(x)*sin(x))/cos(x)^2 4) f=sym(t*sin(t)); int(f,0,pi) ans = pi 已知函数z=sin(xy), 计算 syms x y z; z=sin(x*y); D3z=diff(diff(z,x),y,2) D3z = - 2*x*sin(x*y) - x^2*y*cos(x*y) 3 用数值方法求定积分 fun=inline(x.^2.*sqrt(2.*x.^2+3)); [f,n]=quad(fun,1,5) f = 232.8057 n = 21 fun=inline(x./sin(x.^2)); [f,n]=quadl(fun,pi/4,pi/3) f = 0.3254 n = 18 4 已知数据[x,y]如下表，试求2次拟合多项式f(x)，然后求x=0.05,0.25,0.45,0.65,0.85,1.05各点的函数近似值，并绘出拟合曲线及求得的函数点。 x 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 y -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2 x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2]; a=polyfit(x,y,2) x1=0.05:0.2:1.05; y1=a(1)*x1.^2+a(2)*x1+a(3); yfit=polyval(a,x); plot(x1,y1,or,x,yfit,-g) a = -9.8108 20.1293 -0.0317 5 求方程exsinx=0在区间[-4,-3]上的一个根。 fzero(exp(x)*sin(x),[-4,-3]) ans = -3.1416 6 已知某函数的离散点列如下 X -4 -3 -2 -1 0 1 2 3 4 y 0.97279 3.4234 3.8186 1.8415 0 -0.15853 -0.18141 -2.5766 -7.0272 分别使用hermite插值、三次样条插值、三次方插值估算在x在2.5，1.5，0.5，-0.5，-1.5，-2.5的函数值，并在同一张图上绘制出三条曲线。 x=[-4 -3 -2 -1 0 1 2 3 4]; y=[ 0.97279 3.4234 3.8186 1.8415 0 -0.15853 -0.18141 -2.5766 -7.0272]; x1=-2.5:1:2.5; c1=interp1(x,y,x1,pehip); c2=interp1(x,y,x1,splin); c3=interp1(x,y,x1,cubic); plot(x1,c1,o-r,x1,c2,*--g,x1,c3,:bH) legend(hermite插值,三次样条插值,三次插值) 7 用经典的RK方法计算常微分方程：y’=e-2t-2y 初值条件y(0)=0.1 fun=inline(exp(-2*t)-2*y,t,y); [t,y]=ode45(fun,[0,1],0.1) t = 0 0.0063 0.0126 0.0188 0.0251 0.0501 0.0751 0.1001