Ошибка при вычислении

@sanekkkk · Регистрация: 16.12.2012

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

Код программы

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 all
close all
clc
 
Ea = 0.15;
Eb = 0.15;
Ee = 0.1;
 
 
N = 70;
Llist = 2.^(1:7);
 
for k = 1:length(Llist)
    L = Llist(k);
    
    E = (1-Ea)*Eb + Ea*(1-Eb)
    betta(k) = E^L/((E^L)+(1-E)^L)
    
    h = @(x) -x.*log2(x) - (1-x).*log2(1-x);
    %Ib = 1-h(betta(k));
    %Ib = real(Ib1)
     Ib = 1 + betta(k)*log2(betta(k)) + ( 1 - betta(k))*log2(1-betta(k));
      
    Dw = @(w) ((1-Ea)*(1-Eb)*(1-Ee)+Ea*Eb*Ee).^(L-w).*(((1-Ea)*(1-Eb)*Ee+Ea*Eb*(1-Ee)).^w) +...
        (((1-Ea)*Eb*(1-Ee)+Ea*(1-Eb)*Ee)).^(L-w).*((1-Ea)*Eb*Ee+Ea*(1-Eb)*(1-Ee)).^w;
    
    S = 0;
    for w = 0:L
       % u = Dw(w)/(Dw(w)+Dw(N-w));
        S =   nchoosek(L,w)*Dw(w)*(1-h(Dw(w)/(Dw(w)+Dw(N-w))));
    end
    Ie = 1/(E^L+(1-E)^L) * S;
 
    Rc = 1/L;
    K = log2(L);
    e = zeros(1,K); %заготовка из нулей
    e(1) = E; % начальное значение
    for i = 0:K-1
        Rc = Rc * ( e(i+1)^(2*i) + (1-e(i+1))^(2*i) );
        e(i+2) = e(i+1)^2 + (1-e(i+1))^2 ; 
    end
    Reff(k) = Rc*(Ib-Ie);
    Idiff(k)=(Ib-Ie);
end
Reff(k)
 
 
semilogy(Llist, Idiff )
title ( 'График зависимости Reff от длины сообщения')
xlabel('N');
ylabel('Reff'); 
set(0,'DefaultAxesFontSize',14,'DefaultAxesFontName','Times New Roman');
grid on
 
% figure
% plot(Llist,Ie,'blue')
% %hold on
%plot(Llist,Ib,'red')
%grid on
%figure
%plot(Llist,betta,'blue')
%hold on
% plot(Llist,gamma,'red')
% grid on

Выдается ошибка типа

Warning: Result may not be exact. Coefficient is greater than 9.007199e+15 and is only accurate to 15 digits
> In nchoosek at 92
In iter_orig at 30

Как можно ее исправить?

@VSI · 28.11.2017, 21:28

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

Выдается ошибка типа

Вообще-то это не ошибка, а предупреждение...

Предупреждение. Результат может быть неточным. Коэффициент больше 9.007199e + 15 и имеет точность до 15 цифр

@sanekkkk · 29.11.2017, 14:10 **[ТС]**

ТО есть в данном случае для nchossek не побороть это предупреждение?
Если хочу большие числа?

Добавлено через 16 часов 4 минуты
А как еще вывести на экран например что при L=2 Idiff=
и т.д до L=2^5

@Centurio · 29.11.2017, 17:00

У вас N=70 или всё-таки N=7 в программе?

@sanekkkk · 29.11.2017, 17:22 **[ТС]**

Давайте для ясности N=40 а Llist = 2.^(1:5);

@Centurio · 29.11.2017, 18:25

sanekkkk, у вас код

Matlab M
1
2
3
4
5
S = 0;
    for w = 0:L
       % u = Dw(w)/(Dw(w)+Dw(N-w));
        S =   nchoosek(L,w)*Dw(w)*(1-h(Dw(w)/(Dw(w)+Dw(N-w))));
    end

должен накапливать сумму, а на деле выдаёт на выходе всего лишь значение S при последнем значении w.

@sanekkkk · 29.11.2017, 20:05 **[ТС]**

Код программы

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
clear all
close all
clc
 
