【建模教程】-統(tǒng)計建模與R軟件

----宋停云與您分享--------宋停云與您分享----學(xué)Th學(xué)號學(xué)Th學(xué)號實驗課成績武漢理工大學(xué)學(xué)生實驗報告書實驗課程名稱開課學(xué)院指導(dǎo)老師姓名學(xué)Th姓名學(xué)Th專業(yè)班級201 — 201 學(xué)年 第 學(xué)期----宋停云與您分享--------宋停云與您分享----MonteCarloa=0.05，0.01（1）隨機投點法>MC1=function(n){+k=0;x=runif(n,0,1);y=runif(n,0,sqrt(2))+for(iin1:n){+if(y[i]<sqrt(1+x[i]^2))+k=k+1+}+k/n+}>p1=MC1(100000);[1]0.81134>s1=p1*sqrt(2)>s1[1]1.147408>n1=p1*(1-p1)*(1.96^2)/(0.01^2)>n1[1]5880.237由上可以得到，積分值s1=1.147408，試驗次數(shù)n1=5880.237（2）平均值法----宋停云與您分享--------宋停云與您分享---->MC2=function(n){+x=runif(n)+sum(sqrt(1+x^2))/n+}>p2=MC2(100000)>p2[1]1.147071>s2=p2>s2[1]1.147071>n2=(4/3-(p2^2))*(1.96^2)/(0.01^2)>n2[1]674.6593由上可以得到，積分值s2=1.147071，試驗次數(shù)n2=674.6593O20B33A>##畫出A,B,O三點位置>plot(c(0,0,33),c(0,20,0),xlab="",ylab="")>text(33,0,labels="A",adj=c(0.5,-0.5))>text(0,20,labels="B",adj=c(-.5,0.3))----宋停云與您分享--------宋停云與您分享---->text(0,0,labels="O",adj=c(1.3,0.3))>lines(c(33,0,0),c(0,0,20))x,y>delta_t=0.1;n=201>x=matrix(0,nrow=2,ncol=n);x[,1]=c(33,0)>y=matrix(0,nrow=2,ncol=n);y[,1]=c(0,0)>##畫出相應(yīng)的曲線>for(jin1:(n-1)){+ d=sqrt((x[2,j]-x[1,j])^2+(y[2,j]-y[1,j])^2)+ x[1,j+1]=x[1,j]+2*delta_t*(x[2,j]-x[1,j])/d+ y[1,j+1]=y[1,j]+2*delta_t*(y[2,j]-y[1,j])/d+ x[2,j+1]=0;y[2,j+1]=delta_t*j+}>lines(x[1,],y[1,])>----宋停云與您分享--------宋停云與您分享----0 B251015A0 O0 5 10 15 20 25 3010.2：追逐曲線從圖中可以看出，在兔子到達(dá)洞口前能被狼抓住一個服務(wù)員的售貨亭，顧客到達(dá)時間服從均值20秒，標(biāo)準(zhǔn)差105由題意可得，顧客到達(dá)時間間隔為rnorm(1,20,10)，服務(wù)時間間隔rnorm(1,28.5,5.45)>queue1<-function(lambda,mu,T){----宋停云與您分享--------宋停云與您分享----+ k=0;wt=0;wn=0;ws=0;+ tp=0;nA=0;n=0;t=0+ r=rnorm(1,20,10);tA=r;tD=Inf++ repeat{+ k=k+1;wt[k]=t;wn[k]=n+ if(tA<T){+ ws[k]=min(tA,tD)-t+ if(tA<tD){+ t=tA;n=n+1;nA=nA+1+ r=rnorm(1,20,10);tA<-t+r+ if(n==1){+ r=rnorm(1,28.5,5.45);tD<-t+r+ }+ }else{+ t=tD;n=n-1+ if(n==0){+ tD=Inf+ }else{+ r=rnorm(1,28.5,5.45);tD=t+r+ }+ }+ }else{----宋停云與您分享--------宋停云與您分享----+ ws[k]=if(tD==Inf)0elsetD-t+ if(n>0){+ t=tD;n=n-1+ if(n>0){+ r=rnorm(1,28.5,5.45);tD=t+r+ }+ }else+ tp=1+ }+ if(tp==1)break+ }+ data.frame(Ls=sum(ws*wn)/t,Ws=sum(ws*wn)/nA,+ Pwait=sum(ws[wn>=1])/t)+}>queue1(lambda=1/20,mu=1/28.5,T=12*3600)Ls Ws Pwait1331.97279475.720.999687331.97279475.72售貨亭繁忙的概率為0.99968710----宋停云與您分享--------宋停云與您分享---->queue4<-function(lambda,mu1,mu2,mu3,mu4,T,S=1,K){+ if(K<S)K<-S+ k<-0;wt<-0;wn<-0;ws<-0+ tp<-0;nA<-0;t<-0+ r<-runif(1);tA<--1/lambda*log(r)+ tD<-rep(Inf,S);SS<-rep(0,S+1)++ repeat{+ t1<-if(SS[1]==0)Infelsemin(tD)+ i1<-if(SS[1]==0)1elsewhich.min(tD)+ k<-k+1;wt[k]<-t;wn[k]<-SS[1]+ if(tA<T){+ ws[k]<-min(tA,t1+1)-t+ if(tA<t1+1){+ t<-tA;nA<-nA+1+ r<-runif(1);tA<-t-1/lambda*log(r)+ n<-SS[1];SS[1]<-n+1+ for(iin1:S){+ if(SS[1+i]==0){+ SS[1+i]<-1+ r<-runif(1)----宋停云與您分享--------宋停云與您分享----+ if(i==1)+ {d=runif(1)+ if(0<=d||d<=0.6)+ {tD[i]<-t-1/mu1*log(r)}+ else tD[i]<-t-1/mu2*log(r)}+ if(i==2)+ {e=runif(1)+ if(0<=e||e<=0.6)+ {tD[i]<-t-1/mu3*log(r)}+ else tD[i]<-t-1/mu4*log(r)}++ break+ }+ }+ if(SS[1]==K){+ t1<-min(tD)+ while(tA<t1+1){+ r<-runif(1);tA<-tA-1/lambda*log(r)+ }+ }+ }else{+ t<-t1+1;n<-SS[1];SS[1]<-n-1+ if(n==1){----宋停云與您分享--------宋停云與您分享----+ SS[2:(S+1)]<-0;tD[1:S]<-Inf+ }elseif(n<=S){+ SS[1+i1]<-0;tD[i1]<-Inf+ }else{+ r<-runif(1)+ if(i==1)+ {d=runif(1)+ if(0<=d||d<=0.6)+ {tD[i]<-t-1/mu1*log(r)}+ else tD[i]<-t-1/mu2*log(r)}+ if(i==2)+ {e=runif(1)+ if(0<=e||e<=0.6)+ {tD[i]<-t-1/mu3*log(r)}+ else tD[i]<-t-1/mu4*log(r)}++ }+ }+ }else{+ ws[k]<-if(t1==Inf)0elset1+1-t+ n<-SS[1]+ if(n>0){+ t<-t1+1;SS[1]<-n-1;----宋停云與您分享--------宋停云與您分享----++if(n==1){SS[2:(S+1)]<-0;tD[1:S]<-Inf+++++}elseif(n<=S){SS[1+i1]<-0;tD[i1]<-Inf}else{r<-runif(1)if(i==1)+++++{d=runif(1)if(0<=d||d<=0.6){tD[i]<-t-1/mu1*log(r)}else tD[i]<-t-1/mu2*log(r)}if(i==2)++++++{e=runif(1)if(0<=e||e<=0.6){tD[i]<-t-1/mu3*log(r)}else tD[i]<-t-1/mu4*log(r)}}}else++}tp<-1+ if(tp==1)break+ }+ data.frame(Ls=sum(ws*wn)/t,Ws=sum(ws*wn)/nA,----宋停云與您分享--------宋停云與您分享----+ Plost=sum(ws[wn>=K])/t)+}>queue4(lambda=4,mu1=6,mu2=9,mu3=5,mu4=7.5,T=10*60,S=2,K=6)Ls Ws Plost13.3778331.0076350.1770578從結(jié)果上看一天的運營中，來電顧客的隊長為3.377833人，平均逗留時間為1.007635，理發(fā)店的損失率為0.1770578----宋停云與您分享--------宋停云與您分享----系統(tǒng)建模仿真技術(shù)的歷史現(xiàn)狀和發(fā)展趨勢分析工程133 胡浩 3130212026VVA控制器的仿真軟件，在研究控制策略、控制算法、控制系統(tǒng)的品質(zhì)方面提供了強大的支持。隨著執(zhí)行機構(gòu)技術(shù)的發(fā)展，機、電、液、熱、氣、磁等驅(qū)動技術(shù)的進步，以高可靠性、高精度、高反應(yīng)速度和AMESIM軟件正是能夠提供平臺級仿真技術(shù)的工具。從根據(jù)用戶需求，提供液壓、機械、氣動等設(shè)計分析到復(fù)雜系統(tǒng)的全系統(tǒng)分析，----宋停云與您分享--------宋停云與您分享----到引領(lǐng)協(xié)同仿真技術(shù)的發(fā)展方向，AMESIM的發(fā)展軌跡和方向代表了工程系統(tǒng)仿真技術(shù)的發(fā)展歷程和趨勢。一、系統(tǒng)仿真技術(shù)發(fā)展的現(xiàn)狀工程系統(tǒng)仿真作為虛擬設(shè)計技術(shù)的一部分，與控制仿真、視景仿1、控制器和被控對象的聯(lián)合仿真：MATLAB＋AMESIM，可以覆蓋整個自動控制系統(tǒng)的全部要求。2、被控對象的多學(xué)科、跨專業(yè)的聯(lián)合仿真：AMESIM＋機構(gòu)動力學(xué)＋CFD＋THERMAL＋電磁分析3、實時仿真技術(shù)實時仿真技術(shù)是由仿真軟件與仿真機等半實物仿真系統(tǒng)聯(lián)合實現(xiàn)的，通過物理系統(tǒng)的實時模型來測試成型或者硬件控制器。4、集成進設(shè)計平臺現(xiàn)代研發(fā)制造單位，尤其是設(shè)計研發(fā)和制造一體化的大型單位，PDM/PLMCAE工作的一部分，被要求嵌入流程，與上下游工具配合。----宋停云與您分享--------宋停云與您分享----5、超越仿真技術(shù)本身工程師不必是精通數(shù)值算法和仿真技術(shù)的專家，而只需要關(guān)注自己的專業(yè)對象，其他大量的模型建立、算法選擇和數(shù)據(jù)前后處理等工作都交給軟件自動完成。這一技術(shù)特點極大地提高了仿真的效率，降低了系統(tǒng)仿真技術(shù)的應(yīng)用門檻，避免了因為不了解算法造成的仿真失敗。6、構(gòu)建虛擬產(chǎn)品在通過建立虛擬產(chǎn)品進行開發(fā)和優(yōu)化過程中，關(guān)注以各種特征值為代表的系統(tǒng)性能，實現(xiàn)多方案的快速比較。二、系統(tǒng)仿真技術(shù)的發(fā)展趨勢2．1、屏棄單專業(yè)的仿真單一專業(yè)仿真將退出系統(tǒng)設(shè)計的領(lǐng)域，專注于單一專業(yè)技術(shù)的深2．2、跟隨計算技術(shù)的發(fā)展隨著計算技術(shù)在軟硬件方面的發(fā)展，大型工程軟件系統(tǒng)開始有減在高性能計算方面，將支持包括并行處理、網(wǎng)格計算技術(shù)和高速計算系統(tǒng)等技術(shù)。2．3、平臺化----宋停云與您分享--------宋停云與您分享----要求仿真工具能夠提供建模、運算、數(shù)據(jù)處理（包括二次開發(fā)后的集成和封裝）PDM/PLM2．4、整合和細(xì)分市場整合化：將出現(xiàn)主流的標(biāo)準(zhǔn)工具。其特征是功能涵蓋了現(xiàn)代工業(yè)專業(yè)化：隨著市場需求的細(xì)分，走專業(yè)化道路，將出現(xiàn)極專業(yè)的工具。這些工具將在某些具體的專業(yè)領(lǐng)域提供深入研究的特殊支持，如開發(fā)特殊的庫或模型，專注于具有鮮明行業(yè)特征的技術(shù)，滿足特殊的行業(yè)標(biāo)準(zhǔn)。將出現(xiàn)整合型工具和專業(yè)化工具互補的局面。、智能化將引進更加友好的操作界面，智能化的求解器及模2．6、豐富的二次開發(fā)選項提供源代碼級的二次開發(fā)支持，開放的架構(gòu)滿足不同用戶的專業(yè)開發(fā)要求。----宋停云與您分享--------宋停云與您分享----在強大的工具平臺上，根據(jù)自身的需要，進行二次開發(fā)。這已經(jīng)【】【1】梁炳成,王恒霖,鄭燕紅;用仿真技術(shù)的發(fā)展動向和展望 [J];系統(tǒng)仿真學(xué)報;2001年01期【2】楊明,張冰,馬萍,王子才;VV&A發(fā)展的五大關(guān)鍵問題[J];系統(tǒng)仿真學(xué)報;200311期【