フィルターのクリア

How to obtain numerical value during integration

1 回表示 (過去 30 日間)
okoth ochola
okoth ochola 2023 年 2 月 28 日
コメント済み: Dyuman Joshi 2023 年 2 月 28 日
Hi, I have a slight problem here, I am doing integration based on series of input picked from an array. Have used the following code and have got an output. However, the output am getting is not in numerical form, what can i do to get numerical values? Have pasted part of my code ad the results. I just want to know what I can do to get numerical resulsts as 10, 12 and so forth (numbers picked randomly just for explanation). Please find the attachment below; I need the Power_density_monthly to be array of values, kindy help.
for n=1:c
syms x
expr=(Monthly_parameters(1,n)/Monthly_parameters(2,n))*((x/Monthly_parameters(2,n))^(Monthly_parameters(1,n)-1))*exp(-(x/Monthly_parameters(2,n))^Monthly_parameters(1,n))*density*x^3
Power_density_monthly(n,:) = int(expr,[0 10]);
if n==c
display(Power_density_monthly)
Mp_power_monthly=0.5*density.*(Most_probable_wind_speed.^3)
syms x
expr=(Overal_parameters(1)/Overal_parameters(2))*((x/Overal_parameters(2))^(Overal_parameters(1)-1))*exp(-(x/Overal_parameters(2))^Overal_parameters(1))*density*x^3
Mp_power_Overall=int(expr,[0 10])
end
end
%the reuslts i got aare as follows
Power_density_monthly =
int((588683273498799241*x^3*exp(-((1125899906842624*x)/1214356027258251)^(1054706091853107/562949953421312))*((1125899906842624*x)/1214356027258251)^(491756138431795/562949953421312))/281474976710656000, x, 0, 10)
int((1136847329389978479*x^3*exp(-((562949953421312*x)/1059881698388299)^(888859807670943/281474976710656))*((562949953421312*x)/1059881698388299)^(607384830960287/281474976710656))/562949953421312000, x, 0, 10)
int((2324591596434618903*x^3*exp(-((4503599627370496*x)/5108924510777933)^(8760923315398151/4503599627370496))*((4503599627370496*x)/5108924510777933)^(4257323688027655/4503599627370496))/1125899906842624000, x, 0, 10)
int((374417108231022927*x^3*exp(-((4503599627370496*x)/5360205970664107)^(231329446538139/140737488355328))*((4503599627370496*x)/5360205970664107)^(90591958182811/140737488355328))/225179981368524800, x, 0, 10)
int((972383524037662051*x^3*exp(-((4503599627370496*x)/6882223601054367)^(154273069918933/70368744177664))*((4503599627370496*x)/6882223601054367)^(83904325741269/70368744177664))/562949953421312000, x, 0, 10)
int((304346401206112139*x^3*exp(-((140737488355328*x)/286487278512009)^(643201787191211/281474976710656))*((140737488355328*x)/286487278512009)^(361726810480555/281474976710656))/225179981368524800, x, 0, 10)
int((287711040109225411*x^3*exp(-((140737488355328*x)/420488410245951)^(7139606461783553/2251799813685248))*((140737488355328*x)/420488410245951)^(4887806648098305/2251799813685248))/225179981368524800, x, 0, 10)
int((2140903232335760319*x^3*exp(-((2251799813685248*x)/4266926036145285)^(6738851442847649/2251799813685248))*((2251799813685248*x)/4266926036145285)^(4487051629162401/2251799813685248))/1125899906842624000, x, 0, 10)
int((436526358616317991*x^3*exp(-((1125899906842624*x)/2709338596340557)^(3489855762716625/1125899906842624))*((1125899906842624*x)/2709338596340557)^(2363955855874001/1125899906842624))/281474976710656000, x, 0, 10)
int((2236234975767258893*x^3*exp(-((2251799813685248*x)/3933457441055675)^(6488818419739043/2251799813685248))*((2251799813685248*x)/3933457441055675)^(4237018606053795/2251799813685248))/1125899906842624000, x, 0, 10)
int((159206995598753823*x^3*exp(-((1125899906842624*x)/1816334693374449)^(6826242830561875/2251799813685248))*((1125899906842624*x)/1816334693374449)^(4574443016876627/2251799813685248))/70368744177664000, x, 0, 10)
int((424189302194492611*x^3*exp(-((4503599627370496*x)/6586604795853987)^(1288175588800873/562949953421312))*((4503599627370496*x)/6586604795853987)^(725225635379561/562949953421312))/225179981368524800, x, 0, 10)
  1 件のコメント
Dyuman Joshi
Dyuman Joshi 2023 年 2 月 28 日
If the input data is numerical, use integral (use element-wise operations whereever required)

サインインしてコメントする。

回答 (0 件)

カテゴリ

Help Center および File ExchangeSymbolic Math Toolbox についてさらに検索

タグ

製品


リリース

R2018a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by