☰

Как вывести таблицу с указанными заголовками столбцов в Matlab?

tag-icon • tag-icon matlab
Как вывести таблицу с указанными заголовками столбцов в Matlab?

Мой метод создания вектора из переменных "столбцов" в таблице работает только часть времени. Чем отличается K1 от индекса в следующем коде? Я отлаживаю численный метод, и мне нужны индексы столбцов, X, K1, K2, K, Y. Все работает отлично, пока я не добавлю K1, K2 и 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;
k2=0;
k=0;
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];   Line 22
K2=[K2;k2];
K= [K;k];
end                % end loop
ImprovedEulerTable=table(index,X,K1,K2,K,Y)
end

КОД ВЫЗОВА:

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

БРЕВНО:

>> [index,X,K1,K2,K,Y] = impeulerT(0,1,1,10);
Undefined function or variable 'K1'.

Error in impeulerT (line 22)
K1=[K1;k1];

22  K1=[K1;k1];
>>

решение1

Этот код лучше. Первую строку вывода можно изменить, чтобы она лучше соответствовала тому, что делается вручную. По сути, это означает, что столбцы K1, K2 и K можно сместить на одну строку вверх. Другие столбцы можно оставить прежними. Однако это соответствует Edwards и Penney. Обратите внимание, как здесь определены K1, K2 и K. Ура! MM

Код функции:

function [index,X,Y,K1,K2,K] = 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
% Initialize the lower-case variables
k1=0;
k2=0;
k=0;
% Initialize the upper-case variables
K1=k1;
K2=k2;
K =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,Y,K1,K2,K)
end

Код вызова:

% Improved Euler
[index,X,Y,K1,K2,K] = impeulerTX(0,1,1,10);

решение2

Ошибка заключается в том, что в строке 22 K1=[K1;k1];внутри квадратных скобок находится K1, но ранее он не был определен.

Решение — определить K1=[]; перед циклом for.

Редактировать: то же самое касается всех остальных переменных. Так что этот код работает

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;
k2=0;
k=0;
% Initialize the upper-case variables
K1=[];
K2=[];
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];  % Line 22
  K2=[K2;k2];
  K= [K;k];
end                % end loop
ImprovedEulerTable=table(index,X,K1,K2,K,Y)
end

Связанный контент