Advanced Search

pik1432

245 Reputation

6 Badges

6 years, 122 days
Contact pik1432

MaplePrimes Activity


These are questions asked by pik1432

Difference between 'map' and '~(element-...

pik1432 245 Maple 2020
June 12 2021
3 3

Hello there, 

Would you allow me to ask this one question?

What is the difference between 'map' and '~(element-wise operation)'?

The following worksheet snipped shows an example:



Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/Q20210612.mw .
 

I thought that both operations were the same, but seemingly they are not. 

Best Regards, 

In Kwon Park 

Download Q20210612.mw

How to map an expression to x-axis and a...

pik1432 245 Maple 2019
plot parametric plotting
April 16 2021
1 5

Hello there, 

Would you allow me to ask this question?

What I tried in the worksheet below is to map an expression to x-axis and another expression to y-axis, while both expressions have the same input (or independent variable). However, the resulting plot still shows those two expressions as two separate curves. Would you tell me how to map or assign the expressions to each axis?


 

> 

restart;
> 

with(LinearAlgebra):
> 

Digits := 6:
> 

j := I:
> 

L__adpu := 1.66;
L__aqpu := 1.61;
L__lpu := 0.15;
R__apu := 0.003;
E__tpu := 1.0;
MVA := 555;

1.66

  

1.61

  

.15

  

0.3e-2

  

1.0

  

555

 (1)
> 

A__sat := 0.0312500:
B__sat := 6.93147:
> 

I__tpu_compare := (P__t - j* Q__t) / E__tpu;

-(1.00000*I)*Q__t+1.00000*P__t

 (2)
> 

I__tpu_a := abs(I__tpu_compare);

abs(-(1.00000*I)*Q__t+1.00000*P__t)

 (3)
> 

Phi := arccos(P__t / (I__tpu_a * E__tpu));

arccos(1.00000*P__t/abs(-(1.00000*I)*Q__t+1.00000*P__t))

 (4)
> 

E__a__tilde := E__tpu + (R__apu + j * L__lpu) * I__tpu_compare;

1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t)

 (5)
> 

psi__at := abs(E__a__tilde);

abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))

 (6)
> 

psi__I := A__sat * exp(B__sat * (psi__at - 0.8));

0.312500e-1*exp(6.93147*abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))-5.54518)

 (7)
> 

K__sd := psi__at / (psi__at + psi__I);

abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))/(abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))+0.312500e-1*exp(6.93147*abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))-5.54518))

 (8)
> 

K__sq := K__sd;

abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))/(abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))+0.312500e-1*exp(6.93147*abs(1.0+(0.3e-2+.15*I)*(-(1.00000*I)*Q__t+1.00000*P__t))-5.54518))

 (9)
> 

X__ad := K__sd * L__adpu:
> 

X__d := X__ad + L__lpu:
> 

X__aq := K__sq * L__aqpu:
> 

X__q := X__aq + L__lpu:
> 

delta__i := arctan((X__q * I__tpu_a * cos(Phi) - R__apu * I__tpu_a * sin(Phi))/(E__tpu + R__apu * I__tpu_a * cos(Phi) + X__q * I__tpu_a * sin(Phi))):
> 

e__q := E__tpu * cos(delta__i):
> 

i__d := I__tpu_a * sin(delta__i + Phi):
> 

i__q := I__tpu_a * cos(delta__i + Phi):
> 

i__fd := (e__q + R__apu * i__q + X__d * i__d)/(X__ad):
> 

i_armature := sqrt(i__d^2 + i__q^2):
> 

i_armature2 := subs(P__t = 0, i_armature):
> 

i_armature3 := subs(P__t = 0.5, i_armature):
> 

i_armature4 := subs(P__t = 1.0, i_armature):
> 

plotA := plot(i_armature2, Q__t = -0.6..0.6):
> 

plotB := plot(i_armature3, Q__t = -0.6..0.6):
> 

plotC := plot(i_armature4, Q__t = -0.6..0.6):
> 

plots:-display([plotA, plotB, plotC]);

 
> 

