Интерполяция функции по формулам Лагранжа

@Arseniiy · Регистрация: 03.05.2012

Студворк — интернет-сервис помощи студентам

Имеется код Полинома Лагранжа. Что необходимо поправить чтоб можно было вводить X и Y массивы. и получать вывод график Полинома Лагранжа с узлами интерполяции

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear
function L_n = Lag2 (x,y,xx)
%x-массив координат узлов
%y-массив значений интерполируемой функции
N=length(x); % узнаем число узлов интерполяции (N=n+1)
L_n = zeros(size(xx));
s=0;
for k=1:N
    P=ones(size(xx));
  for i=1:N
      if k~=i
          P=P.*(xx-x(i))./(x(k)-x(i));
          %.* - поэлементное умножение
          %./ - поэлементное деление
      end
  end
    L_n = L_n + y(k)*P;
end
end

@Arseniiy · 04.05.2012, 23:45 **[ТС]**

Понял что программа не работает из за того, что не заданна функция Lag2 которая является функцией ввода для пользователей. Как её организовать? Нужно ввести x и y

@Arseniiy · 09.05.2012, 20:12 **[ТС]**

Программа упростилась.

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
clear
a=0;
b=3.5;
x=[0 0.25 1.25 2.125 3.25]; %x-массив координат узлов
y=[5 4.6 5.7 5.017 4.33]; %y-массив значений интерполируемой функции
xx=linspace(a,b,1000);
N=length(x); % узнаем число узлов интерполяции
L_n = zeros(size(xx)); % создаем нулевой массив значений интерполяционного полинома
% в цикле считаем сумму по узлам
for k=1:N
    % вычисляем произведения
    P=ones(size(xx));
  for i=1:N
      if k~=i
          P=P.*(xx-x(i))./(x(k)-x(i));
          %.* - поэлементное умножение
          %./ - поэлементное деление
      end
  end
    L_n = L_n + y(k)*P; % накапливаем сумму
end
plot(xx,L_n,'r')
hold on
plot(x,y,'bo')

Помогите переделать. Хочется чтобы был диалог с пользователем, а не изменение вручную в коде программы переменных x,y,a,b.

Добавлено через 10 часов 8 минут
Как ввести массив чисел с помощью функции input, или другой функции но чтоб это был массив чисел? Например такой массив [2 3 4 5 6]

@Olzh · 10.04.2014, 06:04

>> M{end+1}=input('введите что-то: ') типо вот те то, что ты хотел...

@Olzh · 21.04.2014, 09:06

хотя так думаю будет проще
x=input('введите значение икс в квадратных скобках ')

Новые блоги и статьи Все статьи Все блоги /
Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .	Фото: Daniel Greenwood kumehtar 13.11.2025

@Arseniiy 1 / 1 / 0 Регистрация: 03.05.2012 Сообщений: 3
	04.05.2012, 23:45 [ТС]
	Понял что программа не работает из за того, что не заданна функция Lag2 которая является функцией ввода для пользователей. Как её организовать? Нужно ввести x и y 0

@Olzh 0 / 0 / 0 Регистрация: 10.04.2014 Сообщений: 2
	10.04.2014, 06:04
	>> M{end+1}=input('введите что-то: ') типо вот те то, что ты хотел... 0

@Olzh 0 / 0 / 0 Регистрация: 10.04.2014 Сообщений: 2
	21.04.2014, 09:06
	хотя так думаю будет проще x=input('введите значение икс в квадратных скобках ') 0