%報童問題(高斯分布)%幾個參數：%賣出報紙掙a元，退回報紙賠b元;%賣報量y,訂報量z，訂報量上限z_up;%觀察期的總天數t,用于積累天數的ti;%總的平均損失量ave_loss，供過于求損失loss1,供小于求損失loss2,ave_loss=loss1+loss2;%最佳訂報量superz,最小損失量ave_loss_min;function [superz,ave_loss_min]=gauss(t,z_up,a,b)   mu=input('請輸入高斯分布的均值mu');                 %輸入高斯分布均值mu   if isempty(mu)       mu=1000;   end   sigma=input('請輸入高斯分布的標準差sigma');         %輸入高斯分布標準差sigma   if isempty(sigma)       sigma=100;   end   y=round(normrnd(mu,sigma,1,t));                   %生成服從高斯分布的隨機賣報量y（矢量）    ave_loss=zeros(1,z_up);                          %各種訂報量情況下的損失統計（z_up個元素的矢量）     for z=1:z_up                                    %外層循環，z_up種情況，各種訂報量         loss1=0;         loss2=0;        for ti=1:t                                   %內層循環，一共實驗多少天          if z>y(ti)                                 %如果供大于求，說明要退報，賠錢               loss1=loss1+b*(z-y(ti));          elseif z<y(ti)                             %如果供小于求，說明還可以賺錢，報紙不夠，少賺了錢               loss2=loss2+a*(y(ti)-z);          else                                       %如果訂報量與賣報量相同，損失不變              loss1=loss1;              loss2=loss2;              end        end         ave_loss(z)=(loss1+loss2)/t;                %求每天的平均損失量     end     ave_loss_min=min(ave_loss)                      %顯示最后的平均最小損失     superz=find(ave_loss==ave_loss_min)             %記錄在給定a和給定b條件下的最佳訂報量          figure(1)                                       %畫給定a和b的時候的最佳訂報量和最小平均損失關系曲線     plot(1:z_up,ave_loss,'b')     grid on     axis([0,1000,0,100])     xlabel('訂報量/份')     ylabel('虧損/元')     title('虧損與訂報量的關系曲線')           figure(2)                                      %每份報紙收益與最佳訂報量的關系曲線     z_a=zeros(1,100);     n1=1;     for a_temp=0:0.01:1                             %收益a由0元遞增至1元         for z=1:z_up                                %外層循環，z_up種情況，各種訂報量           loss1=0;           loss2=0;              for ti=1:t                                %內層循環，一共實驗多少天            if z>y(ti)                               %如果供大于求，說明要退報，賠錢               loss1=loss1+b*(z-y(ti));            elseif z<y(ti)                           %如果供小于求，說明還可以賺錢，報紙不夠，少賺了錢               loss2=loss2+a_temp*(y(ti)-z);            else                                     %如果訂報量與賣報量相同，損失不變              loss1=loss1;              loss2=loss2;            end        end         ave_loss(z)=(loss1+loss2)/t;                %計算平均損失       end     ave_loss_min=min(ave_loss);     z_a(n1)=find(ave_loss==ave_loss_min,1,'first'); %將每種a的最佳訂報量放入z_a矢量      n1=n1+1;     end     a_temp=[0:0.01:1];     plot(a_temp,z_a,'b')     grid on     axis([0,1,0,1000])     xlabel('每份報紙凈賺/元')     ylabel('最佳訂報量/份')     title('每份報紙凈賺與最佳訂報量的關系曲線')          figure(3)                                          %每份報紙賠償與最佳訂報量的關系曲線     z_b=zeros(1,100);     n2=1;     for b_temp=0:0.01:1                                %賠償b由0元遞增至1元         for z=1:z_up                                   %外層循環，z_up種情況，各種訂報量           loss1=0;            loss2=0;              for ti=1:t                                   %內層循環，一共實驗多少天            if z>y(ti)                                  %如果供大于求，說明要退報，賠錢               loss1=loss1+b_temp*(z-y(ti));            elseif z<y(ti)                              %如果供小于求，說明還可以賺錢，報紙不夠，少賺了錢               loss2=loss2+a*(y(ti)-z);            else                                        %如果訂報量與賣報量相同，損失不變              loss1=loss1;              loss2=loss2;            end        end         ave_loss(z)=(loss1+loss2)/t;                   %計算平均損失       end     ave_loss_min=min(ave_loss);     z_b(n2)=find(ave_loss==ave_loss_min,1,'first');    %將每種賠償b的最佳訂報量放入z_b矢量       n2=n2+1;     end     b_temp=[0:0.01:1];     plot(b_temp,z_b,'b')     grid on     axis([0,1,0,1000])     xlabel('每份報紙賠償/元')     ylabel('最佳訂報量/份')     title('每份報紙賠償與最佳訂報量的關系曲線')end