montecarlo模擬誤差分析課程設(shè)計(jì)-文庫(kù)吧資料
2025-01-19 16:36本頁(yè)面
【正文】 is([1,1,0,N])% hist(delta1,M)。)。% % plot(x,MuSigma,39。)。% % plot(x,Mu+Sigma,39。k39。grid on。% % hold on% % %%變換為任意均值和方差的正態(tài)分布%%% %Z1=Sigma*Y1+Mu。)。% % plot(x,1,39。)。% % plot(x,1,39。)。% % hold on% % plot(x,0,39。% % figure(2),% axis([0,N,max(5*Y2),max(5*Y2)])% plot(Y2)。 %%delta_n為誤差分布序列%%作圖%%%%高斯隨機(jī)信號(hào)%%% figure(1),% axis([0,N,max(5*Y1),max(5*Y1)])% plot(Y1)。d_delta=(max(delta)min(delta))./(floor(M)1)。delta2=u2*Y2。Y2=sqrt(2*log(X1)).*sin(2*pi*X2)。X2=rand(1,floor(N))。,39)。 %%產(chǎn)生兩個(gè)在(0,1)上服從均勻分布的,種子為39,每一次都相同的隨機(jī)數(shù)X1和X2%%rand(39。u1=。x=0:1:floor(M)。Mu=1。%%bdclose all。clc。toc。endfigure(2),plot3(c,d,Result2)。for jjj=1:max(size(c))。c=logspace(1,6)。 Result1(j,jj)=shiyan(a(j),b(jj))。for j=1:max(size(a))。warning off。(3). 分別進(jìn)行三維網(wǎng)格作圖和三維曲線作圖，通過觀察曲線獲得最佳的N，M組合。(1). 基于前一個(gè)實(shí)驗(yàn)，構(gòu)建衡量Monte Carlo法和GUM法計(jì)算得到標(biāo)準(zhǔn)不確定度差值的函數(shù)。 自適應(yīng)MCM法在執(zhí)行自適應(yīng)蒙特卡洛方法的基本過程中，蒙特卡洛試驗(yàn)次數(shù)不斷增加，直至所需要的各種結(jié)果達(dá)到統(tǒng)計(jì)意義上的穩(wěn)定。當(dāng)N=M時(shí)，由于對(duì)離散的誤差值統(tǒng)計(jì)運(yùn)算存在舍入誤差導(dǎo)致誤差，此誤差隨著M的增大可消除此項(xiàng)舍入誤差。當(dāng)N改變時(shí)，即當(dāng)N增大時(shí)可使直方圖更為精細(xì)，且此時(shí)不改變直方圖包絡(luò)的基本形狀。根據(jù)以上各圖分析知：當(dāng)N固定的情況下，隨著M值得增大直方圖的平滑性變差，Y軸高度下降。 s=,delta_n_u= N=500000,M=5N。 s=,delta_n_u= N=500000,M=N。 delta_n_u= N=500000,M=N/5。 s=,delta_n_u= N=100000,M=5N。 s=,delta_n_u= N=100000, M=N。（2）當(dāng)更改N與M不同的倍數(shù)關(guān)系時(shí)，可得到不同的計(jì)算結(jié)果，如下各圖所示： N=100000,M=N/5。%%%%%%%%%%%%%%%toc。delta_cancha=deltadelta_mean。 delta_n_u=(delta_n(floor(M/2+iii))delta_n(floor(M/2iii+1)))/6。)。 breakendendplot([delta_n(M/2iii+1),delta_n(M/2+iii)],[H(M/2iii),H(M/2+iii)],39。 area(iii)。for ii=1:M/21area(ii+1)=area(ii)+s_(ii)。s_(i)=s_1(i)+s_2(i)。s_1(i)=d_delta*abs(floor(H(floor(M/2+i))))。%%計(jì)算直方圖的面積%%S=sum(d_delta*abs(H))。grid on。b:39。%%畫包絡(luò)線%%figure(6)HH=envelope(x_,H)。hist(delta,M)。grid on。grid on。grid on。r39。grid on。r39。)。% hold on% plot(x,Mu,39。%%作圖%%%%高斯隨機(jī)信號(hào)%%% subplot(2,2,2)% axis([0,N,6,6])% plot(Z1)。grid on。r39。grid on。r39。grid on。k39。grid on。grid on。
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1