第六章 数学实验.pptVIP

*/16 ? 方程求根方法 ? 函数求极值方法 ? 常微分方程问题 ? 思考题与练习题 代数方程与常微分方程 例6.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 例6.3 球体的吃水深度. 计算半径 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是目标函数的最小值点。 例6.7求一元函数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) x0 = 0.7071 y0 = 0.0711 温室 例6.9 花园靠楼房处有一温室,温室伸入花园 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 ,由牛顿第二定律,得 一阶微分方程主要信息是右端项和初始值: 例7.1 马尔萨斯模型,以1994 年我国人口为12亿为初值，求解常微分方程 N(t)表示人口数量,取人口变化率r =0.015,微分方程 function z=fun1(t,N) z=0.015*N; ode23(fun1,[1994,2020],12) [T,N]=ode23(fun1,[1994,2020],12) 命令窗口 ? 编辑窗口 ? 例2. Logistic模型 创建微分方程右端函数: function z=fun2(t,u) z=