Items tagged with analysis

Feed

Hi,

Although everyone one might know, but still I want to clarify the concept of Sensitivity analysis. Suppose we have an objective function Z in which we have two variables X & Y and two parameters a and b. So in sensitivity analysis, we systematically change the values of parameters ( for eg. changing value of "a" by 20%, then 40%, then 60%, similiarly for parameter "b"). At a time, we only change the value of one parameter and the values of other parameters remain fixed to their original value. Then we would again solve the model ( ie. objective function Z) and calculate the new values of objective function and variables.  Then we try to see the impact of parameter change on the value of objective function Z and on variables X & Y.

The formula for sensitivity analysis for parameter a is  for objective function Z =

 (new values of objective function  - old value of objective function)*100 / (old value of objective function)

Here new value of objective function Z denotes when we use changed value of paramter a (keeping other parameter values fixed to their original value) and old value of objective function denotes value of Z with all parameters (including parameter a) at their original value

Although I can do it by a for loop but it is not an efficient procedure and also time consuming. Can anybody suggest whether any package suitable for sensitivity analysis is available in maple.

Thanks for your attentive reading of problem.

Thanks and Regards,

Nilesh

 

Hello,

I try his code:

 

for i to 6 do u[i, j] := (u[i-1, j]+u[i+1, j]+lambda*(u[i, j-1]+u[i, j+1])+sigma(h, r[i+1])*(u[i+1, j]-u[i-1, j]))/mu end do;


 
Error, too many levels of recursion.

I want to the  results of the analytic equations, for example:

                                                          
 u[1,j]= (u[0, j] + u[2, j] + lambda (u[1, j - 1] + u[1, j + 1])                                                       

      + sigma(h, r[2]) (u[2, j] - u[0, j]))/mu

.

.

.

.
Regards.

Dear,

