ly6常微分,概率导论.ppt

* * * 常微分方程的求解 统计作图 随机试验 设微分方程初值问题： 常微分方程的求解 分为：符号解法和数值解法 符号解法： 命令形式1： dsolve(’eqution’,’var’) 命令形式2： dsolve(’eqution’ , ’cond1,cond2,…’ , ’var’) P158 初始条件 符号解法 dsolve(Dy=y^2,x) 自变量 大写 ?dsolve(Dy=y^2,x) ans = -1/(x-C1) 求解 ?dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) 大写 对应求导阶数 condition ans =-1/3*x^3+125/468+31/468*x^4 “常微分方程初值问题数值解”的提法： 常 微 分 方 程 的 数 值 解 法 假设(1)式的解 存在且唯一 ……..(1) 不求解析解 （无解析解或求解困难） 而在一系列离散点 求 的近似值 通常取等步长h xn=x0+nh 常用的数值方法： 欧拉方法 梯形公式 龙格库塔方法 欧拉方法 基本思路 对(1)在区间[xi,xi+1]上积分 将（*）右边定积分应用于左矩形公式 向前欧拉公式 y1 y2 y3 ……..(1) ……..(*) 将（*）右边定积分应用于不同公式 各种数值方法 欧拉方法 向后欧拉公式 右端点yi+1未知，需迭代求解。常用的方法是：先用向前欧拉法提供初值，然后再利用向后欧拉法迭代。计算公式为： 将（*）右边定积分应用于右矩形公式 梯形公式 将（*）右边定积分应用于梯形公式 龙格库塔方法 向前，向后欧拉公式 用[xi,xi+1]内1个点的导数代替 梯形公式 用[xi,xi+1]内2个点的导数平均值代替 龙格库塔方法的基本思想 在[xi,xi+1]内多取几个点，将它们的导数加权平均代替 ，构造出精度更高的计算公式。 2阶龙格—库塔公式 其中 具有2阶精度。 4阶龙格—库塔公式 其计算精度为4阶。 龙格库塔方法的Matlab实现： [t,y]=ode23(‘fun’,tspan,y0) 其中，fun是定义函数的文件名，该函数fun必须以dy输出量，以t,y为输入量。tspan=[t0 tf]表示积分的起始值和终止值；y0是初始状态列向量。 [t,y]=ode45(‘fun’,tspan,y0) 常 微 分 方 程 的 数 值 解 法 例7-45：用数值积分的方法求解微分方程： 设初始时间 ；终止时间 初始条件 分析：求解 令： （化为一阶微分方程）即原微分方程化为： 写成矩阵形式为 （化为一阶微分方程） 放入函数exf.m中 u function xdot=exf(t,x) u=1-(t.^2)/(2*pi); xdot=[0,1;-1,0]*x+[0 1]*u; 编写函数文件 exf.m clf; t0=0; tf=3*pi; x0t=[0;0]; 定义另外一个函数为主函数 初始和终止时间 初始条件 [t,x]=ode23(exf,[t0,tf],x0t) exf为已定义的子函数 该函数用来记录一阶微分方程 xdot= y=x(:,1); %[t,x]中求出的x是按列排列， 故用ode23求出x后 只要第 一列即为y y2=-1/2*(-2*pi-2+t.^2)/pi-(pi+1)/pi*cos(t); clf, plot(t,y,-, t,y2,o) 解析解为： dsolve(D2y+y=1-t^2 /(2*pi),y(0)=0,Dy(0)=0,t) ans = -1/2*(-2*pi-2+t^2)/pi-(pi+1)/pi*cos(t) 随机试验 古典概率： 事件A发生的概率 在100人的团体中，如果不考虑年龄的差异，研究是否有 两个以上的人生日相同。假设每人的生日在一年365天中 的任意一天是等可能的，那么随机找 个人（365） 问这些人生日各不相同的概率是多少？ 至少有两个人生日相同的概率为多少？ P169 for n=2:100 p0(n)=prod(365:-1:365-n+1)/365^n; p1(n)=1-p0(n); end n=2:100; plot(n,p0,n,p1,--) xlabel(‘人数’),ylabel(‘概率’) legend(‘生日各不相同的概率’,‘至少两人相同的概率’) axis([0 100 -0.1 1.1]),grid on 统计作图 在数据较小、较少的情况下输入 Matlab交互环境 境