X~2分布t分布和f分布地性質(zhì)及近似計算

1、實(shí)用標(biāo)準(zhǔn)文案以下為三個分布的動畫演示程序，直接運(yùn)行：function dandn(action)% % % %X2 分布、 t 分布和 f 分布的性質(zhì)及近似計算% 制作人： 幸瑋 石巖 顧秋容%All Rights Reserved% 2008.12% % % 無輸入?yún)?shù)時進(jìn)行初始程序 ,否則調(diào)界面參數(shù) if nargin1,action= initialize ;else%看是哪種分布data=get(gco, UserData );pdflabel=data(4);精彩文檔實(shí)用標(biāo)準(zhǔn)文案%圖形顯示參數(shù)，調(diào)試時用 switch pdflabelcase 1xstart=0;%x 軸起始點(diǎn)wid

2、th=120;%x 軸范圍xdis=30;%x 軸刻度間距l(xiāng)enth=0.24;%y 軸范圍ydis=0.03;%y 軸刻度間距speed=0.6;%動畫演進(jìn)速率case 2width=6; xdis=3; lenth=0.42; ydis=0.05; speed=0.5;case 3width=2; xdis=0.5; lenth=4;ydis=0.4; speed=0.8;endend%初始化模塊if strcmp(action, initialize )%數(shù)據(jù)初始化%分布函數(shù)選擇pdflabel=0;d=inputdlg( 輸入分布函數(shù)類型 :, 選擇分布函數(shù) ,1, x2 ); if

