四MATLAB求解数学问题市公开课一等奖百校联赛特等奖课件.pptx

（四）MATLAB求解数学问题v数学分析v概率论与数理统计v插值与拟合v优化运算v线性代数v复变函数 10/10/0黄建华制作第1页4.1数学分析v符号方程求解v极限v导数与微分v(重)积分v曲线积分与曲面积分 v空间解析几何与向量代数v级数v微分方程10/10/1黄建华制作第2页4.1.1符号方程求解 主要内容主要内容v 线性方程线性方程v 非线性方程非线性方程10/10/2黄建华制作第3页4.1.1符号方程求解v线性方程线性方程 惯用惯用solve()和和linsolve()函数来处理线性方程问题。函数来处理线性方程问题。详细格式：详细格式：X=solve(方程方程1,.方程方程n,变量变量1,.变量变量m)说明：能够求解方程组，单变量时变量申明能够省略。说明：能够求解方程组，单变量时变量申明能够省略。X=linsolve(A,B)%求解线性方程组求解线性方程组AX=B，返回特解，返回特解X 10/10/3黄建华制作第4页4.1.1符号方程求解v例例4.1.1求解方程：求解方程：命令：命令：clear;x=solve(x2-x-6=0)相当于：clear;syms x;f=x2-x-6;x=solve(f)x=3 -210/10/4黄建华制作第5页4.1.1符号方程求解v例例4.1.2 求解方程组：求解方程组：命令：命令：clear;x,y=solve(x2+y-6=0,y2+x-6=0,x,y)x=2 y=2 -3 -3 1/2-1/2*21(1/2)1/2+1/2*21(1/2)1/2+1/2*21(1/2)1/2-1/2*21(1/2)若将x,y=改用X=，则仅将返回一个解结构。X=x:4x1 sym y:4x1 sym10/10/5黄建华制作第6页4.1.1符号方程求解v例例4.1.3 求解方程组：求解方程组：命令：命令：clear;A=5,0,4,2;1,-1,2,1;4,1,2,0;1,1,1,1;b=3;1;1;0;X=linsolve(A,b)X=1.0000 -1.0000 -1.0000 1.000010/10/6黄建华制作第7页4.1.1符号方程求解v非线性方程非线性方程 惯用惯用fsolve()函数来处理非线性方程问题。函数来处理非线性方程问题。详细格式：详细格式：X=fsolve(fun,x0,options)说明：说明：fun为定义好非线性方程为定义好非线性方程(组组)文件名，其中文件名，其中为调用函数为调用函数符号；符号；x0为求解方程初始向量；为求解方程初始向量；options设置求解过程各种参设置求解过程各种参数，普通采取默认参数数，普通采取默认参数optimset(fsolve),其它参数能够查询其它参数能够查询帮助。帮助。10/10/7黄建华制作第8页4.1.1符号方程求解v例例4.1.4 求解方程组：求解方程组：x0=x(1),x(2)=0.1,0.110/10/8黄建华制作第9页4.1.1符号方程求解首先建立函数文件首先建立函数文件fun.m并保留在默认路径下：并保留在默认路径下：function y=fun(x)y=x(1)-0.5*sin(x(1)-0.3*cos(x(2),.x(2)-0.5*cos(x(1)+0.3*sin(x(2);然后运行命令：然后运行命令：clear;x0=0.1,0.1;x=fsolve(fun,x0,optimset(fsolve)x=0.5414 0.331010/10/9黄建华制作第10页4.1.2极限 主要内容主要内容v单变量函数极限单变量函数极限v多变量函数极限多变量函数极限10/10/10黄建华制作第11页4.1.2极限v单变量函数极限单变量函数极限首先进行首先进行符号变量说明符号变量说明:syms x y t h a 然后定义函数然后定义函数fun，再使用以下命令格式求对应极限，再使用以下命令格式求对应极限:limit(fun,x,a)%求函数求函数fun当当xa时极限时极限limit(fun,a)%默认变量默认变量x或唯一符号变量或唯一符号变量limit(fun)%默认变量默认变量x,且且a=0limit(fun,x,a,right)%右极限右极限 xa+limit(fun,x,a,left)%左极限左极限 xa-10/10/11黄建华制作第12页4.1.2极限v例例4.1.5 举例举例:结果结果 syms x h a f=sin(x)/x;limit(f)1 limit(f,inf)0 limit(x-2)/(x2-4),2)1/4 limit(1/x,x,0,right)inf limit(1/x,x,0,left)-inf limit(sin(x+h)-sin(x)/h,h,0)cos(x)limit(1+a/x)*sin(x),x,a)2*sin(a)10/10/12黄建华制作第13页4.1.2极限v例例4.1.6试求解极限问题试求解极限问题:解：解：syms x a b;f=x*(1+a/x)x*sin(b/x);L=limit(f,x,inf)L=exp(a)*b10/10/13黄建华制作第14页4.1.2极限v多变量函数极限多变量函数极限 假设有二元函数求极限问题：假设有二元函数求极限问题：则能够嵌套使用则能够嵌套使用limit()函数函数:limit(limit（fun,x,x0),y,y0)或或 limit(limit（fun,y,y0),x,x0)假如假如x0或或y0不是确定值，而是另一个变量函数，则次序不是确定值，而是另一个变量函数，则次序不能交换。不能交换。注意：此种使用方法只适合用于极限存在情况。注意：此种使用方法只适合用于极限存在情况。10/10/14黄建华制作第15页4.1.2极限v例例4.1.7 求出二元函数极限值求出二元函数极限值:解：解：syms x y;f=sin(x*y)/x;L=limit(limit(f,x,0),y,0)L1=limit(limit(f,y,0),x,0)L=0 L1=010/10/15黄建华制作第16页4.1.2极限v例例4.1.8 求出二元函数极限值求出二元函数极限值:解解 syms x y k;f=(x2-y2)/(x2+y2);(1)L=limit(limit(f,x,0),y,0)(2)L1=limit(limit(f,y,0),x,0)(3)L2=limit(limit(f,y,k*x),x,0)L=-1 L1=1 L2=-(-1+k2)/(1+k2)10/10/16黄建华制作第17页4.1.3导数和微分 主要内容主要内容v导数和高阶导数导数和高阶导数v高阶混合偏导数高阶混合偏导数v复合函数求导复合函数求导v隐函数求偏导隐函数求偏导v参数方程求导参数方程求导v导数应用导数应用v梯度计算和方向导数梯度计算和方向导数10/10/17黄建华制作第18页4.1.3导数和微分v导数和高阶导数导数和高阶导数首先进行首先进行符号变量说明符号变量说明:syms 然后定义函数然后定义函数f，再使用，再使用以下命令格式求对应导数以下命令格式求对应导数:diff(f)%f对默认变量对默认变量x求一阶导数求一阶导数diff(f,v)%f 对变量对变量v求一阶导数求一阶导数diff(f,n)%f对默认变量对默认变量x 求求n阶导数阶导数diff(f,v,n)%f 对变量对变量v 求求n阶导数阶导数显然，用以上命令能够实现求各阶纯偏导。显然，用以上命令能够实现求各阶纯偏导。微分与导数是共通，只须将求导答案变形一下即可！微分与导数是共通，只须将求导答案变形一下即可！10/10/18黄建华制作第19页4.1.3导数和微分v例例4.1.9 命令：命令：syms a x y;f=x3+y2+a*x*y-3*x+7;f1=diff(f)结果：结果：f1=3*x2+a*y-3 f2=diff(f,y)f2=2*y+a*x f3=diff(f,2)f3=6*x f4=diff(f,y,2)f4=2 10/10/19黄建华制作第20页4.1.3导数和微分v高阶混合偏导数高阶混合偏导数假设有多元函数求偏导问题：假设有多元函数求偏导问题：则能够嵌套使用则能够嵌套使用diff()函数函数:命令格式命令格式:diff(difft(f,x,m),y,n)或或 diff(diff(f,y,n),x,m)10/10/20黄建华制作第21页4.1.3导数和微分v例例4.1.10 求：求：命令：命令：syms x y;z=(x2-2*x)*exp(-x2-y2-x*y);zxzy=diff(diff(z,x),y)zxzy=(2*x-2)*(-2*y-x)*exp(-x2-y2-x*y)-(x2-2*x)*exp(-x2-y2-x*y)+(x2-2*x)*(-2*x-y)*(-2*y-x)*exp(-x2-y2-x*y)10/10/21黄建华制作第22页4.1.3导数和微分化简一下：化简一下：命令：命令：zxzy1=simple(zxzy)zxzy1=exp(-x2-y2-x*y)*(-4*x*y-3*x2+4*y+4*x +5*x3*y+2*x4+2*x2*y2-10*x2*y-4*x3-4*x*y2)10/10/22黄建华制作第23页4.1.3导数和微分v例例4.1.11 求：求：命令：命令：syms x y;f=(x2-2*x)*exp(-x2-y2-x*y);ydx=-diff(f,x)/diff(f,y)ydx=(-(2*x-2)*exp(-x2-y2-x*y)-(x2-2*x)*(-2*x-y)*exp(-x2-y2-x*y)/(x2-2*x)/(-2*y-x)/exp(-x2-y2-x*y)10/10/23黄建华制作第24页4.1.3导数和微分v复合函数求导复合函数求导v例例4.1.12 已知：已知：求：求：命令：命令：syms t x y;t=exp(sin(x);y=sin(exp(t);ydx=diff(y,x)ydx=cos(exp(exp(sin(x)*cos(x)*exp(sin(x)*exp(exp(sin(x)注意不能将注意不能将t=exp(sin(x);y=sin(exp(t);输入次序颠倒，输入次序颠倒，不然结果不然结果0。10/10/24黄建华制作第25页4.1.3导数和微分v隐函数求偏导隐函数求偏导假设有隐函数表示式假设有隐函数表示式f(x1,x2,xn)=0,求偏导问题：求偏导问题：能够使用能够使用diff()函数函数,命令格式命令格式:F=-diff(f,xj)/diff(f,xi)10/10/25黄建华制作第26页4.1.3导数和微分v例例4.1.13 已知已知 求：求：命令：命令：clear;syms x y;F=atan(y/x)-log(sqrt(x2+y2);ydx=-diff(F,y)/diff(F,x)ydx=(-1/x/(1+y2/x2)+1/(x2+y2)*y)/(-y/x2/(1+y2/x2)-1/(x2+y2)*x)simple(ydx)ydx=(x-y)/(y+x)10/10/26黄建华制作第27页4.1.3导数和微分v参数方程求导参数方程求导假设有参数方程表示式假设有参数方程表示式y=f(t),x=g(t),求导数：求导数：能够使用能够使用diff()函数递归调用函数递归调用,命令格式命令格式:dk=diff(dk-1,t)/diff(x,t)其中其中dk-1表示表示k-1阶导数阶导数注意不能用：注意不能用：dk=diff(y,t,k)/diff(x,t,k)10/10/27黄建华制作第28页4.1.3导数和微分v例例4.1.14 已知已知 求：求：命令：命令：syms a b t x y;x=a*cos(t);y=b*sin(t);d1=diff(y,t)/diff(x,t)d1=-b*cos(t)/a/sin(t)pretty(simple(d1)b -a tan(t)10/10/28黄建华制作第29页4.1.3导数和微分d2=diff(d1,t)/diff(x,t)d2=-(b/a+b*cos(t)2/a/sin(t)2)/a/sin(t)pretty(simple(d2)b -2 2 a (-1+cos (t)sin(t)10/10/29黄建华制作第30页4.1.3导数和微分若改用：若改用：xd2=diff(y,t,2)/diff(x,t,2)xd2=b*sin(t)/a/cos(t)pretty(simple(xd2)b tan(t)-a显然不一样。显然不一样。10/10/30黄建华制作第31页4.1.3导数和微分v导数应用导数应用v例例4.1.15 讨论函数讨论函数 极值、单调性和其极值、单调性和其导数函数关系。导数函数关系。命令：命令：clear;syms x y dy d2y;y=x2/(1+x2);dy=simple(diff(y);pretty(dy)x 2 -2 2 (1+x )10/10/31黄建华制作第32页4.1.3导数和微分Px=solve(dy)Px=0 d2y=simplify(diff(y,2)d2y=-2*(-1+3*x2)/(1+x2)3 P2x=solve(d2y)P2x=-1/3*3(1/2)1/3*3(1/2)10/10/32黄建华制作第33页4.1.3导数和微分lims=-5,5;subplot(3,1,1);ezplot(y,lims);hold on;line(0,0,-0.5,1.5),line(-sqrt(3)/3,-sqrt(3)/3,-0.5,1.5);line(sqrt(3)/3,sqrt(3)/3,-0.5,1.5);subplot(3,1,2);ezplot(dy,lims);hold on;line(0,0,-1,1.5);line(-5,5,0,0);%同时绘制横轴同时绘制横轴subplot(3,1,3);ezplot(d2y,-5,5);hold on;line(-5,5,0,0);line(-sqrt(3)/3,-sqrt(3)/3,-1,1.5);line(sqrt(3)/3,sqrt(3)/3,-1,1.5);10/10/33黄建华制作第34页4.1.3导数和微分10/10/34黄建华制作第35页4.1.3导数和微分v梯度计算和方向导数梯度计算和方向导数jacobian(fun,v)%v是求导变量向量，表示是求导变量向量，表示fun对对 v求偏导矩阵即梯度求偏导矩阵即梯度gridient(F)%求求F数值梯度，一维时可用数值梯度，一维时可用diff代替代替dot（jacobian(fun)，v）=jacobian(fun)v%v是某方向单位向量，数量积就是方向导数是某方向单位向量，数量积就是方向导数10/10/35黄建华制作第36页4.1.4积分 主要内容主要内容v不定积分不定积分v定积分与无穷积分定积分与无穷积分v重积分重积分v数值积分数值积分10/10/36黄建华制作第37页4.1.4积分v不定积分不定积分 不定积分问题：不定积分问题：能够使用能够使用int()函数函数:F=int(fun,x)或或 F=int(fun)%当当fun中只有一个自变量中只有一个自变量x，则则x可省可省 最终答案应为：最终答案应为：F(x)+C10/10/37黄建华制作第38页4.1.4积分v例例4.1.16 用用diff()diff()函数求函数求4 4阶导数，再积分，阶导数，再积分，检验是否能够得出一致结果。检验是否能够得出一致结果。命令：命令：syms x;y=sin(x)/(x2+4*x+3);y4=diff(y,4);y0=int(int(int(int(y4);pretty(simple(y0)%对导数积分应该得出原函数对导数积分应该得出原函数 sin(x)-(x+1)(x+3)10/10/38黄建华制作第39页4.1.4积分v例例4.1.17证实：证实：命令：命令：syms a x;f=simple(int(x3*cos(a*x)2,x);f1=x4/8+(x3/(4*a)-3*x/(8*a3)*sin(2*a*x)+.(3*x2/(8*a2)-3/(16*a4)*cos(2*a*x);simple(f-f1)%求两个结果差 -3/16/a4 结果是一个常数，表明答案正确。结果是一个常数，表明答案正确。10/10/39黄建华制作第40页4.1.4积分v例例4.1.18不可积问题：不可积问题：命令：命令：syms x;int(exp(x2)vans=-1/2*i*pi(1/2)*erf(i*x)结果中结果中erf是一个定义函数：表明不可积！10/10/40黄建华制作第41页4.1.4积分v定积分与无穷积分定积分与无穷积分 定积分问题：定积分问题：能够使用能够使用int()函数函数:int(fun,x,a,b)若为无穷积分问题，则只需将命令中若为无穷积分问题，则只需将命令中a(或或b)改为改为-inf(或或inf)即可。如求：即可。如求：用：用：int(fun,x,a,inf)10/10/41黄建华制作第42页4.1.4积分v例例4.1.19不可积问题不可积问题 定积分定积分 可积。可积。命令：命令：syms x;int(exp(-x2/2),x,0,inf)ans=1/2*2(1/2)*pi(1/2)命令：命令：syms x;I=int(exp(-x2/2),x,0,1)I=1/2*erf(1/2*2(1/2)*2(1/2)*pi(1/2)能够使用vpa()函数显示数值：vpa(I,5)=0.8556510/10/42黄建华制作第43页4.1.4积分v例例4.1.20变限积分也可使用定积分求解：变限积分也可使用定积分求解：命令：命令：syms x t;F=int(exp(t),t,2*x,sin(x)F=exp(sin(x)-exp(2*x)命令：命令：Fx=diff(F,x)Fx=cos(x)*exp(sin(x)-2*exp(2*x)10/10/43黄建华制作第44页4.1.4积分v重积分重积分 重积分问题能够先化为累次积分方式再使用重积分问题能够先化为累次积分方式再使用 int()函数嵌套来处理函数嵌套来处理:v例例4.1.21求二重积分：求二重积分：10/10/44黄建华制作第45页4.1.4积分 clear ;syms x y;f=x2*sin(x*y);a1=int(f,x,y,1)a1=(-y2*cos(y)+2*cos(y)+2*y*sin(y)+y4*cos(y2).-2*cos(y2)-2*y2*sin(y2)/y3 I=int(a1,y,0,1)I=-1/2*sin(1)+1/2a2=int(f,y,0,x)a2=-cos(x2)*x+x I1=int(a2,x,0,1)I=-1/2*sin(1)+1/210/10/45黄建华制作第46页4.1.4积分v数值积分数值积分 一元函数数值积分：一元函数数值积分：q=quad(fun,a,b,tol)-采取辛普森计算积分采取辛普森计算积分q=quad8(fun,a,b,tol)-采取采取newton cotes方法计算积分方法计算积分q=quadl(fun,a,b,tol)-采取采取lobatto方法计算方法计算 tol表示绝对误差限，默认表示绝对误差限，默认10-6，a,b是确定值是确定值;fun能够是能够是字符串、内联函数或字符串、内联函数或M函数名。函数名。二重数值积分：二重数值积分：q=dblquad(fun,inmin,inmax,outmin,outmax,tol,method)inmin,inmax是内变量下限和上限，是内变量下限和上限，outmin,outmax是外是外变量下限和上限，只能是常数，即只能计算矩形域上积分。变量下限和上限，只能是常数，即只能计算矩形域上积分。10/10/46黄建华制作第47页4.1.4积分v例例4.1.21(2)求积分：求积分：比较以下三种结果：比较以下三种结果：syms x;i=int(exp(-x2),x,0,1),vpa(i,15)i=1/2*erf(1)*pi(1/2)0.746824132812427f=inline(exp(-x.2,x);q1=quad(f,0,1)q1=0.74682418072642q2=quadl(f,0,1)q2=0.7468241339884510/10/47黄建华制作第48页4.1.5曲线积分与曲面积分 曲线积分与曲面积分总是能够经过类似于重积曲线积分与曲面积分总是能够经过类似于重积分处理方法转化为定积分方式来求解。关键是分处理方法转化为定积分方式来求解。关键是掌握好转化公式和对积分限要求。掌握好转化公式和对积分限要求。10/10/48黄建华制作第49页4.1.5曲线积分与曲面积分v例例4.1.22求对坐标曲线积分：求对坐标曲线积分：C是圆周是圆周x2+y2=ax上半部分顺时针方向。上半部分顺时针方向。命令：命令：clear ;syms x y a;y=sqrt(a*x-x2);ydx=diff(y,x);f=x2+y2+4*x*y*ydx;L=int(f,x,0,a)L=1/6*a310/10/49黄建华制作第50页4.1.6空间解析几何与向量代数 主要内容主要内容v相关向量计算：模、方向余弦和方向角相关向量计算：模、方向余弦和方向角v相关向量计算：数量积、向量积相关向量计算：数量积、向量积v空间曲线和曲面绘制空间曲线和曲面绘制10/10/50黄建华制作第51页4.1.6空间解析几何与向量代数v相关向量计算：模、方向余弦和方向角相关向量计算：模、方向余弦和方向角v例例4.1.23求向量求向量 模、方向余弦和方向角。模、方向余弦和方向角。命令：命令：clear;a=2,-3,5;mo=sqrt(sum(a.2)%模模cx=2/mo;cy=-3/mo;cz=5/mo;c=cx,cy,cz%方向余弦方向余弦ax=acos(cx);ay=acos(cy);az=acos(cz);%方向角方向角A=ax,ay,az*180/pi%将弧度变为角度将弧度变为角度mo=6.1644c=0.3244 -0.4867 0.8111 A=71.0682 119.1216 35.795810/10/51黄建华制作第52页4.1.6空间解析几何与向量代数v相关向量计算：数量积、向量积相关向量计算：数量积、向量积v例例4.1.24求向量求向量 和和 数量积、向量积。数量积、向量积。命令：命令：clear;a=2,-3,1;b=1,-1,3;s=dot(a,b)%数量积等同于数量积等同于a*b ch=a*bs=8 ch=8 xlj1=cross(a,b)%a和和 b向量积向量积 xlj1=-8 -5 1xlj2=cross(b,a)%b和和 a向量积向量积 xlj2=8 5 -110/10/52黄建华制作第53页4.1.6空间解析几何与向量代数v空间曲线和曲面绘制空间曲线和曲面绘制 plot3(X,Y,Z)、mesh()、meshgrid()、surf()、ezmesh()等请参考详细章节和帮助。等请参考详细章节和帮助。v例例4.1.25绘制函数：绘制函数：y=1/xy=1/x围绕围绕y y轴旋转所形成旋转曲面。轴旋转所形成旋转曲面。命令：命令：clear;x=0.2:0.001:0.5;y=1./x;X,Y,Z=cylinder(y,30);%命令 cylinder(x,n)生成.绕母线x旋转曲面，n定义母线分格线条数 mesh(X,Y,Z)10/10/53黄建华制作第54页4.1.7级数 主要内容主要内容v级数求和与审敛级数求和与审敛v泰勒展开泰勒展开v傅立叶展开傅立叶展开 10/10/54黄建华制作第55页4.1.7级数v级数求和与审敛级数求和与审敛 级数求和与审敛实际是同一问题，只要能够求和，自然级数求和与审敛实际是同一问题，只要能够求和，自然收敛。级数求和命令格式：收敛。级数求和命令格式：symsum(fun,变量，起点变量，起点，终点终点)省略变量则对默认变量求和。省略变量则对默认变量求和。v例例4.1.26求求（1）（2）命令：命令：clear;syms n;f1=(2*n-1)/2n ;f2=1/(n*(2*n+1);I1=symsum(f1,n,1,inf)I1=3%收敛收敛I2=symsum(f2,n,1,inf)I2=2-2*log(2)10/10/55黄建华制作第56页4.1.7级数v级数求和与审敛级数求和与审敛v例例4.1.27求：求：命令：命令：clear;syms n m;f1=symsum(1/m,m,1,n);limit(f1-log(n)n,inf)ans=eulergamma%欧拉常数欧拉常数vpa(ans,20)ans=0.5772156649015328606110/10/56黄建华制作第57页4.1.7级数v级数求和与审敛级数求和与审敛v例例4.1.28求求（1）（2）命令：命令：clear;syms n x;f1=sin(x)/n2 ;f2=(-1)(n-1)*xn/n;I1=symsum(f1,n,1,inf)I1=1/6*sin(x)*pi2 I2=symsum(f2,n,1,inf)I2=log(1+x)10/10/57黄建华制作第58页4.1.7级数v 泰勒展开泰勒展开 命令格式：命令格式：taylor(fun,n，变量，变量，a)fun为待展函数；为待展函数；n为展开阶数，缺省是为展开阶数，缺省是6阶；变量为申阶；变量为申明明fun中变量，省略变量则对默认变量展开；中变量，省略变量则对默认变量展开；a为变量求导为变量求导取值点，缺省为取值点，缺省为0，即麦克劳林展开。，即麦克劳林展开。v例例4.1.29将将 展开成幂级数。展开成幂级数。命令：命令：clear;syms x;f=1/(1+x2);taylor(f)ans=1-x2+x4taylor(f,20)ans=1-x2+x4-x6+x8-x10+x12-x14+x16-x1810/10/58黄建华制作第59页4.1.7级数v例例4.1.30将将 展开成展开成(x-1)幂级数。幂级数。命令：命令：clear;syms x;f=1/(x2+4*x+3);taylor(f,10,x,1)ans=7/32-3/32*x+7/128*(x-1)2-15/512*(x-1)3+31/2048*(x-1)4-63/8192*(x-1)5+127/32768*(x-1)6-255/131072*(x-1)7+511/524288*(x-1)8-1023/2097152*(x-1)910/10/59黄建华制作第60页4.1.7级数v 傅立叶展开傅立叶展开 将函数展开为以下格式：将函数展开为以下格式：能够结合能够结合MATLAB积分命令积分命令int()计算：计算：即可进行傅立叶展开。即可进行傅立叶展开。10/10/60黄建华制作第61页4.1.8微分方程 求解微分方程（组）由命令求解微分方程（组）由命令dsolve()完成。格式：完成。格式：dsolve(方程方程1,方程方程n,条件条件1，条件条件m,变量变量1,.,变量变量k)其中方程其中方程i为待解方程；条件为初始状态，缺省则求通解；为待解方程；条件为初始状态，缺省则求通解；变量为微分自变量，缺省为默认。变量为微分自变量，缺省为默认。注意：在输入形式中，注意：在输入形式中，y记为记为Dy，y为为D2y，,y(n)为为Dny。t,x=ode23(方程函数名方程函数名,tspan,x0,选项选项,附加参数附加参数)t,x=ode45(方程函数名方程函数名,tspan,x0,选项选项,附加参数附加参数)分别采取二阶三级和四阶五级分别采取二阶三级和四阶五级RKF方法计算常微分方程数值方法计算常微分方程数值解，解，plot(t,x)为解曲线。为解曲线。10/10/61黄建华制作第62页4.1.8微分方程v例例4.1.31求解微分方程：求解微分方程：(1)(2)命令：命令：clear;syms x y;Y1=dsolve(Dy=1/(x+y),x)Y1=-lambertw(-C1*exp(-1-x)-1-x%lambertw(x)表示一个函数关系表示一个函数关系y*exp(y)=xY2=dsolve(D2y*y-Dy2=0,x)%若若D2y*y改为改为y*D2y会提醒犯错，不符合某种规则会提醒犯错，不符合某种规则 Y2=0 或或exp(C1*x)*C210/10/62黄建华制作第63页4.1.8微分方程v例例4.1.32求解微分方程：求解微分方程：命令：命令：clear;syms x y;f=(1+x2)*D2y=2*x*Dy;c1=y(0)=1,Dy(0)=3;Y=dsolve(f,c1,x)Y=1+3*x+x310/10/63黄建华制作第64页4.1.8微分方程v例例4.1.33求解微分方程：求解微分方程：命令：命令：clear;syms x y;Y=dsolve(D2y-5*Dy+6*y=x*exp(2*x),x)Y=exp(3*x)*C2+exp(2*x)*C1-1/2*x*exp(2*x)*(2+x)10/10/64黄建华制作第65页4.1.8微分方程v例例4.1.34求解微分方程组：求解微分方程组：命令：命令：clear;syms t x y;x,y=dsolve(Dx=x-y,Dy=x+y,t)x=exp(t)*(C1*cos(t)-C2*sin(t)y=exp(t)*(C1*sin(t)+C2*cos(t)10/10/65黄建华制作第66页4.2 概率统计概率统计v随机变量及其分布随机变量及其分布v随机变量函数分布随机变量函数分布v随机变量数字特征随机变量数字特征v参数预计参数预计v假设检验假设检验v方差分析方差分析10/10/66黄建华制作第67页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=hygecdf(x,M,N,K)功效：计算超几何分布累积概率，总共功效：计算超几何分布累积概率，总共M件产品，其中次品件产品，其中次品N 件，抽取件，抽取K件检验，计算发觉次品不多于件检验，计算发觉次品不多于x件概率件概率Fx=P次品数次品数Xx=F(x)命令命令2：x=hygeinv(p,M,N,K)功效：在已知参数功效：在已知参数M、N、K和和p情况下计算随机量情况下计算随机量x，使，使得得p=P0次品数次品数Xx命令命令3：X=hygernd(M,N,K,m,n)功效：在已知参数功效：在已知参数M,N,K情况下产生情况下产生m*n维符合超几何分维符合超几何分布随机数矩阵布随机数矩阵X10/10/67黄建华制作第68页4.2.1 随机变量及其分布随机变量及其分布命令命令4：Px=hygepdf(x,M,N,K)功效：总共功效：总共M件产品，其中次品件产品，其中次品N 件，抽取件，抽取K件检验，计算发件检验，计算发觉恰好觉恰好x件次品概率件次品概率Px=PX=x命令命令5：stairs(x,Px)功效：绘制以功效：绘制以 x为横坐标，为横坐标，Px为纵坐标阶梯平面图；当为纵坐标阶梯平面图；当Px是是分布列分布列(或密度或密度)时，绘制概率密度分布图；当时，绘制概率密度分布图；当Px是累积分是累积分布时，绘制概率分布函数图布时，绘制概率分布函数图注：注：以后碰到命令末尾为：以后碰到命令末尾为：rnd-产生随机数产生随机数X；cdf-产生分布函数产生分布函数F(x)pdf-产生密度函数产生密度函数p(x)或分布列或分布列Px=PX=xinv-计算计算x=F-1(p)p=F(x)10/10/68黄建华制作第69页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=binocdf(x,n,p)功效：计算二项分布累积概率功效：计算二项分布累积概率Fx=PXx=F(x)命令命令2：x=binoinv(y,n,p)功效：计算随机量功效：计算随机量x，使得，使得y=PXx命令命令3：X=binornd(n,p,M,N)功效：产生功效：产生M*N维符合二项分布随机数矩阵维符合二项分布随机数矩阵X命令命令4：Px=binopdf(x,n,p)功效：计算试验中事件恰好发生功效：计算试验中事件恰好发生x次概率次概率10/10/69黄建华制作第70页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=poisscdf(x,lambda)功效：计算累积概率功效：计算累积概率Fx=PXx=F(x)命令命令2：x=poissinv(p,lambda)功效：计算随机量功效：计算随机量x，使得，使得p=PXx命令命令3：X=poissrnd(lambda,M,N)功效：产生功效：产生M*N维随机数矩阵维随机数矩阵X命令命令4：Px=poisspdf(x,lambda)功效：计算概率功效：计算概率Px=PX=x10/10/70黄建华制作第71页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=normcdf(x,mu,sigma)功效：计算累积概率功效：计算累积概率Fx=PXx=F(x)命令命令2：x=norminv(p,mu,sigma)功效：计算随机量功效：计算随机量x，使得，使得p=PXx命令命令3：X=normrnd(mu,sigma,M,N)功效：产生功效：产生M*N维随机数矩阵维随机数矩阵X命令命令4：Px=normpdf(x,mu,sigma)功效：计算分布密度功效：计算分布密度p(x)在在x值值补充：补充：randn()-标准正态分布随机数标准正态分布随机数10/10/71黄建华制作第72页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=expcdf(x,lambda)功效：计算累积概率功效：计算累积概率Fx=PXx=F(x)命令命令2：x=expinv(p,lambda)功效：计算随机量功效：计算随机量x，使得，使得p=PXx命令命令3：X=exprnd(lambda,M,N)功效：产生功效：产生M*N维随机数矩阵维随机数矩阵X命令命令4：Px=exppdf(x,lambda)功效：计算分布密度功效：计算分布密度p(x)在在x值值10/10/72黄建华制作第73页4.2.1 随机变量及其分布随机变量及其分布命令命令1：Fx=unifcdf(x,a,b)功效：计算累积概率功效：计算累积概率Fx=PXx=F(x)命令命令2：x=unifinv(p,a,b)功效：计算随机量功效：计算随机量x，使得，使得p=PXx命令命令3：X=unifrnd(a,b,M,N)功效：产生功效：产生M*N维随机数矩阵维随机数矩阵X命令命令4：Px=unifpdf(x,a,b)功效：计算分布密度功效：计算分布密度p(x)在在x值值补充：补充：rand()-(0,1)均匀分布随机数均匀分布随机数10/10/73黄建华制作第74页4.2.1 随机变量及其分布随机变量及其分布命令命令:gamcdf(x,a,lambda),gaminv(p,a,lambda)gampdf(x,a,lambda),gamrnd(a,lambda,m,n)10/10/74黄建华制作第75页4.2.1 随机变量及其分布随机变量及其分布命令命令:chi2cdf(x,k),chi2inv(p,k),chi2pdf(x,k)chi2rnd(k,m,n)10/10/75黄建华制作第76页4.2.1 随机变量及其分布随机变量及其分布命令命令:tcdf(x,k),tinv(p,k),tpdf(x,k)trnd(k,m,n)10/10/76黄建华制作第77页4.2.1 随机变量及其分布随机变量及其分布命令命令:fcdf(x,p,q),finv(F,p,q),fpdf(x,p,q)frnd(p,q,m,n)10/10/77黄建华制作第78页4.2.1 随机变量及其分布随机变量及其分布程序：程序：clear;px=binopdf(45,100,0.5)%计算计算x=45概率概率 px=0.0485fx=binocdf(45,100,0.5)%