octave-lyh: scripts/control/system/sysrepdemo.m comparison

comparison scripts/control/system/sysrepdemo.m @ 4771:b8105302cfe8

[project @ 2004-02-16 17:45:50 by jwe]

author	jwe
date	Mon, 16 Feb 2004 17:45:50 +0000
parents	65b3519ac3a1
children	4c8a2e4e0717

comparison

equal deleted inserted replaced

-:ef5e598f099b
+:b8105302cfe8
 ch_quit = ch_details+1;
 while(syschoice != ch_quit)
 disp(" ")
 syschoice = menu("Octave System Representation Menu", ...
 "General overview of system representation (DO THIS FIRST)", ...
-"Initialize a system (ss2sys, tf2sys, zp2sys)", ...
+"Initialize a system (ss, tf, zp)", ...
 "Extract data from a system(sys2ss, sys2tf, sys2zp, etc.)", ...
 "Update internal representation (sysupdate)", ...
 "View the internal contents of a system (sysout)", ...
 "Details of internal representation", ...
 "Return to main menu");
 disp("in correcting their programming errors while using the OCST.")
 disp("The details of the internal representation can be seen in ");
 disp(["menu option ",num2str(ch_details)]);
 disp("The data structure used in the OCST is called a \"system data structure.\"");
 disp("A system data structure is contstructed with one of:")
-disp("   fir2sys (FIR transfer function to system)")
+disp("   fir (FIR transfer function to system)")
-disp("   ss2sys (state space matrices to system)")
+disp("   ss (state space matrices to system)")
-disp("   tf2sys (SISO transfer function to system)")
+disp("   tf (SISO transfer function to system)")
-disp("   zp2sys (SISO zero/pole/leading coefficient to system)")
+disp("   zp (SISO zero/pole/leading coefficient to system)")
 disp(" ")
 disp(["These functions are discussed in in menu option ",num2str(ch_init)])
 disp("The data in a system may be extracted using ")
 disp("   sys2fir (FIR transfer function from system")
 disp("   sys2ss (state space matrices from system)")
 disp(" ");
 formopt = 0;
 while(formopt != 4)
 disp("Three data formats may be used to initialize a system:")
 formopt = menu("System data structure initialization menu", ...
-"State space form       (ss2sys)", ...
+"State space form       (ss)", ...
-"Transfer function form (tf2sys)", ...
+"Transfer function form (tf)", ...
-"zero-pole form         (zp2sys)", ...
+"zero-pole form         (zp)", ...
 "Return to System representation menu");
 if(formopt == 1)
 disp("State space representation of a system is based on the usual")
 disp("multi-variable differential equations")
 disp(" ")
 while(ssopt < ssquit)
 ssopt = menu("State space initialization examples", ...
 "Double integrator example", ...
 "Double delay (discrete-time) example", ...
 "Summing junction (D-matrix only) example", ...
-"ss2sys details (help ss2sys)", ...
+"ss details (help ss)", ...
 "return to system initialization menu", ...
 "return to system representation main menu");
 if(ssopt == 1)
 disp("Example: construct a system representation of a")
 disp("double integrator via state-space form")
 run_cmd
 cmd = "b = [0; 1];";
 run_cmd
 cmd = "c = [1, 0];";
 run_cmd
-cmd = "sys = ss2sys(a,b,c);";
+cmd = "sys = ss(a,b,c);";
 run_cmd
 disp("The state space form of the system is seen via sysout:")
 cmd = "sysout(sys)";
 run_cmd
 disp("Notice that the Octave controls  toolbox automatically")
 run_cmd
 cmd = "b = [0; 1];";
 run_cmd
 cmd = "c = [1, 0];";
 run_cmd
-cmd = "sys=ss2sys(a,b,c,[],1e-3);";
+cmd = "sys=ss(a,b,c,[],1e-3);";
 run_cmd
 cmd = "sysout(sys)";
 run_cmd
 disp("Notice that the D matrix was filled in automatically.")
 disp("This is done if D is input as the empty matrix.")
 disp("A summing junction that computes e(t) = r(t) - y(t) may be");
 disp("constructed as follows:");
 disp("First, we set the matrix D:")
 cmd = "D = [1, -1];";
 run_cmd
-disp("ss2sys allows the initialization of signal and state names")
+disp("ss allows the initialization of signal and state names")
 disp("(see option 4), so we initialize these as follows:")
-cmd = "inname = list(\"r(t)\",\"y(t)\");";
+cmd = "inname = {\"r(t)\",\"y(t)\"};";
 run_cmd;
 cmd = "outname = \"e(t)\";";
 run_cmd
 disp("Since the system is continous time and without states,")
-disp("the ss2sys inputs tsam, n, and nz are all zero:")
+disp("the ss inputs tsam, n, and nz are all zero:")
-cmd = "sys = ss2sys([],[],[],D,0,0,0,[],inname,outname);";
+cmd = "sys = ss([],[],[],D,0,0,0,[],inname,outname);";
 run_cmd
 disp("The resulting system is:")
 cmd = "sysout(sys)";
 run_cmd
 disp("A discrete-time summing block can be implemented by setting")
 disp("the sampling time positive:")
