Aakanksha

30 Reputation

4 Badges

9 years, 40 days

MaplePrimes Activity


These are replies submitted by Aakanksha

@Carl Love 

gamma[o] is just a variable. actually i want a the variation of function with respect to gamma[o]

capital GAMMA is the GAMMA function in maple

and summation is over p.

 

I am having maple 9.5

I recently opted maple for my simulations as it can solve special functions easily. I am totaly unaware of its functionalities and just trying to explore. So, my posts may bother you all.

Thanks for all suggestions.

@Carl Love 

Thanks for your solution. It worked. But I want my vertical axis in log axis.

restart:

k:=2.01: m:=2: L1:=1: L2:=3: L3:=5: with(plots):

Warning, the name changecoords has been redefined

a1:=(m*L1)+k-1: a2:=(m*L2)+k-1: a3:=(m*L3)+k-1:

b1:=k-(m*L1): b2:=k-(m*L2): b3:=k-(m*L3):

p1:=((((k*m)/snr)^((a1+1)/2))*MeijerG([[-a1/2,(1-a1)/2],[]],[[b1/2,-b1/2],[-(1+a1)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L1)*GAMMA(k)):

p2:=((((k*m)/snr)^((a2+1)/2))*MeijerG([[-a2/2,(1-a2)/2],[]],[[b2/2,-b2/2],[-(1+a2)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L2)*GAMMA(k)):

p3:=((((k*m)/snr)^((a3+1)/2))*MeijerG([[-a3/2,(1-a3)/2],[]],[[b3/2,-b3/2],[-(1+a3)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L3)*GAMMA(k)):

result1:=[seq([i,evalf(eval(p1,snr=i))],i=1..30)];

result1 := [[1, .1459286401], [2, 0.8869547820e-1], [3, 0.6236755560e-1], [4, 0.4727672818e-1], [5, 0.3756222199e-1], [6, 0.3083117361e-1], [7, 0.2592152085e-1], [8, 0.2220128203e-1], [9, 0.1929787603e-1], [10, 0.1697787638e-1], [11, 0.1508789788e-1], [12, 0.1352319762e-1], [13, 0.1220994115e-1], [14, 0.1109467362e-1], [15, 0.1013780300e-1], [16, 0.9309422223e-2], [17, 0.8586548553e-2], [18, 0.7951250986e-2], [19, 0.7389350534e-2], [20, 0.6889499432e-2], [21, 0.6442517024e-2], [22, 0.6040902793e-2], [23, 0.5678474074e-2], [24, 0.5350092865e-2], [25, 0.5051457327e-2], [26, 0.4778940890e-2], [27, 0.4529466759e-2], [28, 0.4300409197e-2], [29, 0.4089515096e-2], [30, 0.3894841233e-2]]

result2:=[seq([i,evalf(eval(p2,snr=i))],i=1..30)]:

result3:=[seq([i,evalf(eval(p3,snr=i))],i=1..30)]:

plots[pointplot](result1);

 

 

 

 

 


Download ak_BER_9sep.mw

@Markiyan Hirnyk 


restart:

k:=2.01: m:=2: L1:=1: L2:=3: L3:=5: with(plots):

Warning, the name changecoords has been redefined

a1:=(m*L1)+k-1: a2:=(m*L2)+k-1: a3:=(m*L3)+k-1:

b1:=k-(m*L1): b2:=k-(m*L2): b3:=k-(m*L3):

p1:=((((k*m)/snr)^((a1+1)/2))*MeijerG([[-a1/2,(1-a1)/2],[]],[[b1/2,-b1/2],[-(1+a1)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L1)*GAMMA(k)):

p2:=((((k*m)/snr)^((a2+1)/2))*MeijerG([[-a2/2,(1-a2)/2],[]],[[b2/2,-b2/2],[-(1+a2)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L2)*GAMMA(k)):

p3:=((((k*m)/snr)^((a3+1)/2))*MeijerG([[-a3/2,(1-a3)/2],[]],[[b3/2,-b3/2],[-(1+a3)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L3)*GAMMA(k)):

result1:=[seq([i,evalf(eval(p1,snr=i))],i=1..30)];

result1 := [[1, .1459286401], [2, 0.8869547820e-1], [3, 0.6236755560e-1], [4, 0.4727672818e-1], [5, 0.3756222199e-1], [6, 0.3083117361e-1], [7, 0.2592152085e-1], [8, 0.2220128203e-1], [9, 0.1929787603e-1], [10, 0.1697787638e-1], [11, 0.1508789788e-1], [12, 0.1352319762e-1], [13, 0.1220994115e-1], [14, 0.1109467362e-1], [15, 0.1013780300e-1], [16, 0.9309422223e-2], [17, 0.8586548553e-2], [18, 0.7951250986e-2], [19, 0.7389350534e-2], [20, 0.6889499432e-2], [21, 0.6442517024e-2], [22, 0.6040902793e-2], [23, 0.5678474074e-2], [24, 0.5350092865e-2], [25, 0.5051457327e-2], [26, 0.4778940890e-2], [27, 0.4529466759e-2], [28, 0.4300409197e-2], [29, 0.4089515096e-2], [30, 0.3894841233e-2]]

result2:=[seq([i,evalf(eval(p2,snr=i))],i=1..30)]:

result3:=[seq([i,evalf(eval(p3,snr=i))],i=1..30)]:

plots[pointplot](result1,axis[2]=[mode=log]);

Error, (in plot/options2d) unknown or bad argument: axis[2] = [mode = log]

 

 

 

 

 


Download ak_BER_9sep.mw

@Carl Love 

> plots[pointplot(result1)];
plots[INTERFACE_PLOT(POINTS([1., 0.1459286401], [2., 0.08869547820],

[3., 0.06236755560], [4., 0.04727672818], [5., 0.03756222199],

[6., 0.03083117361], [7., 0.02592152085], [8., 0.02220128203],

[9., 0.01929787603], [10., 0.01697787638], [11., 0.01508789788],

[12., 0.01352319762], [13., 0.01220994115], [14., 0.01109467362],

[15., 0.01013780300], [16., 0.009309422223], [17., 0.008586548553],

[18., 0.007951250986], [19., 0.007389350534], [20., 0.006889499432],

[21., 0.006442517024], [22., 0.006040902793], [23., 0.005678474074],

[24., 0.005350092865], [25., 0.005051457327], [26., 0.004778940890],

[27., 0.004529466759], [28., 0.004300409197], [29., 0.004089515096],

[30., 0.003894841233]))]

 

 

 

using plots[pointplot] returned me the above values.

@Carl Love 

How can i use logplot command to simultaneously plot some equations.

That means  I am having three eqs p1, p2 and p3 and I want a logplot of these eqs with respect to i from 0 to 30.

 

@Axel Vogt  I made certain changes to my function and got the same plot as u are having. Thanks for helping. Regards.



k:=1:

m:=1:

a:=k-m:

b:=k+m-1:

z:=(k*m)/snr:

x:=(z^((b+1)/2))/(0.693*GAMMA(k)*GAMMA(m));

x := 1.443001443/snr

c:=x*MeijerG([[-(b+1)/2],[(1-b)/2]],[[a/2,-a/2,-(b+1)/2,-(b+1)/2],[]],z);

c := 1.443001443*MeijerG([[-1], []], [[0, -1, -1], []], 1/snr)/snr

with(plots):

plot(c,snr=0..10);

kernelopts(version);

`Maple 9.50, IBM INTEL NT, Apr 7 2004 Build ID 155251`

 

 

 



Download capacity_kfading_ak.mw

restart: k:=5:

m:=10:

a:=k-m:

b:=k+m-1:

z:=(k*m)/snr:

x:=(z^((b+1)/2))/(0.693*GAMMA(k)*GAMMA(m));

x := 129444.1778*50^(1/2)*(1/snr)^(15/2)

c:=x*MeijerG([[-(b+1)/2],[(1-b)/2]],[[a/2,-a/2,-(b+1)/2,-(b+1)/2],[]],z);

Error, (in MeijerG/normal/parameters) function does not exist: the parameters lists are inconsistent in MeijerG([[-15/2], [-13/2]],[[-15/2, -15/2, -5/2, 5/2], []],<...>)

with(plots):

Warning, the name changecoords has been redefined

plot(c,snr=0..10);

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

Error, empty plot

kernelopts(version);

`Maple 9.50, IBM INTEL NT, Apr 7 2004 Build ID 155251`

 

 

 


Download capacity_kfading_ak_er.mw

@acer 

@acer 

code is working for low values of m but if I choose k=5 and m=10 then error is prompted.

@acer 

plotting with c and exact c gave me the same plot.  however for certain values of k amd m I am getting error for MeijerG function saying "parameters are inconsistent".

@Preben Alsholm 

OK. Thanks for showing me the plot. Will you please check out the same plot with snr defined between the range -5 to 10 and show me the output. 

@Preben Alsholm 

I am having maple 9.5

@Axel Vogt 

I used "with(RealDomain)" to solve my csc function.

However avoiding the use of with(RealDomain) also doesn't overcome the problem.

Thanks.

@gkokovidis 

I changed Digits:=15:  

but changing this doesn't proved to be a solution for my problem.

 

@Carl Love 

Thanks for helping.

I am using maple 9.5. I started working on a new worksheet and it is woking fine. However I encountered new error while plotting my function i.e.  "unable to evaluate the function to numeric values in the region";

which I think is due to negative arguments in MeijerG as evaluation of my function gives

0.62*((1/s)^4.0)*MeijerG([[-3.5,3.0],[]],[[-2.0,2.0],[-4.0]],4.0/s)

and I want to plot this equation for various values of s.

Do I need to specify my special functions using "with" as we do in case of with(plot)  orany other way to remove this error? 

 

1 2 3 Page 2 of 3