如何阻止MATLAB将极小值舍入为0？

Question

我在MATLAB中有一个代码,它可以处理非常小的数字,例如,我的值大约为10 ^ { - 25},但是当MATLAB进行计算时,值本身会舍入为0.注意,我我不是指format显示这些额外的小数,而是将数字本身更改为0.我认为原因是因为默认情况下,MATLAB在小数点后最多使用15位进行计算.如何更改此值以便保留非常小的数字,因为它们在计算中？

编辑:

我的代码如下:

clc;
clear;

format long;

% Import data
P = xlsread('Data.xlsx', 'P');
d = xlsread('Data.xlsx', 'd');
CM = xlsread('Data.xlsx', 'Cov');


Original_PD = P;                       %Store original PD
LM_rows = size(P,1)+1;                           %Expected LM rows
LM_columns = size(P,2);                 %Expected LM columns
LM_FINAL = zeros(LM_rows,LM_columns);    %Dimensions of LM_FINAL

for ii = 1:size(P,2)

P = Original_PD(:,ii);
% c1, c2, ..., cn, c0, f
interval = cell(size(P,1)+2,1);

for i = 1:size(P,1)
    interval{i,1} = NaN(size(P,1),2);
    interval{i,1}(:,1) = -Inf;
    interval{i,1}(:,2) = d;

    interval{i,1}(i,1) = d(i,1);
    interval{i,1}(i,2) = Inf;
end
interval{i+1,1} = [-Inf*ones(size(P,1),1) d];
interval{i+2,1} = [d Inf*ones(size(P,1),1)];

c = NaN(size(interval,1),1);
for i = 1:size(c,1)
    c(i,1) = mvncdf(interval{i,1}(:,1),interval{i,1}(:,2),0,CM);
end

c0 = c(size(P,1)+1,1);
f = c(size(P,1)+2,1);
c = c(1:size(P,1),:);
b0 = exp(1);
b = exp(1)*P;

syms x;

eqn = f*x;
for i = 1:size(P,1)
    eqn = eqn*(c0/c(i,1)*x + (b(i,1)-b0)/c(i,1));
end

eqn = c0*x^(size(P,1)+1) + eqn - b0*x^size(P,1);

x0 = solve(eqn);
x0 = double(x0);

for i = 1:size(x0)
    id(i,1) = isreal(x0(i,1));
end

x0 = x0(id,:);
x0 = x0(x0 > 0,:);

clear x;

for i = 1:size(P,1)
    x(i,:) = (b(i,1) - b0)./(c(i,1)*x0) + c0/c(i,1);
end

% x = [x0 x1 ... xn]
x = [x0'; x];
x = x(:,sum(x <= 0,1) == 0);

% lamda
lamda = -log(x);
LM_FINAL(:,ii) = lamda;
end

问题出在这一步:

 for i = 1:size(P,1)
        x(i,:) = (b(i,1) - b0)./(c(i,1)*x0) + c0/c(i,1);
end

"差异"非常接近0.如何在此步骤中停止此舍入？

例如,当i = 10时,我有以下值:

b_10 = 0.006639735483297
b_0 = 2.71828182845904
c_10 = 0.000190641848119641
c_0 = 0.356210110252579
x_0 = 7.61247930625269

在进行计算后,我们得到:-1868.47805854794 + 1868.47805854794,它产生-2.27373675443232E-12的差值,由MATLAB舍入为0.

编辑2:

这是我的数据文件,用于代码.运行代码后(大约需要一分半钟才能完成运行),11变量x中的行显示为0(即使在双击后检查它的实际值),何时不应该.

Answer 1

您遇到的问题是因为浮点的IEEE标准无法将您的数字与零区分开来,因为它们没有使用足够的位.

看看John D'Errico的大十进制类和可变精度整数算法.另一种选择是使用Java中的Big Integer Class,但如果您不熟悉在MATLAB中使用Java和其他rexternal库,则可能更具挑战性.

你能举例说明你使用1e-25并得到零的计算吗？这是我得到的浮点数small_num和John的高精度浮点数之一,small_hpf在分配它们并乘以时调用pi.

>> small_num = 1e-25

small_num =

   1.0000e-25

>> small_hpf = hpf(1e-25)
small_hpf =
1.000000000000000038494869749191839081371989361591338301396127644e-25
>> small_num * pi

ans =

   3.1416e-25

>> small_hpf * pi
ans =
3.141592653589793236933163473501228686498684350685747717239459106e-25