Аппроксимация рядом Фурье

@hatsumi · Регистрация: 20.01.2016

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

Здравствуйте, есть задание на аппроксимацию исходных данных отрезком ряда Фурье
Исходные данные: массив координат Х и массив значений функции Y

Собственно вот, что получилось:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
function [] = approx ()
clear;
clc;
close all;
%исходные данные
x = 0:1:10;
y = [5 5 7 3 7 9 7 4 3 2 -3];
m = 11; %кол-во точек
T = 11;
n = 5;
%интервал аппроксимации с шагом 0.1
xgrid = 0:0.1:10;
ygrid = zeros(size(xgrid));
%инициализация и расчет
A = zeros(1, 6);
B = zeros(1, 6);
summAi = 0;
summBi = 0;
summ = 0;
 
for i = 0:5
    for k = 0:m-1
        summAi = summAi + y(k+1)*cos(i*2*pi*k/m);
        summBi = summBi + y(k+1)*sin(i*2*pi*k/m);
    end
    A(i+1) = (2/m)*summAi;
    B(i+1) = (2/m)*summBi;
end
 
for k = 1:length(xgrid)
    for i = 2:n
        summ = summ + A(i)*cos(i*2*pi*xgrid(k)/T) + B(i)*sin(i*2*pi*xgrid(k)/T);
    end
    ygrid(k) = A(1)/2 + summ;
end
plot(xgrid, ygrid)
hold on
plot(x, y, 'bo')
end

Но получается далеко не совсем то, что нужно
Подскажите в чем может быть проблема, вроде делал все по формулам..

@Nick07 · 16.10.2017, 22:55

Вам правильный код или найти ошибку?
Если правильный код, то в cftool

@hatsumi · 16.10.2017, 23:58 **[ТС]**

Мне бы хотелось все-таки ошибку в моем коде найти

Добавлено через 1 час 0 минут
Одно понял точно, summ надо обнулять в теле цикла, иначе переменная накапливает все суммы, но вот в чем еще загвоздка не могу понять..

@SSC · 17.10.2017, 07:18

Сообщение от hatsumi

summ надо обнулять в теле цикла, иначе переменная накапливает все суммы, но вот в чем еще загвоздка не могу понять..

Так обнулите суммы, как при определении коэффициентов разложения, так и при расчете результатов.
Кроме этого, коэффициенты разложения Вы помещаете в массивы размерностью от 1 до 6, а при расчете аппроксимирующей кривой используете коэффициенты из диапазона от 2 до 5.
И куда подевалась постоянная составляющая?

@hatsumi · 17.10.2017, 12:01 **[ТС]**

Обнулил сумму в вычислении ряда, но не обнулил суммы в вычислении коэффициентов
Теперь все работает правильно, спасибо!

Новые блоги и статьи Все статьи Все блоги /
Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .	Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .
Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .	модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	16.10.2017, 22:55
	Вам правильный код или найти ошибку? Если правильный код, то в cftool 0

@hatsumi 0 / 0 / 0 Регистрация: 20.01.2016 Сообщений: 9
	16.10.2017, 23:58 [ТС]
	Мне бы хотелось все-таки ошибку в моем коде найти Добавлено через 1 час 0 минут Одно понял точно, summ надо обнулять в теле цикла, иначе переменная накапливает все суммы, но вот в чем еще загвоздка не могу понять.. 0

@hatsumi 0 / 0 / 0 Регистрация: 20.01.2016 Сообщений: 9
	17.10.2017, 12:01 [ТС]
	Обнулил сумму в вычислении ряда, но не обнулил суммы в вычислении коэффициентов Теперь все работает правильно, спасибо! 0