導航:首頁 > 源碼編譯 > 高斯積分演算法matlab

高斯積分演算法matlab

發布時間:2024-04-23 11:00:59

Ⅰ Gauss縐鍒唌atlab鐨勭▼搴

% 楂樻柉-鍕掕╁痙鏁板肩Н鍒

fun=@(x)sin(x).*x.^3; % fun錛氱Н鍒琛ㄨ揪寮
% a,b錛氱Н鍒嗕笂涓嬮檺
a = 0;
b = pi;
tol=1e-8; % tol錛氱Н鍒嗙簿搴︼紝榛樿1e-6

% 璁$畻奼傜Н鑺傜偣
syms x
p=sym2poly(diff((x^2-1)^(n+1),n+1))/(2^n*factorial(n));
tk=roots(p); % 奼傜Н鑺傜偣
% 璁$畻奼傜Н緋繪暟
Ak=zeros(n+1,1);
for i=1:n+1
xkt=tk;
xkt(i)=[];
pn=poly(xkt);
fp=@(x)polyval(pn,x)/polyval(pn,tk(i));
Ak(i)=quadl(fp,-1,1,tol); % 奼傜Н緋繪暟
end
% 縐鍒嗗彉閲忎唬鎹錛屽皢[a,b]鍙樻崲鍒癧-1,1]
xk=(b-a)/2*tk+(b+a)/2;
% 媯楠岀Н鍒嗗嚱鏁癴un鏈夋晥鎬
fun=fcnchk(fun,'vectorize');
% 璁$畻鍙橀噺浠f崲涔嬪悗縐鍒嗗嚱鏁扮殑鍊
fx=fun(xk)*(b-a)/2;

% 璁$畻縐鍒嗗
ql=sum(Ak.*fx)
quadl(fun,0,pi) % 璋冪敤MATLAB鍐呴儴縐鍒嗗嚱鏁版楠

Ⅱ 奼傚姪錛孧atlab楂樻柉鍕掕╁痙奼傚嶆潅縐鍒

楂樻柉涓鍕掕╁痙奼傜Н鍏寮忚緝澶嶆潅,涓嶉傚悎涓鑸璁$畻搴旂敤銆傛枃鐚甗3]鎻愪緵浜嗗彉姝ラ暱楂樻柉涓鍕掕╁痙縐鍒嗙畻娉,鍏跺熀鏈鎬濇兂鏄鍩轟簬澶嶅寲奼傜Н鐨勬濇兂,瀵逛簬杈冮暱鍖洪棿[a,b銆,閫氳繃涓嶆柇鍦頒簩鍒嗗尯闂,姣

Ⅲ matlab瀵歸珮鏂鍑芥暟奼傜Н鍒嗭紝縐鍒嗕笅闄愭槸璐熸棤絀鳳紝涓婇檺鏄痻錛寈鍙栧兼槸浠-255錛255鐨勬暟緇勶紝瑕佹庝箞寮勶紵

搴旇ョ敤璇宸鍑芥暟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相關的資料

熱點內容
求算術平方根演算法設計思想 瀏覽:363
核彈打擊加密文件 瀏覽:820
虛擬機安卓9怎麼安裝toos 瀏覽:578
安卓怎麼下載joox 瀏覽:912
pdf天空 瀏覽:763
php讀取二維數組 瀏覽:348
php編譯安裝參數 瀏覽:278
其實壓力沒那麼大程序員圖片 瀏覽:416
如何查看app內訪問的網頁地址 瀏覽:757
安卓手機信號旁邊的漢字怎麼設置 瀏覽:304
nrf2401單片機 瀏覽:713
清除電腦文件夾垃圾的方法 瀏覽:226
天河程序員 瀏覽:192
成都程序員公積金 瀏覽:768
程序員為什麼叫程序猿 瀏覽:484
加西貝拉壓縮機價格 瀏覽:788
海信聚好看如何用u盤安裝app 瀏覽:71
加密狗怎麼寫的 瀏覽:560
安卓手機如何能調最大聲音 瀏覽:668
編程開發工具大全 瀏覽:572