Matlab求解積分上限是自變量的計(jì)算方法

數(shù)學(xué)中經(jīng)常會(huì)遇到積分上限是自變量的積分，例如計(jì)算大氣透過(guò)率的積分：
$$\Gamma (z)=exp(?{\int }_0^z\mathrm{A}(t)dt).$$
但Matlab的數(shù)值積分中沒(méi)有直接計(jì)算此類(lèi)積分的函數(shù)以調(diào)用，因此，可采用兩步實(shí)現(xiàn)：
1、定義積分核函數(shù)：

funA = @(t) A(t);

2、定義積分函數(shù)：

funGamma = @(z) integral(funA, 0, z);

之后可用for循環(huán)對(duì)每個(gè)自變量求出積分，并作圖。

h = (0:0.1:30); %海拔 Gamma = zeros(1,length(h)); %預(yù)分配內(nèi)存 for i = 1:1:length(h)Gamma(i) = exp(-funGamma(h(i))); end figure(1) plot(Gamma,h);

此為一個(gè)簡(jiǎn)單易于實(shí)現(xiàn)的方法，共大家參考。

總結(jié)

以上是生活随笔為你收集整理的2020-09-27的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。