麻烦的“双重”数值数据类型

20视图(30天)

显示旧的评论

弗朗西斯科·法夫里约5小时前

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1976269-trouble-with-the-double-numeric-data-type

评论道: 弗朗西斯科·法夫里9分钟前

我有下面的代码,提取从一个“r = 1:长度(μ)”周期不正常工作,我们计算结果为r = 24:

清晰;

信谊x t

r = 24;

deltamu = 0.5;

σ= 5;

λ= 1;

Cstar = 7;

μ= (4:deltamu: 10);

P = 0(1、长度(μ));

TP = 0(1、长度(μ));

TN = 0(1、长度(μ));

s = 1;

sigmae = 5;

exp_t = exp (- t ^ 2);

达峰时间= +正;

tmin =((μ(r) +λ。*σ。^ 2 x)。/ (sqrt(2) *σ));

误差补函数=(2. /√(π)。* int (exp_t t tmin,达峰时间);

f_lambda = (lambda. / 2)。* exp ((lambda. / 2) *(2。*μ(r) +λ。*σ^ 2 - 2。* x))。*误差补函数;

P (s, r) =双(int (f_lambda, x, Cstar +正));

phi_e = 1. / (sigmae。*√(2。* pi))。* exp (- (1/2)。* ((x-Cstar)。/ sigmae)。^ 2);

φ= int (phi_e x,负无穷,x);

fxP = f_lambda。*φ;

TP (s, r) =双(int (fxP, x, Cstar +正));

fxN1 = f_lambda;

fxN2 = f_lambda。*φ;

IntN1 =双(int (fxN1 x,负无穷,Cstar));

IntN2 =双(int (fxN2 x,负无穷,Cstar));

IntN = IntN1-IntN2;

TN (s, r) =双(IntN)

它的工作原理,计算向量TN的24日组件

TN =

列1到20

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

列21到29

0.2809 0 0 0 0 0 0 0 0

然而,令人惊奇的是,如果我们设定而不是r = 25日,我们获得以下信息:

“错误使用symengine

不能包含剩余的符号函数调用表达式转换成数组的两倍。参数必须表达评估数量。信谊误差/双(第729行)

Xstr = mupadmex (symobj:双,年代。年代,0);

2的评论
显示1年长的评论藏1年长的评论

Dyuman Joshi 1分钟前

并不是所有的积分可以通过MATLAB象征性地评估(或其他符号动力学)。在这种情况下可以使用 vpaintegral 数值近似积分的值-

                                信谊x t
                               
                                r = 25;
                               
                                deltamu = 0.5;
                               
                                σ= 5;
                               
                                λ= 1;
                               
                                Cstar = 7;
                               
                                μ= (4:deltamu: 10);
                               
                                P = 0(1、长度(μ));
                               
                                TP = 0(1、长度(μ));
                               
                                TN = 0(1、长度(μ));
                               
                                s = 1;
                               
                                sigmae = 5;
                               
                                exp_t = exp (- t ^ 2);
                               
                                达峰时间= +正;
                               
                                tmin =((μ(r) +λ。*σ。^ 2 x)。/ (sqrt(2) *σ));
                               
                                误差补函数=(2. /√(π)。* int (exp_t t tmin,达峰时间);
                               
                                f_lambda = (lambda. / 2)。* exp ((lambda. / 2) *(2。*μ(r) +λ。*σ^ 2 - 2。* x))。*误差补函数;
                               
                                双(vpaintegral P (s, r) = (f_lambda, x, Cstar +正));
                               
                                phi_e = 1. / (sigmae。*√(2。* pi))。* exp (- (1/2)。* ((x-Cstar)。/ sigmae)。^ 2);
                               
                                φ= int (phi_e x,负无穷,x);
                               
                                fxP = f_lambda。*φ;
                               
                                双(vpaintegral TP (s, r) = (fxP, x, Cstar +正));
                               
                                fxN1 = f_lambda;
                               
                                fxN2 = f_lambda。*φ;
                               
                                IntN1 =双(vpaintegral (fxN1 x,负无穷,Cstar));
                               
                                IntN2 =双(vpaintegral (fxN2 x,负无穷,Cstar))
                               
                                  IntN2 = 0.1042

                                IntN = IntN1-IntN2;
                               
                                TN (s, r) = IntN;
                               
                                disp (TN)
                               
                                  列1到19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0列20到29日0 0 0 0 0 0 0 0 0 0.2809

弗朗西斯科·法夫里 9分钟前

Dyuma,非常感谢你的好的建议。效果好,解决了我的问题!它仍然是为什么的神秘事务r = 24 - > 25防止正确的计算

登录置评。

在回答这个问题。

答案(1)

VBBV 10分钟前

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1976269-trouble-with-the-double-numeric-data-type answer_1247794

编辑:VBBV 3分钟前

以下版本的代码运行没有错误。

                              清晰;
                             
                              信谊x t
                             
                              % r = 24;
                             
                              deltamu = 0.5;
                             
                              σ= 5;
                             
                              λ= 1;
                             
                              Cstar = 7;
                             
                              μ= (4:deltamu: 10)
                             
                              P = 0(1、长度(μ));
                             
                              TP = 0(1、长度(μ));
                             
                              TN = 0(1、长度(μ));
                             
                              s = 1;
                             
                              sigmae = 5;
                             
                              exp_t = exp (- t ^ 2);
                             
                              达峰时间= +正;
                             
                              为r = 1:长度(亩)
                             
                              tmin =((μ(r) +λ。*(σ^ 2)- x)。/ (sqrt(2) *σ));
                             
                              误差补函数=(2. /√(π)。* int (exp_t t tmin,达峰时间);
                             
                              f_lambda = (lambda. / 2)。* exp ((lambda. / 2) *(2。*μ(r) +λ。*(σ^ 2)- x))。*误差补函数;
                             
                              P (s, r) =双(int (f_lambda, x, Cstar +正));
                             
                              phi_e = 1. / (sigmae。*√(2。* pi))。* exp (- (1/2)。* ((x-Cstar)。/ sigmae)。^ 2);
                             
                              φ= int (phi_e x,负无穷,x);
                             
                              fxP = f_lambda。*φ;
                             
                              TP (s, r) =双(int (fxP, x, Cstar +正));
                             
                              fxN1 = f_lambda;
                             
                              fxN2 = f_lambda。*φ;
                             
                              IntN1 =双(int (fxN1 x,负无穷,Cstar));
                             
                              IntN2 =双(int (fxN2 x,负无穷,Cstar));
                             
                              IntN = IntN1-IntN2;
                             
                              TN (s, r) =双(IntN)
                             
                              结束

如果我理解正确的话,你的意思是写这行代码

                              % - - - - - - - - - - - - - - - - - - - - - - - > > - - - - - - - > >
                             
                              tmin =((μ(r) +λ。*(σ^ 2)- x)。/ (sqrt(2) *σ));
                             
                              % - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - > > - - - - - - - > >
                             
                              f_lambda = (lambda. / 2)。* exp ((lambda. / 2) *(2。*μ(r) +λ。*(σ^ 2)- x))。*误差补函数;