Вопросы с меткой [численные-методы]
Численные методы — методы решения математических задач в численном виде Представление как исходных данных в задаче, так и её решения — в виде числа или набора чисел.
60
вопросов
1
голос
0ответов
15 показов
W-функция Ламберта
При вычислении уравнения 2x - exp(-x) = 2, получилось, что x = (1 / 2exp(x)) + 1. Потом решил проверить в Wolfram, а там решение в виде W-функции Ламберта, я немного не понимаю как мы приходим к виду ...
2
голоса
0ответов
86 показов
Метод Рунге-Кутта и метод Милна(поиск шага интегрирования)
Мне нужно написать метод Рунге Кутта(4 порядка) и метод Милна для решения задачи Коши, сами методы несложные, а вот подбор шага интегрирования поставил меня в тупик, я так понял, что нужно найти ...
1
голос
1ответ
56 показов
Экстраполяция синуса и косинуса
Только что закончилось одно мини-соревнование по машинному обучению, и сейчас я хочу узнать, как решается одна задача оттуда.
Горизонтальные координаты парашютиста от времени имеют следующий вид: x = ...
0
голосов
1ответ
36 показов
не могу вывести корень уравнения, используя метод половинного деления
Возможно я где-то туплю. Но на данный момент я не понимаю...Решаю эту задачу уже третий час. Моя программа видит только одно: интервал. Не корень.
Задача: Написать программу которая находит решение ...
1
голос
0ответов
49 показов
Неточность в нахождении минимума функции в R
Имеется функция:
Holt <- function(data_table, coef)
{
alpha <- coef[1]
beta <- coef[2]
ldt = length(data_table)
ut = bt = numeric(ldt)
for(i in (1:(ldt)))
{
if(i == 1)
{
...
7
голосов
3ответа
362 показа
Какой алгоритм использовать для решения системы с приведенной матрицей коэффициентов?
Какой алгоритм поиска решения (x1, x2, ...) использовать для системы с матрицей коэффициентов как на приведенной картинке?
-1
голос
1ответ
60 показов
Применение метода Вегстейна
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
double f(double x)
{
return (x*x)*log10(x)-3.8;
}
double psi(double x, double *x1)
{
double h=(*x1-x);
*x1=x;
...
0
голосов
0ответов
28 показов
Не работает реализация Рунге-Кутты 4 порядка MATHLAB
В первой итерации рассчитываются шаги k1,k2,k3,k4, но как только я добавляю последний шаг формулы метода РК X(i+1) = X(i) + 1/6.*(k1X + 2.*k2X + 2.*k3X + k4X) все ломается и в значениях k появляется ...
0
голосов
1ответ
59 показов
Как конвертировать c/c++ заголовочные .h файлы для OpenBLAS/LAPACK в c# библиотеку
Начнем с того, что lapack.h содержит 13 717 строк и описание 1 828 функций. Вручную заменить каждое с/с++ описание
#define LAPACK_cbbcsd LAPACK_GLOBAL(cbbcsd,CBBCSD)
void LAPACK_cbbcsd(
char const*...
1
голос
1ответ
48 показов
Как при помощи x64 версии OpenBLAS/LAPACK решать стационарные уравнения теплопроводности?
Итак - x64-версия библиотеки OpenBLAS.dll прекрасно подключилась в x64 c#-проект (Подключение фортрановской LAPACK x64 в x64 проект на C#). Как же теперь с её помощью решать стационарные (и даже ...
4
голоса
2ответа
580 показов
Почему не стоит пользоваться функцией Pow для целых чисел?
p.s. Если что я говорю про C++ функцию Pow.
Мне не раз приходилось встречать рекомендации по типу "не пользуйтесь функцией Pow". В связи с чем многие дают эту рекомендацию?
void MyPow(int&...
2
голоса
1ответ
84 показа
Как учитывать погрешность при вычислениях значения функции
Задается x и задается некая погрешность, например, 0,001. Моя погрешность вычисления функции не должна превышать заданную. Как мне учитывать это при вычислениях?
#include <iostream>
int main()...
0
голосов
1ответ
430 показов
Сумма ряда в С++
Задача. Сумма ряда от 1 до n 1/x^i
В консоль выводит ноль, не могу найти ошибку. Объясните, пожалуйста
#include <iostream>
using namespace std;
void main()
{
setlocale(LC_ALL, "ru");
...
1
голос
1ответ
361 показ
matlab метод Эйлера для решения диф уравнения
Реализовал метод Эйлера, который работает для дифференциального уравнения номер 1, а для второго выдает ошибку, помогите, пожалуйста!
clear all
clc
fun='t.^2*cos(t)+2/t';
fun=inline(fun);
t(1)=0.1;...
0
голосов
0ответов
81 показ
Производная от многочлена Ньютона с++
есть код такой вот код , приближаю функцию многочленном Лагранжа , нужна производная от такого многочлена в точке с = (2pi/2)+0.2 как можно ее задать в с++?
#include <iostream>
#include <...
0
голосов
0ответов
57 показов
Степенной метод , нахождение второго по величине собственного значения
Добрый день(вечер)всем .
Есть у меня код "степенной метод" по поиску максимального по модулю собственного значения и собственых векторов матрицы .
Хочу доработать его до поиска второго по величине а ...
1
голос
1ответ
135 показов
Как реализовать алгоритм трилинейной интерполяции для 3-х переменных?
у меня есть код, на с++, реализации билинейной интерполяции для двух переменных
class point{
public:
double x,y,f;
void init(double x1,double y1,double f1){
x=x1;
y=y1;
...
1
голос
1ответ
46 показов
Подскажите, пожалуйста, книги про численные методы или математическим моделям с примерами кода на с++/python
Потребовалось реализовать математическую модель какого-либо физического явления.
Подскажите, пожалуйста, литературу (желательно с примером кода на Python / C++), если таковая существует, или ссылки ...
2
голоса
1ответ
388 показов
Метод Рунге-Кутта 4-го порядка для CДУ C#
Нужно решить систему из двух дифференциальных уравнений методом Рунге-Кутты 4-го порядка с коррекцией шага по правилу Рунге.
Так как мне нужно корректировать шаг,то я хотела написать метод rk4,как и ...
2
голоса
1ответ
56 показов
Почему после i = 5 выводится одинаковый результат?
Недавно начал изучать C, не могу понять, почему после i = 5 выводится одинаковый результат?
#include <stdio.h>
#include <math.h>
#define SIZE 100
int main()
{
int x = 2;
...
-3
голоса
1ответ
37 показов
Помогите понять как решать такого типа задачи по численным методам
При каких a и b сходится метод Гаусса-Зейделя для системы с матрицей
A=(a,b,0; b,a,0; 0,0,a) в норме ||.||{1} ({1} - нижний индекс 1)
5
голосов
0ответов
59 показов
Конкретный пример использования собственных векторов
Хотелось бы увидеть пример использования в той же физике на конкретном примере.
В моем представлении это должно выглядеть как-то так: у нас есть n физических объектов с определенными свойствами, на ...
1
голос
0ответов
85 показов
Линейная интерполяция для рассчета траектории движения обьекта на canvas
Пытаюсь написать свою игру на canvas, но никак не могу нормально рассчитать траекторию полета снарядов при выстреле. По задумке - снаряд летит с фиксированной скоростью по вектору от текущих координат ...
6
голосов
1ответ
925 показов
Численное решение задачи теплопроводности
Необходимо решить задачу теплопроводности на отрезке
При решении использовала явную схему
N = 10 # максимальное число шагов по х
K = 10 # максимальное число шагов по t
l = 1 # значение х на правой ...
1
голос
1ответ
730 показов
Поиск корня функции методом дихотомии
Выполняю задание курсовой работы по Численным Методам и Моделированию на ЭВМ: найти корень функции и отобразить на графике функции методами итераций, дихотомии и Ньютона. Для данного задания (и ...
1
голос
1ответ
70 показов
Решить нелинейное уравнение методом хорд
Есть уравнение x^2 = sin(x). Нужно написать программу,которая находит корни на указанном промежутке. Использовать нужно эту формулу -
Вот мой код
float e = 0.0001;// погрешность
float x0 = -...
1
голос
2ответа
112 показов
Как уменьшить ошибку интерполяции на границах интервала
Необходимо выполнить интерполяцию данных кубическим сплайном. Все вроде хорошо, но полиномы на концах имеют уж слишком большую ошибку. Что в таком случае обычно делают?
Пример интерполяции синуса и ...
2
голоса
1ответ
393 показа
В чём отличие между аналитическим методом рядов для решения задачи Коши для ОДУ и численным методом разложения в ряд Тейлора?
Я так понимаю, общая формула (7.7) представляет собой аналитический метод, тогда как выглядит численный?
Или то, когда я решаю конкретную задачу на каком-то примере - считаю производные до нужного ...
-7
голосов
1ответ
84 показа
Ряд Тейлора, написал как знаю, кому не трудно гляньте плз [закрыт]
#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <math.h>
double f(long n, long x) {
long tt = (2*n)+1;
long b = 1, c;
double res;
res = pow(x, (2*n+1));
...
-1
голос
1ответ
112 показов
Вычислить корни уравнения ln(x)+x-0.5=0 методом Ньютона с заданной точностью eps [закрыт]
Накидал кода, но правильного ответа нету. Где может быть ошибка?)
#include "pch.h"
#include <iostream>
#include <cmath>
using namespace std;
double fx, f, fa, x, e;
int main()
{
...
-2
голоса
2ответа
46 показов
Не происходит расчет матриц определенного типа
#include<iostream>
#include <math.h>
#include <stdlib.h>
using namespace std;
int main() {
int i, j, n, m;
//создаем массив
cout << "vvedite chislo yravneni: ";
...
1
голос
0ответов
241 показ
Численное решение уравнения колебания волны методом сеток С++
Согласно заданию, мне необходимо решить уравнение колебания волны методом сеток. Для метода был написан код:
#include <iostream>
#include <cmath>
using namespace std;
double (*p) (...
1
голос
1ответ
94 показа
Метод прямоугольника не параллелится
Здравствуйте пишу простой алгоритм численного метода (прямоугольника). Распараллелил эту функцию с помощью OpenMP. Но время выполнения параллельного алгоритма больше, чем выполнение с одним поток. В ...
1
голос
1ответ
1k показов
Метод покоординатного спуска (метод Гаусса - Зейделя) C#
Есть алгоритм, который нужно реализовать. Основа для кода готова, смотрите ниже. Проблема заключается в шагах 3, 8. Как понять эти шаги, что используется в C# для нахождения argmin? z равняется 0.5
...
0
голосов
1ответ
167 показов
Алгоритм регрессии для логистичекой функции (4PL или 5PL)
Необходим алгоритм регрессии (нахождения коэффициентов) для 4 или 5 параметровой логистической функции по методу наименьших квадратов.
F(x) = D+(A-D)/((1+(x/C)^B)^E)
Даю уточнения. Алгоритм нужен ...
7
голосов
1ответ
458 показов
Алгоритм для гравитационной задачи N тел
К своему стыду очень слаб в численных методах.
Суть задачи в том, чтобы смоделировать движение тел с течением времени под воздействием силы гравитации друг на друга в количестве >= 3.
В данный момент ...
1
голос
1ответ
151 показ
Как написать подсчет Определителя Больших Матриц методом Гаусса с выбором глобального главного элемента?
Нужно написать подсчет определителя матрицы, с помощью метода Гаусса с выбором глобального главного элемента.
Все алгоритмы которые я нахожу, это просто метод Гаусса. А мне нужен частный случай с ...
4
голоса
1ответ
89 показов
Как называется по-русски распределение, плотность которого 1/x?
Ричард Хэмминг опубликовал в 1970 году статью "On the Distribution of Numbers". В этой статье он решил вопрос о том, как выглядит распределение мантисс всех представимых чисел в арифметике с плавающей ...
-3
голоса
1ответ
2k показов
Вычисление функции c помощью ряда Тейлора
Вещественная функция может быть представлена рядом Тейлора, радиус сходимости которого определяется положением особых точек функции. Любой конечный отрезок ряда Тейлора является полиномом, поэтому ...
3
голоса
1ответ
597 показов
Интерактивные Графики в Python: как реализовать аналог функции Manipulate из Mathematica
Добрый день! Недавно начал изучать Python, решая задачу по численным методам. Хотел узнать, можно ли как-то в Python реализовать аналог функции Manipulate как в Mathematica, чтобы была возможность ...
1
голос
1ответ
658 показов
Проблема вывода в Python: не выводит сообщение, которое написано в самом конце функции
Добрый день!
Пишу программу, которая бы решала нелинейное уравнение методом Ньютона(касательных). Проблема не в реализации данного метода, а в выводе данных.
Почему то не выводит сообщение, которое ...
2
голоса
1ответ
5k показов
Таблицы в Python
Добрый день! Недавно начал изучать Python, и возникла необходимость, загрузить данные, которые вычисляются в функциях в таблицу.
Приведу конкретно мой случай: Написал программу которая вычисляет ...
6
голосов
3ответа
3k показов
Хранение разреженной матрицы в разреженном строчном формате (CSR или CRS)
Здравствуйте.
Подскажите, пожалуйста, алгоритм упаковки разреженной матрицы в разреженном строчном формате (Compressed Sparse Rows или Compressed Row Storage) и способ доступа к элементам упакованной ...
5
голосов
1ответ
349 показов
Поиск корней полинома
Как найти корни полинома, заданного через boost::math::tools::polynomial<double> на промежутке [0,1] используя тот же boost? Чтение документации не особо помогло в решении.
-2
голоса
1ответ
107 показов
Небольшая программа на C++ [закрыт]
Доброго времени суток всем!
Давно учил С++, уже его подзабыл, но сейчас необходимо решить задание. Буду благодарен любой помощи!
Задание:
Есть нелинейное уравнение cos(2/x) - 2sin(1/x) + 1/x = 0 , ...
2
голоса
1ответ
2k показов
Численные методы (метод хорд)
Было задание найти корни уравнения методом хорд. Код выглядит следующим образом:
#include <iostream>
#include <cmath>
using namespace std;
double f(double x)
{
return pow(cos(x), 3) ...
2
голоса
1ответ
310 показов
Найти погрешность функции
Оценить абсолютную и относительную погрешность вычисления величины F,
если для x1, x2, x3 погрешность Δ = 0.007
Как я понимаю, нужно найти абсолютную погрешность от возведения в степень
...
3
голоса
1ответ
658 показов
Какие существуют численные методы для метода наименьших квадратов?
У меня имеется следующий код в Matlab:
fileID = fopen('input.txt','r');
formatSpec = '%f';
y = fscanf(fileID,formatSpec);
step = 0.1;
x0 = step;
xn = length(y)*0.1;
x = x0:step:xn;
fitfunc = 'a + ...
0
голосов
0ответов
54 показа
Вопрос о модификациях метода Эйлера
Доброго времени суток! Помогите пожалуйста с ответом на один вопрос. Уверен тут есть люди знающие... Я надеюсь что вопрос не удалят, так как он не о программировании.
В общем суть вопроса: всем ...
1
голос
1ответ
73 показа
Волновое уравнение - проблема реализации
Задача: решить уравнение колебания струны методом конечных разностей.
Условие устойчивости:
dt = dx/(2*A);
Первое начальное условие (функция натяжения струны):
Uprev = UAllocate(n);
(UAllocate() - ...