-cmd = "sys = ss2sys([],[],[],D,1e-3,0,0,[],inname,outname);";
+cmd = "sys = ss([],[],[],D,1e-3,0,0,[],inname,outname);";
 run_cmd
 disp("The resulting system is:")
 cmd = "sysout(sys)";
 run_cmd
 prompt
 elseif(ssopt == 4)
-help ss2sys
+help ss
 disp(" ")
 disp(" ")
 disp("Notice that state-space form allows a single system to have")
 disp("both continuous and discrete-time states and to have both continuous")
 disp("and discrete-time outputs.  Since it's fairly easy to make an")
 while(tfopt < 5)
 tfopt = menu("Transfer function initialization menu", ...
 "Continuous time initialization" , ...
 "Discrete time initialization" , ...
 "User specified signal names" , ...
-"tf2sys details (help tf2sys)", ...
+"tf details (help tf)", ...
 "Return to system initialization menu", ...
 "Return to system representation main menu");
 if(tfopt == 1) # continuous time
 disp("A transfer function is represented by vectors of the")
 disp("coefficients of the numerator and denominator polynomials");
 disp("is generated by the following commands:")
 cmd = "num = [5, -1]";
 run_cmd
 cmd = "denom = [1, -2, 6]";
 run_cmd
-cmd = "sys = tf2sys(num,denom);";
+cmd = "sys = tf(num,denom);";
 run_cmd
 disp("alternatively, the system can be generated in a single command:");
-cmd = "sys = tf2sys([5, -1], [1, -2, 6]);";
+cmd = "sys = tf([5, -1], [1, -2, 6]);";
 run_cmd
 disp("Notice the output of sys: it is an Octave data structure.")
 disp("The details of its member variables are explained under")
 disp("System Representation Menu option 5 (the details of system form)")
 disp(" ");
 elseif(tfopt == 2) # discrete time
 disp("A transfer function is represented by vectors of the")
 disp("coefficients of the numerator and denominator polynomials");
 disp("Discrete-time transfer functions require ")
 disp("the additional parameter of a sampling period:")
-cmd = "sys=tf2sys([5, -1], [1, 2, -6], 1e-3);";
+cmd = "sys=tf([5, -1], [1, 2, -6], 1e-3);";
 run_cmd
 cmd = "sysout(sys)";
 run_cmd
 disp("The OCST recognizes discrete-time transfer functions and")
 disp("accordingly prints them with the frequency domain variable z.");
 disp("names to the signals in the other examples.  We may initialize a transfer");
 disp("function with user-specified names as follows: Consider a simple ")
 disp("double-integrator model of aircraft roll dynamics with ")
 disp("input \"aileron angle\" and output \"theta\".  A ")
 disp("system for this model is generated by the command")
-cmd = "aircraft=tf2sys(1, [1, 0, 0], 0,\"aileron angle\",\"theta\");";          run_cmd
+cmd = "aircraft=tf(1, [1, 0, 0], 0,\"aileron angle\",\"theta\");";          run_cmd
 disp("The sampling  time parameter 0 indicates that the system")
 disp("is continuous time.  A positive sampling time indicates a")
 disp("discrete-time system (or sampled data system).")
 cmd = "sysout(aircraft)";
 run_cmd
 disp("These signal names are used in OCST plots and design functions.");
 disp("(Run the frequency response demo to see an example of the use of ");
 disp("signal names in plots.)")
 prompt
 elseif(tfopt == 4) # help
-help  tf2sys
+help  tf
 prompt
 elseif(tfopt == 6) # return to main menu
 formopt = 4;
 endif
 endwhile
 while(zpopt < 5)
 zpopt = menu("Zero-pole initialization menu", ...
 "Continuous time initialization" , ...
 "Discrete time initialization" , ...
 "User specified signal names" , ...
-"zp2sys details (help zp2sys)", ...
+"zp details (help zp)", ...
 "Return to system initialization menu", ...
 "Return to system representation main menu");
 if(zpopt == 1) # continuous time
 disp("A zero-pole form representation of a system includes vectors")
 disp("of the system poles and zeros and a scalar leading coefficient.");
 run_cmd
 cmd = "denom = [1, -2, 6]";
 run_cmd
 cmd = "k = 5";
 run_cmd
-cmd = "sys = zp2sys(num,denom,k);";
+cmd = "sys = zp(num,denom,k);";
 run_cmd
 disp("alternatively, the system can be generated in a single command:");
