Решение

@Tovarisch · 14.01.2015, 22:59 **[ТС]**

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

большое спасибо, буду пробовать. надеюсь, Вы не возненавидите меня, если у меня снова возникнут к Вам вопросы

Зосима · 14.01.2015, 23:46

Не по теме:

ни в коем случае :jokingly:

@Tovarisch · 16.01.2015, 14:39 **[ТС]**

здравствуйте)
1. ввела я в help flag, и выбило мне colormap, что-то я не понимаю, куда ap_fl вставлять(
2. и еще подскажите, пожалуйста, где в Вашей программе введена

Matlab M
1
fi=1.05;%радианы

и как можно поменять на 5, 10 градусов?

Не по теме:

простите, что донимаю Вас, просто больше не к кому обратиться(

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

1. не, я имею ввиду самим ввести новую переменную, в зависимости от которой будем считать W одним или иным способом

см. строки: 36-42-47

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
65
66
67
68
function main
clc
global s ap_fl
 
ap_fl = 1; % установить 1 если сигнал апериодический
 
s = 1; % амплитуды сигнала
Fi0 = [0; pi; pi]; % нач. приближение фаз сигнала
Sigm0 = Sigma(Fi0) % нач. значение сигма
 
% тут должен быть цикл с итерациями
% для нахождения максимума
% но я пока использую встроенную
% функцию матлаба ;)
fff = @(x) 1./Sigma(x); % обратная функция
Fi = fminsearch( fff, Fi0) % фазы, при которых обратная ф-ция достигает минимума
Sm = Sigma(Fi) % макс. значение сигма
 
end
 
 
%*** Вспомогательные функции ***
%*******************************
function u = calcS(F) % функция расчета комплексного сигнала
global s
F = F(:); % делаем столбец
s = s(:); % делаем столбец
u = s.*exp(1i*F); % умножаем амплитуды на экспонету
end
 
function w = calcW(F) % функция расчета фильтра
global Rn ap_fl
n = length(F); 
S = calcS(F); % вычисляем сигнал
R = 0;
if ap_fl==1 % если сигнал апериодичен
    B = [zeros(n-1,1); S; zeros(n-1,1)]; % дополняем массив нулями
    for i = 1:2*n-1 % сдвигаем окно
        ss = B(i:i+n-1); % берем кусочек массива из окна
        R = R + ss*ss';
    end
else  % если сигнал периодичен
    for k=1:length(F)
        ss = circshift(S,k-1);
        R = R + ss*ss';
    end
end
N0 = 5; % шум
Rn = R + N0*eye(size(R));
w = Rn\S;
end
 
function Z = Sigma(F) % функция расчета сигмы
global Rn
S = calcS(F);
W = calcW(F);
Z = abs(W'*S)^2/(W'*Rn*W);
end
 
function g = gradS(F) % функция расчета градиента
h = 1e-2; % приращение
g = zeros(size(F)); % заготовка вектора градиента
for i = 1:length(F) % для каждой компоненты фазы
    df = zeros(size(F)); % заготовка массива приращения из нулей
    df(i) = h; % приращение i-й компоненты
    g(i) = (Sigma(F+df)-Sigma(F-df))/(2*h); % классическая формула производной
end
end

2. я посмотрел в твоей первой программе, что такое fi и понял, что это приращение для расчета градиента

в моем случае это параметр h (строка 61) и его нужно брать чем меньше тем лучше

чтобы перевести градусы в радианы, нужно их умножить на pi и разделить на 180

@Tovarisch · 16.01.2015, 15:39 **[ТС]**

Благодарю Вас)) еще раз извините, за то, что отнимаю Ваше время

@Nick07 · 16.01.2015, 16:31

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

чтобы перевести градусы в радианы, нужно их умножить на pi и разделить на 180

а можно попросить МАТЛАБ это сделать:

degtorad
Convert angles from degrees to radians
http://www.mathworks.com/help/... torad.html

radtodeg
Convert angles from radians to degrees
http://www.mathworks.com/help/... todeg.html

@Tovarisch · 16.01.2015, 16:43 **[ТС]**

