Handling a very big difference between numbers(ratio)

Question

Vadim Patrick Nave 2022-8-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1789470-handling-a-very-big-difference-between-numbers-ratio

编辑： Infinite_king 2023-12-4

Hello,

I was working on a code that demands operating with statistics of long numbers (such as 2000000) vs the errors(such as 0.0003) so different equations that demands both of the numbers are equal 0 or Inf, how can I improve my situation?

I tried to use vpa, it helps but it's not a magic pill and the problems appers later in the code.

Is there anyway to operate with such numbers in terms of division without getting results as Inf or 0?

Thank you and a blessed week,

Vadim

5 个评论
显示 3更早的评论隐藏 3更早的评论

Steven Lord 2022-8-29

I was working on a code that demands operating with statistics of long numbers (such as 2000000) vs the errors(such as 0.0003) so different equations that demands both of the numbers are equal 0 or Inf, how can I improve my situation?

It's likely to be difficult if not impossible to offer any specific suggestions without seeing the specific equations you're using.

Vadim Patrick Nave 2022-8-29

在 MATLAB Online 中打开

Thank you Steven, you are right. Sorry.

While dy is a very small number

x is very large number

% Fit line analytically:
matrix = [ sum(x.^2./dy.^2)       sum(x./dy.^2);
    sum(x./dy.^2)          sum(1./dy.^2)];
y_sigma = y./dy.^2;
free_vector = [sum(x.*y_sigma);
    sum(1.*y_sigma)];
solution=matrix\free_vector;           
a=solution(1);b=solution(2);
err=sqrt(diag(inv(matrix)));
da=err(1);db=err(2);
% if there are errors in x, co-add them and re-fit
if nargin>10
    dy = sqrt((a*dx).^2+dy.^2);
    matrix = [ sum(x.^2./dy.^2)       sum(x./dy.^2);
        sum(x./dy.^2)          sum(1./dy.^2)];
    y_sigma = y./dy.^2;
    free_vector = [sum(x.*y_sigma);
        sum(1.*y_sigma)];
    solution=matrix\free_vector;           
    a=solution(1);b=solution(2);
    err=sqrt(diag(inv(matrix)));
    da=err(1);db=err(2);
end

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Infinite_king 2023-12-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1789470-handling-a-very-big-difference-between-numbers-ratio#answer_1365029

编辑：Infinite_king 2023-12-4

在 MATLAB Online 中打开

Hi Vadim Patrick Nave,

I understand that you want to work with very large and very small numbers and perform arithmetic operations on them without running into ‘Inf’ or ‘nan’ values.

I suggest you to use ‘sym’ function which was available in ‘Symbolic Math Toolbox’. First convert the number to a symbolic number or matrix to symbolic matrix. Then you can perform simple arithmetic operations and finally you can use ‘double’ function to convert the answer to double.

Refer below code snippet,

% let x be a matrix of numbers
x = rand(5);
disp(x);
    0.6448    0.1331    0.2477    0.0948    0.9077
    0.0157    0.2160    0.5678    0.6035    0.4680
    0.1102    0.4070    0.8393    0.2991    0.9093
    0.2817    0.2833    0.9535    0.0625    0.5194
    0.2113    0.5092    0.7433    0.0416    0.2276
% now convert the matrix to symbolic matrix
x_sym = sym(x);
% now perform simple arithmetic operations
% op 1
% op 2
% for example, addition
x_sym = x_sym + 5;
% now convert the values to double
% make sure the numbers are within range of double
res = double(x_sym);
disp(res);
    5.6448    5.1331    5.2477    5.0948    5.9077
    5.0157    5.2160    5.5678    5.6035    5.4680
    5.1102    5.4070    5.8393    5.2991    5.9093
    5.2817    5.2833    5.9535    5.0625    5.5194
    5.2113    5.5092    5.7433    5.0416    5.2276

For more information on how to use ‘sym’ function and ‘Symbolic Math Toolbox’, please refer the following MATLAB documentations,

Hope this is helpful.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Walter Roberson 2023-12-4

编辑：Walter Roberson 2023-12-4

在 MATLAB Online 中打开

and final you can use ‘eval’ function to convert the answer to double.

No you cannot. There is no documented meaning for eval() of a symbolic expression.

syms x

f = piecewise(1 <= x & x <= 5, x, x^2)

f =

char(f)

ans = 'piecewise(x in Dom::Interval([1], [5]), x, symtrue, x^2)'

eval(f)

Error using sym/eval
Invalid expression. Check for missing multiplication operator, missing or unbalanced delimiters, or other syntax error. To construct matrices, use brackets instead of parentheses.

Symbolic expressions displayed to the user are in a user-friendly-ish language, not in a programming language. They are not expressions in MATLAB, and by the time they reach the user, they are not expressions in the internal MuPAD programming language. They can use constructs that are not accepted by MATLAB, such as the above in-fix notation use of the "in" operator, and such as the :: construction . The order of parameters they show might not be the same as the order that a numeric function of the same name expects, and might not be the same as the order that a symbolic function of the same name expects, and might not be the same as the order that the internal MuPAD programming language expects.

Walter Roberson 2023-12-4

在 MATLAB Online 中打开

x = 5;
x_sym = sym(x);
x_sym = x_sym * 5;
x_sym = x_sym + 5;
x_sym = x_sym / 5;
% how to evaluate x_sym, assuming the resulting value is within ranage of
% double.
double(x_sym)
ans = 6

No eval() needed.

Handling a very big difference between numbers(ratio)

5 个评论
显示 3更早的评论隐藏 3更早的评论

回答（1 个）

6 个评论
显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Handling a very big difference between numbers(ratio)

5 个评论 显示 3更早的评论隐藏 3更早的评论

回答（1 个）

6 个评论 显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

5 个评论
显示 3更早的评论隐藏 3更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论