二維隨機(jī)變量數(shù)值特征及其可視化

sh_y_f_1.m【clc,cleark=200; %計(jì)算機(jī)模擬次數(shù)z=1.4166; %期望的理論值(單位：萬(wàn)元)xx=zeros(1,k);yy=zeros(1,k);zz=zeros(1,k);figure(1)hndl=plot([00],[0,k/2],':','erasemode','xor');%實(shí)時(shí)動(dòng)畫(huà)要素1holdonaxis([0,2,0,1.5*k]);theAxs=axis;%取出當(dāng)前坐標(biāo)軸，使得后繼圖形保持與當(dāng)前坐標(biāo)相同尺度set(gca,'ytick',[0k/3k/21.5*k],'yticklabel','0|0|0|200');plot([theAxs(1),theAxs(2)],[k/3,k/3],'r:','linewidth',2);%畫(huà)線標(biāo)識(shí)期望理論值plot([z,z],[theAxs(3),theAxs(4)],'r-','linewidth',3); %畫(huà)線標(biāo)識(shí)利潤(rùn)期望fori=1:k %進(jìn)行k次模擬x=unifrnd(10,20);%隨機(jī)產(chǎn)生服從均勻分布U(10,20)的生產(chǎn)量y=unifrnd(10,20);%隨機(jī)產(chǎn)生服從均勻分布U(10，20)的需求量xx(i)=x;yy(i)=y;% %計(jì)算利潤(rùn) ify<=xzz(i)=0.1*y;elsezz(i)=0.05*(x+y);end% % plot([0zz(i)],[i+k/2i+k/2],'linewidth',1); %將點(diǎn)投影到x軸,畫(huà)相應(yīng)的投影示意線plot(zz(i),k/3,'o'); %畫(huà)投點(diǎn)坐標(biāo)的示意點(diǎn)set(hndl,'Xdata',[zz(i)zz(i)],'Ydata',[k/3,i+k/2]);%用當(dāng)前數(shù)據(jù)替換圖形句柄的X和Y數(shù)據(jù)drawnowifi<10 %if塊語(yǔ)句，前幾次模擬慢顯示pause(0.5)endxlabel('試驗(yàn)重復(fù)次數(shù)n')ylabel('每次試驗(yàn)的利潤(rùn)')endhist(zz)figure(2)plot(xx,yy,'r*')xlabelC試驗(yàn)?zāi)M的進(jìn)貨量')ylabel('試驗(yàn)?zāi)M的需求量')holdonplot([1020201010],[1010202010],'m-')axis([822822])】sh_y_f_2.m【function[E_heiE_hongp]=sh_y_f_6clc,clearnum_hei=0;num_hong=0;N=input('請(qǐng)輸入試驗(yàn)次數(shù)N');NN=N;p=zeros(3,4);%分布列whileN>0% 畫(huà)出盒子 axis([018018])holdonplot([516],[33],'LineWidth',2,'color','b');plot([55],[313],'LineWidth',2,'color','b');plot([5,16],[1313],'LineWidth',2,'color','b');plot([1616],[313],'LineWidth',2,'color','b');% 固定求的位置即大小 alpha=0:pi/200:2*pi;r=0.5;xx=[6789101112131415];%第一次抽球前球所在位置的球心橫坐標(biāo)yy=[487594710126]; %第一次抽球前球所在位置的球心縱坐標(biāo)ii=1:10;% 模擬抽球試驗(yàn)% 第一次抽球pp1=fix(10*unifrnd(0,1))+1;%抽出的是第ppi個(gè)球，對(duì)應(yīng)與第xx(ppi)所表示顏色的球[num_heinum_hong]=chouqiu(num_hei,num_hong,xx,yy,r,alpha,ii,ppi);ii=[i:ppi-ippi+i:i0]; %剩余的球的編號(hào)% 第二次抽球 ppii=fix(9*unifrnd(0,i))+i;pp2二ii(pp11); %抽出的是第pp2個(gè)球，對(duì)應(yīng)與第xx(pp2)所表示顏色的球pp=sort([ppipp2]);[num_heinum_hong]=chouqiu(num_hei,num_hong,xx,yy,r,alpha,ii,pp2);ii=[i:pp(i)-ipp(i)+i:pp(2)-ipp(2)+i:i0] ;%剩余的球的編號(hào)% 第三次抽球 pp22=fix(8*unifrnd(0,1))+1;pp3二ii(pp22); %抽出的是第pp2個(gè)球，對(duì)應(yīng)與第xx(pp2)所表示顏色的球pp=sort([pp1pp2pp3]);[num_heinum_hong]=chouqiu(num_hei,num_hong,xx,yy,r,alpha,ii,pp3);ii=[1:pp(1)-1pp(1)+1:pp(2)-1pp(2)+1:pp(3)-1pp(3)+1:10];%剩余的球的編號(hào)N=N-1;clfifnum_hei==0&&num_hong==2p(1,3)=p(1,3)+1;elseifnum_hei==0&&num_hong==3p(1,4)=p(1,4)+1;elseifnum_hei==1&&num_hong==1p(2,2)=p(2,2)+1;elseifnum_hei==1&&num_hong==2p(2,3)=p(2,3)+1;elseifnum_hei==2&&num_hong==0p(3,1)=p(3,1)+1;elseifnum_hei==2&&num_hong==1p(3,2)=p(3,2)+1;enddisp('第k此模擬結(jié)果：’)disp(['黑球數(shù)為num_hei=',num2str(num_hei)])disp(['紅球數(shù)為num_hong=',num2str(num_hong)])num_hei=0;num_hong=0;endk二sum(p'); %N次試驗(yàn)中出現(xiàn)黑球數(shù)分別為012的次數(shù)r=sum(p);%N次試驗(yàn)中出現(xiàn)紅球數(shù)分別為0123的次數(shù)E_hei=(0*k(1)+1*k(2)+2*k(3))/NN;E_hong=(0*r(1)+1*r(2)+2*r(3)+3*r(4))/NN;p=p/NN;bar3(p,'detached')title('bar3以參數(shù)detached繪制條形圖')% 子函數(shù) function[num_heinum_hong]=chouqiu(num_hei,num_hong,xx,yy,r,alpha,ii,pp)n=length(ii);fori=1:nx=r*cos(alpha)+xx(ii(i));y=r*sin(alpha)+yy(ii(i));holdonplot(x,y); %半徑為r的圓，xO,yO為原點(diǎn)坐標(biāo)axisequalifii(i)==1||ii(i)==5%定義該球?yàn)楹谇騠ill(x,y,'k')elseifii(i)==10%定義該球?yàn)榘浊騟lseendfill(x,y,'w')fill(x,y,'r')%定義該球?yàn)榧t球endhndl=plot(x,y,'k','EraseMode','XOR');ifpp==1||pp==5hndl1=fill(x,y,'k');num_hei=num_hei+1;%定義黑球的函數(shù)句柄elseifpp==10elsehndl2=fill(x,y,'w');%定義白球的函數(shù)句柄endhndl3=fill(x,y,'r');num_hong=num_hong+1;%定義紅球的函數(shù)句柄fort=0:0.5:(16-yy(pp))%可以調(diào)節(jié)速度0.5可變xx0=xx(pp);yy0=yy(pp)+t;x=r*cos(alpha)+xx0;y=r*sin(alpha)+yy0;set(gca,'UserData',hndl);hndl=get(gca,'UserData');set(hndl,'XData',x);set(hndl,'YData',y);ifpp==1||pp==5 %抽出的是黑球set(gca,'UserData',hndl1);hndl1=get(gca,'UserData');set(hndl1,'XData',x);set(hndl1,'YData',y);elseifpp==10 %抽出的是白球set(gca,'UserData',hndl2);hndl2=get(gca,'Use