спасибо, буду знать!)

@Nick07 · 16.01.2015, 17:17

Как предлагает уважаемый ЗОСИМА, видимо, методически правильнее, ибо нагляднее, но мне проще через операторы МАТЛАБ (не надо держать в голове формулу пересчета).

@Tovarisch · 17.01.2015, 16:50 **[ТС]**

1. а зачем это делать через Матлаб, если эту переменную я ввожу 1 раз на протяжении всей программы?

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

h = 1e-2;

2. очень странно, что для м-последовательности для периодического и апериодического случаев, N=3, 7,15 сигма возрастает, а для кода Баркера, N=5, 7, 13- не возрастает

Добавлено через 20 часов 47 минут
Здравствуйте. Подскажите, пожалуйста, почему для кода Баркера, N=5, 7, 13 - сигма не возрастает?????

Добавлено через 1 час 9 минут
Даже если убирать 2 крайних боковых лепестка сигма все равно не растет(

Зосима · 17.01.2015, 19:30

Солнце, нет нужды печалиться, если для некоторых случаев сигма не увеличивается!

Смотри, как ищется максимум: берем начальную точку, смотрим на знак производной/градиента (т.е. узнаем уменьшается функция или возрастает) соответственно берем следущую точку и так двигаемся, пока не нацдем максимум... НО возможен вариант, когда ты везунчик и начальная точка оказывается в точке максимума!

тогда как бы ты не меняла точку, функция будет убывать

думается мне, что это аккурат твой случай

Кстати, если не секрет, какой физический смысл этой сигмы?

@Tovarisch · 17.01.2015, 21:27 **[ТС]**

Ох, благодарю Вас, надеюсь это мой случай

Конечно это не секрет, сигма это сигнал/помеха+шум, по крайней мере, так говорил преподаватель...
Подскажите пожалуйста:
1. Почему меняя значение в строке

Matlab M
1
h = 1e-2;

значение сигмы не меняется?

2. Не могу понять где записано в Вашей программе значение лямбда, и где происходит деление на 2, при не возрастании сигмы?

Зосима · 18.01.2015, 13:13

Скорее всего так и есть, ибо код Баркера оптимальный код с точки щрения мощности боковых лепестков АКФ

Приращение h у меня используется в функции расчета градиента gradS, которую я пока не использовал и до цикла и лямбда у меня пока руки не дошли.

до завтра потерпит?

Вопрос: какие условия выхода из цикла? (обычно задают некторое малое число-точность и продолжают вычисления, пока разность функций (или норма градиента) не станет меньше этого значения.

Добавлено через 1 час 29 минут
Набросал цикл, но не проверял: может не работать или зависать

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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
function main
clc
global s ap_fl
 
ap_fl = 1; % установить 1 если сигнал апериодический
 
s = 1; % амплитуды сигнала
Fi0 = [0; pi; pi]; % нач. знвчение фаз сигнала
Sigm0 = Sigma(Fi0) % нач. значение сигма
 
eps = 1e-6; % точность
Sigm1 = 10*eps+Sigm0; % начальное значение сигма1, чтобы вошло в цикл
 
while abs(Sigm1-Sigm0)>eps*
** *Sigm0 = Sigma(Fi0);*
** *lambda = Sigm0/sum(gradS(Fi0).^2); % формула 17
** *Fi1 = Fi0 + lamda*gradS(Fi0);*
** *Sigm1 = Sigma(Fi1);
** *while Sigm1 < Sigm0*
** * * *lambda = lambda/2;
** * * *Fi1 = Fi0 + lamda*gradS(Fi0);*
** * * Sigm1 = Sigma(Fi1);
** *end
** *Fi0 = Fi1;
end
Fi1 % оптимальные фазы
Sigma(Fi1) % макс. значение сигма
 
end
 
 
%*** Вспомогательные функции ***
%*******************************
function u = calcS(F) % функция расчета комплексного сигнала
global s
F = F(:); % делаем столбец
s = s(:); % делаем столбец
u = s.*exp(1i*F); % умножаем амплитуды на экспонету
end
 
function w = calcW(F) % функция расчета фильтра
global Rn ap_fl
n = length(F);*
S = calcS(F); % вычисляем сигнал
R = 0;
if ap_fl==1 % если сигнал апериодичен
** *B = [zeros(n-1,1); S; zeros(n-1,1)]; % дополняем массив нулями
** *for i = 1:2*n-1 % сдвигаем окно
** * * *ss = B(i:i+n-1); % берем кусочек массива из окна
** * * *R = R + ss*ss';
** *end
else *% если сигнал периодичен
** *for k=1:length(F)
** * * *ss = circshift(S,k-1);
** * * *R = R + ss*ss';
** *end
end
N0 = 5; % шум
Rn = R + N0*eye(size(R));
w = Rn\S;
end
 
function Z = Sigma(F) % функция расчета сигмы
global Rn
S = calcS(F);
W = calcW(F);
Z = abs(W'*S)^2/(W'*Rn*W);
end
 
function g = gradS(F) % функция расчета градиента
h = 0.1*pi/180; % приращение 0.1 градус
g = zeros(size(F)); % заготовка вектора градиента
for i = 2:length(F) % для каждой компоненты фазы, кроме первой
** *df = zeros(size(F)); % заготовка массива приращения из нулей
** *df(i) = h; % приращение i-й компоненты
** *g(i) = (Sigma(F+df)-Sigma(F-df))/(2*h); % классическая формула производной
end
end

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

Не по теме:

*вот Ъ! вместо пробелов звездочки повтыкались! :(

@Tovarisch · 18.01.2015, 18:06 **[ТС]**

Здравствуйте. Благодарю Вас за помощь, потерпит

На счет условия выхода из цикла, мне известно, только то, что я ранее писала в условии, как только встречусь с преподавателем, уточню.
В программе, к сожалению, выбивается ошибка:
Warning: Matrix is singular, close to singular or badly scaled. Results may be
inaccurate. RCOND = NaN.
ans =

NaN + NaNi
не могу понять почему так происходит??

Не по теме:

Простите меня за наглость:gsad:

Зосима · 19.01.2015, 16:40

а, вот еще хитрый момент заметил: когда идет формула для сигма оно обозначается как (формулы 14, 15):

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma \left(W_0,S_0 \right) = \Sigma_0(S_0) = \frac{\left|W_0'S_0 \right|^2}{W_0'RW_0};$

а дальше идет !внимание!

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma_0 \left(\Phi_0 \right)$

т.е. я правильно понимаю, что тут вместо S₀ подставляется фаза Ф₀, т.е. получаем формулу:

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma_0(\Phi _0) = \frac{\left|W_0'\Phi_0 \right|^2}{W_0'RW_0};$

Что скажешь?

Зосима · 19.01.2015, 18:12

Вот такую штуку набросал, где используется $https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma _0 \left(\Phi _0 \right)$

что-то максимизирует, но я сам слабо понимаю что к чему

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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
function main
clc
global s ap_fl
 
ap_fl = 1; % установить 1 если сигнал апериодический
 
S = [+1 +1 +1 -1 +1]; % код Баркера 5
 
s = abs(S); % амплитуды сигнала
fi = angle(S); % фазы сигналов
Fi0 = [0;diff(fi(:))]; %fi(:); % разности фаз
Fi1 = Fi0;
Sigm0 = Sigma(calcS(Fi0)) % нач. значение сигма
 
eps = 1e-6; % точность
Sigm1 = 10*eps+Sigm0; % начальное значение сигма1, чтобы вошло в цикл
 
while abs(Sigm1-Sigm0)>eps & norm(gradS(Fi0))>eps
    Sigm0 = Sigma(calcS(Fi0));
    Gr = gradS(calcS(Fi0));
    lambda = Sigma(Fi0)/norm(Gr); % формула 17
    Fi1 = Fi0 + lambda*gradS(calcS(Fi0));
    Sigm1 = Sigma(calcS(Fi1))
    while Sigm1 < Sigm0 & norm(gradS(Fi0))>eps
        lambda = lambda/2;
        Fi1 = Fi0 + lambda*gradS(calcS(Fi0));
        Sigm1 = Sigma(calcS(Fi1));
    end
    Fi0 = Fi1;
end
 
disp('Оптимальные фазы :')
Fi1 % оптимальные фазы
disp('Макс. отн. С/Ш :')
Sigma(calcS(Fi1)) % макс. значение сигма
end
 
 
%*** Вспомогательные функции ***
%*******************************
function u = calcS(F) % функция расчета комплексного сигнала
global s
F = F(:); % делаем столбец
s = s(:); % делаем столбец
u = s.*exp(1i*F); % умножаем амплитуды на экспонету
end
 
function w = calcW(F) % функция расчета фильтра
global Rn ap_fl
n = length(F);
S = calcS(F); % вычисляем сигнал
R = 0;
if ap_fl==1 % если сигнал апериодичен
    B = [zeros(n-1,1); S; zeros(n-1,1)]; % дополняем массив нулями
    for i = 1:2*n-1 % сдвигаем окно
        ss = B(i:i+n-1); % берем кусочек массива из окна
        R = R + ss*ss';
    end
else % если сигнал периодичен
    for k=1:length(F)
        ss = circshift(S,k-1);
        R = R + ss*ss';
    end
end
N0 = 5; % шум
Rn = R + N0*eye(size(R));
w = Rn\S;
end
 
function Z = Sigma(X) % функция расчета сигмы
global Rn
W = calcW(calcS(X));
Z = abs(W'*X)^2/(W'*Rn*W);
end
 
function g = gradS(F) % функция расчета градиента
h = 0.0001*pi/180; % приращение 0.1 градус
g = zeros(size(F)); % заготовка вектора градиента
for i = 2:length(F) % для каждой компоненты фазы, кроме первой
    df = zeros(size(F)); % заготовка массива приращения из нулей
    df(i) = h; % приращение i-й компоненты
    g(i) = (Sigma(F+df)-Sigma(F-df))/(2*h); % классическая формула производной
end
end

Еще хочу у тебя спросить за фильтр: когда мы считаем $https://www.cyberforum.ru/cgi-bin/latex.cgi?\Sigma \left(\Phi \right)$ , то каким берется W ? от Ф или от S(Ф) ?

@Tovarisch · 19.01.2015, 19:27 **[ТС]**

Здравствуйте

1. Я так понимаю, что сигма(Ф), она так отмечена, потому что рассчитывается каждый раз, после добавления фазы = 60 градусов, и затем добавления градиента. А сигма(S) обозначается так, потому что рассчитывается для сигнала без добавления фаз.

Не по теме:

это мои глупые предположения :confusion:

2. Изначально я брала фильтр от первого сигнала, 1 -1 -1, и вставляла каждый раз для подсчета новой сигмы, но преподавателю не понравился такой вариант, и он решил, что нужно каждый раз высчитывать фильтр для нового сигнала, т.е. к тому, к которому прибавляла 60 градусов и градиенты, и считать сигму, каждый раз с новым фильтром.

Не по теме:

надеюсь, для Вас хоть что-то прояснилось:gsmile:

3. Я использовала сигналы м-последовательности, например N=7 :-1 1 -1 1 1 1 -1, и для N=15.А также коды Баркера для N=5, 13.

Не по теме:

К сожалению, не могу добиться по Вашей программе, тех результатов, которые получила по своей допотопной программе. Преподаватель говорил, что при периодическом сигнале 1 -1 -1 максимальное значение, которого может достичь сигма = 0.3750, что я и получила подставляя данные каждый раз вручную, в свою программу.:gsad:

Зосима · 20.01.2015, 10:43

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

каждый раз, после добавления фазы = 60 градусов

где это ты такое нашла про 660 градусов?

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

он решил, что нужно каждый раз высчитывать фильтр для нового сигнала

в моих вариантах так и происходит

меняем фазу - меняется сигнал и фильтр

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

Преподаватель говорил, что при периодическом сигнале 1 -1 -1 максимальное значение, которого может достичь сигма = 0.3750

Хорошо, будем опираться на этот результат

вот рабочая программа:

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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
function main
clc
global s ap_fl
 
ap_fl = 1; % установить 1 если сигнал апериодический
 
S = [1 -1 -1]; 
 
s = abs(S); % амплитуды сигнала
fi = angle(S); % фазы сигналов
%Fi0 = [0;diff(fi(:))];  % разности фаз
Fi0 = fi(:);
Fi1 = Fi0;
Sigm0 = Sigma(Fi0); % нач. значение сигма
 
eps = 1e-6; % точность сигма
Geps = 1e-14; % точность градиента
Sigm1 = 10*eps+Sigm0; % начальное значение сигма1, чтобы вошло в цикл
 
k = 1; % множитель lambda
while abs(Sigm1-Sigm0)>eps & norm(gradS(Fi0))>Geps
    Sigm0 = Sigma(Fi0)
    lambda = Sigma(Fi0)/norm(gradS(Fi0)); % формула 17
    Fi1 = Fi0 + k*lambda*gradS(Fi0);
    Sigm1 = Sigma(Fi1);
    while Sigm1 < Sigm0 & norm(gradS(Fi0))>Geps
        k = k/2;
        Fi1 = Fi0 + k*lambda*gradS(Fi0);
        Sigm1 = Sigma(Fi1);
    end
    Fi0 = Fi1;
end
 
disp('Оптимальные фазы :')
Fi1 % оптимальные фазы
disp('Макс. отн. С/Ш :')
Sigma(Fi1) % макс. значение сигма
end
 
 
%*** Вспомогательные функции ***
%*******************************
function u = calcS(F) % функция расчета комплексного сигнала
global s
F = F(:); % делаем столбец
s = s(:); % делаем столбец
u = s.*exp(1i*F); % умножаем амплитуды на экспонету
end
 
function w = calcW(F) % функция расчета фильтра
global Rn ap_fl
n = length(F);
S = calcS(F); % вычисляем сигнал
R = 0;
if ap_fl==1 % если сигнал апериодичен
    B = [zeros(n-1,1); S; zeros(n-1,1)]; % дополняем массив нулями
    for i = 1:2*n-1 % сдвигаем окно
        ss = B(i:i+n-1); % берем кусочек массива из окна
        R = R + ss*ss';
    end
else % если сигнал периодичен
    for k=1:length(F)
        ss = circshift(S,k-1);
        R = R + ss*ss';
    end
end
N0 = 5; % шум
Rn = R + N0*eye(size(R));
w = Rn\S;
end
 
function Z = Sigma(Fi) % функция расчета сигмы
global Rn
S = calcS(Fi);
W = calcW(Fi);
Z = abs(W'*S)^2/(W'*Rn*W);
end
 
function g = gradS(F) % функция расчета градиента
h = 0.0001*pi/180; % приращение 0.1 градус
g = zeros(size(F)); % заготовка вектора градиента
for i = 1:length(F) % для каждой компоненты фазы, кроме первой
    df = zeros(size(F)); % заготовка массива приращения из нулей
    df(i) = h; % приращение i-й компоненты
    g(i) = (Sigma(F+df)-Sigma(F-df))/(2*h); % классическая формула производной
end
end

@Tovarisch · 20.01.2015, 12:32 **[ТС]**

Здравствуйте)
1. Преподаватель сначала хотел, чтоб было 5, 10, 15 градусов, но потом решил что нужно добавлять 60, посмотреть, что будет получаться, а потом взять снова 5 градусов. Зачем это нужно, мне не понятно.

2. Но я опять же не могу понять, почему по Вашей данной программе для периодического сигнала 1 -1 -1, сигма=0.3519 ( Может кроме 1 в 5 строчке программы еще что-то нужно менять?

Не по теме:

Огромное спасибо Вам за помощь! Извините, что донимаю Вас своими глупыми объяснениями и вопросами :gsorry:

Зосима · 20.01.2015, 12:40

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

Преподаватель сначала хотел, чтоб было 5, 10, 15 градусов

что именно "5, 10, 15 градусов"? куда их добавлять?

мне кажется, что 0.35 довольно близко к 0.37, т.е. результаты правдоподобные

Любопытно узнать как у преподавателя получилось это значение? может шум N0 был другой?

Не по теме:

мне самому интересно решать трудные задания ;)

@Tovarisch · 20.01.2015, 13:10 **[ТС]**

1. Допустим 5 градусов для сигнала 1 -1 -1, нужно было добавлять к второй фазе, при этом вычислить новую сигму, а затем к третей фазе и заново вычислить сигму, и тоже самое проделывать с 10, 15 градусами.Т.е. по моей старой программе добавлять нужно в строчки 18 и 33.

2. К сожалению, при сравнении сигм, даже увеличение сотой играет роль, так сказал преподаватель)

3. Почему-то преподаватель взял сигнал 1 120градусов 120, шум=5, посчитал сигма, которое = 0.3750, и сказал, что это максимум, чего может достигнуть сигма при периодическом сигнале1 -1 -1. Т.е. я так поняла,что он заведомо знает какие макс. значения сигмы должны быть для периодических сигналов. Поэтому его очень интересует, сигма при апериодических случаев.

Новые блоги и статьи Все статьи Все блоги /
Был там один разговор по поводу свободы в материальном мире. kumehtar 19.08.2026 Суть: рассматривается живое существо, оказавшееся внутри довольно странной системы (этого мира) и пытающееся обустроить в ней свой кусок пространства. Жизнь действительно предъявляет каждому. . .	Когда логика программы не спасает от человеческих ошибок Maks 18.08.2026 В последнее время всё чаще и чаще сталкиваюсь с таким явлением, как абсолютная невнимательность (или глупость) пользователей. Проявляется это чаще всего на работе в коллективе. Допустим, человек с. . .	Лето уходит kumehtar 17.08.2026	Мысли в слух kumehtar 17.08.2026 Забавно, насколько сейчас стала доступна информация. Например о магии, духовном развитии, медитациях, и других подобных направлениях, ранее зачастую тайных, передаваемых от учителя к ученику. Хотя. . .
Перемещение строк из ТЧ в другой документ с учетом текущего пробега Maks 17.08.2026 Реализация из решения ниже выполнена на примере нетипового документа "Автозапчасти", с ТЧ "Шины". За основу взят алгоритм отсюда: https:/ / www. cyberforum. ru/ blogs/ 359708/ 10838. html Задача: . . .	Саморегулирующийся социальный контракт для сервера cross-section. Hrethgir 14.08.2026 С кодом конечно таких глубоких размышлений пока не было, впрочем я уже привык к алгоритмизации. Суть предмета записи: снова в диалоге с нейросетью (я взял пока себе ник для учётки админа - Rector). . . .	Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	14.01.2015, 22:59 [ТС]
	большое спасибо, буду пробовать. надеюсь, Вы не возненавидите меня, если у меня снова возникнут к Вам вопросы 0

Зосима
	14.01.2015, 23:46
	Не по теме: ни в коем случае :jokingly: 0

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	16.01.2015, 15:39 [ТС]
	Благодарю Вас)) еще раз извините, за то, что отнимаю Ваше время 0

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	16.01.2015, 16:43 [ТС]
	спасибо, буду знать!) 0

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	16.01.2015, 17:17
	Как предлагает уважаемый ЗОСИМА, видимо, методически правильнее, ибо нагляднее, но мне проще через операторы МАТЛАБ (не надо держать в голове формулу пересчета). 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	17.01.2015, 19:30
	Солнце, нет нужды печалиться, если для некоторых случаев сигма не увеличивается! Смотри, как ищется максимум: берем начальную точку, смотрим на знак производной/градиента (т.е. узнаем уменьшается функция или возрастает) соответственно берем следущую точку и так двигаемся, пока не нацдем максимум... НО возможен вариант, когда ты везунчик и начальная точка оказывается в точке максимума! тогда как бы ты не меняла точку, функция будет убывать думается мне, что это аккурат твой случай Кстати, если не секрет, какой физический смысл этой сигмы? 1

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	18.01.2015, 18:06 [ТС]
	Здравствуйте. Благодарю Вас за помощь, потерпит На счет условия выхода из цикла, мне известно, только то, что я ранее писала в условии, как только встречусь с преподавателем, уточню. В программе, к сожалению, выбивается ошибка: Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN. ans = NaN + NaNi не могу понять почему так происходит?? Не по теме: Простите меня за наглость:gsad: 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	19.01.2015, 16:40
	а, вот еще хитрый момент заметил: когда идет формула для сигма оно обозначается как (формулы 14, 15): $https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma \left(W_0,S_0 \right) = \Sigma_0(S_0) = \frac{\left\|W_0'S_0 \right\|^2}{W_0'RW_0};$ а дальше идет !внимание! $https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma_0 \left(\Phi_0 \right)$ т.е. я правильно понимаю, что тут вместо S₀ подставляется фаза Ф₀, т.е. получаем формулу: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\huge \Sigma_0(\Phi _0) = \frac{\left\|W_0'\Phi_0 \right\|^2}{W_0'RW_0};$ Что скажешь? 1

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	19.01.2015, 19:27 [ТС]
	Здравствуйте 1. Я так понимаю, что сигма(Ф), она так отмечена, потому что рассчитывается каждый раз, после добавления фазы = 60 градусов, и затем добавления градиента. А сигма(S) обозначается так, потому что рассчитывается для сигнала без добавления фаз. Не по теме: это мои глупые предположения :confusion: 2. Изначально я брала фильтр от первого сигнала, 1 -1 -1, и вставляла каждый раз для подсчета новой сигмы, но преподавателю не понравился такой вариант, и он решил, что нужно каждый раз высчитывать фильтр для нового сигнала, т.е. к тому, к которому прибавляла 60 градусов и градиенты, и считать сигму, каждый раз с новым фильтром. Не по теме: надеюсь, для Вас хоть что-то прояснилось:gsmile: 3. Я использовала сигналы м-последовательности, например N=7 :-1 1 -1 1 1 1 -1, и для N=15.А также коды Баркера для N=5, 13. Не по теме: К сожалению, не могу добиться по Вашей программе, тех результатов, которые получила по своей допотопной программе. Преподаватель говорил, что при периодическом сигнале 1 -1 -1 максимальное значение, которого может достичь сигма = 0.3750, что я и получила подставляя данные каждый раз вручную, в свою программу.:gsad: 0

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	20.01.2015, 12:32 [ТС]
	Здравствуйте) 1. Преподаватель сначала хотел, чтоб было 5, 10, 15 градусов, но потом решил что нужно добавлять 60, посмотреть, что будет получаться, а потом взять снова 5 градусов. Зачем это нужно, мне не понятно. 2. Но я опять же не могу понять, почему по Вашей данной программе для периодического сигнала 1 -1 -1, сигма=0.3519 ( Может кроме 1 в 5 строчке программы еще что-то нужно менять? Не по теме: Огромное спасибо Вам за помощь! Извините, что донимаю Вас своими глупыми объяснениями и вопросами :gsorry: 0

@Tovarisch 0 / 0 / 0 Регистрация: 13.01.2015 Сообщений: 38
	20.01.2015, 13:10 [ТС]
	1. Допустим 5 градусов для сигнала 1 -1 -1, нужно было добавлять к второй фазе, при этом вычислить новую сигму, а затем к третей фазе и заново вычислить сигму, и тоже самое проделывать с 10, 15 градусами.Т.е. по моей старой программе добавлять нужно в строчки 18 и 33. 2. К сожалению, при сравнении сигм, даже увеличение сотой играет роль, так сказал преподаватель) 3. Почему-то преподаватель взял сигнал 1 120градусов 120, шум=5, посчитал сигма, которое = 0.3750, и сказал, что это максимум, чего может достигнуть сигма при периодическом сигнале1 -1 -1. Т.е. я так поняла,что он заведомо знает какие макс. значения сигмы должны быть для периодических сигналов. Поэтому его очень интересует, сигма при апериодических случаев. 0