Поиск оптимальной последовательности

@Юрий124 · Регистрация: 01.06.2013

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

Здравствуйте, мне нужно дописать готовую рабочую программу. Суть заключается в том, что сейчас программа задает одну какую-то произвольную случайную последовательность, производит ее обработку путем постепенной заменой 1 на -1, если происходит улучшение, и возвращается к исходному если такого не происходит, затем строит автокорреляционную функцию оптимальной последовательности, которая получилась после последовательной замены 1 на -1 и обратно, и исходной первоначальной последовательности.
Теперь нужно усовершенствовать программу так, чтобы она задавала не одну такую случайную последовательность а 100. При этом проводила все те же операции что и в первом случае, затем из всех полученных 100 оптимальных последовательностей строила по среднеариметическому закону одну автокорреляционную функцию, и еще самую оптимальную из всех 100. Вот программа

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
clear; clc;
n = 32;
 
A = ones(1,n).*(-1).^randi([1 10],1,n);
 
Aopt = A;
 
for i=1:n
    Aopt_new = Aopt;
    Aopt_new(i) = Aopt(i) * -1;
    
    B1 = xcorr(Aopt); 
    B2 = xcorr(Aopt_new);
    
    maxB1 = max(abs([ B1(1:31) B1(33:end) ]));
    maxB2 = max(abs([ B2(1:31) B2(33:end) ]));
    
    if maxB2<maxB1
        fprintf('Поменен знак на %d позиции\n', i);
        Aopt = Aopt_new;
    end
    
end
 
 
 
figure;
 
subplot(221); plot(xcorr(A)); title('Кросс-корреляция исходного вектора A');
subplot(222); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
A
Aopt

@sergsh · 24.09.2013, 14:58

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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
clear; clc;
%%
m=10;
nn=32;
Aopt100=zeros(nn,m);
Aopt_100=zeros(nn,m);
%
for j=1:m
n = nn;
A = ones(1,n).*(-1).^randi([1 10],1,n);
Aopt = A;
for i=1:n
    Aopt_new = Aopt;
    Aopt_new(i) = Aopt(i) * -1;
    B1 = xcorr(Aopt); 
    B2 = xcorr(Aopt_new);
    maxB1 = max(abs([ B1(1:31) B1(33:end) ]));
    maxB2 = max(abs([ B2(1:31) B2(33:end) ]));
    if maxB2<maxB1
        fprintf('Поменен знак на %d позиции\n', i);
        Aopt = Aopt_new;
    end
