Дописать программу

@Yurok007 · Регистрация: 27.04.2010

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

Здравствуйте, мне нужно исправить программу

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
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

так, чтобы случайная последовательность 1 и -1 задавалась случайным образом ( как сейчас у меня ), но так чтобы она один раз ее задала и сохранила, чтобы потом при закрытии программы и открытии мне выдавалась та же самая случайная последовательность, что и последний раз! Затем сделать так чтобы таких последовательностей было уже 100 и точно так же сохранила их.

Зосима · 23.10.2013, 15:10

lomt, ой! точно! бошка дырявая!

@Yurok007 · 29.10.2013, 11:57 **[ТС]**

а как теперь сделать, чтобы вместо макс. бок. лепестка и его номера определялся минимальный боковой лепесток и его номер и также выводился график. И еще как сделать чтобы теперь на экран 100 графиков не выводились, а выводился только один для минимального бок. лепестка????

Зосима · 29.10.2013, 14:01

Вот, причесал

что не нужно - закоментировал.

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
clear; clc; close all
 
m=100;
nn=32;
Aopt100=zeros(nn,m);
Aopt_100=zeros(nn,m);
AB = load('Adata.mat');
A = AB.AA;
for j=1:m
    n = nn;
    Aopt = A(j,:);
    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
    Q(j) = maxB1;
    %
    Aopt100(:,j)=Aopt;
    %figure
    %subplot(121); plot(xcorr(A(j,:))); title('Кросс-корреляция исходного вектора A');
    %subplot(122); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
    % A(:,j)
    % Aopt
end
 
[q1 id1] = max(Q); % поределяем макс. бок лепесток и его номер
[q2 id2] = min(Q); % поределяем макс. бок лепесток и его номер
figure
subplot(121)
plot(abs(xcorr(Aopt100(:,id1)))), hold on
plot([0 2*size(Aopt100,1)],[q1 q1],':r')
title(['Макс. боковой лепесток, r_{max}=',num2str(q1)])
xlim([0 2*size(Aopt100,1)])
 
subplot(122)
plot(abs(xcorr(Aopt100(:,id2)))), hold on
plot([0 2*size(Aopt100,1)],[q2 q2],':r')
title(['Мин. боковой лепесток r_{min}=',num2str(q2)])
xlim([0 2*size(Aopt100,1)])
 
figure
plot(1:length(Q),Q,'-ob')
xlabel('Номер последовательности')
ylabel('Макс. боковой лепесток')
hold on
plot([1 length(Q)],[mean(Q) mean(Q)],':r')
title(['Средний уровень бок.лепестка, r_{cp}=',num2str(mean(Q))])

@Yurok007 · 10.11.2013, 14:30 **[ТС]**

Сообщение от Зосима

Вот, причесал

что не нужно - закоментировал.

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
clear; clc; close all
 
m=100;
nn=32;
Aopt100=zeros(nn,m);
Aopt_100=zeros(nn,m);
AB = load('Adata.mat');
A = AB.AA;
for j=1:m
    n = nn;
    Aopt = A(j,:);
    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
    Q(j) = maxB1;
    %
    Aopt100(:,j)=Aopt;
    %figure
    %subplot(121); plot(xcorr(A(j,:))); title('Кросс-корреляция исходного вектора A');
    %subplot(122); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
    % A(:,j)
    % Aopt
end
 
[q1 id1] = max(Q); % поределяем макс. бок лепесток и его номер
[q2 id2] = min(Q); % поределяем макс. бок лепесток и его номер
figure
subplot(121)
plot(abs(xcorr(Aopt100(:,id1)))), hold on
plot([0 2*size(Aopt100,1)],[q1 q1],':r')
title(['Макс. боковой лепесток, r_{max}=',num2str(q1)])
xlim([0 2*size(Aopt100,1)])
 
