uapaaa第12講-計(jì)算機(jī)模擬-資料下載頁(yè)
2025-08-04 10:07本頁(yè)面
【正文】 mfunction lp c= lp co ns t( x) % 約束條件if 3*x(1)+x (2)10= amp。 3* x(1)+x(2 )10=0 .5 % 約束條件的誤差為? lpc=1。else lpc=0。end% function [sol,r1,r2]=randlp(a,b,n) %隨機(jī)模擬解非線性規(guī)劃 debug=1。 a=0。 %試驗(yàn)點(diǎn)下界 b=10。 %試驗(yàn)點(diǎn)上界 n=1000。 %試驗(yàn)點(diǎn)個(gè)數(shù) r1=unifrnd(a,b,n,1)。 %n?1階的［ a,b］均勻分布隨機(jī)數(shù)矩陣 r2=unifrnd(a,b,n,1)。 sol=[r1(1) r2(1)]。 z0=inf。 for i=1:n x1=r1(i)。 x2=r2(i)。 lpc=lpconst([x1 x2])； if lpc==1 z=mylp([x1 x2])； if zz0 z0=z； sol=[x1 x2]； end end end To Matlab(randlp) 返回 實(shí)驗(yàn)作業(yè) 2. 某報(bào)童以每份 元的價(jià)格買進(jìn)報(bào)紙 , 以 元的價(jià)格出售 . 根據(jù)長(zhǎng)期統(tǒng)計(jì) , 報(bào)紙每天的銷售量及百分率為 銷售量 200 210 220 230 240 250百分率 已知當(dāng)天銷售不出去的報(bào)紙，將以每份 元的價(jià)格退還報(bào)社 . 試用模擬方法確定報(bào)童每天買進(jìn)報(bào)紙數(shù)量 , 使報(bào)童的平均總收入為最大 ?編一個(gè)福利彩票電腦選號(hào)的程序。 4. 某設(shè)備上安裝有四只型號(hào)規(guī)格完全相同的電子管，已知電子管壽命為 10002022小時(shí)之間的均勻分布。當(dāng)電子管損壞時(shí)有兩種維修方案，一是每次更換損壞的那一只；二是當(dāng)其中一只損壞時(shí)四只同時(shí)更換。已知更換時(shí)間為換一只時(shí)需1小時(shí)， 4只同時(shí)換為 2小時(shí)。更換時(shí)機(jī)器因停止運(yùn)轉(zhuǎn)每小時(shí)的損失為 20元，又每只電子管價(jià)格 10元，試用模擬方法決定哪一個(gè)方案經(jīng)濟(jì)合理？ 5. 導(dǎo)彈追蹤問題： 設(shè)位于坐標(biāo)原點(diǎn)的甲艦向位于 x軸上點(diǎn) A(1, 0)處的乙艦發(fā)射導(dǎo)彈，導(dǎo)彈頭始終對(duì)準(zhǔn)乙艦 .如果乙艦以最大的速度(是常數(shù) )沿平行于 y軸的直線行駛，導(dǎo)彈的速度是 5，模擬導(dǎo)彈運(yùn)行的軌跡 .又乙艦行駛多遠(yuǎn)時(shí)，導(dǎo)彈將它擊中？ 初始化 :i=0,k1=0,k2=0,k3=0 i=i+1 骰子 點(diǎn)數(shù) ? k1=k1+1 k2=k2+1 k3=k3+1 k1=k1+1 i＜ 20? E=(k2+k3)/20 E1=0*k1/20+1*k2/20+2*k3/20 停止 硬幣正面 ? Y N N Y 1， 2， 3 4， 5 6 投擲硬幣的計(jì)算機(jī)模擬 產(chǎn)生服從 U（ 0， 1）的隨機(jī)數(shù) R1 將區(qū)間 [0， 1]兩等分： 若 ，則對(duì)應(yīng)硬幣正面 若 ，則對(duì)應(yīng)硬幣反面 1 ?? R 1 ?? R擲骰子的計(jì)算機(jī)模擬 產(chǎn)生服從 U（ 0， 1）的隨機(jī)數(shù) R2 將區(qū)間 [0， 1]六等份： 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 1 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 2 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 3 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 4 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 5 若 ，則對(duì)應(yīng) 骰子點(diǎn)數(shù)為 6 6102 ?? R62612 ?? R63622 ?? R65642 ?? R64632 ?? R165 2 ?? R初始化 :i=0,k1=0,k2=0,k3=0 i=i+1 R2=? k1=k1+1 k2=k2+1 k3=k3+1 k1=k1+1 i＜ 20? E=(k2+k3)/20 E1=0*k1/20+1*k2/20+2*k3/20 停止 R1= Y N N Y R23/6 其它 R25/6 To Matlab(liti1)
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1