四MATLAB求解數(shù)學(xué)問題市公開課一等獎百校聯(lián)賽特等獎?wù)n件

（四）MATLAB求解數(shù)學(xué)問題數(shù)學(xué)分析概率論與數(shù)理統(tǒng)計(jì)插值與擬合優(yōu)化運(yùn)算線性代數(shù)復(fù)變函數(shù)10/10/0黃建華制作第1頁4.1數(shù)學(xué)分析符號方程求解極限導(dǎo)數(shù)與微分(重)積分曲線積分與曲面積分空間解析幾何與向量代數(shù)級數(shù)微分方程10/10/1黃建華制作第2頁4.1.1符號方程求解

主要內(nèi)容

線性方程

非線性方程10/10/2黃建華制作第3頁4.1.1符號方程求解線性方程

慣用solve()和linsolve()函數(shù)來處理線性方程問題。詳細(xì)格式：X=solve('方程1',….'方程n','變量1',...'變量m')說明：能夠求解方程組，單變量時(shí)變量申明能夠省略。

X=linsolve(A,B)%求解線性方程組AX=B，返回特解X

10/10/3黃建華制作第4頁4.1.1符號方程求解例4.1.1求解方程：命令：>>clear;x=solve('x^2-x-6=0')相當(dāng)于：clear;symsx;f=x^2-x-6;x=solve(f)x=

3-210/10/4黃建華制作第5頁4.1.1符號方程求解例4.1.2求解方程組：命令：>>clear;[x,y]=solve('x^2+y-6=0','y^2+x-6=0','x','y')x=2y=2-3-31/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=

，則僅將返回一個(gè)解結(jié)構(gòu)。

X=x:[4x1sym]y:[4x1sym]10/10/5黃建華制作第6頁4.1.1符號方程求解例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.00001.000010/10/6黃建華制作第7頁4.1.1符號方程求解非線性方程

慣用fsolve()函數(shù)來處理非線性方程問題。詳細(xì)格式：X=fsolve(@fun,x0,options)說明：fun為定義好非線性方程(組)文件名，其中@為調(diào)用函數(shù)符號；x0為求解方程初始向量；options設(shè)置求解過程各種參數(shù)，普通采取默認(rèn)參數(shù)optimset('fsolve'),其它參數(shù)能夠查詢幫助。

10/10/7黃建華制作第8頁4.1.1符號方程求解例4.1.4求解方程組：x0=[x(1),x(2)]=[0.1,0.1]10/10/8黃建華制作第9頁4.1.1符號方程求解首先建立函數(shù)文件fun.m并保留在默認(rèn)路徑下：functiony=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))];然后運(yùn)行命令：>>clear;x0=[0.1,0.1];>>x=fsolve(@fun,x0,optimset('fsolve'))x=0.54140.331010/10/9黃建華制作第10頁4.1.2極限

主要內(nèi)容單變量函數(shù)極限多變量函數(shù)極限10/10/10黃建華制作第11頁4.1.2極限單變量函數(shù)極限首先進(jìn)行符號變量說明:symsxytha然后定義函數(shù)fun，再使用以下命令格式求對應(yīng)極限:limit(fun,x,a)%求函數(shù)fun當(dāng)xa時(shí)極限limit(fun,a)%默認(rèn)變量x或唯一符號變量limit(fun)%默認(rèn)變量x,且a=0limit(fun,x,a,'right')%右極限xa+limit(fun,x,a,'left')%左極限xa-10/10/11黃建華制作第12頁4.1.2極限例4.1.5舉例:結(jié)果>>

symsxha>>

f=sin(x)/x;>>limit(f)1>>limit(f,inf)0>>limit((x-2)/(x^2-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極限例4.1.6試求解極限問題:解：>>symsxab;f=x*(1+a/x)^x*sin(b/x);

>>L=limit(f,x,inf)L=exp(a)*b10/10/13黃建華制作第14頁4.1.2極限多變量函數(shù)極限假設(shè)有二元函數(shù)求極限問題：則能夠嵌套使用limit()函數(shù):

limit(limit（fun,x,x0),y,y0)或

limit(limit（fun,y,y0),x,x0)假如x0或y0不是確定值，而是另一個(gè)變量函數(shù)，則次序不能交換。注意：此種使用方法只適合用于極限存在情況。10/10/14黃建華制作第15頁4.1.2極限例4.1.7求出二元函數(shù)極限值:解：>>symsxy;f=sin(x*y)/x;

>>L=limit(limit(f,x,0),y,0)

>>L1=limit(limit(f,y,0),x,0)L=0L1=010/10/15黃建華制作第16頁4.1.2極限例4.1.8求出二元函數(shù)極限值:解>>symsxyk;f=(x^2-y^2)/(x^2+y^2);

(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=-1L1=1L2=-(-1+k^2)/(1+k^2)10/10/16黃建華制作第17頁4.1.3導(dǎo)數(shù)和微分

主要內(nèi)容導(dǎo)數(shù)和高階導(dǎo)數(shù)高階混合偏導(dǎo)數(shù)復(fù)合函數(shù)求導(dǎo)隱函數(shù)求偏導(dǎo)參數(shù)方程求導(dǎo)導(dǎo)數(shù)應(yīng)用梯度計(jì)算和方向?qū)?shù)10/10/17黃建華制作第18頁4.1.3導(dǎo)數(shù)和微分導(dǎo)數(shù)和高階導(dǎo)數(shù)首先進(jìn)行符號變量說明:syms然后定義函數(shù)f，再使用以下命令格式求對應(yīng)導(dǎo)數(shù):diff(f)

%f對默認(rèn)變量x求一階導(dǎo)數(shù)diff(f,v)

%f對變量v求一階導(dǎo)數(shù)diff(f,n)

%f對默認(rèn)變量x求n階導(dǎo)數(shù)diff(f,v,n)

%f對變量v求n階導(dǎo)數(shù)顯然，用以上命令能夠?qū)崿F(xiàn)求各階純偏導(dǎo)。微分與導(dǎo)數(shù)是共通，只須將求導(dǎo)答案變形一下即可！10/10/18黃建華制作第19頁4.1.3導(dǎo)數(shù)和微分例4.1.9

