montecarlo模擬誤差分析課程設(shè)計(參考版)
2025-01-16 16:36本頁面
【正文】 22。適合解決概率密度函數(shù)非對稱的問題(2) 最短包含區(qū)間與GUM法同時計算不確定度的時候，哪一個略顯保守。程序運(yùn)行結(jié)果：Figure 1Figure 2Figure 3Figure 4計算結(jié)果：delta_n_u2 =，s = 實驗需要討論的問題：(1) 最短包含區(qū)間MCM法與普通MCM法有何區(qū)別。delta_cancha=deltadelta_mean。grid on。r*39。 for r=1:M/2 if area2(r+q)area2(r)=area2(jj+q)area2(jj)。for j=1:M1 area2(j+1)=area2(j)+s2(j)。s2(iiii)=d_delta*abs((H(iiii)))。 %置信概率q=P*M。 delta_n_u1=(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+1),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))))。%%計算直方圖的面積%%S=sum(d_delta*abs(H))。grid on。b:39。 %%畫包絡(luò)線%%figure(4)HH=envelope(x_,H)。grid on。figure(3)。figure(2)axis([1,1,0,N])hist(Y2,M)。 %%delta_n為誤差分布序列%% 畫直方圖figure(1)axis([1,1,0,N])hist(Y1,M)。d_delta=(max(delta)min(delta))/(M1)。delta2=u2*Y2。Y2=sqrt(2*log(X1)).*cos(2*pi*X2)。X2=rand(1,N)。,31)。%%rand(39。u1=。x=0:1:M。 %%設(shè)定參數(shù)值%%%%隨機(jī)信號點(diǎn)數(shù)N，均值1，標(biāo)準(zhǔn)差%%N=10^6。clc。(1). 先確定q(P=,M=1000000,q=PM=10000)(2). 重新計算包絡(luò)線下的面積(不是對稱的時候)(3). 根據(jù)算法：，計算r*(4). r*對應(yīng)的區(qū)間長度極為最短包含區(qū)間實驗示例程序：%%x為均勻分布，正態(tài)分布，反正弦分布%%y=sin(x)為何種分布tic。 基于最短包含區(qū)間的MCM法如果PDF不對稱，可采用最短包含區(qū)間。程序運(yùn)行結(jié)果：Figure 1Figure 2實驗需要討論的問題：如何根據(jù)三維網(wǎng)格曲線和三維曲線獲得最佳的N，M組合。s=sqrt((sum(delta_cancha.^2))./(floor(N)1))。 %%理論計算標(biāo)準(zhǔn)不確定度%%delta_mean=mean(delta)。grid on。ro39。if (area(iii)*S)=0。end %% %的直方圖面積for iii=1:floor(M./2)。 area(1)=s_(1)。s_2(i)=d_delta.*abs(floor(H(floor(M./2i+1))))。%% 計算直方圖的面積%%%%s_1表示正向直方圖的每一個單元的面積%%s_2表示反向直方圖的每一個單元的面積%%s_表示正反向兩兩對稱每一對單元的面積%%area表示以中心為對稱軸的累加面積i=1:1:floor(M./2)。% hold on。)。% plot(delta_n,HH,39。grid on。% %%合成誤差幅度直方圖%%% figure(5)% axis([1,1,0,N]) H=hist(delta,floor(M))。% %%正態(tài)分布誤差2幅度直方圖%%% figure(4)% axis([1,1,0,N])% hist(delta2,M)。% % hold on% %%正態(tài)分布誤差1幅度直方圖%%% figure(3)% ax
高考資料相關(guān)推薦
鄂ICP備17016276號-1