How can I shrink the time of getting the results please?

Question

Avan Al-Saffar 2015-5-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/220318-how-can-i-shrink-the-time-of-getting-the-results-please

评论： Avan Al-Saffar 2015-6-9

I run the following code for 13 days until I got the results so could anyone advise me if there is any truck that I can use to shrink the time please? Because I need to run the code for different values!

function FisherN01analyticalsmallomega
y0=1;
a=1
c=50;
b = 2.5 
t =(0:0.1:10000);
 y = logisticOscilanalytical(t,c,y0,b,a);
Y = @(T) interp1(t,y,T);
g = @(t)  ( (b.*sin(b.*t) + cos(b.*t).* (2.*Y(t)./c)).^2./(a.^2.*sin(b.*t).^3 .* (Y(t) - (Y(t).^2).^2))) );
t= (0:10:10000);
N = length(t);
JT = zeros(1,N);
       for i = 1:N
           JT = integral(g,0,t(1,i),'ArrayValued',true);
       end
           figure(1)
           plot(t,JT,'-','LineWidth',2)
1;
% function y = logisticOscilanalytical(t, c, y0, b, a)
% ee = exp( (a/ b) * (1- cos( b*t)));
% y = (-c.* y0.* ee)./ ( (y0 - c) - y0.* ee);
% end

2 个评论
显示无隐藏无

Jan 2015-5-28

I've formatted the code, because it was not readable. Please use the "{} Code" button. Do not post code which is commented out - this confuses the readers.

Avan Al-Saffar 2015-5-28

Many thanks and sorry.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mike Hosea 2015-6-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/220318-how-can-i-shrink-the-time-of-getting-the-results-please#answer_181203

编辑：Mike Hosea 2015-6-1

在 MATLAB Online 中打开

That code never gave you results because it isn't working--there is a syntax error. Specific suggestions are therefore hard to provide. Some more or less generic suggestions would be:

Debug your code using a smaller t range.
Avoid the 'ArrayValued' option unless your function is really array-valued. Using that option with a scalar-valued problem avoids the need to make your integrand handle array inputs, but it slows INTEGRAL significantly. If only compensating for an integrand that can't handle array inputs (and ONLY if the integrand can't handle array inputs), consider using ARRAYFUN, e.g. integrate

g1 = @(t)arrayfun(g,t);

If you're going to use an interpolated function (this will affect accuracy negatively), supply the abscissas as 'Waypoints'. In this example you'll need to save the original t values (they are overwritten in the supplied code).

q = integral(g1,a,b,'Waypoints',t);

If you want to generate a table of results, consider only integrating between the tabulated points and adding to an accumulator variable as you go, i.e.

JT = zeros(1,N);
s = 0;
for i = 2:N
    q = integral(g,t(1,i-1),t(1,i));
    s = s + q;
    JT(i) = s; 
    % The above could be JT(i) = JT(i-1) + q; to do without s.
end

It looks like your integrand has a singularity when sin(b*t) == 0. You'll have to integrate between those values. Unfortunately, they can't just be 'Waypoints', and there is no guarantee that INTEGRAL will be able to handle the singularity.