Как вывести таблицу в Matlab, изменив оператор функции?

Question

Моя техника создания вектора из «столбцовых» переменных в таблице работает только часть времени.Чем отличается K1 от индекса в следующем коде?Я отлаживаю числовой метод, и мне нужен индекс столбцов, X, K1, K2, K, Y. Все работает нормально, пока я не добавлю K1, K2 и K?Как это сделать?

ММ

Что нужно для правильной инициализации K1, k1, K2, k2, K и k.Исправленный код приведен ниже.

ФУНКЦИОНАЛЬНЫЙ КОД:

function [index,X,K1,K2,K,Y] = impeulerT(x,y,x1,n)
% modified version of Improved Euler method found in
% Elementary Differential Equations by Edwards and Penney 
X=x;               % initial x
Y=y;               % initial y
x1 = x1;           % final x
n = n;             % number of subintervals
h = (x1-x)/n;      % step size
index = 0;         % initialize index
k1=0; K1=k1;       % initialize k1
k2=0; K2=k2;       % initialize k2
k=0; K=k;          % initialize k
for i=1:n;         % begin loop
k1=f(x,y);         % first slope
k2=f(x+h,y+h*k1);  % second slope
k=(k1+k2)/2;       % average slope
x=x+h;             % new x
y=y+h*k;           % new y
X=[X;x];           % update x-column       
Y=[Y;y];           % update y-column
index = [index;i]; % update index-column
K1=[K1;k1];        % update K1 column
K2=[K2;k2];        % update K2 column
K= [K;k];          % update K column
end                % end loop
ImprovedEulerTable=table(index,X,K1,K2,K,Y)
clear
end

КОД ВЫЗОВА:

[index,X,K1,K2,K,Y] = impeulerT(0,1,1,10);

ЖУРНАЛ:

Output argument "index" (and maybe others) not
assigned during call to "impeulerT".

Answer 1

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

Чтобы исправить свой код, просто удалите или прокомментируйте оператор clear после создания таблицы:

function [index,X,K1,K2,K,Y] = impeulerT(x,y,x1,n)
% modified version of Improved Euler method found in
% Elementary Differential Equations by Edwards and Penney 
    X=x;               % initial x
    Y=y;               % initial y
    x1 = x1;           % final x
    n = n;             % number of subintervals
    h = (x1-x)/n;      % step size
    index = 0;         % initialize index
    k1=0; K1=k1;       % initialize k1
    k2=0; K2=k2;       % initialize k2
    k=0; K=k;          % initialize k
    for i=1:n;         % begin loop
        k1=(x/y);         % first slope
        k2=(x+h)/(y+h*k1);  % second slope
        k=(k1+k2)/2;       % average slope
        x=x+h;             % new x
        y=y+h*k;           % new y
        X=[X;x];           % update x-column       
        Y=[Y;y];           % update y-column
        index = [index;i]; % update index-column
        K1=[K1;k1];        % update K1 column
        K2=[K2;k2];        % update K2 column
        K= [K;k];          % update K column
    end                % end loop
    ImprovedEulerTable=table(index,X,K1,K2,K,Y);
end

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