eye=ones(1,30); eye=eye.';

x=[eye x_1 x_2 x_4];

[b,bint,r,rint,stats]=regress(y,x,0.05) b =

17.3244 1.3070 -3.6956 0.3486

bint =

5.7282 0.6829 -7.4989 0.0379 r =

-0.0988 -0.0795 -0.1195 -0.0441 0.4660 -0.0133 0.2912 0.2735 -0.2351 0.1031 -0.4033 0.1747 0.0400 -0.1504 0.1284 0.1637 -0.0527 -0.1907 -0.0870 -0.0165 -0.1292 -0.3002 -0.2933 -0.1679 -0.2177 0.1116 0.3035 0.0693 0.2474 0.2270 rint =

-0.5270 -0.5309 -0.5106 -0.4731 0.0813 28.9206 1.9311 0.1077 0.6594 0.3294 0.3718 0.2716 0.3848 0.8507

第 41 页 共 46 页

-0.4609 0.4343 -0.1374 0.7197 -0.0870 0.6340 -0.5960 0.1258 -0.3280 0.5341 -0.8190 0.0125 -0.2618 0.6112 -0.4032 0.4832 -0.5933 0.2925 -0.3207 0.5775 -0.2841 0.6116 -0.4830 0.3776 -0.6248 0.2434 -0.5348 0.3609 -0.4423 0.4092 -0.5609 0.3024 -0.7181 0.1177 -0.7243 0.1377 -0.5548 0.2190 -0.6449 0.2095 -0.2994 0.5226 -0.1037 0.7106 -0.3714 0.5099 -0.1807 0.6755 -0.1890 0.6430

stats =

0.9054 82.9409 0.0000 0.0490 参数 参数估计值 // 17.3244 // 1.3070 // -3.6956 // 0.3486 R//=0.9054 F=82.9409 p<0.0001 s//=0.0490 y=17.3244+1.3070x1-3.6956x2+0.3486

参数置信区间 [5.7282，28.9206] [0.6829，1.9311] [-7.4989，0.1077 ] [ 0.0379，0.6594] 九、计算机模拟

（一）实验目的

1.掌握计算机模拟的内容。

2.学会利用软件求解计算机模拟问题。

（二）实验课时

9.1.当天生产的产品必须售出，否则就会变质。该产品单位成本为 2.5 元，单位产品售价为5元。企业为避免存货过多而造成损失，拟从以下两种库存方案中选出一个较优的方案：

第 42 页 共 46 页

方案甲：按前1天的销售量作为当天的库存量；

方案乙：按前2天的平均销售量作为当天的库存量。 假定市场对该产品的每天需求量是一个随机变量，但从以往的统计分析得知它服从正态分布N(135 , 22.42)。

function [TL1,TL2]=kucum(T,S1,S21,S22)TL1=0;TL2=0;N=1;while N<TQ1=S1;Q2=(S21+S22)/2;D=normrnd(135,22.4);if D<Q1; S3=Q1;else S3=D;endif D<Q2 S4=Q2;elseS4=D;endL1=5*S3-2.5*Q1;L2=5*S4-2.5*Q2;TL1=TL1+L1;TL2=TL2+L2;N=N+1;S1

=S3;S22=S21;S21=S4;end TL1<

TL2，所以方案乙较优

9.2.某报童以每份0.03元的价格买进报纸，以0.05元的价格出售．根据长期统计，报纸每天的销售量及百分率为

200 210 220 230 240 250 销售量

0.10 0.20 0.40 0.15 0.10 0.05 百分率

第 43 页 共 46 页

已知当天销售不出去的报纸，将以每份0.02元的价格退还报社．试用模拟方法确定报童每天买进报纸数量，使报童的平均总收入为最大？

主文件main.m：

BUYMIN=200; % 每天的最小购买量 BUYMAX=250; % 每天的最大购买量

SIMUDAY=1.0e+5; % 模拟时间 sell_amount=200:10:250; % 销售量 percentage=[0.1 0.3 0.7 0.85 0.95 1]; % 百分率 buy_amount=0; ave_profit=0;

for loop_buy=BUYMIN:BUYMAX sum_profit=0;

for loop_day=1:SIMUDAY

index=find(percentage>=rand); % 产生随机数，用于决定当天的销售量 sum_profit=sum_profit+GetProfit(loop_buy,sell_amount(index(1))); end

buy_amount=[buy_amount,loop_buy]; % 循环嵌套 ave_profit=[ave_profit,sum_profit/SIMUDAY]; % 循环嵌套 end

buy_amount(1)=[]; % 第一个元素置空 ave_profit(1)=[];

[val,id]=max(ave_profit) % 显示最大平均收入val

buy=buy_amount(id) % 显示在平均收入最大情况下的每天的购买量buy xlabel='每天的购买量'; ylabel='平均利润';

plot(buy_amount,ave_profit,'*:'); % 函数GetProfit.m代码： function re=GetProfit(a,b)

if a

else % 供过于求：报童购买量大于销售量 re=b*(0.05-0.03)+(a-b)*(0.02-0.03); end 运行结果：

第 44 页 共 46 页