end
%%
Aopt100(:,j)=Aopt;
figure(j);
subplot(221); plot(xcorr(A)); title('Кросс-корреляция исходного вектора A');
subplot(222); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
% A
% Aopt
end
%%
 Amean=mean(Aopt100');
 %%
 figure(m+1), plot(xcorr(Amean) );
 %%
 for j=1:m
n = nn;
A = Aopt100(:,j);
Aopt = A;
for i=1:n
    Aopt_new = Aopt;
    Aopt_new(i) = Aopt(i) * -1;
    B1 = xcorr(Aopt); 
    B2 = xcorr(Aopt_new);
    maxB1 = max(abs([ B1(1:31) B1(33:end) ]));
    maxB2 = max(abs([ B2(1:31) B2(33:end) ]));
    if maxB2<maxB1
        fprintf('Поменен знак на %d позиции\n', i);
        Aopt = Aopt_new;
    end
end
%%
Aopt_100(:,j)=Aopt;
% figure(j+20);
% subplot(221); plot(xcorr(A)); title('Кросс-корреляция исходного вектора A');
% subplot(222); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
% A
% Aopt
end
%%
 A_mean=mean(Aopt_100');
 %%
 figure(m+2), plot(xcorr(A_mean) );

Примерно так, но возможно что я вас не очень понял - что нужно ...
Проверьте

параметр m - у меня стоит 10, если нужно 100 поставьте

Matlab M
1
m=100;

@Юрий124 · 24.09.2013, 21:23 **[ТС]**

Спасибо, но немного не то......когда программа считает среднее арифметическое всех боковых лепестков она не должна складывать значение бокового лепестка с центральным, значение центрального лепестка должно остаться постоянным: в данном случае здесь 32, в зависимости от того как мы задаем nn. Боковые лепестки берутся не все, а только у каждого случая первый, т.е следующее число в последовательности за 32. как считалось в моей программе. А оптимальная автокорреляционная функция выбирает самый наилучший из 100 оптимальных автокорреляционных функций, т.е тот в котором боковой лепесток после оптимизации получился самый маленький по модулю. И в итоге в конце я должен получить 2 этих графика: один который построен путем среднеарифметического преобразования, а второй путем выбора из 100 оптимальных автокорреляционных функций самого наилучшего.

@sergsh · 25.09.2013, 12:35

Сообщение от Юрий124

Спасибо, но немного не то......когда программа считает среднее арифметическое всех боковых лепестков она не должна складывать значение бокового лепестка с центральным, значение центрального лепестка должно остаться постоянным: в данном случае здесь 32, в зависимости от того как мы задаем nn. Боковые лепестки берутся не все, а только у каждого случая первый, т.е следующее число в последовательности за 32. как считалось в моей программе. А оптимальная автокорреляционная функция выбирает самый наилучший из 100 оптимальных автокорреляционных функций, т.е тот в котором боковой лепесток после оптимизации получился самый маленький по модулю. И в итоге в конце я должен получить 2 этих графика: один который построен путем среднеарифметического преобразования, а второй путем выбора из 100 оптимальных автокорреляционных функций самого наилучшего.

К сожалению мне ничего не понятно в ваших пояснениях.
Если возможно напишите проще
Или может быть вам поможет кто-то еще ...

@Юрий124 · 25.09.2013, 15:03 **[ТС]**

сейчас моя программа получает рандомно определенную последовательность из 1 и -1, и строит автокорреляционную функцию её, затем происходит последовательное изменение знака в этой последовательности и программа уже строит автокорреляционную функцию оптимальную где боковой лепесток стал меньше. А теперь мне нужно чтобы эта программа выдавала 100 таких последовательностей, производила там ту же самую операцию оптимизации и в конце чтобы выдавала 2 графика: первый, где берутся боковые лепестки складываются и делятся на 100 и по результатам этих значений строился график; и второй - из всех 100 автокорреляционных функций оптимальных уже, строился график с самым маленьким боковым лепестком.

Боковой лепесток это - после оператора xcorr выдается последовательность, допустим 2 3 9 -6 32 -6 9 3 2 вот мне нужно уменьшить по модулю только число -6. Только их программа потом складывает и делит на 100. Другие числа не трогаются.

Новые блоги и статьи Все статьи Все блоги /
Контроль заполнения и очистка дат в зависимости от значения перечислений 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. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .	Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ

@Юрий124 0 / 0 / 0 Регистрация: 01.06.2013 Сообщений: 13
	24.09.2013, 21:23 [ТС]
	Спасибо, но немного не то......когда программа считает среднее арифметическое всех боковых лепестков она не должна складывать значение бокового лепестка с центральным, значение центрального лепестка должно остаться постоянным: в данном случае здесь 32, в зависимости от того как мы задаем nn. Боковые лепестки берутся не все, а только у каждого случая первый, т.е следующее число в последовательности за 32. как считалось в моей программе. А оптимальная автокорреляционная функция выбирает самый наилучший из 100 оптимальных автокорреляционных функций, т.е тот в котором боковой лепесток после оптимизации получился самый маленький по модулю. И в итоге в конце я должен получить 2 этих графика: один который построен путем среднеарифметического преобразования, а второй путем выбора из 100 оптимальных автокорреляционных функций самого наилучшего. 0

@Юрий124 0 / 0 / 0 Регистрация: 01.06.2013 Сообщений: 13
	25.09.2013, 15:03 [ТС]
	сейчас моя программа получает рандомно определенную последовательность из 1 и -1, и строит автокорреляционную функцию её, затем происходит последовательное изменение знака в этой последовательности и программа уже строит автокорреляционную функцию оптимальную где боковой лепесток стал меньше. А теперь мне нужно чтобы эта программа выдавала 100 таких последовательностей, производила там ту же самую операцию оптимизации и в конце чтобы выдавала 2 графика: первый, где берутся боковые лепестки складываются и делятся на 100 и по результатам этих значений строился график; и второй - из всех 100 автокорреляционных функций оптимальных уже, строился график с самым маленьким боковым лепестком. Боковой лепесток это - после оператора xcorr выдается последовательность, допустим 2 3 9 -6 32 -6 9 3 2 вот мне нужно уменьшить по модулю только число -6. Только их программа потом складывает и делит на 100. Другие числа не трогаются. 0