Ea = 0.15;
Eb = 0.15;
Ee = 0.1;
  
N = 40;
Llist = 2.^(1:6);
 
for k = 1:length(Llist)
    L = Llist(k);
    
    E = (1-Ea)*Eb + Ea*(1-Eb)
    betta(k) = E^L/((E^L)+(1-E)^L)
    
    h = @(x) -x.*log2(x) - (1-x).*log2(1-x);
    %Ib = 1-h(betta(k));
    %Ib = real(Ib1)
     Ib = 1 + betta(k)*log2(betta(k)) + ( 1 - betta(k))*log2(1-betta(k));
      
    Dw = @(w) ((1-Ea)*(1-Eb)*(1-Ee)+Ea*Eb*Ee).^(L-w).*(((1-Ea)*(1-Eb)*Ee+Ea*Eb*(1-Ee)).^w) +...
        (((1-Ea)*Eb*(1-Ee)+Ea*(1-Eb)*Ee)).^(L-w).*((1-Ea)*Eb*Ee+Ea*(1-Eb)*(1-Ee)).^w;
    
    S = 0;
    for w = 0:L
       % u = Dw(w)/(Dw(w)+Dw(N-w));
        S = S + nchoosek(L,w)*Dw(w)*(1-h(Dw(w)/(Dw(w)+Dw(N-w))));
    end
    Ie = 1/(E^L+(1-E)^L) * S;
 
    Rc = 1/L;
    K = log2(L);
    e = zeros(1,K); %заготовка из нулей
    e(1) = E; % начальное значение
    for i = 0:K-1
        Rc = Rc *( ((e(i+1)^2)^i) + (1-((e(i+1))^2)^i) );
        e(i+2) = e(i+1)^2 + (1-e(i+1))^2 ; 
    end
    Reff(k) = Rc*(Ib-Ie);
    Idiff(k)=(Ib-Ie);
end
Reff(k)
  
semilogy(Llist, Idiff )
title ( 'График зависимости Reff от длины сообщения')
xlabel('N');
ylabel('Reff'); 
set(0,'DefaultAxesFontSize',14,'DefaultAxesFontName','Times New Roman');
grid on

График должнен быть возрастающий и стремиться к 1(фигура во вложении)
Где-то есть ошибка,но не могу ее найти.
Считал все по формулам вроде-формулы во вложении.

График должен иметь примерно следующий вид(получил его если из вышеприведенного кода из строки 30 убрать h)

Формула для h в том же вложении под графиком.

@Centurio · 29.11.2017, 20:09

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

График должнен быть возрастающий и стремиться к 1(фигура во вложении)

График чего? Фигуры во вложении я не наблюдаю. Там несколько формул, причём не все, формулы для h, например, нет.
А. выложили, теперь вижу.
Половина формул неизвестно где используется. Как вычисляется Reff?

@sanekkkk · 29.11.2017, 20:14 **[ТС]**

Прошу прощения.
Вот такое должно быть в итоге(хотя бы примерно такой вид)
формула для H во вложении и ПРИМЕРНО ЖЕЛАЕМЫЙ ГРАФИК.
Собственно нужен график semilogy(Llist, Idiff )
Reff(k) = Rc*(Ib-Ie); эта формула верна
Глаза замылились,не могу найти ошибку
Формулы вроде перепроверил

@faraonCar · 29.11.2017, 22:02

Так поставьте везде Ib(k) и Rc(i), а то у Вас онb получают только одно значение, последнее, для k=length(Llist) и для i = K-1.

Добавлено через 1 минуту
Ну и S(w) правьте, как Вам сказали выше.

@sanekkkk · 29.11.2017, 22:31 **[ТС]**

Про S поправил
А про Ib(k) и Rc(i) не понял может .
Вставил ошибку выдает .

Код программы

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
clear all
close all
clc
 
Ea = 0.15;
Eb = 0.15;
Ee = 0.1;
  
N = 40;
Llist = 2.^(1:6);
 
