Changeset 6a02de in git
- Timestamp:
- Aug 11, 2004, 11:59:59 AM (19 years ago)
- Branches:
- (u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
- Children:
- 9433a26de2fd01924b6e4703e7e8faac77f9e821
- Parents:
- 597f38ebceb55e2871e4816fe483390649592d80
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/control.lib
r597f38 r6a02de 1 version="$Id: control.lib,v 1.1 2 2004-08-10 14:47:02plural Exp $";1 version="$Id: control.lib,v 1.13 2004-08-11 09:59:59 plural Exp $"; 2 2 category="Miscellaneous"; 3 3 info=" … … 21 21 AUXILIARY PROCEDURES: 22 22 ncdetection(ring r); computes an ideal, presenting an involution map on non-comm algebra r; 23 involution(m, map theta); applies the involution, presented by theta to m of typpoly, vector, ideal, module;24 declare(string NameOfRing, string Variables[,string Parameters, string Ordering]); defines the ring, optional parametes are a string of parameters and a sting ofordering,23 involution(m, map theta); applies the involution, presented by theta to m of type poly, vector, ideal, module; 24 declare(string NameOfRing, Variables[,string Parameters, Ordering]); defines the ring, optional parametes are strings of parameters and ordering, 25 25 view(); Well-formatted output of lists, modules and matrixes 26 26 "; … … 386 386 example 387 387 {"EXAMPLE:";echo = 2; 388 //-------------------------------389 ring r=0,(s1,s2,s3),dp;390 module R=[0,-s3,s2],391 [s3,0,-s1];392 R=transpose(R);393 view( R );394 view( control(R) );395 396 module Rc=[0,-s3,s2],397 [s3,0,-s1],398 [-s2,s1,0];399 Rc=transpose(Rc);400 view( Rc );401 view( control(Rc) );402 kill r;403 //----------------------------------404 //reflector antenna405 ring A = (0, K1, K2, Te, Kp, Kc),(Dt, delta), dp;406 module R;407 R = [Dt, -K1, 0, 0, 0, 0, 0, 0, 0],408 [0, Dt+K2/Te, 0, 0, 0, 0, -Kp/Te*delta, -Kc/Te*delta, -Kc/Te*delta],409 [0, 0, Dt, -K1, 0, 0, 0, 0, 0],410 [0, 0, 0, Dt+K2/Te, 0, 0, -Kc/Te*delta, -Kp/Te*delta, -Kc/Te*delta],411 [0, 0, 0, 0, Dt, -K1, 0, 0, 0],412 [0, 0, 0, 0, 0, Dt+K2/Te, -Kc/Te*delta, -Kc/Te*delta, -Kp/Te*delta];413 R=transpose(R);414 view(R);415 view(control(R));416 kill A;417 //----------------------------------418 //Flexible Rod419 ring A = 0,(D1, delta), dp;420 module R;421 R = [D1, -D1*delta, -1],422 [2*D1*delta, -D1-D1*delta^2, 0];423 R=transpose(R);424 view(R);425 view(control(R));426 kill A;427 //-------------------------------------428 //TwoPendula429 ring r=(0,m1,m2,M,g,L1,L2),Dt,dp;430 module R;431 R = [m1*L1*Dt^2, m2*L2*Dt^2, -1, (M+m1+m2)*Dt^2],432 [m1*L1^2*Dt^2-m1*L1*g, 0, 0, m1*L1*Dt^2],433 [0, m2*L2^2*Dt^2-m2*L2*g, 0, m2*L2*Dt^2];434 R=transpose(R);435 view(R);436 view(control(R));437 kill r;438 //----------------------------------------439 388 //Wind Tunnel 440 389 ring A = (0,a, omega, zeta, k),(D1, delta),dp; … … 446 395 view(R); 447 396 view(control(R)); 448 kill A;449 //-------------------------------------------450 397 451 398 }; … … 539 486 example 540 487 {"EXAMPLE:"; echo = 2; 488 //Cauchy 541 489 ring r=0,(s1,s2,s3,s4),dp; 542 490 module R= [s1,-s2], … … 572 520 example 573 521 {"EXAMPLE:"; echo = 2; 522 //Cauchy 574 523 ring r=0,(s1,s2,s3,s4),dp; 575 524 module R= [s1,-s2], … … 585 534 // 586 535 //Some example rings with defined systems 587 //-------------------------------------------------------------------------- 536 //---------------------------------------------------------------------------- 588 537 //autonomy: 589 538 // … … 597 546 return(@r); 598 547 }; 599 //------------------------------------------------------- 548 //---------------------------------------------------------------------------- 600 549 proc exAut2() 601 550 { … … 609 558 }; 610 559 611 //----------------------------------- 560 //---------------------------------------------------------------------------- 612 561 proc exCauchy() 613 562 { … … 619 568 return(@r); 620 569 }; 621 //----------------------------------- 570 //---------------------------------------------------------------------------- 622 571 proc exCauchy2() 623 572 { … … 631 580 return(@r); 632 581 }; 633 //---------------------------------------- 582 //---------------------------------------------------------------------------- 634 583 proc exEvasLecture() 635 584 { … … 641 590 return(@r); 642 591 }; 643 //--------------------------------------------------------------------------- 592 //---------------------------------------------------------------------------- 644 593 //control 645 594 // 595 proc ex1() 596 { 597 ring @r=0,(s1,s2,s3),dp; 598 module R=[0,-s3,s2], 599 [s3,0,-s1]; 600 R=transpose(R); 601 export R; 602 return(@r); 603 }; 604 //---------------------------------------------------------------------------- 605 proc ex2() 606 { 607 ring @r=0,(s1,s2,s3),dp; 608 module R=[0,-s3,s2], 609 [s3,0,-s1], 610 [-s2,s1,0]; 611 R=transpose(R); 612 export R; 613 return(@r); 614 }; 615 //---------------------------------------------------------------------------- 646 616 proc exAntenna() 647 617 { … … 660 630 }; 661 631 662 //---------------------------------------------------------------------------- -----------------632 //---------------------------------------------------------------------------- 663 633 664 634 proc exEinstein() … … 722 692 //-------------------------------------------------------------------------------------------- 723 693 724 proc ex()725 {726 };727 694 728 695 //---------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.