用matlab求高斯函式積分,用matlab求高斯函式積分

2022-12-07 08:06:01 字數 3540 閱讀 4030

1樓:匿名使用者

積分割槽間呢?

%%%**********

a1=0.3826;

b1=452.1;

c1=9.185;

a2=0.5569;

b2=455;

c2=18.23;

a3=-0.03431;

b3=497.9;

c3=11.48;

a4=0.2741;

b4=554.2;

c4=85.93;

%積分上下限

x0=1;

x1=3;

%高斯積分點以及權係數

gx=[-0.9061799,-0.5384693,0,0.5384693,0.9061799];

gweight=[0.2369269,0.4786287,0.5688889,0.4786287,0.2369269];

intf=0;

for i=1:5

x=(x0+x1)/2+(x1-x0)/2*gx(i);

intf=intf+gweight(i)*(a1*exp(-((x-b1)/c1).^2)+a2*exp(-((x-b2)/c2).^2)+a3*exp(-((x-b3)/c3).

^2)+a4*exp(-((x-b4)/c4).^2));

endintf=intf*(x1-x0)/2;

intf

2樓:年紅旭萬愷

就是要多峰值的高斯函式吧?下面的這個行不行?

ezplot('0.5*exp(-(x+4)^2)+exp(-(x/2)^2)+0.7*exp(-(x-5)^2)')

3樓:高壓學徒

x=1:del:n,做成離散函式。積分=int(x);del越小越精確,n在del確定的基礎上覆蓋f的範圍就行,做個plot(f)看看確保積分的範圍足夠。

如何用matlab計算正態分佈概率分佈函式的積分? 5

4樓:小小小小顰

具體抄操作步驟如下:

1、首先,襲提出問題,然後使用baimatlab計算下圖du中的積分問

zhi題,如下圖所示,dao然後進入下一步。

2、其次,完成上述步驟後,開啟matlab軟體,並按照以下**清除工作區,如下圖所示,然後進入下一步。

3、接著,完成上述步驟後,在matlab中定義符號變數並定義函式,**如下圖所示,然後進入下一步。

4、然後,完成上述步驟後,可以獲得新定義的函式,**如下圖紅框標註所示,然後進入下一步。

5、最後,完成上述步驟後,要計算積分的值,請使用以下**:fz=z;int1=int(fz*c,fa,0,pi)int2=int(int1,theta,0,pi*2),如下圖所示。這樣,問題就解決了。

5樓:匿名使用者

y = cdf('norm' ,x,a,b);

'norm' (normal distribution)%正態分佈復x就是你要求的從制負無窮到x的積分

a 為平均值

b 為標準差

例如,計算均值為0 標準差為1 的分佈,從負無窮到 1 的積分n=cdf('normal',1,0,1)n =

0.84134

matlab對高斯函式求積分,積分下限是負無窮,上限是x,x取值是從-255:255的陣列,要怎麼弄?

6樓:匿名使用者

應該用誤差函式erf來求。

1、首先,積分上下限:

∫(-∞,x)應分成∫(-∞,0)+∫(0,x)=-∫(0,-∞)+∫(0,x)

2、被積變數t應作變換:

t1=t/σ → t=σ*t1

相應的積分限x變為x/σ

3、係數:

dt=σ*dt1,σ和原係數分母中的σ約分,餘下1/√(2π),與erf函式的係數對照,應該乘以1/(2√2)

綜上,原表示式的計算如下(σ的取值自定):

x=-255:255;

sigma=100;

f=1/(2*(2)^0.5) * ( erf(x/sigma)-erf(-inf) )

matlab 兩個高斯函式卷積後求積分

高斯函式積分 30

7樓:

分享一種解法,利用高斯分佈/正態分佈密度函式的性質和伽瑪函式【γ(α)】求解。設a=[1/(δ√π)]^(1/2)、積分(1)、(2)、(3)、(4)式分別用i1、i2、i3、i4表示。

∵x~n(μ,δ²),其密度函式f(x)=(1/√2)a²e^[-(x-μ)²/(2δ²)],∴e(x)=∫(-∞,∞)xf(x)=μ,d(x)=e(x²)-[e(x)]²=δ²,∴e(x²)=∫(-∞,∞)x²f(x)=μ²+δ²。

∴對i1,易得i1=a(√2)/a²e(x)=0; 對i2,易得i2=a(√2)/a²e(x²)=δ²(√2)/a; 對i3,易得i3=a(√2)/a²e(xⁿ)。利用被積函式「xⁿf(x)」 的奇偶性質,n為奇數時,i3=0、n為偶數時,i3=[2a(δ√2)^(n+1)]γ((n+1)/2)。

對i4,∵x²/(2δ²)+ikx=((x+iδ²k)²/(2δ²)+δ²k²/2,∴i4=[(√2)/a]e^(-δ²k²/2)。

供參考。

8樓:為了生活奔波

一類具有最高的代數精度的內插型求積公式(表2)。求積公式(2)含有2(m+1)個自由引數(xj和aj),恰當選擇這些引數,能使公式(2)的代數精度達到2m+1。高斯求積理論中的一個基本定理斷言:

只要把結點x0,x1,…,xm取為區間[α,b]上關於權函式ω(x)的m+1次正交多項式的零點,內插型求積公式(2)即達到最高代數精度2m+1。這裡[α,b] 可以是有限或無限區間,ω(x)為取正值的權函式。許多有關數值積分的論著都列舉出各種高斯型公式的結點和係數的數值。

可以證明:對每個連續函式,當結點個數趨於無窮時,高斯型公式所給出的近似值序列收斂到相應積分的精確值,而牛頓-科茨公式則不具有這種性質。高維數值積分的主要方法有蒙特卡羅法、代數方法和數論方法。

高斯函式與三角函式乘積的積分如何用matlab求? 10

9樓:匿名使用者

x=1:0.1:10;

y=exp(-x.^2).*cos(x) ;

trapz(x,y)

執行結果

ans =     0.034609

matlab高斯積分題,求大神幫助 20

10樓:匿名使用者

你好!function y=fun(x)

y=2*x;

endclc

clear all

a=0;

b=5;

%兩點高斯積分即可滿足精度要求

s=(fun((b+a)/2-(b-a)/2*0.5773502692)*1...

+fun((b+a)/2+(b-a)/2*0.5773502692)*1)*(b-a)/2;

fprintf('積分值是:%f',s)

matlab求定積分解析解,matlab求一個定積分解析解

這個求不出來的。warning explicit integral could not be found.e的表示式的積分找不到。只能通過數值計算的方法求出近似值,這個近似值也保持了足夠的精度 dx 0.01 x 0 dx 2 y exp exp x x s dx cumtrapz y 計算區間內曲...

matlab求積分的數值解,MATLAB求積分的數值解

可以積分出來的,只是好像你的r和 怎麼是反的 clear all clc syms x r 0.3136 beta 0.579 y r 2 pi beta 0.5 x 1.5 exp x r 2 2 beta x s int x exp 3 x y,x,0,inf s double s s 0.08...

matlab畫函式影象問題,用matlab怎麼畫基本的函式影象

x 10 0.001 10 plot x,sin x x 10 0.001 10 plot x,exp x 隱函式有點麻煩,你需要把具體的函式給我,我幫你畫 y 163380 1.60000 t.2 266667 t.2 t.2 t.2 1 1 t 這一句改成這樣就成了。如果看不清楚,可以把所有都改...