3、strcmp(d, x2 )pdflabel=1;elseif strcmp(d, t )pdflabel=2;精彩文檔實(shí)用標(biāo)準(zhǔn)文案elseif strcmp(d, f )pdflabel=3;end%不符合的輸入報錯if pdflabel=0errordlg( 分布函數(shù)只能選擇 x2,t 與 f) returnend%整體用戶界面建立，圖形顯示窗口嵌入 figNumber=figure( .Name ,Changing of the X2 destribution function NumberTitle ,off , .Visible ,off );axes( .Units ,normal

4、ized , .Visible ,off , .Position ,0.06 0.065 0.70 0.9);%顯示初始現(xiàn)實(shí)畫面switch pdflabelcase 1%初始顯示參數(shù)xstart=0; %x 軸起始點(diǎn) width=120;%x 軸范圍xdis=30;%x 軸刻度間距l(xiāng)enth=0.24;%y 軸范圍ydis=0.03; %y 軸刻度間距%畫背景圖plot(xstart xstart+width,lenth lenth,b );hold onplot(xstart xstart+width,lenth lenth,r );plot(xstart xstart+width,len

5、th lenth,black);legend( X2 分布 , 正態(tài)分布 );精彩文檔實(shí)用標(biāo)準(zhǔn)文案axis(xstart xstart+width 0 lenth);grid onset(gca, XTick ,xstart:xdis:xstart+width); set(gca, YTick ,0:ydis:lenth);hold off%分布函數(shù)動畫初始參數(shù)m=0;n=1;case 2width=6;xdis=3;lenth=0.42;ydis=0.05;b );r- ); black );plot(-width -width width,0 lenth lenth, hold onplo

6、t(-width -width width,0 lenth lenth, plot(-width -width width,0 lenth lenth, legend( T 分布 ,正態(tài)分布 );grid onaxis(-width width 0 lenth); set(gca, XTick ,-width:xdis:width); set(gca, YTick ,0:ydis:lenth);hold off m=0;n=2;otherwise width=2; xdis=0.5; lenth=4;ydis=0.4;plot(width,0, b ); hold on精彩文檔實(shí)用標(biāo)準(zhǔn)文案pl

7、ot(width,0, r );legend( f 分布 ,正態(tài)分布 );axis(0 width 0 lenth);grid onset(gca, XTick ,0:xdis:width);set(gca, YTick ,0:ydis:lenth);hold off%f 分布輸入 nn=inputdlg( 輸入 n 的值:,獲取 n ,1,5 ); n=str2double(n1);if n=4errordlg( 要使正態(tài)逼近式成立 n 必須大于 4!) returnendm=1;end%按鈕顯示位置top=0.95;left=0.80;btnWid=0.15;btnHt=0.08;spac

8、ing=0.02;%“ Start ”按鈕btnNumber=1;labelStr= Start ;callbackStr= dandn(start) ;%顯示按鈕btnPos=left top-btnHt btnWid btnHt;uicontrol( .Style ,pushbutton , .Units ,normalized , .精彩文檔實(shí)用標(biāo)準(zhǔn)文案Position ,btnPos, .String ,labelStr, .Callback ,callbackStr, .UserData ,btnNumber m n pdflabel);%調(diào)用時的句柄參數(shù)傳遞%關(guān)閉按鈕yPos=to

9、p-(btnHt+spacing);uicontrol( .Style ,push , .Units ,normalized , .Position ,left yPos-btnHt btnWid btnHt,.String ,Close , .Callback ,close(gcf) );set(figNumber, Visible ,on ); % 顯示初始圖像%“ Start ”按鈕執(zhí)行程序 elseif strcmp(action, start )%參數(shù)初始化 data=get(gco, UserData );method=data(1); m=data(2); n=data(3);%

10、“ Stop ”按鈕設(shè)置，點(diǎn)擊開始后替換“ Start ”按鈕 btnNumber=2;labelStr= Stop ; callbackStr= dandn(stop) ;top=0.95;left=0.80; btnWid=0.15;btnHt=0.08;精彩文檔實(shí)用標(biāo)準(zhǔn)文案btnPos=left top-btnHt btnWid btnHt;stophld=uicontrol( .Style ,pushbutton , .Units ,normalized , .Position ,btnPos, .String ,labelStr, .Callback ,callbackStr, .U

11、serData ,btnNumber m n pdflabel);%畫圖程序switch pdflabel case 1while (method=1)%當(dāng)未點(diǎn)擊“ Stop ”則繼續(xù)執(zhí)行%畫 x2 分布密度函數(shù) a=n/2; b=2;x=xstart:0.1:xstart+width;y=gampdf(x,a,b); plot(x,y, b ,EraseMode ,normal ); hold on%畫近似正太分布密度函數(shù) y1=normpdf(x,n,(2*n)0.5); plot(x,y1, r- );%圖形參數(shù)調(diào)節(jié) axis(xstart xstart+width 0 lenth);

12、set(gca, XTick ,xstart:xdis:xstart+width);set(gca, YTick ,0:ydis:lenth); grid on ;legend( X2 分布 , 正態(tài)分布 ); %顯示當(dāng)前 ntemp=num2str(floor(n); str=strcat( n= ,temp);text(92,0.2,str); %顯示當(dāng)前最大間距 dis=num2str(max(abs(y-y1); str=strcat( dis= ,dis);精彩文檔實(shí)用標(biāo)準(zhǔn)文案text(92,0.187,str);%動畫指令drawnow;hold off ;n=n+speed;%看

13、 stop 是否按下data=get(gco, UserData ); method=data(1);%更新傳遞參數(shù)set(stophld, Userdata ,btnNumber m n pdflabel); endcase 2while (method=1)%畫 t 分布密度函數(shù) x=-width:0.1:width;y=tpdf(x,n); plot(x,y, EraseMode ,normal );hold on %畫近似正太分布密度函數(shù) y1=normpdf(x,0,(n/(n-2)0.5); plot(x,y1, r- );%圖形參數(shù)調(diào)節(jié) axis(-width width 0 l

14、enth); set(gca, XTick ,-width:xdis:width);set(gca, YTick ,0:ydis:lenth);grid on ;legend( T 分布 ,正態(tài)分布 );%顯示當(dāng)前 n temp=num2str(floor(n); z=strcat( n= ,temp);text(3.3,0.335,z);%顯示當(dāng)前最大間距 dis=num2str(max(abs(y-y1); str=strcat( dis= ,dis); text(3.3,0.315,str);%動畫指令 drawnow;精彩文檔實(shí)用標(biāo)準(zhǔn)文案hold off ; n=n+speed;%看

15、stop 是否按下 data=get(gco, UserData ); method=data(1); %更新傳遞參數(shù) set(stophld, Userdata ,btnNumber m n pdflabel); end otherwisewhile (method=1)%畫 f 分布密度函數(shù) x=0:0.01:width;y=fpdf(x,m,n); plot(x,y, EraseMode ,normal ); hold on%畫近似正太分布密度函數(shù) y1=normpdf(x,n/(n-2),(2*n2*(m+n-2)/m/(n-2)2/(n-4)0.5); plot(x,y1, r);%

16、圖形參數(shù)調(diào)節(jié) axis(0 width 0 lenth); set(gca, XTick ,0:xdis:width); set(gca, YTick ,0:ydis:lenth); grid on ;legend( f 分布 ,正態(tài)分布 );%顯示當(dāng)前 m temp1=num2str(floor(m); str=strcat( m= ,temp1); text(1.55,3.32,str);%顯示當(dāng)前 n temp2=num2str(floor(n); str=strcat( n= ,temp2); text(1.78,3.32,str);%顯示當(dāng)前最大間距 dis=num2str(max(

17、abs(y-y1); str=strcat( dis= ,dis); text(1.55,3,str);%動畫指令精彩文檔實(shí)用標(biāo)準(zhǔn)文案drawnow;hold off ;m=m+speed;%看 stop 是否按下data=get(gco, UserData ); method=data(1);%更新傳遞參數(shù)set(stophld, Userdata ,btnNumber m n pdflabel); endend%“Close ”按鈕執(zhí)行程序 elseif strcmp(action, stop )data=get(gco, UserData );% 讀原有傳遞數(shù)據(jù)top=0.95;left

18、=0.80; btnWid=0.15;btnHt=0.08;%用“ continue ”按鈕替代“ Stop btnNumber=1;labelStr= Continue ; callbackStr= dandn(continue) ;data(1)=btnNumber; % 更新按鈕標(biāo)識btnPos=left top-btnHt btnWid btnHt; uicontrol( .Style ,pushbutton , .Units ,normalized , . Position ,btnPos, .String ,labelStr, .Callback ,callbackStr, .精彩

19、文檔實(shí)用標(biāo)準(zhǔn)文案UserData ,data);%“ continue ”按鈕執(zhí)行程序elseif strcmp(action, continue )%讀取傳遞數(shù)據(jù)data=get(gco, UserData );method=data(1);m=data(2);n=data(3);%用“ Stop ”按鈕替代“ continue ”btnNumber=2;labelStr= Stop ;callbackStr= dandn(stop) ;top=0.95;left=0.80;btnWid=0.15;btnHt=0.08;btnPos=left top-btnHt btnWid btnHt;c

20、onthld=uicontrol( .Style ,pushbutton , .Units ,normalized , .Position ,btnPos, .String ,labelStr, .Callback ,callbackStr, .UserData ,btnNumber m n pdflabel);%動畫switch pdflabelcase 1while (method=1)精彩文檔實(shí)用標(biāo)準(zhǔn)文案%畫系 x2 分布密度函數(shù)a=n/2; b=2;% 畫 x2 分布密度函數(shù)x=xstart:0.1:xstart+width;y=gampdf(x,a,b);plot(x,y, b ,E

21、raseMode ,normal ); hold on%畫近似正太分布密度函數(shù) y1=normpdf(x,n,(2*n)0.5); plot(x,y1, r- );%圖形參數(shù)調(diào)節(jié)axis(xstart xstart+width 0 lenth); set(gca, XTick ,xstart:xdis:xstart+width);set(gca, YTick ,0:ydis:lenth);grid on ;legend( X2 分布 , 正態(tài)分布 );%顯示當(dāng)前 n temp=num2str(floor(n); str=strcat( n= ,temp);text(92,0.2,str); %

22、顯示當(dāng)前最大間距 dis=num2str(max(abs(y-y1);str=strcat( dis= ,dis); text(92,0.187,str);%動畫指令drawnow;hold off n=n+speed;%看 stop 是否按下data=get(gco, UserData ); method=data(1);%更新傳遞參數(shù)set(conthld, UserData ,btnNumber m n pdflabel); endcase 2while (method=1)x=-width:0.1:width;y=tpdf(x,n);plot(x,y, EraseMode ,norma

23、l );精彩文檔實(shí)用標(biāo)準(zhǔn)文案hold on %畫近似正太分布密度函數(shù) y1=normpdf(x,0,(n/(n-2)0.5); plot(x,y1, r- );%圖形參數(shù)調(diào)節(jié) axis(-width width 0 lenth); set(gca, XTick ,-width:xdis:width);set(gca, YTick ,0:ydis:lenth); grid on ;legend( T 分布 ,正態(tài)分布 ); %顯示當(dāng)前 n temp=num2str(floor(n); z=strcat( n= ,temp); text(3.3,0.335,z);%顯示當(dāng)前最大間距 dis=num2str(max(abs(y-y1); str=strcat( dis= ,dis); text(3.3,0.315,str);%動畫指令drawnow;hold offn=n+speed;%看 stop 是否按下 data=get(