《X~分布t分布和f分布的性质及近似计算》由会员分享，可在线阅读，更多相关《X~分布t分布和f分布的性质及近似计算（11页珍藏版）》请在金锄头文库上搜索。
1、以下为三个分布的动画演示程序，直接运行：function dandn(action)% % X2 分布、t 分布和 f 分布的性质及近似计算% 制作人： 幸玮 石岩 顾秋容% All Rights Reserved% 2008.12% %无输入参数时进行初始程序,否则调界面参数 if nargin1,action=initialize; else%看是哪种分布data=get(gco,UserData);pdflabel=data(4);%图形显示参数，调试时用switch pdflabelcase 1xstart=0; %x 轴起始点width=120; %x 轴范围xdis=30; %x
2、轴刻度间距lenth=0.24; %y 轴范围ydis=0.03; %y 轴刻度间距speed=0.6; %动画演进速率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;end end%初始化模块 if strcmp(action,initialize)%数据初始化%分布函数选择pdflabel=0;d=inputdlg(输入分布函数类型:,选择分布函数,1,x2);if strcmp(d,x2)pdflabel=1;elseif strc
3、mp(d,t)pdflabel=2;elseif strcmp(d,f)pdflabel=3;end%不符合的输入报错if pdflabel=0errordlg(分布函数只能选择 x2,t 与 f)returnend %整体用户界面建立，图形显示窗口嵌入figNumber=figure( .Name,Changing of the X2 destribution function, .NumberTitle,off, .Visible,off);axes( .Units,normalized, .Visible,off, .Position,0.06 0.065 0.70 0.9);%显示初始
4、现实画面switch pdflabelcase 1%初始显示参数xstart=0; %x 轴起始点width=120; %x 轴范围xdis=30; %x 轴刻度间距lenth=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,lenth lenth,black);legend(X2 分布,正态分布);axis(xstart xstart+width
5、0 lenth);grid onset(gca,XTick,xstart:xdis:xstart+width);set(gca,YTick,0:ydis:lenth);hold off%分布函数动画初始参数m=0;n=1;case 2width=6;xdis=3;lenth=0.42;ydis=0.05;plot(-width -width width,0 lenth lenth,b);hold onplot(-width -width width,0 lenth lenth,r-);plot(-width -width width,0 lenth lenth,black);legend(T
6、分布,正态分布);grid onaxis(-width width 0 lenth);set(gca,XTick,-width:xdis:width);set(gca,YTick,0:ydis:lenth);hold offm=0;n=2;otherwisewidth=2;xdis=0.5;lenth=4;ydis=0.4;plot(width,0,b);hold onplot(width,0,r);legend(f 分布,正态分布);axis(0 width 0 lenth);grid onset(gca,XTick,0:xdis:width);set(gca,YTick,0:ydis:le
7、nth);hold off%f 分布输入 nn=inputdlg(输入 n 的值:,获取 n,1,5);n=str2double(n1);if n=4errordlg(要使正态逼近式成立 n 必须大于 4!)returnendm=1;end%按钮显示位置top=0.95;left=0.80;btnWid=0.15;btnHt=0.08;spacing=0.02;%“Start”按钮btnNumber=1;labelStr=Start;callbackStr=dandn(start);%显示按钮btnPos=left top-btnHt btnWid btnHt;uicontrol( .Styl
8、e,pushbutton, .Units,normalized, .Position,btnPos, .String,labelStr, .Callback,callbackStr, .UserData,btnNumber m n pdflabel); %调用时的句柄参数传递%关闭按钮yPos=top-(btnHt+spacing);uicontrol( .Style,push, .Units,normalized, .Position,left yPos-btnHt btnWid btnHt, .String,Close, .Callback,close(gcf);set(figNumber
9、,Visible,on); %显示初始图像 %“Start”按钮执行程序 elseif strcmp(action,start)%参数初始化data=get(gco,UserData);method=data(1);m=data(2);n=data(3);%“Stop”按钮设置，点击开始后替换“Start”按钮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;stophld=uicontrol
10、( .Style,pushbutton, .Units,normalized, .Position,btnPos, .String,labelStr, .Callback,callbackStr, .UserData,btnNumber m n pdflabel);%画图程序switch pdflabelcase 1while(method=1) %当未点击“Stop”则继续执行%画 x2 分布密度函数a=n/2; b=2; x=xstart:0.1:xstart+width;y=gampdf(x,a,b);plot(x,y,b,EraseMode,normal);hold on%画近似正太分
11、布密度函数y1=normpdf(x,n,(2*n)0.5);plot(x,y1,r-);%图形参数调节axis(xstart xstart+width 0 lenth);set(gca,XTick,xstart:xdis:xstart+width);set(gca,YTick,0:ydis:lenth);grid on;legend(X2 分布,正态分布);%显示当前 ntemp=num2str(floor(n);str=strcat(n=,temp);text(92,0.2,str);%显示当前最大间距dis=num2str(max(abs(y-y1);str=strcat(dis=,dis
12、);text(92,0.187,str);%动画指令drawnow;hold off;n=n+speed;%看 stop 是否按下data=get(gco,UserData);method=data(1);%更新传递参数set(stophld,Userdata,btnNumber m n pdflabel);endcase 2while(method=1)%画 t 分布密度函数x=-width:0.1:width;y=tpdf(x,n);plot(x,y,EraseMode,normal);hold on%画近似正太分布密度函数y1=normpdf(x,0,(n/(n-2)0.5);plot(
13、x,y1,r-);%图形参数调节axis(-width width 0 lenth);set(gca,XTick,-width:xdis:width);set(gca,YTick,0:ydis:lenth);grid on;legend(T 分布,正态分布);%显示当前 ntemp=num2str(floor(n);z=strcat(n=,temp);text(3.3,0.335,z);%显示当前最大间距dis=num2str(max(abs(y-y1);str=strcat(dis=,dis);text(3.3,0.315,str);%动画指令drawnow;hold off;n=n+spe
14、ed;%看 stop 是否按下data=get(gco,UserData);method=data(1);%更新传递参数set(stophld,Userdata,btnNumber m n pdflabel);endotherwisewhile(method=1)%画 f 分布密度函数x=0:0.01:width;y=fpdf(x,m,n);plot(x,y,EraseMode,normal);hold on%画近似正太分布密度函数y1=normpdf(x,n/(n-2),(2*n2*(m+n-2)/m/(n-2)2/(n-4)0.5);plot(x,y1,r);%图形参数调节axis(0 width 0 lenth);set(gca,XTick,0:xdis:width);set(gca,YTick,0:ydis:lenth);grid on;legend(f 分布,正态分布);%显示当前
