Iは、確率密度関数{P（Z、PHI）を}プロット及び添付のEQに示すように、それを統合する必要がある必要がある確率密度関数

Question

との統合。＃1 enter image description here AFとVfは定数である 、 φは角度であり、 はz（距離、数値、小数でも可）です。 誰かが私をMATLABで案内できますか？これらの方程式を書くにはどうしたらいいですか？あなたの関数をintregateする

Answer 1

、あなたはM-ファイルまたは無名関数

F = @（Z、PHI）P（Z、PHI）* pを（Z、PHI）

作成する必要がありますどちらかあなたはPとPを同様に構築します。次に、fを2回積分するためにode45などの数値積分器の1つを使用する必要があります.fを1回、phiを1回積分します。

Answer 2

私が正しく理解すれば、-lf/2とlf/2の間の一様確率分布に、正弦波の第1四半期のような別の確率分布を乗算します。結果の確率分布を知りたい。

基本的にlf/2> pi/2の場合、同じ分布になります。正弦分布は完全に一様分布の内部にある。 （lf/2）<（pi/2）の場合、サイン分布の一部分の一様分布が切り詰められます。確率分散を分割して、積分が1になるように分割したいとします。それは確率分布のままでなければならない。

sin（x）の積分はcos（x）です。だから、あなたは（1-COS（LF/2））によりdevideその場合

以下

にそれがより見えるようにするスクリプトです：

lf=2; 
xx = linspace(-lf,lf,1E4); 
p1 = (xx>-lf/2&xx<lf/2)*(1/lf); 
p2 = zeros(size(xx)); 
p2(xx>0&xx<pi/2) = sin(xx(xx>0&xx<pi/2)); 

p3 = p2.*p1.*lf; 
if lf<pi 
    p3 = p3./(1-cos(lf/2)); 
end 
plot(xx,p1,xx,p2,xx,p3) 
legend({'uniform distribution','sine','result'}) 
%integrals (actually Riemann sums): 
sum(p1.*(xx(2)-xx(1))) 
sum(p2.*(xx(2)-xx(1))) 
sum(p3.*(xx(2)-xx(1)))