Увеличьте MATLAB®скорость путем сокращения точности вычислений. Уменьшайте точность при помощи арифметики переменной точности, обеспеченной vpa
и digits
функции в Symbolic Math Toolbox™. Когда вы уменьшаете точность, вы получаете эффективность путем сокращения точности. Для получения дополнительной информации смотрите, Выбирают Numeric or Symbolic Arithmetic.
Например, находя Дзета-функцию Римана большого матричного C
занимает много времени. Во-первых, инициализируйте C
.
[X,Y] = meshgrid((0:0.0025:.75),(5:-0.05:0)); C = X + Y*i;
Затем найдите время потраченным, чтобы вычислить zeta(C)
.
tic zeta(C); toc
Elapsed time is 340.204407 seconds.
Теперь повторите эту операцию с более низкой точностью при помощи vpa
. Во-первых, измените точность, используемую vpa
к более низкой точности 10
цифры при помощи digits
. Затем используйте vpa
уменьшать точность C
и найдите zeta(C)
снова. Операция значительно быстрее.
digits(10) vpaC = vpa(C); tic zeta(vpaC); toc
Elapsed time is 113.792543 seconds.
Примечание
vpa
выход является символьным. Чтобы использовать символьный выход с функцией MATLAB, которая не принимает символьные значения, преобразуйте символьные значения в двойную точность при помощи double
.
Для больших матриц различие во время вычисления может быть четным более значительное. Например, рассмотрите 1001
- 301
матричный C
.
[X,Y] = meshgrid((0:0.0025:.75),(5:-0.005:0)); C = X + Y*i;
Выполнение zeta(vpa(C))
с 10-разрядной точностью занимает 15 минут, при выполнении zeta(C)
берет в три раза более долго.
digits(10) vpaC = vpa(C); tic zeta(vpaC); toc
Elapsed time is 886.035806 seconds.
tic zeta(C); toc
Elapsed time is 2441.991572 seconds.
Примечание
Если вы хотите увеличить точность, смотрите Точность Увеличения Числовых Вычислений.