命令：>>symsaxy;f=x^3+y^2+a*x*y-3*x+7;>>f1=diff(f)結(jié)果：f1=3*x^2+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導(dǎo)數(shù)和微分高階混合偏導(dǎo)數(shù)假設(shè)有多元函數(shù)求偏導(dǎo)問題：則能夠嵌套使用diff()函數(shù):命令格式:diff(difft(f,x,m),y,n)或

diff(diff(f,y,n),x,m)10/10/20黃建華制作第21頁4.1.3導(dǎo)數(shù)和微分例4.1.10

求：命令：>>symsxy;z=(x^2-2*x)*exp(-x^2-y^2-x*y);

>>zxzy=diff(diff(z,x),y)

zxzy=(2*x-2)*(-2*y-x)*exp(-x^2-y^2-x*y)-(x^2-2*x)*exp(-x^2-y^2-x*y)+(x^2-2*x)*(-2*x-y)*(-2*y-x)*exp(-x^2-y^2-x*y)

10/10/21黃建華制作第22頁4.1.3導(dǎo)數(shù)和微分化簡一下：命令：>>zxzy1=simple(zxzy)

zxzy1=exp(-x^2-y^2-x*y)*(-4*x*y-3*x^2+4*y+4*x+5*x^3*y+2*x^4+2*x^2*y^2-10*x^2*y-4*x^3-4*x*y^2)10/10/22黃建華制作第23頁4.1.3導(dǎo)數(shù)和微分例4.1.11

求：命令：>>symsxy;f=(x^2-2*x)*exp(-x^2-y^2-x*y);

>>ydx=-diff(f,x)/diff(f,y)

ydx=(-(2*x-2)*exp(-x^2-y^2-x*y)-(x^2-2*x)*(-2*x-y)*exp(-x^2-y^2-x*y))/(x^2-2*x)/(-2*y-x)/exp(-x^2-y^2-x*y)10/10/23黃建華制作第24頁4.1.3導(dǎo)數(shù)和微分復(fù)合函數(shù)求導(dǎo)例4.1.12已知：求：命令：>>symstxy;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));輸入次序顛倒，不然結(jié)果0。10/10/24黃建華制作第25頁4.1.3導(dǎo)數(shù)和微分隱函數(shù)求偏導(dǎo)假設(shè)有隱函數(shù)表示式f(x1,x2,…,xn)=0,求偏導(dǎo)問題：

能夠使用diff()函數(shù),命令格式:F=-diff(f,xj)/diff(f,xi)10/10/25黃建華制作第26頁4.1.3導(dǎo)數(shù)和微分例4.1.13

已知

求：命令：>>clear;symsxy;F=atan(y/x)-log(sqrt(x^2+y^2);>>ydx=-diff(F,y)/diff(F,x)ydx=(-1/x/(1+y^2/x^2)+1/(x^2+y^2)*y)/(-y/x^2/(1+y^2/x^2)-1/(x^2+y^2)*x)

>>simple(ydx)ydx=(x-y)/(y+x)10/10/26黃建華制作第27頁4.1.3導(dǎo)數(shù)和微分參數(shù)方程求導(dǎo)假設(shè)有參數(shù)方程表示式y(tǒng)=f(t),x=g(t),求導(dǎo)數(shù)：

能夠使用diff()函數(shù)遞歸調(diào)用,命令格式:dk=diff(dk-1,t)/diff(x,t)其中dk-1表示k-1階導(dǎo)數(shù)注意不能用：dk=diff(y,t,k)/diff(x,t,k)10/10/27黃建華制作第28頁4.1.3導(dǎo)數(shù)和微分例4.1.14

已知

求：命令：>>symsabtxy;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---------atan(t)10/10/28黃建華制作第29頁4.1.3導(dǎo)數(shù)和微分>>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------------------------22a(-1+cos(t))sin(t)10/10/29黃建華制作第30頁4.1.3導(dǎo)數(shù)和微分若改用：>>xd2=diff(y,t,2)/diff(x,t,2)xd2=b*sin(t)/a/cos(t)>>pretty(simple(xd2))

btan(t)--------a顯然不一樣。10/10/30黃建華制作第31頁4.1.3導(dǎo)數(shù)和微分導(dǎo)數(shù)應(yīng)用例4.1.15

討論函數(shù)

極值、單調(diào)性和其導(dǎo)數(shù)函數(shù)關(guān)系。命令：>>clear;symsxydyd2y;y=x^2/(1+x^2);dy=simple(diff(y));pretty(dy)

x2---------22(1+x)

10/10/31黃建華制作第32頁4.1.3導(dǎo)數(shù)和微分>>Px=solve(dy)

Px=0>>d2y=simplify(diff(y,2))d2y=-2*(-1+3*x^2)/(1+x^2)^3

>>

P2x=solve(d2y)P2x=-1/3*3^(1/2)1/3*3^(1/2)10/10/32黃建華制作第33頁4.1.3導(dǎo)數(shù)和微分>>lims=[-5,5];subplot(3,1,1);ezplot(y,lims);holdon;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);holdon;line([0,0],[-1,1.5]);line([-5,5],[0,0]);%同時(shí)繪制橫軸subplot(3,1,3);ezplot(d2y,[-5,5]);holdon;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導(dǎo)數(shù)和微分10/10/34黃建華制作第35頁4.1.3導(dǎo)數(shù)和微分梯度計(jì)算和方向?qū)?shù)jacobian(fun,v)%v是求導(dǎo)變量向量，表示fun對v求偏導(dǎo)矩陣即梯度gridient(F)%求F數(shù)值梯度，一維時(shí)可用diff代替dot（jacobian(fun)，v）=jacobian(fun)·v

%v是某方向單位向量，數(shù)量積就是方向?qū)?shù)10/10/35黃建華制作第36頁4.1.4積分主要內(nèi)容不定積分定積分與無窮積分重積分?jǐn)?shù)值積分10/10/36黃建華制作第37頁4.1.4積分不定積分不定積分問題：能夠使用int()函數(shù):F=int(fun,x)或F=int(fun)%當(dāng)fun中只有一個(gè)自變量x，則x可省最終答案應(yīng)為：F(x)+C10/10/37黃建華制作第38頁4.1.4積分例4.1.16用diff()函數(shù)求4階導(dǎo)數(shù)，再積分，