subplot(122)
plot(abs(xcorr(Aopt100(:,id2)))), hold on
plot([0 2*size(Aopt100,1)],[q2 q2],':r')
title(['Мин. боковой лепесток r_{min}=',num2str(q2)])
xlim([0 2*size(Aopt100,1)])
 
figure
plot(1:length(Q),Q,'-ob')
xlabel('Номер последовательности')
ylabel('Макс. боковой лепесток')
hold on
plot([1 length(Q)],[mean(Q) mean(Q)],':r')
title(['Средний уровень бок.лепестка, r_{cp}=',num2str(mean(Q))])

Только что заметил прежнюю ошибку, как только я раскомментирую %figure
%subplot(121); plot(xcorr(A(j,

)); title('Кросс-корреляция исходного вектора A');
%subplot(122); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
% A(:,j)
% Aopt
эти строчки снова выдается ошибка ??? Index exceeds matrix dimensions.

Error in ==> AKF5 at 31
A(:,j)

Добавлено через 8 минут

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

Только что заметил прежнюю ошибку, как только я раскомментирую %figure
%subplot(121); plot(xcorr(A(j,

)); title('Кросс-корреляция исходного вектора A');
%subplot(122); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
% A(:,j)
% Aopt
эти строчки снова выдается ошибка ??? Index exceeds matrix dimensions.

Error in ==> AKF5 at 31
A(:,j)

заметил такую вещь, что график, который выводится для минимального бокового лепестка ( самый последний) отличается от графика который находится по адресу [q2 id2]. Хоть и значения минимального бокового лепестка одинаковые, но графики совершенно разные, можно это как то исправить?

@Yurok007 · 14.11.2013, 15:52 **[ТС]**

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

Только что заметил прежнюю ошибку, как только я раскомментирую %figure
%subplot(121); plot(xcorr(A(j,

)); title('Кросс-корреляция исходного вектора A');
%subplot(122); plot(xcorr(Aopt)); title('Кросс-корреляция оптим. вектора A');
% A(:,j)
% Aopt
эти строчки снова выдается ошибка ??? Index exceeds matrix dimensions.

Error in ==> AKF5 at 31
A(:,j)

Добавлено через 8 минут

заметил такую вещь, что график, который выводится для минимального бокового лепестка ( самый последний) отличается от графика который находится по адресу [q2 id2]. Хоть и значения минимального бокового лепестка одинаковые, но графики совершенно разные, можно это как то исправить?

похоже исправить никак нельзя((((((((((((

Новые блоги и статьи Все статьи Все блоги /
Запустил конкурс "тем и промптов для текстовых квестов созданных почти чисто ИИ" Adler 06.10.2026 Всем привет! За последние три-четыре дня я создал более 16 текстовых квестовых игр используя преимущественно по одному запросу к ИИ на игру. Мне так понравилось смотреть все ветки/ сцены во всех. . .	ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .	Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Nekobox - outbounds[0].transport: unknown transport type: raw damix 01.10.2026 Фикс ошибки Правым кликом по серверу -> отладочная информация -> edit Заменить "net": "raw", на "net": "tcp", Нажать кнопку reload.
Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .	Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .	Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .	Запрет дублирования строк в табличной части Maks 13.09.2026 Реализация из решения ниже выполнена на нетиповом справочнике "Нормы ТО" с табличной часть "Виды ТО", разработанного в КА2, со следующими реквизитами: - ВидТО (СправочникСсылка. ВидыТО); - ВидГСМ. . .

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	23.10.2013, 15:10
	lomt, ой! точно! бошка дырявая! 0

@Yurok007 0 / 0 / 0 Регистрация: 27.04.2010 Сообщений: 16
	29.10.2013, 11:57 [ТС]
	а как теперь сделать, чтобы вместо макс. бок. лепестка и его номера определялся минимальный боковой лепесток и его номер и также выводился график. И еще как сделать чтобы теперь на экран 100 графиков не выводились, а выводился только один для минимального бок. лепестка???? 0