I am trying to use some artificial data to see if my code is working.. but there is a error for the division and subtraction part.. See below the Code…
function Sa = trial(lambdaMax,lambda,T)t = 0;I = 0;Sa = [];u = {10,2,11,4,5,6};t = t - log(u)/lambdaMax;while t < T u = {10,2,11,4,5,6}; if (u <= lambda(t)/lambdaMax) I = I+1; Sa(I) = t; end u = {10,2,11,4,5,6}; t = t - log(u)/lambdaMax;endlambdaMax=50;T=20;lambda =@(Sa) lambdaMax*cos(Sa); Sa = trial(lambdaMax,lambda,T);figurehold on%plot(Sa,lambda(Sa),'*')
xlabel('t')ylabel ('cos(x)')X = linspace(min(Sa),max(Sa),10); Y = pchip(Sa,lambda(Sa),X);plot(X,Y)line(repmat(Sa,2,1),repmat([0;1],1,length(Sa)),'color','r' )
Thanks all in advance 🙂
Best Answer