Eigenvalue problem associated with free surface waves

Javier Mora Torres

2021 1 2

0 个回答

1 次查看（30 天）

0 个投票

Hello,

I need to solve non homogeneous first order differential equation related to Rayleigh Waves.

How can I perform numerical integration in ode23 or ode45? I need recommendations. Thank you.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ritankar Sahu 2021-1-3

在 MATLAB Online 中打开

Hi,

I don't know if you are familier with the 'ode45' or not. So I am sharing a part of my one project I think you will get some idea from it. If afterwards you need some more, I can find for you.

Good Luck!!

% Initial Conditions
A=300;
B=1;
C=500;
[TOUT, YOUT] = ode45(@model, [0 1000], [A;B;C]); % that [0 1000] is the limit of simulation
figure("name", "A")
plot(TOUT, YOUT(:,3));
grid on;
title('Solution of State Space model with ODE45 for A');
function myFunc = model(TOUT, YOUT)
Fa=1.5;
Fb=0.5;
Fd=2;
Ta=300;
Tb=Ta;
H=10.^6;
c=5000;
h=500;
k=0.2*(1-exp(1-(X(3)/173)));
myFunc = [Fa+Fb-Fd; (1 - X(2)) / X(1) * Fa - X(2) / X(1) * Fb - X(2) * k; 
    (Ta - X(3)) / X(1) * Fa + (Tb - X(3)) / X(1) * Fb + H / (c * X(1)) - h / c * X(2) * k]; % Equations Matrix
end