for k = 1:length(Llist)
    L = Llist(k);
    
    E = (1-Ea)*Eb + Ea*(1-Eb)
    betta(k) = E^L/((E^L)+(1-E)^L)
    
    h = @(x) -x.*log2(x) - (1-x).*log2(1-x);
    %Ib = 1-h(betta(k));
    %Ib = real(Ib1)
     Ib(k) = 1 + betta(k)*log2(betta(k)) + ( 1 - betta(k))*log2(1-betta(k));
      
    Dw = @(w) ((1-Ea)*(1-Eb)*(1-Ee)+Ea*Eb*Ee).^(L-w).*(((1-Ea)*(1-Eb)*Ee+Ea*Eb*(1-Ee)).^w) +...
        (((1-Ea)*Eb*(1-Ee)+Ea*(1-Eb)*Ee)).^(L-w).*((1-Ea)*Eb*Ee+Ea*(1-Eb)*(1-Ee)).^w;
    
    S = 0;
    for w = 0:L
        u = Dw(w)/(Dw(w)+Dw(N-w));
        S = S + nchoosek(L,w)*Dw(w)*(1-h(u));
    end
    Ie(w) = 1/(E^L+(1-E)^L) * S;
 
    Rc = 1/L;
    K = log2(L);
    e = zeros(1,K); %заготовка из нулей
    e(1) = E; % начальное значение
    for i = 0:K-1
        Rc(i) = Rc *( ((e(i+1)^2)^i) + (1-((e(i+1))^2)^i) );
        e(i+2) = e(i+1)^2 + (1-e(i+1))^2 ; 
    end
    Reff(k) = Rc*(Ib-Ie);
    Idiff(k)=(Ib-Ie);
end
Reff(k)
  
semilogy(Llist, Idiff )
title ( 'График зависимости Reff от длины сообщения')
xlabel('N');
ylabel('Reff'); 
set(0,'DefaultAxesFontSize',14,'DefaultAxesFontName','Times New Roman');
grid on

Добавлено через 31 секунду
Скорректируете может быть?
Добавлено через 11 минут
Достаточно данных предоставил?

@faraonCar · 30.11.2017, 08:04

Посмотрите внимательно все циклы и исправьте их.

@sanekkkk · 30.11.2017, 08:28 **[ТС]**

Так я вам прислал исправленный вариант,но повторюсь ,что выдаёт ошибки

Добавлено через 2 минуты
21 ,31 и 38 строки надо поправить ?

@faraonCar · 30.11.2017, 10:37

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

@sanekkkk · 30.11.2017, 20:05 **[ТС]**

Код программы

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
clear all
close all
clc
 
Ea = 0.15;
Eb = 0.15;
Ee = 0.1;
  
N = 40;
Llist = 2.^(1:6);
 
for k = 1:length(Llist)
    L = Llist(k);
    
    E = (1-Ea)*Eb + Ea*(1-Eb)
    betta(k) = E^L/((E^L)+(1-E)^L)
    
    h = @(x) -x.*log2(x) - (1-x).*log2(1-x);
    %Ib = 1-h(betta(k));
    %Ib = real(Ib1)
     Ib(k) = 1 + betta(k)*log2(betta(k)) + ( 1 - betta(k))*log2(1-betta(k));
      
    Dw = @(w) ((1-Ea)*(1-Eb)*(1-Ee)+Ea*Eb*Ee).^(L-w).*(((1-Ea)*(1-Eb)*Ee+Ea*Eb*(1-Ee)).^w) +...
        (((1-Ea)*Eb*(1-Ee)+Ea*(1-Eb)*Ee)).^(L-w).*((1-Ea)*Eb*Ee+Ea*(1-Eb)*(1-Ee)).^w;
    
    S = 0;
    for w = 0:L
        u = Dw(w)/(Dw(w)+Dw(N-w));
        S = S + nchoosek(L,w)*Dw(w)*(1-h(u));
    end
    Ie(k) = 1/(E^L+(1-E)^L) * S;
 
    Rc = 1/L;
    K = log2(L);
    e = zeros(1,K); %заготовка из нулей
    e(1) = E; % начальное значение
    for i = 0:K-1
        Rc(i+1) = Rc *( ((e(i+1)^2)^i) + (1-((e(i+1))^2)^i) );
        e(i+2) = e(i+1)^2 + (1-e(i+1))^2 ; 
    end
    Reff(k) = Rc*(Ib(k)-Ie(k));
    Idiff(k)=(Ib(k)-Ie(k));