i__fd2 := subs(P__t = 0, i__fd):
> 

i__fd3 := subs(P__t = 0.5, i__fd):
> 

i__fd4 := subs(P__t = 1.0, i__fd):
> 

 

> 

plotD := plot(i__fd2, Q__t = -0.6..0.6):
> 

plotE := plot(i__fd3, Q__t = -0.6..0.6):
> 

plotF := plot(i__fd4, Q__t = -0.6..0.6):
> 

plots:-display([plotD, plotE, plotF]);

 
> 

plot([i__fd2, i_armature2, Q__t = -0.6..0.6]); #failed, two expressions are plotted seperately.

 
> 

 

> 

 


Best Regards, 

In Kwon Park 

Download Q20210416.mw

How to treat differentiation as an opera...

pik1432 245 Maple 2019
differentiation functional-operator
April 07 2021
1 1

Hello there, 

Would you tell me how to treat differentiation as an operator?

Here is a simple example:

(The content of the worksheet is supposed to be here, but that did not work)

Download Q20210407.mw

The 'coeff' operation shows 'psi__d0' as the only coefficient of 'Delta__delta(t)', while what I wanted to see is the 'Desired' expression in the worksheet, '(psi__d0 + psi__q0 / omega__0 * D)' as the coefficients. Here, 'D' represents the differentiation. 

Thank you, 

In Kwon Park 

RE: rewriting an expression with substit...

pik1432 245 Maple 2019
substitution duplicate-question
March 16 2021
2 5

Hello there, 

Would you allow me to ask this (perhaps simple) question?

My goal is to express an equation as 'desired', but with no success with algsubs()/subs()/simplify(). 

Would you please show me the correct way?

 

> 

restart:
> 

PowerBalanceEq := 0 = e1(t) * i1(t) + e2(t) * i2(t) + e3(t) * i3(t);

0 = e1(t)*i1(t)+e2(t)*i2(t)+e3(t)*i3(t)

 (1)
> 

eq_i1 := i1(t) = solve(PowerBalanceEq, i1(t));

i1(t) = -(e2(t)*i2(t)+e3(t)*i3(t))/e1(t)

 (2)
> 

n21eq := n21 = e2(t) / e1(t);

n21 = e2(t)/e1(t)

 (3)
> 

eq_i2 := algsubs(n21eq, eq_i1);

i1(t) = -(e2(t)*i2(t)+e3(t)*i3(t))/e1(t)

 (4)
> 

eq_i3 := subs(n21eq, eq_i1);

i1(t) = -(e2(t)*i2(t)+e3(t)*i3(t))/e1(t)

 (5)
> 

eq_i4 := simplify(eq_i1, {e2(t) / e1(t) = n21});

i1(t) = (-i2(t)*n21*e1(t)-e3(t)*i3(t))/e1(t)

 (6)
> 

desired := i1(t) = -n21*i2(t) - e3(t)*i3(t)/e1(t);

i1(t) = -n21*i2(t)-e3(t)*i3(t)/e1(t)

 (7)
> 

 


Best Regards, 

In Kwon Park 

Download Q20210316.mw

RE: how to check 'arctan(omega*L/R)' is ...

pik1432 245 Maple 2020
assumptions evalb verification
March 08 2021
2 2

Hello there, 

I would think that those two expressions are same as long as L::positive, omega::positive, R::positive, but when I tried in Maple, I could not get that result. The following worksheet shows my attempts. Would you please tell me how to get the answer, true, instead of false from the expression comparison?


 

> 

restart:
> 

arctan_expression := arctan(omega*L/R);

arctan(omega*L/R)

 (1)
> 

arctan_expression2 := arctan(R, L*omega);

arctan(R, L*omega)

 (2)
> 

evalb(arctan_expression = arctan_expression2) assuming L::positive, omega::positive, R::positive;

false

 (3)
> 

is(arctan_expression = arctan_expression2);

false

 (4)
> 

 


Best Regards, 

In Kwon Park 

Download Q20210308.mw

 

2 3 4 5 6 7 8 Last Page 4 of 12