檢驗(yàn)是否能夠得出一致結(jié)果。命令：>>symsx;y=sin(x)/(x^2+4*x+3);y4=diff(y,4);y0=int(int(int(int(y4))));>>pretty(simple(y0))%對導(dǎo)數(shù)積分應(yīng)該得出原函數(shù)

sin(x)---------------(x+1)(x+3)10/10/38黃建華制作第39頁4.1.4積分例4.1.17證實(shí)：命令：>>symsax;f=simple(int(x^3*cos(a*x)^2,x));f1=x^4/8+(x^3/(4*a)-3*x/(8*a^3))*sin(2*a*x)+...(3*x^2/(8*a^2)-3/(16*a^4))*cos(2*a*x);simple(f-f1)%求兩個(gè)結(jié)果差

-3/16/a^4結(jié)果是一個(gè)常數(shù)，表明答案正確。10/10/39黃建華制作第40頁4.1.4積分例4.1.18不可積問題：命令：>>symsx;int(exp(x^2))ans=-1/2*i*pi^(1/2)*erf(i*x)結(jié)果中erf是一個(gè)定義函數(shù)：表明不可積！10/10/40黃建華制作第41頁4.1.4積分定積分與無窮積分定積分問題：能夠使用int()函數(shù):int(fun,x,a,b)若為無窮積分問題，則只需將命令中a(或b)改為-inf(或inf)即可。如求：用：int(fun,x,a,inf)10/10/41黃建華制作第42頁4.1.4積分例4.1.19不可積問題定積分可積。命令：>>symsx;int(exp(-x^2/2),x,0,inf)ans=1/2*2^(1/2)*pi^(1/2)命令：>>symsx;I=int(exp(-x^2/2),x,0,1)I=1/2*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/2)能夠使用vpa()函數(shù)顯示數(shù)值：vpa(I,5)=0.8556510/10/42黃建華制作第43頁4.1.4積分例4.1.20變限積分也可使用定積分求解：命令：>>symsxt;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積分重積分重積分問題能夠先化為累次積分方式再使用int()函數(shù)嵌套來處理:例4.1.21求二重積分：10/10/44黃建華制作第45頁4.1.4積分>>clear;symsxy;f=x^2*sin(x*y);a1=int(f,x,y,1)a1=(-y^2*cos(y)+2*cos(y)+2*y*sin(y)+y^4*cos(y^2)...-2*cos(y^2)-2*y^2*sin(y^2))/y^3>>I=int(a1,y,0,1)I=-1/2*sin(1)+1/2>>a2=int(f,y,0,x)

a2=-cos(x^2)*x+x>>I1=int(a2,x,0,1)I=-1/2*sin(1)+1/210/10/45黃建華制作第46頁4.1.4積分?jǐn)?shù)值積分

一元函數(shù)數(shù)值積分：q=quad(fun,a,b,tol)----采取辛普森計(jì)算積分q=quad8(fun,a,b,tol)----采取newtoncotes方法計(jì)算積分q=quadl(fun,a,b,tol)-----采取lobatto方法計(jì)算tol表示絕對誤差限，默認(rèn)10-6，a,b是確定值;fun能夠是字符串、內(nèi)聯(lián)函數(shù)或M函數(shù)名。二重?cái)?shù)值積分：q=dblquad(fun,inmin,inmax,outmin,outmax,tol,method)

inmin,inmax是內(nèi)變量下限和上限，outmin,outmax是外變量下限和上限，只能是常數(shù)，即只能計(jì)算矩形域上積分。10/10/46黃建華制作第47頁4.1.4積分例4.1.21(2)求積分：比較以下三種結(jié)果：>>symsx;i=int('exp(-x^2)','x',0,1),vpa(i,15)

i=1/2*erf(1)*pi^(1/2)0.746824132812427>>f=inline('exp(-x.^2','x');q1=quad(f,0,1)q1=0.74682418072642>>q2=quadl(f,0,1)q2=0.7468241339884510/10/47黃建華制作第48頁4.1.5曲線積分與曲面積分

曲線積分與曲面積分總是能夠經(jīng)過類似于重積分處理方法轉(zhuǎn)化為定積分方式來求解。關(guān)鍵是掌握好轉(zhuǎn)化公式和對積分限要求。10/10/48黃建華制作第49頁4.1.5曲線積分與曲面積分例4.1.22求對坐標(biāo)曲線積分：C是圓周x2+y2=ax上半部分順時(shí)針方向。命令：>>clear;symsxya;y=sqrt(a*x-x^2);ydx=diff(y,x);f=x^2+y^2+4*x*y*ydx;L=int(f,x,0,a)L=1/6*a^310/10/49黃建華制作第50頁4.1.6空間解析幾何與向量代數(shù)

主要內(nèi)容相關(guān)向量計(jì)算：模、方向余弦和方向角相關(guān)向量計(jì)算：數(shù)量積、向量積空間曲線和曲面繪制10/10/50黃建華制作第51頁4.1.6空間解析幾何與向量代數(shù)相關(guān)向量計(jì)算：模、方向余弦和方向角例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%將弧度變?yōu)榻嵌萴o=6.1644c=0.3244-0.48670.8111A=71.0682119.121635.795810/10/51黃建華制作第52頁4.1.6空間解析幾何與向量代數(shù)相關(guān)向量計(jì)算：數(shù)量積、向量積例4.1.24求向量和數(shù)量積、向量積。命令：>>clear;a=[2,-3,1];b=[1,-1,3];s=dot(a,b)%數(shù)量積等同于a*b’ch=a*b's=8ch=8

xlj1=cross(a,b)%a和b向量積xlj1=-8-51xlj2=cross(b,a)%b和a向量積xlj2=85-110/10/52黃建華制作第53頁4.1.6空間解析幾何與向量代數(shù)空間曲線和曲面繪制

plot3(X,Y,Z)、mesh()、meshgrid()、surf()、ezmesh()等請參考詳細(xì)章節(jié)和幫助。例4.1.25繪制函數(shù)：y=1/x圍繞y軸旋轉(zhuǎn)所形成旋轉(zhuǎn)曲面。命令：>>clear;x=0.2:0.001:0.5;y=1./x;

>>[X,Y,Z]=cylinder(y,30);%命令cylinder(x,n)生成...繞母線x旋轉(zhuǎn)曲面，n定義母線分格線條數(shù)>>mesh(X,Y,Z)10/10/53黃建華制作第54頁4.1.7級數(shù)主要內(nèi)容級數(shù)求和與審斂泰勒展開傅立葉展開

