Странная ошибка в Матлабе

@Rykova · Регистрация: 28.11.2015

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

Всем здравствуйте! Столкнулась с очень странным результатом работы Матлаба.

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
    angle_max = 15;
    angle_rad = angle_max*pi/180;
    x2 = 100;
    X = linspace (0, tan(angle_rad)*x2, 50);
    theta = atan (X ./x2);
    k_x = k * sin(theta);
    j = 1;
    [Ti,Wi] = Points(N,0,2*pi);
    E_rez = zeros(size(X));
 
    
    for x1 = k_x
        p_new = @(theta) sqrt(2*(R^2-x1*R*cos(theta) - x2*R*sin(theta)));
        K_new = @(theta) -1i*k*p_new(theta).*(besselj(1,k*p_new(theta))+1i*bessely(1,k*p_new(theta)))/8;
        M = zeros(1,N+1);
        for i = 1:N+1
                Value = K_new(Ti(i));
                if (isinf(Value) | isnan(Value))
                    Value = 1/(2*R);
                end
                M(i) = Value;
                
        end
 
        E_rez(j)  = M.*Solve*(Wi');
        j = j+1;
 
    end

Solve - вектор размера N+1

Ошибка заключается в том, что при больших N массив E_rez становится нулевым, но если обратиться к элементу массива, к примеру, E_rez(50), то он выдает ненулевое значение. Не понимаю, как такое возможно. Кто-нибудь с таким встречался? Помогите пожалуйста.

@S_el · 28.11.2015, 12:35

Rykova, приведите полный пример кода демонстрирующий проблему.

@SSC · 30.11.2015, 07:19

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

for x1 = k_x

Внутри цикла x1 и k_x не изменяются, да и присваивание в качестве параметра цикла неудачное решение

Новые блоги и статьи Все статьи Все блоги /
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .

@S_el 2444 / 1842 / 406 Регистрация: 15.12.2013 Сообщений: 8,243
	28.11.2015, 12:35
	Rykova, приведите полный пример кода демонстрирующий проблему. 0