Olá, estou tentando fazer uma tabela dentro do Matlab, mas não está saindo com os resultados corretos.
Rotina no Matlab:
alpha=0.008;
gamma=3.3;
g=9.80665; % aceleração gravitacional
delta_1=0.07; % se f<=fp
delta_2=0.09; % se f>fp
Tp=2:16; % período de pico do espectro [s] -> variando de 2 a 16 segundos
fp=1./Tp; % frequência do pico do espectro de onda [Hz]
T=2:16; % período da onda
f=1./T; % frequência da onda
nrows = 14; %16-2 do período de pico do espectro
ncols = 14;
E = ones(nrows,ncols);
for fp = 1:ncols % coluna
for f = 1:nrows % linha
if f<=fp
E(f,fp) = (alpha*(g^2)*((2*pi)^(-4))*(f^(-5)))*exp((-5/4)*((f/fp)^(-4)))*gamma*exp((-1/2)*((f-fp)/(delta_1*fp))^2);
elseif f>fp
E(f,fp) = (alpha*(g^2)*((2*pi)^(-4))*(f^(-5)))*exp((-5/4)*((f/fp)^(-4)))*gamma*exp((-1/2)*((f-fp)/(delta_2*fp))^2);
else
E(f,fp) = 0;
end
end
end
E
Em anexo tem uma figura que é a realidade (como saiu) e outra a expectativa (como era para ser).
Desde já, agradeço pela ajuda