10/10/54黃建華制作第55頁4.1.7級數(shù)級數(shù)求和與審斂

級數(shù)求和與審斂實(shí)際是同一問題，只要能夠求和，自然收斂。級數(shù)求和命令格式：

symsum(fun,變量，起點(diǎn)，終點(diǎn))

省略變量則對默認(rèn)變量求和。例4.1.26求（1）（2）命令：>>clear;symsn;f1=(2*n-1)/2^n;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級數(shù)級數(shù)求和與審斂例4.1.27求：命令：>>clear;symsnm;f1=symsum(1/m,m,1,n);limit(f1-log(n)n,inf)

ans=eulergamma%歐拉常數(shù)>>vpa(ans,20)ans=0.5772156649015328606110/10/56黃建華制作第57頁4.1.7級數(shù)級數(shù)求和與審斂例4.1.28求（1）（2）命令：>>clear;symsnx;f1=sin(x)/n^2;f2=(-1)^(n-1)*x^n/n;>>I1=symsum(f1,n,1,inf)

I1=1/6*sin(x)*pi^2>>I2=symsum(f2,n,1,inf)

I2=log(1+x)10/10/57黃建華制作第58頁4.1.7級數(shù)泰勒展開命令格式：taylor(fun,n，變量，a)

fun為待展函數(shù)；n為展開階數(shù)，缺省是6階；變量為申明fun中變量，省略變量則對默認(rèn)變量展開；a為變量求導(dǎo)取值點(diǎn)，缺省為0，即麥克勞林展開。例4.1.29將展開成冪級數(shù)。命令：>>clear;symsx;f=1/(1+x^2);>>taylor(f)ans=1-x^2+x^4>>taylor(f,20)ans=1-x^2+x^4-x^6+x^8-x^10+x^12-x^14+x^16-x^1810/10/58黃建華制作第59頁4.1.7級數(shù)例4.1.30將展開成(x-1)冪級數(shù)。命令：>>clear;symsx;f=1/(x^2+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級數(shù)傅立葉展開將函數(shù)展開為以下格式：能夠結(jié)合MATLAB積分命令int()計(jì)算：即可進(jìn)行傅立葉展開。10/10/60黃建華制作第61頁4.1.8微分方程求解微分方程（組）由命令dsolve()完成。格式：

dsolve(('方程1',…,'方程n','條件1'，…，'條件m','變量1',...,'變量k')其中方程i為待解方程；條件為初始狀態(tài)，缺省則求通解；變量為微分自變量，缺省為默認(rèn)。注意：在輸入形式中，y'記為Dy，y''為D2y，…,y(n)為Dny。[t,x]=ode23(方程函數(shù)名,tspan,x0,選項(xiàng),附加參數(shù))[t,x]=ode45(方程函數(shù)名,tspan,x0,選項(xiàng),附加參數(shù))分別采取二階三級和四階五級RKF方法計(jì)算常微分方程數(shù)值解，plot(t,x)為解曲線。10/10/61黃建華制作第62頁4.1.8微分方程例4.1.31求解微分方程：(1)(2)命令：>>clear;symsxy;Y1=dsolve('Dy=1/(x+y)','x')Y1=-lambertw(-C1*exp(-1-x))-1-x%lambertw(x)表示一個(gè)函數(shù)關(guān)系y*exp(y)=x>>Y2=dsolve('D2y*y-Dy^2=0','x')%若D2y*y改為y*D2y會提醒犯錯(cuò)，不符合某種規(guī)則

Y2=0或exp(C1*x)*C210/10/62黃建華制作第63頁4.1.8微分方程例4.1.32求解微分方程：命令：>>clear;symsxy;f='(1+x^2)*D2y=2*x*Dy';c1='y(0)=1,Dy(0)=3';Y=dsolve(f,c1,'x')Y=1+3*x+x^310/10/63黃建華制作第64頁4.1.8微分方程例4.1.33求解微分方程：命令：>>clear;symsxy;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微分方程例4.1.34求解微分方程組：命令：>>clear;symstxy;[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概率統(tǒng)計(jì)隨機(jī)變量及其分布隨機(jī)變量函數(shù)分布隨機(jī)變量數(shù)字特征參數(shù)預(yù)計(jì)假設(shè)檢驗(yàn)方差分析10/10/66黃建華制作第67頁4.2.1隨機(jī)變量及其分布超幾何分布H(n,M,N)命令1：Fx=hygecdf(x,M,N,K)功效：計(jì)算超幾何分布累積概率，總共M件產(chǎn)品，其中次品N件，抽取K件檢驗(yàn)，計(jì)算發(fā)覺次品不多于x件概率Fx=P{次品數(shù)X≤x}=F(x)命令2：x=hygeinv(p,M,N,K)功效：在已知參數(shù)M、N、K和p情況下計(jì)算隨機(jī)量x，使得p=P{0≤次品數(shù)X≤x}命令3：X=hygernd(M,N,K,m,n)功效：在已知參數(shù)M,N,K情況下產(chǎn)生m*n維符合超幾何分布隨機(jī)數(shù)矩陣X10/10/67黃建華制作第68頁4.2.1隨機(jī)變量及其分布命令4：Px=hygepdf(x,M,N,K)功效：總共M件產(chǎn)品，其中次品N件，抽取K件檢驗(yàn)，計(jì)算發(fā)覺恰好x件次品概率Px=P{X=x}命令5：stairs(x,Px)功效：繪制以x為橫坐標(biāo)，Px為縱坐標(biāo)階梯平面圖；當(dāng)Px是分布列(或密度)時(shí)，繪制概率密度分布圖；當(dāng)Px是累積分布時(shí)，繪制概率分布函數(shù)圖注：以后碰到命令末尾為：rnd----產(chǎn)生隨機(jī)數(shù)X；cdf----產(chǎn)生分布函數(shù)F(x)pdf----產(chǎn)生密度函數(shù)p(x)或分布列Px=P{X=x}inv----計(jì)算x=F-1(p)→p=F

(x)10/10/68黃建華制作第69頁4.2.1隨機(jī)變量及其分布二項(xiàng)分布B(n,p)命令1：Fx=binocdf(x,n,p)功效：計(jì)算二項(xiàng)分布累積概率Fx=P{X≤x}=F(x)命令2：x=binoinv(y,n,p)功效：計(jì)算隨機(jī)量x，使得y=P{X≤x}命令3：X=binornd(n,p,M,N)功效：產(chǎn)生M*N維符合二項(xiàng)分布隨機(jī)數(shù)矩陣X命令4：Px=binopdf(x,n,p)功效：計(jì)算試驗(yàn)中事件恰好發(fā)生x次概率10/10/69黃建華制作第70頁4.2.1隨機(jī)變量及其分布泊松分布X~P(λ)命令1：Fx=poisscdf(x,lambda)功效：計(jì)算累積概率Fx=P{X≤x}=F(x)命令2：x=poissinv(p,lambda)功效：計(jì)算隨機(jī)量x，使得p=P{X≤x}命令3：X=poissrnd(lambda,M,N)功效：產(chǎn)生M*N維隨機(jī)數(shù)矩陣X命令4：Px=poisspdf(x,lambda)功效：計(jì)算概率Px=P{X=x}10/10/70黃建華制作第71頁4.2.1隨機(jī)變量及其分布正態(tài)分布X~N(μ,σ2)命令1：Fx=normcdf(x,mu,sigma)功效：計(jì)算累積概率Fx=P{X≤x}=F(x)命令2：x=norminv(p,mu,sigma)功效：計(jì)算隨機(jī)量x，使得p=P{X≤x}命令3：X=normrnd(mu,sigma,M,N)功效：產(chǎn)生M*N維隨機(jī)數(shù)矩陣X命令4：Px=normpdf(x,mu,sigma)功效：計(jì)算分布密度p(x)在x值補(bǔ)充：randn()---標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)10/10/71黃建華制作第72頁4.2.1隨機(jī)變量及其分布指數(shù)分布X~exp(λ)命令1：Fx=expcdf(x,lambda)功效：計(jì)算累積概率Fx=P{X≤x}=F(x)命令2：x=expinv(p,lambda)功效：計(jì)算隨機(jī)量x，使得p=P{X≤x}命令3：X=exprnd(lambda,M,N)功效：產(chǎn)生M*N維隨機(jī)數(shù)矩陣X命令4：Px=exppdf(x,lambda)功效：計(jì)算分布密度p(x)在x值10/10/72黃建華制作第73頁4.2.1隨機(jī)變量及其分布均勻分布X~U(a,b)命令1：Fx=unifcdf(x,a,b)功效：計(jì)算累積概率Fx=P{X≤x}=F(x)命令2：x=unifinv(p,a,b)功效：計(jì)算隨機(jī)量x，使得p=P{X≤x}命令3：X=unifrnd(a,b,M,N)功效：產(chǎn)生M*N維隨機(jī)數(shù)矩陣X命令4：Px=unifpdf(x,a,b)功效：計(jì)算分布密度p(x)在x值補(bǔ)充：rand()---(0,1)均勻分布隨機(jī)數(shù)10/10/73黃建華制作第74頁4.2.1隨機(jī)變量及其分布

Γ分布命令:gamcdf(x,a,lambda),gaminv(p,a,lambda)gampdf(x,a,lambda),gamrnd(a,lambda,m,n)10/10/74黃建華制作第75頁4.2.1隨機(jī)變量及其分布

Χ2分布命令:chi2cdf(x,k),chi2inv(p,k),chi2pdf(x,k)chi2rnd(k,m,n)10/10/75黃建華制作第76頁4.2.1隨機(jī)變量及其分布

T分布命令:tcdf(x,k),tinv(p,k),tpdf(x,k)trnd(k,m,n)10/10/76黃建華制作第77頁4.2.1隨機(jī)變量及其分布

F分布命令:fcdf(x,p,q),finv(F,p,q),fpdf(x,p,q)frnd(p,q,m,n)10/10/77黃建華制作第78頁4.2.1隨機(jī)變量及其分布

例4.2.1某人向空中拋硬幣100次，落下為正面概率為0.5。這100次中正面向上次數(shù)記為X：(1)試計(jì)算x=45概率和x≤45概率；(2)繪制分布函數(shù)圖象和分布列圖象。程序：》clear;px=binopdf(45,100,0.5)%計(jì)算x=45概率px=0.0485fx=binocdf(45,100,0.5)%計(jì)算x≤45概率fx=0.1841》x=1:100;p1=binocdf(x,100,0.5);plot(x,p1,'+');title('分布函數(shù)圖')10/10/78黃建華制作第79頁4.2.1隨機(jī)變量及其分布p2=binopdf(x,100,0.5);plot(x,p2,'*r');title('概率分布圖')10/10/79黃建華制作第80頁4.2.1隨機(jī)變量及其分布

例4.2.2設(shè)X~N(2,0.25)(1)求概率P{1<X<2.5}；(2)繪制分布函數(shù)圖象和分布密度圖象;(3)畫出區(qū)間[1.5,1.9]上分布密度曲線下方區(qū)域。程序:(1)p=normcdf(2.5,2,0.5)-normcdf(1,2,0.5)

p=0.8186(2)x=0:0.1:4;px=normpdf(x,2,0.5);fx=normcdf(x,2,0.5);plot(x,px,'+b');holdon;plot(x,fx,'*r');legend('正態(tài)分布函數(shù)','正態(tài)分布密度');(3)specs=[1.5,1.9];pp=normspec(specs,2,0.5)10/10/80黃建華制作第81頁4.2.1隨機(jī)變量及其分布10/10/81黃建華制作第82頁4.2.2隨機(jī)變量函數(shù)分布依據(jù)概率統(tǒng)計(jì)教材中定理：假如已知隨機(jī)變量X密度fX(x),隨機(jī)變量函數(shù)Y=g(X)單調(diào)，則Y密度函數(shù)為：fY(x)=fX(h(y))|h'(y)|,其中x=h(y)是y=g(x)反函數(shù)。假如y=g(x)不單調(diào)，則將定義域分成若干單調(diào)區(qū)間進(jìn)行討論。也可利用：據(jù)此意思，計(jì)算隨機(jī)變量函數(shù)分布相當(dāng)于編程10/10/82黃建華制作第83頁

例4.2.3設(shè)隨機(jī)變量X服從均勻分布U[0,1]，求Y=eX分布。程序：clear;x=solve('y=exp(x)')

x=log(y)dy=diff(x,'y')

dy=1/yfy=1*abs(dy)

fy=1/|y|注：取值區(qū)域需要自己確定，用積分求法作為練習(xí)！4.2.2隨機(jī)變量函數(shù)分布10/10/83黃建華制作第84頁4.2.3隨機(jī)變量數(shù)字特征隨機(jī)變量數(shù)學(xué)期望1.數(shù)組平均值---Y=mean(X)功效：當(dāng)X為向量時(shí)，輸出一個(gè)平均數(shù)；當(dāng)X為矩陣時(shí)，輸出為行向量，對應(yīng)于矩陣每列平均值；所以計(jì)算矩陣全部數(shù)平均值，應(yīng)用嵌套：mean(mean(X))或m=mean(X(:))與這類似有：求和(sum),最大(max),最小(min)等2.離散型隨機(jī)變量期望----EX=sum(X.*P)功效：計(jì)算隨機(jī)值向量X與對應(yīng)概率向量P乘積之和3.連續(xù)型隨機(jī)變量期望----EX=int(x*fx,x,a,b)功效：用積分計(jì)算期望10/10/84黃建華制作第85頁4.2.3隨機(jī)變量數(shù)字特征

例4.2.4設(shè)隨機(jī)變量X分布列，求期望。程序：clear;x=[-1,0,2,3];p=[1/8,1/4,3/8,1/4];EX=sum(x.*p)

1.3750X-1023P1/81/43/81/410/10/85黃建華制作第86頁4.2.3隨機(jī)變量數(shù)字特征例4.2.5設(shè)隨機(jī)變量X分布密度為：且EX=3/5，求常數(shù)a,b值。程序：clear;symsabx;fx=a+b*x^2;EX=int(x*fx,x,0,1)EX=1/4*b+1/2*a

F=int(fx,x,0,1)F=a+1/3*b

f1=EX-3/5;f2=f-1;[a,b]=solve(f1,f2)

a=3/5,b=6/510/10/86黃建華制作第87頁4.2.3隨機(jī)變量數(shù)字特征

例4.2.6設(shè)隨機(jī)變量X分布密度為：求隨機(jī)變量Y=|X|期望。程序：clear;symsx;fx1=0.5*exp(x);fx2=0.5*exp(-x);EY=int(-x*fx1,x,-inf,0)+int(x*fx2,x,0,inf)

EY=110/10/87黃建華制作第88頁4.2.3隨機(jī)變量數(shù)字特征隨機(jī)變量方差1.統(tǒng)計(jì)數(shù)據(jù)方差---D=var(X,1)功效：當(dāng)X為向量時(shí)，輸出一個(gè)標(biāo)量；當(dāng)X為矩陣時(shí)，輸出為行向量，對應(yīng)于矩陣每列方差值；所以計(jì)算矩陣全部數(shù)方差值，應(yīng)用嵌套：var(var(X))缺省1，計(jì)算：不然計(jì)算：2.統(tǒng)計(jì)數(shù)據(jù)標(biāo)準(zhǔn)差---S=std(X,1)功效：使用方法和1解釋同上3.普通隨機(jī)變量方差----DX=E(X2)-(EX)2功效：用積分或級數(shù)編程計(jì)算10/10/88黃建華制作第89頁4.2.3隨機(jī)變量數(shù)字特征

例4.2.7設(shè)隨機(jī)變量X分布密度為：求隨機(jī)變量X期望和方差。程序：clear;symsx;fx=2/pi*(cos(x))^2;EX=int(x*fx,x,-pi/2,pi/2)E2X=int(x^2*fx,x,-pi/2,pi/2)DX=E2X-EX^2

10/10/89黃建華制作第90頁4.2.3隨機(jī)變量數(shù)字特征常見分布期望和方差1.二項(xiàng)分布---[E,D]=binostat(n,p)說明:n,p能夠是標(biāo)量,向量,矩陣,則E,D是對應(yīng)標(biāo)量,向量,矩陣2.超幾何分布---[E,D]=hygestat(M,N,K)3.泊松分布---[E,D]=poissstat(lambda)4.均勻分布---[E,D]=unifstat(a,b)5.指數(shù)分布---[E,D]=expstat(lambda)6.正態(tài)分布---[E,D]=normstat(mu,sigma)其它：gamstat(),tstat(),fstat(),chi2stat()等等10/10/90黃建華制作第91頁4.2.3隨機(jī)變量數(shù)字特征協(xié)方差與相關(guān)系數(shù)計(jì)算1.隨機(jī)變量協(xié)方差---cov(X,Y)=E(X-EX)(Y-EY)2.隨機(jī)變量相關(guān)系數(shù)---ρ=cov(X,Y)/sqrt(DX*DY)3.統(tǒng)計(jì)數(shù)據(jù)協(xié)方差cov(X)-----當(dāng)X為向量時(shí),cov(X)=var(X);當(dāng)X為矩陣時(shí),結(jié)果為X協(xié)方差矩陣.對角線是X每列方差,Xij為X第i列和第j列協(xié)方差值。cov(X,Y)-----計(jì)算向量X和Y協(xié)方差值4.統(tǒng)計(jì)數(shù)據(jù)相關(guān)系數(shù)corrcoef(X),corrcoef(X,Y)-----說明與使用方法與cov()相同10/10/91黃建華制作第92頁4.2.3隨機(jī)變量數(shù)字特征矩計(jì)算1.隨機(jī)變量k階中心矩---Bk=moment(X,k)2.隨機(jī)變量k階原點(diǎn)矩---Ak=sum(X.^k)/length(X)10/10/92黃建華制作第93頁4.2.4參數(shù)預(yù)計(jì)慣用分布參數(shù)預(yù)計(jì)1.正態(tài)分布參數(shù)預(yù)計(jì)格式：[muhat,sigmahat,muci,sigmaci]=normfit(X,alpha)功效：數(shù)組X服從正態(tài)分布,給定顯著水平alpha,缺省時(shí)為0.05,前二項(xiàng)給出點(diǎn)預(yù)計(jì)，后二項(xiàng)給出區(qū)間預(yù)計(jì)。X為矩陣時(shí),針對列進(jìn)行計(jì)算。2.二項(xiàng)分布參數(shù)預(yù)計(jì)(n重已知,p未知)格式：[phat,puci]=binofit(X,n,alpha)3.泊松分布參數(shù)預(yù)計(jì)格式：[lbdhat,lbdci]=poissfit(X,alpha)4.均勻分布參數(shù)預(yù)計(jì)格式：[ahat,bhat,aci,bci]=unifit(X,alpha)10/10/93黃建華制作第94頁4.2.4參數(shù)預(yù)計(jì)5.指數(shù)分布參數(shù)預(yù)計(jì)格式：[lbdhat,lbdci]=expfit(X,alpha)6.通用命令mle()格式：[輸出參數(shù)項(xiàng)]=mle('分布函數(shù)名',X,alpha[,N])說明：分布函數(shù)名有：bino(二項(xiàng)),geo(幾何),hyge(超幾何)poiss(泊松),uinf(均勻),unid(離散均勻),exp(指數(shù))norm(正態(tài)),t(T分布),f(F分布),beta(貝塔),gam(伽嗎)N當(dāng)二項(xiàng)時(shí)需要，其它沒有。10/10/94黃建華制作第95頁

例4.2.8設(shè)生成一組均值為15,方差為2.52正態(tài)分布隨機(jī)數(shù)據(jù)，然后對這組數(shù)據(jù)進(jìn)行置信度97%參數(shù)預(yù)計(jì)。程序：clear;w=normrnd(15,2.5,50,1);或w=15+2.5*randn(50,1);alpha=0.03;[mh,sh,mc,sc]=normfit(w,alpha)運(yùn)行一次：mh=15.1076sh=2.4038mc=14.3478~15.8674sc=1.9709~3.07034.2.4參數(shù)預(yù)計(jì)10/10/95黃建華制作第96頁

例4.2.9設(shè)從一大批產(chǎn)品中抽取100個(gè)產(chǎn)品，經(jīng)檢驗(yàn)知有60個(gè)一級品，求這批產(chǎn)品一級品率(置信度95%)。程序：clear;alpha=0.05;N=100;X=60;[Ph,Pc]=mle('bino',X,alpha,N)運(yùn)行一次：Ph=0.6000

Pc=0.4972~0.6967

4.2.4參數(shù)預(yù)計(jì)10/10/96黃建華制作第97頁4.2.5假設(shè)檢驗(yàn)單正態(tài)總體均值假設(shè)檢驗(yàn)1.方差已知(u檢驗(yàn)或z檢驗(yàn))格式：[H,P,ci,Zval]=ztest(X,Mu,sigma,alpha,tail)功效：對正態(tài)分布總體采樣X進(jìn)行Z檢驗(yàn)，判斷采樣均值在已知標(biāo)準(zhǔn)差sigma下是否等于假設(shè)值Mu；給定顯著水平alpha,缺省時(shí)為0.05；tail是假設(shè)備選項(xiàng)(即備擇假設(shè))，有三個(gè)值：tail=0是默認(rèn)值，可省略,說明備選項(xiàng)為"均值不等于M"；tail=1,說明備選項(xiàng)為"均值大于M"；tail=-1,說明備選項(xiàng)為"均值小于M"。H=0說明接收原假設(shè)，H=1拒絕原假設(shè)；P為假設(shè)成立概率，P值非常小時(shí)對假設(shè)置疑；ci給出均值置信；Zval給出統(tǒng)計(jì)量值。10/10/97黃建華制作第98頁4.2.5假設(shè)檢驗(yàn)

例4.2.10某面粉廠包裝車間包裝面粉，每袋面粉重量服從正態(tài)分布，機(jī)器正常運(yùn)轉(zhuǎn)時(shí)每袋面粉重量均值為50kg，標(biāo)準(zhǔn)差1。某日隨機(jī)抽取了9袋，重量分別為：49.7，50.6，51.8，52.4，49.8，51.1，52，51.5，51.2機(jī)器運(yùn)轉(zhuǎn)是否正常？程序：clear;x=[49.7,50.6,51.8,52.4,49.8,51.1,52,51.5,51.2];sigma=1;mu=50;[h,p,ci,z]=ztest(x,mu,sigma)結(jié)果：h=1%拒絕原假設(shè)即認(rèn)為機(jī)器不正常p=7.6083e-004%p=0.00076083很小，對假設(shè)置疑ci=50.468951.7755%均值偏高z=3.366710/10/98黃建華制作第99頁4.2.5假設(shè)檢驗(yàn)單正態(tài)總體均值假設(shè)檢驗(yàn)2.方差未知(t檢驗(yàn))格式：[H,P,ci,stats]=ttest(X,Mu,alpha,tail)功效：對正態(tài)分布總體采樣X進(jìn)行t檢驗(yàn)，對H,Mu,alpha,tail,P,ci解釋同上；stats是個(gè)結(jié)構(gòu)，包含三個(gè)元素：tstat(統(tǒng)計(jì)值)、df(自由度)和sd(樣本標(biāo)準(zhǔn)差)。

例4.2.11某燈泡廠出廠標(biāo)準(zhǔn)是壽命不少于小時(shí)，現(xiàn)隨機(jī)從該廠生產(chǎn)一批燈泡中抽取了20只，壽命分別為：1558，1627，2101，1786，1921，1843，1655，16751935，1573，2023，1968，1606，1751，1511，12472076，1685，1905，1881假設(shè)燈泡壽命服從正態(tài)分布問這批燈泡是否到達(dá)了出廠標(biāo)準(zhǔn)？(a=0.01)10/10/99黃建華制作第100頁4.2.5假設(shè)檢驗(yàn)原假設(shè)H0：x≥備擇假設(shè)H1：x<程序：clear;x=[1558,1627,2101,1786,1921,1843,1655,1675,1935,1573,2023,1968,1606,1751,1511,1247,2076,1685,1905,1881];alpha=0.01;mu=;[h,p,ci,stats]=ttest(x,mu,alpha,-1)結(jié)果：h=1%拒絕原假設(shè)即認(rèn)為不符合出廠標(biāo)準(zhǔn)p=5.9824e-005%p很小，對假設(shè)置疑ci=1.0e+003*-Inf1.8895%均值偏低stats=tstat:-4.8176df:19sd:216.897310/10/100黃建華制作第101頁4.2.5假設(shè)檢驗(yàn)雙正態(tài)總體均值假設(shè)檢驗(yàn)比較兩個(gè)方差相等正態(tài)總體均值是否相等(T檢驗(yàn))格式：[H,P,ci,stats]=ttest2(X,Y,alpha,tail)功效：對兩個(gè)正態(tài)分布總體采樣X、Y進(jìn)行T檢驗(yàn)，對H,P,alpha解釋同上；tail是假設(shè)備選項(xiàng)(即備擇假設(shè))，有三個(gè)值：tail=0是默認(rèn)值，可省略,說明備選項(xiàng)為"均值不相等"；tail=1,說明備選項(xiàng)為"X均值大于Y均值"；tail=-1,說明備選項(xiàng)為"X均值小于Y均值"。ci給出均值差置信區(qū)間；stats是個(gè)結(jié)構(gòu)，包含三個(gè)元素：tstat(統(tǒng)計(jì)值)、df(自由度)和sd(標(biāo)準(zhǔn)差Sw)。10/10/101黃建華制作第102頁4.2.5假設(shè)檢驗(yàn)程序：clear;x=[2461,2404,2407,2439,2394,2401,2543,2463,2392,2458];y=[2496,2485,2538,2596,2556,2582,2494,2528,2537,2492];alpha=0.01;[h,p,ci,st]=ttest2(x,y,alpha,-1)結(jié)果：h=1%拒絕原假設(shè)即認(rèn)為壽命未提升p=6.3361e-005%p很小，對假設(shè)置疑ci=-Inf-44.6944st=tstat:-4.8567df:18sd:43.3705

例4.2.12某燈泡廠在采取一項(xiàng)新工藝前后，分別抽取了10只進(jìn)行壽命試驗(yàn)，壽命分別為：舊燈泡：2461,2404,2407,2439,2394,2401,2543,2463,2392,2458新燈泡：2496,2485,2538,2596,2556,2582,2494,2528,2537,2492假設(shè)燈泡壽命服從正態(tài)分布，能否定為采取新工藝后，燈泡壽命提升了？(a=0.01)10/10/102黃建華制作第103頁4.2.5假設(shè)檢驗(yàn)兩個(gè)總體一致性假設(shè)檢驗(yàn)比較兩個(gè)不知道確切分布總體均值是否相等格式：[P,H,stats]=ranksum(X,Y,alpha)功效：對兩個(gè)總體采樣X、Y進(jìn)行檢驗(yàn)，對H,P,alpha解釋同上；stats是個(gè)結(jié)構(gòu)，包含二個(gè)元素：zval(均值差正態(tài)統(tǒng)計(jì)值)和ranksum(統(tǒng)計(jì)秩和值)。10/10/103黃建華制作第104頁4.2.5假設(shè)檢驗(yàn)程序：clear;x=[33.592,33.862,33.751,33.673,33.847,...33.778,33.631,33.911,33.785,33.928];y=[34.221,33.947,33.856,34.039,34.000,33.924,34.125,...34.273,33.968,33.923];alpha=0.05;[p,h,st]=ranksum(x,y,alpha)結(jié)果：p=7.6854e-004%p很小，對假設(shè)置疑h=1%拒絕原假設(shè)即認(rèn)為直徑?jīng)]有顯著不一樣st=zval:-3.3639ranksum:60