end
Reff(k)
  
semilogy(Llist, Idiff )
title ( 'График зависимости Reff от длины сообщения')
xlabel('N');
ylabel('Reff'); 
set(0,'DefaultAxesFontSize',14,'DefaultAxesFontName','Times New Roman');
grid on

Все равно ошибку выдает.
Где неверно сделал?
ругается на строку Reff(k) = Rc*(Ib(k)-Ie(k));

@faraonCar · 30.11.2017, 20:20

Строки с ошибкой логики

Matlab M
1
2
3
4
5
6
7
8
9
    for w = 0:L
        u = Dw(w)/(Dw(w)+Dw(N-w));     должно быть u(w)???
        S = S + nchoosek(L,w)*Dw(w)*(1-h(u)); должно быть S(w)=S(какой значение тут???? предыдущее наверно)+...
    end
 
    e(k) = 1/(E^L+(1-E)^L) * S; какой значение принимает S наверняка массив, а цикла нет
 
 
      Reff(k) = Rc*(Ib(k)-Ie(k)); значение Rc постоянно 1/L или то, что получено в цикле Rc(i+1)

@Centurio · 01.12.2017, 21:40

sanekkkk, в формуле для Rc произведение вычисляется до K-1. Чему равно K?

@sanekkkk · 02.12.2017, 11:24 **[ТС]**

ActionScript 3
1
 строка 34 K = log2(L);

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

Код программы

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
clear all
close all
clc
 
Ea = 0.15;
Eb = 0.15;
Ee = 0.1;
  
N = 40;
Llist = 2.^(1:6);
 
for k = 1:length(Llist)
    L = Llist(k);
    
    E = (1-Ea)*Eb + Ea*(1-Eb)
    betta(k) = E^L/((E^L)+(1-E)^L)
    
    h = @(x) -x.*log2(x) - (1-x).*log2(1-x);
    %Ib = 1-h(betta(k));
    %Ib = real(Ib1)
     Ib(k) = 1 + betta(k)*log2(betta(k)) + ( 1 - betta(k))*log2(1-betta(k));
      
    Dw = @(w) ((1-Ea)*(1-Eb)*(1-Ee)+Ea*Eb*Ee).^(L-w).*(((1-Ea)*(1-Eb)*Ee+Ea*Eb*(1-Ee)).^w) +...
        (((1-Ea)*Eb*(1-Ee)+Ea*(1-Eb)*Ee)).^(L-w).*((1-Ea)*Eb*Ee+Ea*(1-Eb)*(1-Ee)).^w;
    
    S = 0;
    for w = 0:L
        u(w) = Dw(w)/(Dw(w)+Dw(N-w));
        S = S + nchoosek(L,w)*Dw(w)*(1-h(u(w)));%тут не знаю как верно оформить-накопление суммы нужно
    end
    Ie(k) = 1/(E^L+(1-E)^L) * S;
 
    Rc = 1/L;
    K = log2(L);
    e = zeros(1,K); %заготовка из нулей
    e(1) = E; % начальное значение
    for i = 0:K-1
        Rc(i+1) = Rc *( ((e(i+1)^2)^i) + (1-((e(i+1))^2)^i) );
        e(i+2) = e(i+1)^2 + (1-e(i+1))^2 ; 
    end
    Reff(k) = Rc(i+1)*(Ib(k)-Ie(k));
    Idiff(k)=(Ib(k)-Ie(k));
end
Reff(k)
  
semilogy(Llist, Idiff )
title ( 'График зависимости Reff от длины сообщения')
xlabel('N');
ylabel('Reff'); 
set(0,'DefaultAxesFontSize',14,'DefaultAxesFontName','Times New Roman');
grid on

Matlab M
1
u = Dw(w)/(Dw(w)+Dw(N-w));

должно быть u(w)???- вроде поправил

Matlab M
1
S = S + nchoosek(L,w)*Dw(w)*(1-h(u));

должно быть S(w)=S -тут накопление суммы,но не знаю как оформить это

