Numerical integration with many parameters.

Question

JACINTA ONWUKA 2019-6-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/465484-numerical-integration-with-many-parameters

编辑： Torsten 2019-6-4

采纳的回答： Torsten

I will like to have the values of n2 depending on 10 values of Mycp. please i need help with this my code below.

L=1;

T=100;

r=0.03;

I1=0.5;

p=0.005;

epsilon=0.5;

beta=0.1;

rho=0.5;

Mycp=0:1:10;

delta=1-Mycp/100

tau=(1/(beta*(L+delta*p)))*log((L*(I1+delta*p))/(delta*p*(L-I1)));

t05 =(1/(beta*(L+delta*p)))*log((L*(0.05*L+delta*p))/(delta*p*(L-0.05*L)));

I2= (L*delta*p*(exp (beta*(L+delta*p)*t)-1)) ./ (L + delta*p* exp(beta*(L+delta*p)*t));

fun2=@(t,rho,I2,r)rho*I2*exp(-r*t);

n2= integral(@(t)fun2(t,rho,I2,r),t05,tau, 'ArrayValued',1)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Torsten 2019-6-4

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/465484-numerical-integration-with-many-parameters#answer_377825

编辑：Torsten 2019-6-4

function main
  L = 1;
  T = 100;
  r = 0.03;
  I1 = 0.5;
  p = 0.005;
  epsilon = 0.5;
  beta = 0.1;
  rho = 0.5;
  Mycp = 0:1:10;
  n2 = zeros(numel(Mycp),1);
  
  for i = 1:numel(Mycp)
    delta = 1-Mycp(i)/100;
    tau = (1/(beta*(L+delta*p)))*log((L*(I1+delta*p))/(delta*p*(L-I1)));
    t05 =(1/(beta*(L+delta*p)))*log((L*(0.05*L+delta*p))/(delta*p*(L-0.05*L))); 
    I2= @(t)(L*delta*p*(exp (beta*(L+delta*p)*t)-1)) ./ (L + delta*p* exp(beta*(L+delta*p)*t));
    fun2 = @(t)rho*I2(t).*exp(-r*t);
    n2(i) = integral(fun2,t05,tau);
  end
  
  plot(Mycp,n2)
end