例4.2.13兩臺機(jī)床加工同一個(gè)軸，抽樣測量產(chǎn)品直徑(mm):機(jī)床甲：33.592,33.862,33.751,33.673,33.847,33.778,33.631,33.911,33.785,33.928機(jī)床乙：34.221,33.947,33.856,34.039,34.000,33.924,34.125,34.273,33.968,33.923在a=0.05下能否定為兩臺機(jī)床加工直徑?jīng)]有顯著不一樣？10/10/104黃建華制作第105頁4.2.5假設(shè)檢驗(yàn)兩個(gè)樣本含有相同連續(xù)分布假設(shè)檢驗(yàn)檢驗(yàn)兩個(gè)樣本是否含有相同連續(xù)分布格式：[

H

,P,ksstat]=kstest2(X,Y,alpha,tail)功效：對兩個(gè)總體采樣X、Y進(jìn)行檢驗(yàn)，對H,P,alpha解釋同上；tail是假設(shè)備選項(xiàng)(即備擇假設(shè))，有三個(gè)值：tail=0是默認(rèn)值，可省略,說明備選項(xiàng)為"不相等"；tail=1,說明備選項(xiàng)為"大于"；tail=-1,說明備選項(xiàng)為"小于"。ksstat表示測試統(tǒng)計(jì)量值。10/10/105黃建華制作第106頁4.2.5假設(shè)檢驗(yàn)程序：clear;x=randn(1,10);y=randn(1,10)+x;[h,p,kst]=kstest2(x,y)結(jié)果：h=0%接收原假設(shè)即認(rèn)為兩樣本含有相同類型連續(xù)分布p=0.6751%表示假設(shè)成立概率為0.6751kst=0.3000