Matlab M
1
e(k) = 1/(E^L+(1-E)^L) * S;

какой значение принимает S наверняка массив, а цикла нет- есть же цикл по k???

Matlab M
1
Reff(k) = Rc*(Ib(k)-Ie(k));

значение Rc постоянно 1/L или то, что получено в цикле Rc(i+1) -НУЖНО ПОЛУЧАТЬ ЗНАЧЕНИЯ ИЗ ЦИКЛА

По поводу Rc я сам запутался если честно.
Давайте вместе посмотрим формулы и поймем может быть?ВО вложении формулы

@Centurio · 02.12.2017, 12:12

sanekkkk, пока скажу вот что:
При степени двойки, равной 5 и более, при вычислении получаются очень маленькие числа, что в дальнейшем приводит к возникновению в ответе неопределённости NaN (берётся логарифм от нуля). Поэтому корректно программа будет считать только до 2⁴.

@sanekkkk · 02.12.2017, 12:33 **[ТС]**

Хорошо,пусть до 24.Но давайте сделаем ее работу корректной.

Новые блоги и статьи Все статьи Все блоги /
Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это дополнительная запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .

@sanekkkk 1 / 1 / 0 Регистрация: 16.12.2012 Сообщений: 202 Записей в блоге: 3
	29.11.2017, 14:10 [ТС]
	ТО есть в данном случае для nchossek не побороть это предупреждение? Если хочу большие числа? Добавлено через 16 часов 4 минуты А как еще вывести на экран например что при L=2 Idiff= и т.д до L=2^5 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	29.11.2017, 17:00
	У вас N=70 или всё-таки N=7 в программе? 0

@sanekkkk 1 / 1 / 0 Регистрация: 16.12.2012 Сообщений: 202 Записей в блоге: 3
	29.11.2017, 17:22 [ТС]
	Давайте для ясности N=40 а Llist = 2.^(1:5); 0

@sanekkkk 1 / 1 / 0 Регистрация: 16.12.2012 Сообщений: 202 Записей в блоге: 3
	29.11.2017, 20:14 [ТС]
	Прошу прощения. Вот такое должно быть в итоге(хотя бы примерно такой вид) формула для H во вложении и ПРИМЕРНО ЖЕЛАЕМЫЙ ГРАФИК. Собственно нужен график semilogy(Llist, Idiff ) Reff(k) = Rc*(Ib-Ie); эта формула верна Глаза замылились,не могу найти ошибку Формулы вроде перепроверил Миниатюры 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	29.11.2017, 22:02
	Так поставьте везде Ib(k) и Rc(i), а то у Вас онb получают только одно значение, последнее, для k=length(Llist) и для i = K-1. Добавлено через 1 минуту Ну и S(w) правьте, как Вам сказали выше. 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	30.11.2017, 08:04
	Посмотрите внимательно все циклы и исправьте их. 0

@sanekkkk 1 / 1 / 0 Регистрация: 16.12.2012 Сообщений: 202 Записей в блоге: 3
	30.11.2017, 08:28 [ТС]
	Так я вам прислал исправленный вариант,но повторюсь ,что выдаёт ошибки Добавлено через 2 минуты 21 ,31 и 38 строки надо поправить ? 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	30.11.2017, 10:37
	Просмотрите все циклы внимательно и поправьте, пишу с телефона, листинг кода полностью не отображается, но то что я видел до этого во многих циклах одна из та же ошибка. Переменная получает последнее значение цикла. 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	01.12.2017, 21:40
	sanekkkk, в формуле для Rc произведение вычисляется до K-1. Чему равно K? 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	02.12.2017, 12:12
	sanekkkk, пока скажу вот что: При степени двойки, равной 5 и более, при вычислении получаются очень маленькие числа, что в дальнейшем приводит к возникновению в ответе неопределённости NaN (берётся логарифм от нуля). Поэтому корректно программа будет считать только до 2⁴. 0

@sanekkkk 1 / 1 / 0 Регистрация: 16.12.2012 Сообщений: 202 Записей в блоге: 3
	02.12.2017, 12:33 [ТС]
	Хорошо,пусть до 24.Но давайте сделаем ее работу корректной. 0