数学实验七代数方程与常微分方程.pptVIP

? 方程求根方法 ? 函数求极值方法 ? 常微分方程问题 ? 思考题与练习题 代数方程与常微分方程 例1求解 3 次方程 x3 + 1 = 0 。 求多项式根(零点)方法: R= roots(P) P是多项式 P(x) = a1xn + a2 xn-1 + …… + an x + an+1 系数[a1,a2,···,an+1],R为多项式全部零点。 求数值解 P=[1,0,0,1]; R=roots(P) R = -1.0000 0.5000 + 0.8660i 0.5000 - 0.8660i 求符号解 sym x; solve(x^3+1=0) ans = [ -1] [ 1/2-1/2*i*3^(1/2)] [ 1/2+1/2*i*3^(1/2)] 一、方程求根方法 相关命令：poly, polyval, fzero X=1:7;P=poly(X) 1 -28 322 -1960 6769 -13132 13068 -5040 roots(P) ans = 7.0000 6.0000 5.0000 4.0000 3.0000 2.0000 1.0000 x=1:0.1:7; y=polyval(P,x); plot(x,y,X,zeros(1,7),o) 多项式求根方法 p=[1 -30 0 2552]; roots(p) ans = 26.3146 11.8615 -8.1761 r x 例2 球体的吃水深度. 计算半径 r =10 cm的球体,密度 ? =0.638.浸入水深度 x = ? 解：重量 体积 ? x3 –30x2 +2552 =0 求函数零点方法 fun=inline(x.^3-30*x.^2+2552); x=fzero(fun,10) x = 11.8615 求一元函数最小值方法 Xmin=fminbnd(fun，x1，x2) fun是目标函数,[x1，x2]是最小值点搜索区间,Xmin是目标函数的最小值点。 例3求一元函数f(x) = 0.5 – x exp(– x2)在区间[0，2]内的最小值，并绘出函数图形标出最小值点。 fun=inline(0.5-x.*exp(-x.^2)); fplot(fun,[0,2]),hold on [x0,y0]=fminbnd(fun,0,2) plot(x0,y0,o) ;hold off x0 = 0.7071 y0 = 0.0711 二、函数求极值方法 温室 例4 花园靠楼房处有一温室,温室伸入花园 2 米,高3米.温室上方是楼房窗台,要将梯子从花园地上放靠在楼房墙上不损坏温室，用 7 米长的梯子是否可行？ 解:设梯子长度为L, 梯子与地面的 夹角为 数学模型: L=inline(2./cos(alpha)+3./sin(alpha)) [x,Lmin]=fminbnd(L,0.8,0.9) x = 0.8528 Lmin = 7.0235 MATLAB求常微分方程初值问题 数值方法是先创建函数文件,用以描述微分方程右端二元函数,然后用ode23()求出数值解 引例 炮弹在飞行过程中,空气阻力与飞行速度v的平方成正比,如果初始速度v0 ,由牛顿第二定律,得 一阶微分方程主要信息是右端项和初始值: 三、常微分方程问题 常微分方程组初值问题 一阶常微分方程组初值问题数值求解方法 [T,y] = ode23( F ,Tspan,y0) 其中, F是函数文件, 表示 微分方程右端函数 Tspan = [t0 Tfinal] —— 求解区域; y0 —— 初始条件 注: 函数F(t,y) 必须返回列向量. 数值解 y 的每一行对应于列向量T中的每一行数据 例5 马尔萨斯模型,以