I have a perfectly working when all parameters are known (figure 1), however I want to perform a sensitivity analysis by derivating the code if one parameter is unknown. Because of multiple possible answers and because of the complexity of the formula, I cannot run this script and get solutions. Any ideas how I can this calculation lighter so it is able to run? Values should be real and positive (so 1 or 2 solutions are the only one I'm interested in)

Any ideas, how I can make this code runnable? (file is below)

I'm stuck on this for a while now :/ So I hope someone will be able to help me

Many thanks in advance!l

Question.mw

Figure 1: [URL=http://s1240.photobucket.com/user/laggstar/media/Parameter%20f%20known.png.html][IMG]http://i1240.photobucket.com/albums/gg494/laggstar/Parameter%20f%20known.png[/IMG][/URL]

 

Figure 2: [URL=http://s1240.photobucket.com/user/laggstar/media/Parameter%20f%20unknown.png.html][IMG]http://i1240.photobucket.com/albums/gg494/laggstar/Parameter%20f%20unknown.png[/IMG][/URL]

I'm working on my master thesis and trying to solve this equation on maple, after I developed my equation and solving it for N=2,3 and then trying to solve it for N=4 and 5 it took so much long time(more than two days without get any answer), so i'm asking you now to help me reducing the time for solving this complicated equation. 

Shklovskii_N_spheres_with_tails.mw

a := .4; lB := .714; R := 3.2; Z := 24; L := 88; `ℓ` := 2*Z/eta; eta := 1; b := .7; N := 4; Lx = (-N*l+L)/N; r = Lx+2*R

.4

 

.714

 

3.2

 

24

 

88

 

48

 

1

 

.7

 

3

 

Lx = 88/3-l

 

r = Lx+6.4

(1)

F := N*(l*ln(R/a)-l*ln(l/R)+(l-`ℓ`)^2/(2*R))+((N-1)*((L-l)/(N-1)))*ln((L-l)/((N-1)*a))+2*(L-l)*ln((L-l)/((N-1)*a))/(N-1)+sum(((N-q)*(l-`ℓ`)*(l-`ℓ`))/(q*(2*R+l)+(q-1)*(L-l)/(N-1)), q = 1 .. N)+sum((2*N-2*q)*(l-`ℓ`)*ln(((L-l)/(N-1)+(2*q-1)*R+(q-1)*(L-l)/(N-1))/((2*q-1)*R+(q-1)*(L-l)/(N-1))), q = 1 .. N-1)+sum((2*(l-`ℓ`))*ln(((L-l)/(N-1)+(2*q-1)*R+(q-1)*(L-l)/(N-1))/((2*q-1)*R+(q-1)*(L-l)/(N-1))), q = 1 .. N)+sum((N-q-1)*((2*(L-l)/(N-1)+2*q*R+(q-1)*(L-l)/(N-1))*ln(2*(L-l)/(N-1)+2*q*R+(q-1)*(L-l)/(N-1))-2*((L-l)/(N-1)-2*q*R+(q-1)*(L-l)/(N-1))*ln((L-l)/(N-1)-2*q*R+(q-1)*(L-l)/(N-1))+(2*q*R+(q-1)*(L-l)/(N-1))*ln(2*q*R+(q-1)*(L-l)/(N-1))), q = 1 .. N-1)+sum(2*((2*(L-l)/(N-1)+2*q*R+(q-1)*(L-l)/(N-1))*ln(2*(L-l)/(N-1)+2*q*R+(q-1)*(L-l)/(N-1))-2*((L-l)/(N-1)-2*q*R+(q-1)*(L-l)/(N-1))*ln((L-l)/(N-1)-2*q*R+(q-1)*(L-l)/(N-1))+(2*q*R+(q-1)*(L-l)/(N-1))*ln(2*q*R+(q-1)*(L-l)/(N-1))), q = 1 .. N-1)+(2*(L-l)/(N-1)+2*N*R+(N-1)*((L-l)/(N-1)))*ln(2*(L-l)/(N-1)+2*N*R+(N-1)*((L-l)/(N-1)))-2*((L-l)/(N-1)+2*N*R+(N-1)*((L-l)/(N-1)))*ln((L-l)/(N-1)+2*N*R+(N-1)*((L-l)/(N-1)))+(2*N*R+(N-1)*((L-l)/(N-1)))*ln(2*N*R+(N-1)*((L-l)/(N-1)))

6.238324626*l-3*l*ln(.3125000000*l)+.4687500000*(l-48)^2+2*(88-l)*ln(110.0000000-1.250000000*l)+2.*(l-48.)^2/(6.400000000+l)+(l-48.)^2/(56.80000000+1.500000000*l)+6.*(l-48.)*ln(14.75000000-.1562500000*l)+4.*(l-48.)*ln((97.60000000-1.*l)/(53.60000000-.5000000000*l))+2.*(l-48.)*ln((148.-1.500000000*l)/(104.-1.*l))+3.*(94.40000000-1.*l)*ln(94.40000000-1.*l)-6.*(37.60000000-.5000000000*l)*ln(37.60000000-.5000000000*l)+35.64092141+2.*(144.8000000-1.500000000*l)*ln(144.8000000-1.500000000*l)-4.*(75.20000000-1.*l)*ln(75.20000000-1.*l)+2.*(56.80000000-.5000000000*l)*ln(56.80000000-.5000000000*l)+(195.2-2*l)*ln(195.2-2*l)-2*(151.2-(3/2)*l)*ln(151.2-(3/2)*l)+(107.2-l)*ln(107.2-l)

(2)

s := diff(F, l)

.9375000000*l-3*ln(.3125000000*l)+4.*(l-48.)/(6.400000000+l)-2.*(l-48.)^2/(6.400000000+l)^2-41.76167537-2*ln(110.0000000-1.250000000*l)-2.500000000*(88-l)/(110.0000000-1.250000000*l)-1.500000000*(l-48.)^2/(56.80000000+1.500000000*l)^2-.9375000000*(l-48.)/(14.75000000-.1562500000*l)-3.000000000*ln(144.8000000-1.500000000*l)+4.*ln(75.20000000-1.*l)-1.000000000*ln(56.80000000-.5000000000*l)+2*(l-48.)/(56.80000000+1.500000000*l)+2.*(l-48.)*(-1.500000000/(104.-1.*l)+1.*(148.-1.500000000*l)/(104.-1.*l)^2)*(104.-1.*l)/(148.-1.500000000*l)+6.*ln(14.75000000-.1562500000*l)+4.*ln((97.60000000-1.*l)/(53.60000000-.5000000000*l))+2.*ln((148.-1.500000000*l)/(104.-1.*l))+4.*(l-48.)*(-1./(53.60000000-.5000000000*l)+.5000000000*(97.60000000-1.*l)/(53.60000000-.5000000000*l)^2)*(53.60000000-.5000000000*l)/(97.60000000-1.*l)-3.*ln(94.40000000-1.*l)+3.000000000*ln(37.60000000-.5000000000*l)-2*ln(195.2-2*l)+3*ln(151.2-(3/2)*l)-ln(107.2-l)

(3)

solve(s, l)

104.4945628

(4)

NULL



Download Shklovskii_N_spheres_with_tails.mwShklovskii_N_spheres_with_tails.mw

 

I have the following expression 

lambda1-2.9881355172134833*lambda2+16.044692204509636*lambda1^2-114.31727386075951*lambda1*lambda2+202.36051939112247*lambda2^2

 

I want to get the expression in a form in which all the coeffs are fixed to 7 decimal places. 

 

I have tried with 

evalf[7](lambda1-2.9881355172134833*lambda2+16.044692204509636*lambda1^2-114.31727386075951*lambda1*lambda2+202.36051939112247*lambda2^2)

Please help me plot the graphs of the following equations:-

eq1 := diff(E[X(t)], t) = lambda-delta*E[X(t)]-beta*E[X(t)*V(t)];
 eq2 := diff(E[Y(t)], t) = -delta*E[Y(t)]+beta*exp(-rho*tau)*E[X(t)*V(t)];

eq3 := diff(E[V(t)], t) = gamma+N*kappa*E[Y(t)]-mu*E[V(t)]-beta*E[X(t)*V(t)]
NB:- E[X(t)]= the expected value of variable X at time t. 
Please help me know how to input this...

I am trying to model a multibody system using euler lagrange equation. There are 3 links in the system attached with each other by revolute joints. The first link has one end fixed with the revolute joint.

Length of links= l1,l2,l3

Angle with x axis=Theta1,Theta2,Theta3

Mass of links=m1,m2,m3

Moment of Inertia = M1,M2,M3

I am not able to find the solution of ODE system comprising of 3 ODEs which I get after solving Euler Lagrange Eq for...

 


Dear Sir,

I attached an example file to this email via speedyshare, please check it out, also please see the following pasted warning message from that file.

http://www.speedyshare.com/files/29834055/DSsolve_range_example.mw

My analysis gives a warning. Though I have results that look reliable, it is nasty to receive this kind of warning...

Actually what this warning means?
Does it affect the results badly?
How about reliability...

Hi all

I need to find the fourier coefiicients.  In the following Maple program I get the fourier transform of the signal which is fine, BUT actually I am interested in finding A0, An and Bn of the signal denoted by Plot(f(t),t=0..2*Pi);

step:= 2*Pi/2000;
   start_point1:=Pi/3+Pi/6;
   duration1:=Pi/3;
   start_point2:=4*Pi/3-Pi/6;
   duration2:=Pi/3;
   f:= t -> piecewise(t >= start_point1...

Hi:

I asked Maple to evaluate (using "evalf") an infinite sum and it says (very quickly) that the answer is identically 0 (not 10^(-7)).

When I numerically  evaluate a large number of terms of the sum, the answer gets smaller as the number of terms increases, so it seems that Maple is probably correct. The fact that the answer is immediate suggest that Maple is doing some analysis rather than adding a large number of numerical values.

I'd like to...

I'm trying to plot phase analysis for the Friedmann's equation, however when running the protocol and trying to plot these phase analysis, it comes up with error that rho wasn't allowed in the set of options!

Here's the code i'm having difficulty with!?!

> SFsys := proc (lambda, W, ics, T, v, u) local sys, flatF, rho0line, plot0, plotDSp, plotDSm, plotE, flatF_Plot, IC, df1, df2; 
sys := eval(DSys, {Lambda = lambda, gamma = W}); 

 

Hey!

I'm really new on using maple 14 and I need help please.

I want to either connect maple 14 software with an oscilloscope for data acquisition or import the data captured by the oscilloscope to maple 14 and manipulate this data with the software for dynamic system analysis so I was wondering if anyone knows a programming code or a command to do this please =( I'm stuck

 

Thanks for your attention...

Hi to you all, I'm brand new to this forum and I hope you could help me.

I wrote a Maple14 worksheet which perform data analysis, reading information from a set of file (usually 50 files) coming from an acquisition device. They are text file but the acquisition software produce them in the form:

namefile-0Ep.irr, namefile-1Ep.irr .... namefile-49Ep.irr

where "namefile" is the name I choose when i start the acquisition.

In order to perform...

I can't seem to find anything on cluster anlaysis (kmean etc...) in Maple's help.  I thought graph theory might contain something close but not really.  I haven't tried but I suppose it shouldn't be too hard to create a few procedures for such things?  Can someone provide a few examples?

Of the 3 M's Maple seems to have left this one out?

Page 1 of 1