导出微分方程.pptVIP

一个引理 设集合E0={g(x)?C1 │g(a) =g(b)=0} 如果在[a,b]连续函数 f(x)满足 那么 f (x) ? 0 对?g (x) ?E0 ，总有 另一种方法－变分法 A B x y c 设曲线为 满足 y(0)=0, y(c)=H 在曲线上P(x,y）处质点速度为 又设从A到P的弧长为s,则 从而质点沿曲线由A到B需时间 那么我们的问题成为 求某个 使得 引进集合 显然若 是最速曲线函数,则 于是函数 在 取得最小值 故得 设集合 那么对 依复合函数求导法 注意第二项 为了计算 ,记 上式乘以 这里 于是导出 注意从降线定义可知 故 1）可求解析解 解法 2）也可以用数值方法，例如欧拉法求解 得到方程为 由于在原点y = 0 ,可改写方程 ? 求解析解提示： 解析解 function cycloid(G,H,n) if nargin==2 %两个参数则默认n为100 n=100; end g=9.8; h=H/n; minc=0;maxc=1/sqrt(2*g*h*n); x=0;y=0; while abs(G-x)1e-4 x=0; c=(minc+maxc)/2; %二分法求c值 for j=1:n y=j*h; v=sqrt(2*g*y); x=x+c*v*