rat

Рациональное дробное приближение

Описание

пример

R = rat(X) возвращает рациональное дробное приближение X к в допуске по умолчанию, 1e-6*norm(X(:),1). Приближение является символьным массивом, содержащим усеченное непрерывное дробное разложение.

пример

R = rat(X,tol) аппроксимирует X к в допуске, tol.

пример

[N,D] = rat(___) возвращает два массива, N и D, таким образом, что N./D аппроксимирует X, использование любого из вышеупомянутых синтаксисов.

Примеры

свернуть все

Аппроксимируйте значение π использование рационального представления количества pi.

Математическое количество π не рациональное число, но количество pi это аппроксимирует его, рациональное число, поскольку все числа с плавающей запятой рациональны.

Найдите рациональное представление pi.

format rat
pi
ans = 
     355/113   

Получившееся выражение является вектором символов. Также можно использовать rats(pi) получить тот же ответ.

Используйте rat видеть непрерывное дробное разложение pi.

R = rat(pi)
R = 
'3 + 1/(7 + 1/(16))'

Результатом является приближение непрерывным дробным разложением. Если вы рассматриваете первые два члена разложения, вы получаете приближение 3+17=227, который только соглашается с pi с 2 десятичными знаками в дробной части.

Однако, если вы рассматриваете все три члена, распечатанные rat, можно восстановить значение 355/113, который соглашается с pi с 6 десятичным знакам в дробной части.

3+17+116=355113

Задайте допуск к дополнительной точности в приближении.

R = rat(pi,1e-7)
R = 
'3 + 1/(7 + 1/(16 + 1/(-294)))'

Получившееся приближение, 104348/33215, соглашается с pi к 9 десятичным числам.

Создайте матрицу 4 на 4.

format short;
X = hilb(4)
X = 4×4

    1.0000    0.5000    0.3333    0.2500
    0.5000    0.3333    0.2500    0.2000
    0.3333    0.2500    0.2000    0.1667
    0.2500    0.2000    0.1667    0.1429

Выразите элементы X как отношения маленьких целых чисел с помощью rat.

[N,D] = rat(X)
N = 4×4

     1     1     1     1
     1     1     1     1
     1     1     1     1
     1     1     1     1

D = 4×4

     1     2     3     4
     2     3     4     5
     3     4     5     6
     4     5     6     7

Эти две матрицы, N и D, аппроксимированный X с N./D.

Просмотрите элементы X как отношения с помощью format rat.

format rat
X
X = 
       1              1/2            1/3            1/4     
       1/2            1/3            1/4            1/5     
       1/3            1/4            1/5            1/6     
       1/4            1/5            1/6            1/7     

В этой форме это ясно тот N содержит числители каждой дроби и D содержит знаменатели.

Входные параметры

свернуть все

Входной массив в виде числового массива класса single или double.

Типы данных: single | double
Поддержка комплексного числа: Да

Допуск в виде скаляра. N и D аппроксимированный X, таким образом, что N./D - X < tol. Допуском по умолчанию является 1e-6*norm(X(:),1).

Выходные аргументы

свернуть все

Непрерывная дробь, возвращенная как символьный массив с m строки, где m число элементов в X. Точность рационального приближения через непрерывные дроби увеличивается с количеством членов разложения.

Числитель, возвращенный как числовой массив. N./D аппроксимирует X.

Знаменатель, возвращенный как числовой массив. N./D аппроксимирует X.

Алгоритмы

Даже при том, что все числа с плавающей запятой являются рациональными числами, иногда желательно аппроксимировать их простыми рациональными числами, которые являются дробями, числитель которых и знаменатель являются маленькими целыми числами. Рациональные приближения сгенерированы путем усечения расширений непрерывной дроби.

rat функция аппроксимирует каждый элемент X непрерывной дробью формы

ND=D1+1D2+1+1Dk.

D s получен путем повторного собирания целой части и затем взятия обратной величины дробной части. Точность приближения увеличивается экспоненциально с количеством членов разложения и хуже когда X = sqrt(2). Для X = sqrt(2) , ошибка с k условия о 2.68*(.173)^k, таким образом, каждое дополнительное условие увеличивает точность меньше чем на одну десятичную цифру. Требуется 21 условие, чтобы получить полную точность с плавающей точкой.

Расширенные возможности

Смотрите также

|

Представлено до R2006a