montecarlo模擬誤差分析課程設(shè)計(jì)-資料下載頁(yè)
2025-06-04 15:04本頁(yè)面
【正文】 r*對(duì)應(yīng)的區(qū)間長(zhǎng)度極為最短包含區(qū)間 實(shí)驗(yàn) 示例程序： %%x為均勻分布，正態(tài)分布，反正弦分布 %%y=sin(x)為何種分布 tic。 clear。clc。close all。 %%設(shè)定參數(shù)值 %% %%隨機(jī)信號(hào)點(diǎn)數(shù) N，均值 1，標(biāo)準(zhǔn)差 %% N=10^6。 M=N/10。 x=0:1:M。 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為誤差分布序列 %% 畫(huà)直方圖 19 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)。 %%畫(huà)包絡(luò)線 %% figure(4) HH=envelope(x_,H)。 plot(delta_n,HH,39。b:39。)。grid on。 hold on。 %%計(jì)算直方圖的面積 %% S=sum(d_delta*abs(H))。 %% 計(jì)算直方圖的面積 %% %%s_1表示正向直方圖的每一個(gè)單元的面 積 %%s_2表示反向直方圖的每一個(gè)單元的面積 %%s_表示正反向兩兩對(duì)稱每一對(duì)單元的面積 %%area表示以中心為對(duì)稱軸的累加面積 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/21 area(ii+1)=area(ii)+s_(ii)。 end %% 計(jì)算 %的直方圖面積 for iii=1:M/2。 area(iii)。 if (area(iii)*S)=0。 break end end plot([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。 20 %% 計(jì)算最短包含 區(qū)間面積 P=。 %置信概率 q=P*M。 % 參數(shù) q iiii=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 end end plot([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 %%理論計(jì)算標(biāo)準(zhǔn)不確定度 %% delta_mean=mean(delta)。 delta_cancha=deltadelta_mean。 s=sqrt((sum(delta_cancha.^2))/(N1)) %%%%%%%%%%%%%%% toc。 程序運(yùn)行結(jié)果 ： Figure 1 21 Figure 2 Figure 3 Figure 4 22 計(jì)算結(jié)果 ： delta_n_u2 =， s = 實(shí)驗(yàn)需要討論的問(wèn)題： (1) 最短包含區(qū)間 MCM 法與普通 MCM 法有何區(qū)別。最短包含區(qū)間 MCM 適合解決怎樣的問(wèn)題。 適合解決概率密度函數(shù)非對(duì)稱的問(wèn)題 (2) 最短包含區(qū)間與 GUM法同時(shí)計(jì)算不確定度的時(shí)候，哪一個(gè)略顯保守。 由計(jì)算得 delta_n_u2 =， s = ，顯然，在評(píng)定不確定度的 時(shí)候，最短包含區(qū)間略顯保守。
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1