二维随机变量数值特征及其可视化

sh_y_f_1.m【clc,cleark=200; %计算机模拟次数z=1.4166; %期望的理论值(单位：万元)xx=zeros(1,k);yy=zeros(1,k);zz=zeros(1,k);figure(1)hndl=plot([00],[0,k/2],':','erasemode','xor');%实时动画要素1holdonaxis([0,2,0,1.5*k]);theAxs=axis;%取出当前坐标轴，使得后继图形保持与当前坐标相同尺度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);%画线标识期望理论值plot([z,z],[theAxs(3),theAxs(4)],'r-','linewidth',3); %画线标识利润期望fori=1:k %进行k次模拟x=unifrnd(10,20);%随机产生服从均匀分布U(10,20)的生产量y=unifrnd(10,20);%随机产生服从均匀分布U(10，20)的需求量xx(i)=x;yy(i)=y;% %计算利润 ify<=xzz(i)=0.1*y;elsezz(i)=0.05*(x+y);end% % plot([0zz(i)],[i+k/2i+k/2],'linewidth',1); %将点投影到x轴,画相应的投影示意线plot(zz(i),k/3,'o'); %画投点坐标的示意点set(hndl,'Xdata',[zz(i)zz(i)],'Ydata',[k/3,i+k/2]);%用当前数据替换图形句柄的X和Y数据drawnowifi<10 %if块语句，前几次模拟慢显示pause(0.5)endxlabel('试验重复次数n')ylabel('每次试验的利润')endhist(zz)figure(2)plot(xx,yy,'r*')xlabelC试验模拟的进货量')ylabel('试验模拟的需求量')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('请输入试验次数N');NN=N;p=zeros(3,4);%分布列whileN>0% 画出盒子 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];%第一次抽球前球所在位置的球心横坐标yy=[487594710126]; %第一次抽球前球所在位置的球心纵坐标ii=1:10;% 模拟抽球试验% 第一次抽球pp1=fix(10*unifrnd(0,1))+1;%抽出的是第ppi个球，对应与第xx(ppi)所表示颜色的球[num_heinum_hong]=chouqiu(num_hei,num_hong,xx,yy,r,alpha,ii,ppi);ii=[i:ppi-ippi+i:i0]; %剩余的球的编号% 第二次抽球 ppii=fix(9*unifrnd(0,i))+i;pp2二ii(pp11); %抽出的是第pp2个球，对应与第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] ;%剩余的球的编号% 第三次抽球 pp22=fix(8*unifrnd(0,1))+1;pp3二ii(pp22); %抽出的是第pp2个球，对应与第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];%剩余的球的编号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此模拟结果：’)disp(['黑球数为num_hei=',num2str(num_hei)])disp(['红球数为num_hong=',num2str(num_hong)])num_hei=0;num_hong=0;endk二sum(p'); %N次试验中出现黑球数分别为012的次数r=sum(p);%N次试验中出现红球数分别为0123的次数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以参数detached绘制条形图')% 子函数 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为原点坐标axisequalifii(i)==1||ii(i)==5%定义该球为黑球fill(x,y,'k')elseifii(i)==10%定义该球为白球elseendfill(x,y,'w')fill(x,y,'r')%定义该球为红球endhndl=plot(x,y,'k','EraseMode','XOR');ifpp==1||pp==5hndl1=fill(x,y,'k');num_hei=num_hei+1;%定义黑球的函数句柄elseifpp==10elsehndl2=fill(x,y,'w');%定义白球的函数句柄endhndl3=fill(x,y,'r');num_hong=num_hong+1;%定义红球的函数句柄fort=0:0.5:(16-yy(pp))%可以调节速度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