Перевести код из MatCad в MatLab

@jekkon · Регистрация: 26.05.2020

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

Коллеги, пожалуйста помогите перевести код из MatCad в MatLab или укажите где у меня ошибка. Код Маткада во вложении
https://drive.google.cоm/оpen?... XltKAHJdqd
мой код матлаб

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
function [Z_k] = vysota(I, X, Y, U, g)
 
Z_k=zeros(length(X),length(Y));
 
    a1=zeros(1,I);
    a2=zeros(1,I);
    a3=zeros(1,I);
    a4=zeros(1,I);
    H=zeros(1,I);
    phi=zeros(1,I);
    
sigma_v = U^2/((10)*sqrt(2)*(pi)*g); % СКО волнения
k_0 = g/U^2;                % граничное волновое число
        
for n=1:I  
    disp(n)
        a1(n) = rand;
        a2(n) = rand;
        a3(n) = rand;
        a4(n) = rand;
        phi(n)=2*pi*a4(n);                                                  %начальные фазы;                  
        H(n) =(sigma_v*sqrt(-2*log(a1(n))))./sqrt(I);                       % амплитуды парциальных плоских волн
        k = (k_0*[sin(2*pi*a3(n))  cos(2*pi*a3(n))])/a2(n);                 % волновые векторы  
        
    for i=1:length(X)
         
        for j=1:length(Y)
            Z_k(i,j)= (H(n)*cos(k*[X(i) Y(j)].'+ phi(n)));              % высота точки в координате  XY
        end
      
    end
 
%     A=zeros(i,j);
%     Z=A+Z_k(:,:);
%     Z_k1=Z;
%     B=Z_k1+Z_k;
end
 
return

@jekkon · 26.05.2020, 17:10 **[ТС]**

Всем спасибо! Решил сам! Пример оставляю для сообщества:

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
function [Z_k] = vysota(I, X, Y, U, g)
 
Z_k=zeros(length(X),length(Y));
 
    a1=zeros(1,I);
    a2=zeros(1,I);
    a3=zeros(1,I);
    a4=zeros(1,I);
    H=zeros(1,I);
    phi=zeros(1,I);
    
sigma_v = U^2/(sqrt(10)*sqrt(2)*sqrt(pi)*g); % СКО волнения
disp(sigma_v)
k_0 = g/U^2;                % граничное волновое число
    A=zeros(length(X),length(Y));    
for n=1:I  
    disp(n)
        
        a1(n) = rand;
        a2(n) = rand;
        a3(n) = rand;
        a4(n) = rand;
        phi=2*pi*a4(n);                                                  %начальные фазы;                  
        H =(sigma_v/sqrt(I))*sqrt(-2*log(a1(n)));                       % амплитуды парциальных плоских волн
    disp(H)
        k = (k_0*[sin(2*pi*a3(n))  cos(2*pi*a3(n))])/a2(n);                 % волновые векторы  
    disp(k)
    
    for i=1:length(X)
         
        for j=1:length(Y)
            
            Z_k(i,j)= (H*cos(k*[X(i) Y(j)]'+ phi));              % высота точки в координате  XY
        
        end
        
      
    end
 
     
     A=A+Z_k(:,:);
end
 
Z_k=A;
return

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию 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 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта 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 На первой гифке отладочные линии отключены, а на второй включены:. . .