montecarlo模擬誤差分析課程設計-資料下載頁
2025-01-13 16:36本頁面
【正文】 x_=[1:M]。u1=。u2=。%%rand(39。state39。,31)。X1=rand(1,N)。X2=rand(1,N)。Y1=sin(X1)。Y2=sqrt(2*log(X1)).*cos(2*pi*X2)。delta1=u1*Y1。delta2=u2*Y2。delta=delta1+delta2。d_delta=(max(delta)min(delta))/(M1)。 %%d_delta為誤差分布的間距delta_n=[min(delta):d_delta:max(delta)]。 %%delta_n為誤差分布序列%% 畫直方圖figure(1)axis([1,1,0,N])hist(Y1,M)。grid on。figure(2)axis([1,1,0,N])hist(Y2,M)。grid on。figure(3)。axis([1,1,0,N])hist(delta,M)。grid on。H=hist(delta,M)。 %%畫包絡線%%figure(4)HH=envelope(x_,H)。plot(delta_n,HH,39。b:39。)。grid on。hold on。%%計算直方圖的面積%%S=sum(d_delta*abs(H))。 %% 計算直方圖的面積%%%%s_1表示正向直方圖的每一個單元的面積%%s_2表示反向直方圖的每一個單元的面積%%s_表示正反向兩兩對稱每一對單元的面積%%area表示以中心為對稱軸的累加面積i=1:1:M/2。s_1(i)=d_delta*abs(floor(H(floor(M/2+i))))。s_2(i)=d_delta*abs(floor(H(floor(M/2i+1))))。s_(i)=s_1(i)+s_2(i)。 area(1)=s_(1)。for ii=1:M/21area(ii+1)=area(ii)+s_(ii)。end %% %的直方圖面積for iii=1:M/2。 area(iii)。if (area(iii)*S)=0。 breakendendplot([delta_n(M/2iii+1),delta_n(M/2+iii)],[H(M/2iii+1),H(M/2+iii)],39。ro39。)。grid on。 delta_n_u1=(delta_n(floor(M/2+iii))delta_n(floor(M/2iii+1)))/6。 %% 計算最短包含區(qū)間面積P=。 %置信概率q=P*M。 % 參數(shù)qiiii=1:M。s2(iiii)=d_delta*abs((H(iiii)))。area2(1)=s2(1)。for j=1:M1 area2(j+1)=area2(j)+s2(j)。end%%%%%%for jj=1:Mq area2(jj)。 for r=1:M/2 if area2(r+q)area2(r)=area2(jj+q)area2(jj)。 break end endendplot([delta_n(r),delta_n(r+q)],[H(r),H(r+q)],39。r*39。)。grid on。 delta_n_u2=(delta_n(floor(r+q))delta_n(floor(r)))/6%%理論計算標準不確定度%%delta_mean=mean(delta)。delta_cancha=deltadelta_mean。s=sqrt((sum(delta_cancha.^2))/(N1))%%%%%%%%%%%%%%%toc。程序運行結(jié)果：Figure 1Figure 2Figure 3Figure 4計算結(jié)果：delta_n_u2 =，s = 實驗需要討論的問題：(1) 最短包含區(qū)間MCM法與普通MCM法有何區(qū)別。最短包含區(qū)間MCM適合解決怎樣的問題。適合解決概率密度函數(shù)非對稱的問題(2) 最短包含區(qū)間與GUM法同時計算不確定度的時候，哪一個略顯保守。由計算得delta_n_u2 =，s = ，顯然，在評定不確定度的時候，最短包含區(qū)間略顯保守。22
高考資料相關推薦
鄂ICP備17016276號-1