MATLAB的答案

不能使用Matlab集成功能

1视图(30天)
大卫餐桌
大卫餐桌 2023年2月9日
评论道: 大卫餐桌2023年2月9日
我想做一个基本的集成的公式4不同的值。做一些阅读后,显然我需要创建一个函数处理,但我没有得到任何的成功。我已经把我的下面。我可能会错过更多。
L = 0:100;
L_bar = 25;
sigma_d = (0.25, 0.5, 0.75, 1.0);
日志(L_bar) -μ= 0.5 * (sigma_d)。^ 2;
L_tilda = exp(μ);
% eta_r = p * exp (r / L);
%我想整合埃塔从零到无穷大,P变化
%变量给4种不同的情节。
r = 0:100;
%我只是测试下,看看一个循环运行之前
%绘制任何
P = 0(1、长度(L));
有趣= @ (r) p * exp (r / L);
为2 = 1:元素个数(sigma_d)
P = 1. / (L *√(2 *π)* sigma_d (ii)。* exp(日志(L / L_tilda (ii)。^ 2 / (2 * sigma_d (2) ^ 2));
eta_r =积分(有趣,0,正);
结束

在回答这个问题。

接受的答案

Torsten
Torsten 2023年2月9日
跑:
L_bar = 25;
sigma_d = (0.25, 0.5, 0.75, 1.0);
日志(L_bar) -μ= 0.5 * (sigma_d)。^ 2;
L_tilda = exp(μ);
r = 0:100;
为2 = 1:元素个数(sigma_d)
P = @(左)/ 1. (L * sqrt(2 *π)* sigma_d (ii)。* exp(日志(L / L_tilda (ii)。^ 2 / (2 * sigma_d (2) ^ 2));
有趣= @(左,右)P (L)。* exp (r / L);
eta_r(二世:)=积分(@ (L)乐趣(左,右),0,正无穷,“ArrayValued”,真正的);
结束
情节(r, eta_r)
网格在
4评论
显示 3年长的评论
大卫餐桌
大卫餐桌 2023年2月9日
完美的。非常感谢你

登录置评。

答案(1)

Dyuman Joshi
Dyuman Joshi 2023年2月9日
跑:
您必须显式地改变函数句柄改变它的定义
P =π;
有趣= @ (x) P * x;
P = 3;
%的乐趣(3)不等于3 * 3 = 9
有趣的(3)
ans = 9.4248
L = 0:100;
L_bar = 25;
sigma_d = (0.25, 0.5, 0.75, 1.0);
日志(L_bar) -μ= 0.5 * (sigma_d)。^ 2;
L_tilda = exp(μ);
% eta_r = p * exp (r / L);
%我想整合埃塔从零到无穷大,P变化
%变量给4种不同的情节。
r = 0:100;
%我只是测试下,看看一个循环运行之前
%绘制任何
为2 = 1:元素个数(sigma_d)
P = 1. / (L *√(2 *π)* sigma_d (ii)。* exp(日志(L / L_tilda (ii)。^ 2 / (2 * sigma_d (2) ^ 2));
有趣= @ (x) p * exp (x / L);
eta_r =积分(有趣,0,正无穷,“ArrayValued”,真正的)
结束
警告:正或NaN值。
eta_r = 1×101
南0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0006 0.0030 0.0109 0.0307 0.0718 0.1437 0.2535 0.4022 0.5842 0.7870 0.9943 1.1886 1.3546 1.4810 1.5614 1.5946 1.5834 1.5336 1.4531 1.3500 1.2326
警告:正或NaN值。
eta_r = 1×101
南0.0000 0.0000 0.0003 0.0023 0.0097 0.0269 0.0572 0.1019 0.1598 0.2281 0.3028 0.3800 0.4560 0.5276 0.5924 0.6491 0.6965 0.7345 0.7630 0.7827 0.7940 0.7979 0.7951 0.7867 0.7733 0.7560 0.7354 0.7122 0.6871
警告:正或NaN值。
eta_r = 1×101
南0.0002 0.0060 0.0263 0.0626 0.1109 0.1656 0.2219 0.2764 0.3268 0.3717 0.4106 0.4433 0.4701 0.4914 0.5076 0.5192 0.5268 0.5309 0.5319 0.5303 0.5265 0.5209 0.5137 0.5053 0.4958 0.4855 0.4746 0.4632 0.4514
警告:正或NaN值。
eta_r = 1×101
南0.0099 0.0513 0.1074 0.1642 0.2156 0.2596 0.2959 0.3252 0.3482 0.3658 0.3789 0.3882 0.3942 0.3977 0.3989 0.3984 0.3963 0.3931 0.3889 0.3839 0.3783 0.3722 0.3658 0.3590 0.3521 0.3450 0.3378 0.3305 0.3233
因为你是除以0 (L)的第一个元素,你会得到一个NaN,因此积分的警告()解决。
3评论
显示特殊照顾年长的评论
大卫餐桌
大卫餐桌 2023年2月9日
Aplogies
所以我的P的定义
P = 1. / (L *√(2 *π)* sigma_d)。* exp(日志(L / L_tilda)。^ 2 / (2 * sigma_d ^ 2))
所以L无所不能在P和L定义为一个长度。我想我也许可以这样一个固定的常数L = 100,保持r = 0:100
我要计算积分
我改变L_tilda和σ4 * 4的值不同的情节。
希望这更清楚:)

登录置评。

类别

找到更多的在编程在帮助中心和文件交换

标签

s manbetx 845


释放

R2022b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的