例4.2.1410/10/106黃建華制作第107頁4.2.5假設(shè)檢驗(yàn)正態(tài)分布假設(shè)檢驗(yàn)檢驗(yàn)樣本是否含有某種連續(xù)分布命令1：[H,P,jbstat,cv]=jbtest(X,alpha)功效：對采樣X進(jìn)行檢驗(yàn)是否服從正態(tài)分布，對H,P,alpha解釋同上；jbstat表示測試統(tǒng)計(jì)量值;cv為是否拒絕假設(shè)臨界值。適合大樣本命令2：[H,P,lstat,cv]=lillietest(X,alpha)功效：對采樣X進(jìn)行檢驗(yàn)是否服從正態(tài)分布，對H,P,alpha解釋同上；jbstat表示測試統(tǒng)計(jì)量值;cv為是否拒絕假設(shè)臨界值。適合小樣本10/10/107黃建華制作第108頁4.2.5假設(shè)檢驗(yàn)程序：clear;m1=ones(1,11)*2.55;m2=ones(1,12)*2.65;m3=ones(1,17)*2.75;m4=ones(1,19)*2.85;m5=ones(1,26)*2.95;m6=ones(1,24)*3.05;m7=ones(1,22)*3.15;m8=ones(1,19)*3.25;m9=ones(1,13)*3.35;M=[m1,m2,m3,m4,m5,m6,m7,m8,m9];[h,p,lst,cv]