matlab 实验六1数值微积分与常微分方程求解.doc

实验六 数值微积分与常微分方程求解 一、实验目的 1. 掌握微分与积分的数值计算方法 2. 掌握常微分方程的数值求解方法 二、实验内容和步骤 1. 求函数指定点的数值导数 2. 求下列函数的导数 （1），求 （2）求 3. 用数值方法求定积分 （1） （2） (3) (4) 4. 求微分方程的数值解 （1） （2） 三、程序清单及运行结果 1. 程序清单： for x=1:1:3 fx=[x,x^2,x^3;1,2*x,3*x^2;0,2,6*x]; diff(fx) end 运行结果： 2.（1）程序清单： x=0:0.05:1; y=x.^10+10.^x+1./log10(x); dy=diff(y)/0.05 运行结果： （2）程序清单： x=0:0.01:1; y=log(1+x); dy=diff(y,2)/0.01; dy(99) 运行结果： 3.(1)程序清单： fx=inline(x.^2.*sqrt(2*x.*x+3)); quad(fx,1,5) 运行结果： (2)程序清单： fx=inline(x./sin(x).^2); quad(fx,pi/4,pi/3) 运行结果： （3）程序清单： fx=inline(abs(cos(x+y))); dblquad(fx,0,pi,0,pi) 运行结果： （4）程序清单： syms x y; fx=x*y; int(int(fx,y^2,y+2),-1,2); 运行结果： 4.（1）程序清单： 建立函数文件funst.m。 function yp=funst(t,y) yp=-(1.2+sin(10*t))*y; 求解微分方程。 t0=0; tf=5; y0=1; [t,y]=ode23(@funst,[t0,tf],y0) 运行结果： t = 0 0.0667 0.1375 0.2003 0.2695 0.3528 0.4362 0.5033 0.5663 0.6369 0.6913 0.7457 0.8081 0.8738 0.9591 1.0277 1.0963 1.1600 1.2246 1.3082 1.3714 1.4347 1.5001 1.5842 1.6530 1.7219 1.7858 1.8501 1.9319 1.9953 2.0587 2.1236 2.2053 2.2745 2.3438 2.4080 2.4719 2.5501 2.6140 2.6779 2.7419 2.8193 2.9046 2.9685 3.0323 3.0959 3.1721 3.2364 3.3007 3.3642 3.4395 3.5328 3.5965 3.6602 3.7238 3.7998 3.8642 3.9286 3.9920 4.0671 4.1408 4.2144 4.2800 4.3432 4.4159 4.4812 4.5465 4.6094 4.6812 4.7567 4.8322 4.8990 4.9620 5.0000  y = 1.0000 0.9035 0.7822 0.6823 0.5984 0.5402 0.5182 0.5106 0.4976 0.4656 0.4280 0.3844 0.3355 0.2936 0.2595 0.2469 0.2421 0.2382 0.2290 0.2054 0.1819 0.1585 0.1387 0.1225 0.1164 0.1140 0.1122 0.1080