☰

Como gerar uma tabela com colunas de cabeçalho especificadas no Matlab?

tag-icon • tag-icon matlab
Como gerar uma tabela com colunas de cabeçalho especificadas no Matlab?

Minha técnica de criar um vetor de variáveis ​​de "coluna" em uma tabela só funciona parte do tempo. O que há de diferente entre K1 e índice no código a seguir? Estou depurando um método numérico e preciso do índice das colunas X, K1, K2, K, Y. Tudo funciona bem até eu adicionar K1, K2 e K?

CÓDIGO DE FUNÇÃO:

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

CÓDIGO DE CHAMADA:

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

REGISTRO:

>> [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];
>>

Responder1

Este código é melhor. A primeira linha da saída pode ser alterada para corresponder melhor ao que é feito manualmente. Basicamente, o que isso significa é que as colunas K1, K2 e K podem ser deslocadas uma linha para cima. Outras colunas podem ser deixadas iguais. No entanto, isso corresponde a Edwards e Penney. Observe como K1, K2 e K são definidos aqui. Saúde! MILÍMETROS

Código de função:

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

Código de chamada:

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

Responder2

O erro está no fato de que na linha 22 K1=[K1;k1];tem K1 entre colchetes, mas não foi definido antes.

A solução é definir K1=[]; antes do loop for.

Editar: o mesmo vale para todas as outras variáveis. Então esse código funciona

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

informação relacionada