Cause of: "Error: Inputs must be floats, namely single or double."
13 views (last 30 days)
Show older comments
I am trying to solve these time-dependent second order differential equations using ode45. My differential equations are
Coef1, coef2, coef3 are coefficients defined as arrays of real numbers. is also an array of real numbers. γ,ω,Ω are constants.
% ω in above differential equation
omega = 2262.72;
% Ω in above differential equation
Omega = 148.328;
% γ in above differential equation
gamma = 1.519;
tau = 0.1;
topt = 3;
tTHz = 5;
ts = linspace(0,10,1e5);
alpha = 191.7789;
beta = 0.82481;
m2 = 3.745755e-5;
erfnc = -sqrt(2*pi)*exp(-2*tTHz*1i-1/2)/32 * ( exp(4*tTHz*1i) + erf( real(sqrt(2).*ts - sqrt(2)*tTHz + sqrt(2)*1i/2) ) - exp(4*tTHz*1i) + erf( real(-sqrt(2).*ts + sqrt(2)*tTHz + sqrt(2)*1i/2) ) + 2*exp(2*tTHz*1i + 1/2) + 2*erf( real(sqrt(2).*ts - sqrt(2)*tTHz ) )*exp(2*tTHz*1i + 1/2)+1 ) ;
ts1 = ts;
coef1 = real( ( -sqrt(2)*exp(-2*tTHz*1i-1/2)/32 .* ( 2*sqrt(2)*exp( -(sqrt(2).*ts - sqrt(2)*tTHz + sqrt(2)*1i/2).^2 ) + 2*sqrt(2)*exp(4*tTHz*1i).*exp( -(-sqrt(2).*ts + sqrt(2)*tTHz + sqrt(2)*1i/2).^2 ) + 4*sqrt(2)*exp(2*tTHz*1i+1/2).*exp(-(sqrt(2).*ts-sqrt(2)*tTHz).^2))) ./ (erfnc + m2) );
ts2 = ts;
coef2 = real( alpha ./ (erfnc + m2) );
ts3 = ts;
coef3 = real( beta ./ (erfnc + m2) );
%Eopt is the last term in my differential equation expressed above.
Eopt = sin(omega.*ts).*exp(-4*log(2).*(ts-topt).^2 / tau.^2)
To solve the ODE I follwed the instruction described here https://www.mathworks.com/help/matlab/ref/ode45.html under ODE with Time-Dependent Terms. I wrote a function named "myode" and saved it in the same folder as rest of code.
function g = myode(t, X, omega, Omega, gamma, Eopt, ts1, coef1, ts2, coef2, ts3, coef3)
coef1 = interp1(ts1,coef1,t); % Interpolate the data set (ts1,coef1) at time t
coef2 = interp1(ts2,coef2,t); % Interpolate the data set (ts2,coef2) at time t
coef3 = interp1(ts3,coef3,t); % Interpolate the data set (ts3,coef3) at time t
g = {X(2) ; -coef1 *X(2) - coef2 * X(1) + coef3 * X(3) - gamma * X(2) + Eopt ; X(4) ; -omega^2 * X(3) + Omega^2 * X(1) - gamma * X(4)} ;
Now I solve for and .
% Time interval
tspan = 0 : 1/Fs :15;
% Initial condition
conds = [0;0;0;0];
[t,X] = ode45( @(t,X) myode(t, X, omega, Omega, gamma, Eopt, ts1, coef1, ts2, coef2, ts3, coef3), tspan, conds);
When I run this code I get erros saying
Error using odearugements
Inputs must be floats, namely single or double.
Error in ode45
odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
Please help me why I get these errors and how to fix it. If there is a better method to solve these differential equations please let me know.
0 Comments
Accepted Answer
Walter Roberson
on 4 Nov 2019
It is not valid to return a cell array from your ode function.
You are trying to return a series of values of different sizes. You must instead return exactly the same number of values as you have initial conditions (4 in your case)
I suspect you should be interpolating Eopt by time, not just coef1, coef2, coef3 . All of those appear to be using the same time base, so it is not clear why you are complicating matters by using different variables for the time bases.
More Answers (1)
Fabio Freschi
on 4 Nov 2019
What I can see is
1) Fs is missing. In any case, tspan shpuld be a 2 entries array with initial and final values
tspan = [0,15];
2) inside myode the output must be an double array, not a cell array.
g = [X(2) ; -coef1 *X(2) - coef2 * X(1) + coef3 * X(3) - gamma * X(2) + Eopt ; X(4) ; -omega^2 * X(3) + Omega^2 * X(1) - gamma * X(4)] ;
The problem is that if I change the curly brackets with square brackets, another error appears, because Eopt is a vector and not a scalar. Could you please check how you can change your Eopt?
1 Comment
Walter Roberson
on 4 Nov 2019
It is perfectly fine to pass a vector of length greater than 2 as tspan: it tells ode45 to return values at exactly those times, whereas with a tspan that is a vector of length 2, it would return whatever time values it felt like.
See Also
Categories
Find more on Error Functions in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!