-cmd = "sys = zp2sys([5, -1],[1, -2, 6],5);";
+cmd = "sys = zp([5, -1],[1, -2, 6],5);";
 run_cmd
 disp("Notice the output of sys: it is an Octave data structure.")
 disp("The details of its member variables are explained under")
 disp("System Representation Menu option 5 (the details of system form)")
 disp(" ");
 elseif(zpopt == 2) # discrete time
 disp("A zero-pole form representation of a system includes vectors")
 disp("of the system poles and zeros and a scalar leading coefficient.");
 disp(" ")
 disp("Discrete-time systems require the additional parameter of a sampling period:")
-cmd = "sys=zp2sys([5, -1],[1, 2, -6],5,1e-3);";
+cmd = "sys=zp([5, -1],[1, 2, -6],5,1e-3);";
 run_cmd
 cmd = "sysout(sys)";
 run_cmd
 disp("The OCST recognizes discrete-time transfer functions and")
 disp("accordingly prints them with the frequency domain variable z.");
 disp("names to the signals in the other examples.  We may initialize a transfer");
 disp("function with user-specified names as follows: Consider a simple ")
 disp("double-integrator model of aircraft roll dynamics with ")
 disp("input \"aileron angle\" and output \"theta\".  A ")
 disp("system for this model is generated by the command")
-cmd = "aircraft=zp2sys([],[0, 0],1,0,\"aileron angle\",\"theta\");";            run_cmd
+cmd = "aircraft=zp([],[0, 0],1,0,\"aileron angle\",\"theta\");";            run_cmd
 disp("The sampling  time parameter 0 indicates that the system")
 disp("is continuous time.  A positive sampling time indicates a")
 disp("discrete-time system (or sampled data system).")
 cmd = "sysout(aircraft)";
 run_cmd
 disp("These signal names are used in OCST plots and design functions.");
 disp("(Run the frequency response demo to see an example of the use of ");
 disp("signal names in plots.)")
 prompt
 elseif(zpopt == 4) # help
-help  zp2sys
+help  zp
 prompt
 elseif(zpopt == 6) # return to main menu
 formopt = 4;
 endif
 endwhile
 endif
 endwhile
 elseif(syschoice == ch_extract)  # extract system information
 disp("Extract information from a system data structure in a selected format:")
-disp("The actions of operations ss2sys, tf2sys, and zp2sys are reversed by")
+disp("The actions of operations ss, tf, and zp are reversed by")
 disp("respective functions sys2ss, sys2tf, and sys2zp.  The latter two");
 disp("functions are applicable only to SISO systems.")
 formopt = 0;
 while(formopt != 8)
 formopt = menu("Extract system information", ...
 help sys2tf
 elseif(formopt == 3)
 help sys2zp
 elseif(formopt == 4)
 help sysgetsignals
-cmd="sys=ss2sys(rand(4),rand(4,2),rand(3,4));";
+cmd="sys=ss(rand(4),rand(4,2),rand(3,4));";
 run_cmd
 printf("Example: All signals names can be extracted by\n");
 cmd = "[Ast,Ain,Aout,Ayd] = sysgetsignals(sys)";
 run_cmd
 printf("Example: Input signal names can be extracted as\n");
 cmd = "Ain = sysgetsignals(sys,\"in\")";
 run_cmd
 printf("Example: The name of output signal 2 can be extracted as\n");
 cmd = "Aout = sysgetsignals(sys,\"out\",2)";
 run_cmd
-printf("\nNotice that Aout is returned as a list; the signal name\n");
+printf("\nNotice that Aout is returned as a cell array; the signal name\n");
 printf("itself is obtained by specifying the input parameter strflg\n");
 cmd = "Aout = sysgetsignals(sys,\"out\",2,1)";
 run_cmd
 prompt
 cmd = "help syssetsignals";
 prompt
 endwhile
 elseif(syschoice== ch_update)
 disp("The OCST system data structure format will store a system in the same format")
 disp("as that with which it was initialized.  For example, consider the following:")
-cmd = "sys=zp2sys([1, 2],[3, 4, 5],6)";
+cmd = "sys=zp([1, 2],[3, 4, 5],6)";
 run_cmd
 disp(" ")
 disp("Notice the internal variables in the structure include zer, pol, and k,")
 disp("the required variables for zero-pole form.  We can update the system")
 disp("to include state-space form as follows:")
 help sysupdate
 prompt
 elseif(syschoice == ch_view)
 disp("The sysout command can be used to view a system in any desired format.")
 disp("For example, consider the system created as follows:")
-cmd = "aircraft=zp2sys(1,[0, 0],1,0,\"aileron angle\",\"theta\");";             run_cmd
+cmd = "aircraft=zp(1,[0, 0],1,0,\"aileron angle\",\"theta\");";             run_cmd
 disp("The system may be viewed in its default format (zero-pole) as follows")
 cmd = "sysout(aircraft)";
 run_cmd
 disp(" ")
 disp("The system may be viewed in state-space or transfer function form as well:")

Mercurial > hg > octave-lyh

comparison scripts/control/system/sysrepdemo.m @ 4771:b8105302cfe8