From bug-octave-request at bevo dot che dot wisc dot edu Wed Jun 7 09:46:39 2000 Subject: c2d matched pole/zero equivalent. From: Ben Sapp To: bug-octave at bevo dot che dot wisc dot edu Date: Wed, 7 Jun 2000 02:57:28 -0600 --Boundary-=_nWlrBbmQBhCDarzOwKkYHIDdqSCD Content-Type: text/plain Content-Transfer-Encoding: 8bit Bug report for Octave from the CVS about a week ago configured for RedHat Linux Description: ----------- The c2d function does not have the ability to use matched pole/zero equivalent for conversion. Fix: --- Apply the attached patch. The patch is not perfect in that it does not work on systems that are not SISO or purely continuous. The books I have do not discuss how to handle MIMO systems. If someone can suggest a place to read up on how to do this for MIMO and discrete/continuous systems I would do it. Configuration (please do not edit this section): ----------------------------------------------- uname output: Linux pcben 2.2.12-20 #1 Mon Sep 27 10:40:35 EDT 1999 i686 unknown configure opts: %config_opts% Fortran compiler: g77 FFLAGS: -O F2C: F2CFLAGS: FLIBS: -lg2c -lm -L/usr/local/lib/gcc-lib/i686-pc-linux-gnu/2.95.2 -L/usr/local/lib -lm CPPFLAGS: INCFLAGS: -I. -I. -I./liboctave -I./src -I./libcruft/misc -I./glob -I./glob C compiler: gcc, version 2.95.2 19991024 (release) CFLAGS: -g -O2 -Wall CPICFLAG: -fPIC C++ compiler: c++, version 2.95.2 19991024 (release) CXXFLAGS: %OCTAVE_CONF_BUG_CXXFLAGS% CXXPICFLAG: -fPIC LDFLAGS: -g LIBFLAGS: -L. RLD_FLAG: -Xlinker -rpath -Xlinker /usr/local/lib/octave-2.1.30 TERMLIBS: -lncurses LIBS: LEXLIB: LIBPLPLOT: LIBDLFCN: LIBGLOB: ./glob/libglob.a DEFS: -DOCTAVE_SOURCE=1 -DSEPCHAR=':' -DSEPCHAR_STR=":" -DUSE_READLINE=1 -DBOUNDS_CHECKING=1 -D__NO_MATH_INLINES=1 -DCXX_NEW_FRIEND_TEMPLATE_DECL=1 -DHAVE_LIBM=1 -DF77_APPEND_UNDERSCORE=1 -DSIZEOF_SHORT=2 -DSIZEOF_INT=4 -DSIZEOF_LONG=4 -DHAVE_ALLOCA_H=1 -DHAVE_ALLOCA=1 -DNPOS=std::string::npos -DSTDC_HEADERS=1 -DHAVE_DIRENT_H=1 -DTIME_WITH_SYS_TIME=1 -DHAVE_SYS_WAIT_H=1 -DHAVE_ASSERT_H=1 -DHAVE_CURSES_H=1 -DHAVE_DLFCN_H=1 -DHAVE_FCNTL_H=1 -DHAVE_FLOAT_H=1 -DHAVE_FNMATCH_H=1 -DHAVE_GLOB_H=1 -DHAVE_GRP_H=1 -DHAVE_LIMITS_H=1 -DHAVE_MEMORY_H=1 -DHAVE_NCURSES_H=1 -DHAVE_POLL_H=1 -DHAVE_PWD_H=1 -DHAVE_SGTTY_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_SYS_IOCTL_H=1 -DHAVE_SYS_PARAM_H=1 -DHAVE_SYS_POLL_H=1 -DHAVE_SYS_RESOURCE_H=1 -DHAVE_SYS_SELECT_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_SYS_TIME_H=1 -DHAVE_SYS_TIMES_H=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_UTSNAME_H=1 -DHAVE_TERMCAP_H=1 -DHAVE_TERMIO_H=1 -DHAVE_UNISTD_H=1 -DHAVE_VARARGS_H=1 -DHAVE_ATEXIT=1 -DHAVE_BCOPY=1 -DHAVE_BZERO=1 -DHAVE_DUP2=1 -DHAVE_ENDGRENT=1 -DHAVE_ENDPWENT=1 -DHAVE_EXECVP=1 -DHAVE_FCNTL=1 -DHAVE_FORK=1 -DHAVE_GETCWD=1 -DHAVE_GETEGID=1 -DHAVE_GETEUID=1 -DHAVE_GETGID=1 -DHAVE_GETGRENT=1 -DHAVE_GETGRGID=1 -DHAVE_GETGRNAM=1 -DHAVE_GETHOSTNAME=1 -DHAVE_GETPGRP=1 -DHAVE_GETPID=1 -DHAVE_GETPPID=1 -DHAVE_GETPWENT=1 -DHAVE_GETPWNAM=1 -DHAVE_GETPWUID=1 -DHAVE_GETTIMEOFDAY=1 -DHAVE_GETUID=1 -DHAVE_GETWD=1 -DHAVE_LOCALTIME_R=1 -DHAVE_LSTAT=1 -DHAVE_MEMMOVE=1 -DHAVE_MKDIR=1 -DHAVE_MKFIFO=1 -DHAVE_ON_EXIT=1 -DHAVE_PIPE=1 -DHAVE_POLL=1 -DHAVE_PUTENV=1 -DHAVE_RENAME=1 -DHAVE_RINDEX=1 -DHAVE_RMDIR=1 -DHAVE_SELECT=1 -DHAVE_SETGRENT=1 -DHAVE_SETPWENT=1 -DHAVE_SETVBUF=1 -DHAVE_SIGACTION=1 -DHAVE_SIGPENDING=1 -DHAVE_SIGPROCMASK=1 -DHAVE_SIGSUSPEND=1 -DHAVE_STAT=1 -DHAVE_STRCASECMP=1 -DHAVE_STRDUP=1 -DHAVE_STRERROR=1 -DHAVE_STRFTIME=1 -DHAVE_STRNCASECMP=1 -DHAVE_STRPTIME=1 -DHAVE_TEMPNAM=1 -DHAVE_UMASK=1 -DHAVE_UNLINK=1 -DHAVE_USLEEP=1 -DHAVE_VFPRINTF=1 -DHAVE_VSPRINTF=1 -DHAVE_VSNPRINTF=1 -DHAVE_WAITPID=1 -DHAVE_LIBDL=1 -DHAVE_DLOPEN=1 -DHAVE_DLSYM=1 -DHAVE_DLERROR=1 -DHAVE_DLCLOSE=1 -DWITH_DL=1 -DWITH_DYNAMIC_LINKING=1 -DHAVE_TIMEVAL=1 -DHAVE_FINITE=1 -DHAVE_ISNAN=1 -DHAVE_ISINF=1 -DHAVE_ACOSH=1 -DHAVE_ASINH=1 -DHAVE_ATANH=1 -DHAVE_ERF=1 -DHAVE_ERFC=1 -DHAVE_ST_BLKSIZE=1 -DHAVE_ST_BLOCKS=1 -DHAVE_ST_RDEV=1 -DHAVE_TM_ZONE=1 -DHAVE_GR_PASSWD=1 -DEXCEPTION_IN_MATH=1 -DRETSIGTYPE=void -DSYS_SIGLIST_DECLARED=1 -DHAVE_SYS_SIGLIST=1 -DHAVE_POSIX_SIGNALS=1 -DHAVE_GETRUSAGE=1 -DHAVE_TIMES=1 -DGNUPLOT_HAS_MULTIPLOT=1 -DGNUPLOT_HAS_FRAMES=1 User-preferences (please do not edit this section): -------------------------------------------------- EDITOR = "emacs" EXEC_PATH = ":/usr/local/bin:/bin:/usr/bin:/sbin:/usr/sbin:/usr/local/sbin:/usr/X11R6/bin:/home/ben/bin:/usr/local/matlab5/bin:/usr/local/mathematica/bin:/usr/local/StarOffice/bin" IMAGEPATH = ".:/usr/local/share/octave/2.1.30/imagelib//" INFO_FILE = "/usr/local/info/octave.info" INFO_PROGRAM = "info" LOADPATH = ":/home/ben/octave/Linux:/home/ben/octave/scripts:" PAGER = "less -e -P'-- less ?pB(%pB\\%):--. (f)orward, (b)ack, (q)uit$'" PS1 = "\\s:\\#> " PS2 = "> " PS4 = "+ " automatic_replot = 1 beep_on_error = 0 completion_append_char = " " default_eval_print_flag = 1 # default_global_variable_value = default_return_value = [] default_save_format = "ascii" define_all_return_values = 0 do_fortran_indexing = 0 echo_executing_commands = 0 empty_list_elements_ok = "warn" fixed_point_format = 0 gnuplot_binary = "gnuplot" gnuplot_command_end = "\n" gnuplot_command_plot = "pl" gnuplot_command_replot = "rep" gnuplot_command_splot = "sp" gnuplot_command_title = "t" gnuplot_command_using = "u" gnuplot_command_with = "w" gnuplot_has_frames = 1 gnuplot_has_multiplot = 1 history_file = "/home/ben/.octave_hist" history_size = 1024 ignore_function_time_stamp = "system" implicit_num_to_str_ok = 0 implicit_str_to_num_ok = 0 initialize_global_variables = 0 max_recursion_depth = 256 ok_to_lose_imaginary_part = "warn" output_max_field_width = 10 output_precision = 5 page_output_immediately = 0 page_screen_output = 1 prefer_column_vectors = 1 print_answer_id_name = 1 print_empty_dimensions = 1 print_rhs_assign_val = 0 propagate_empty_matrices = 1 resize_on_range_error = 1 return_last_computed_value = 0 save_precision = 15 saving_history = 1 silent_functions = 0 split_long_rows = 1 string_fill_char = " " struct_levels_to_print = 2 suppress_verbose_help_message = 0 treat_neg_dim_as_zero = 0 warn_assign_as_truth_value = 1 warn_divide_by_zero = 1 warn_function_name_clash = 1 warn_future_time_stamp = 1 warn_missing_semicolon = 0 warn_variable_switch_label = 0 whitespace_in_literal_matrix = "traditional" --Boundary-=_nWlrBbmQBhCDarzOwKkYHIDdqSCD Content-Type: text/english; name="c2d.m.patch" Content-Transfer-Encoding: base64 Content-Disposition: attachment; filename="c2d.m.patch" SW5kZXg6IGMyZC5tCj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09 PT09PT09PT09PT09PT09PT09PT09PT0KUkNTIGZpbGU6IC9jdnMvb2N0YXZlL3NjcmlwdHMvY29u dHJvbC9zeXN0ZW0vYzJkLm0sdgpyZXRyaWV2aW5nIHJldmlzaW9uIDEuMwpkaWZmIC11IC1yMS4z IGMyZC5tCi0tLSBjMmQubQkyMDAwLzAxLzMxIDA3OjQwOjU3CTEuMworKysgYzJkLm0JMjAwMC8w Ni8wNyAxNDozMTozNApAQCAtMzMsNyArMzMsNiBAQAogIyMgdXNlIHRoZSBtYXRyaXggZXhwb25l bnRpYWwgKGRlZmF1bHQpCiAjIyBAaXRlbSAiYmkiCiAjIyB1c2UgdGhlIGJpbGluZWFyIHRyYW5z Zm9ybWF0aW9uCi0jIyBAZW5kIHRhYmxlCiAjIyBAZXhhbXBsZQogIyMgICAgIDIoei0xKQogIyMg cyA9IC0tLS0tCkBAIC00MSw2ICs0MCwxMyBAQAogIyMgQGVuZCBleGFtcGxlCiAjIyBGSVhNRTog VGhpcyBvcHRpb24gZXhpdHMgd2l0aCBhbiBlcnJvciBpZiBAdmFye3N5c30gaXMgbm90IHB1cmVs eQogIyMgY29udGludW91cy4gKFRoZSBAY29kZXtleH0gb3B0aW9uIGNhbiBoYW5kbGUgbWl4ZWQg c3lzdGVtcy4pCisjIyBAaXRlbSAibWF0Y2hlZCIKKyMjIHVzZSB0aGUgbWF0Y2hlZCBwb2xlL3pl cm8gZXF1aXZhbGVudCB0cmFuc2Zvcm1hdGlvbiAKKyMjIEBleGFtcGxlCisjIyBAZW5kIGV4YW1w bGUgICAgICAKKyMjIEZJWE1FOiBJRiB0aGUgc3lzdGVtIGlzIG5vdCBTSVNPIG9yIHB1cmVseSBj b250aW51b3VzIHRoaXMgb3B0aW9uIAorIyMgZXhpdHMgd2l0aCBhbiBlcnJvcgorIyMgQGVuZCB0 YWJsZQogIyMgQGl0ZW0gdAogIyMgc2FtcGxpbmcgdGltZTsgcmVxdWlyZWQgaWYgc3lzIGlzIHB1 cmVseSBjb250aW51b3VzLgogIyMKQEAgLTcxLDcgKzc3LDcgQEAKICMjIENyZWF0ZWQ6IE9jdG9i ZXIgMTk5MwogIyMgVXBkYXRlZCBieSBKb2huIEluZ3JhbSBmb3Igc3lzdGVtIGRhdGEgc3RydWN0 dXJlIEF1Z3VzdCAxOTk2CiAKLWZ1bmN0aW9uIGRzeXMgPSBjMmQgKHN5cywgb3B0LCBUKQorZnVu Y3Rpb24gZHN5cyA9IG15YzJkIChzeXMsIG9wdCwgVCkKICAgIyMgcGFyc2UgaW5wdXQgYXJndW1l bnRzCiAgIGlmKG5hcmdpbiA8IDEgfCBuYXJnaW4gPiAzKQogICAgIHVzYWdlKCJkc3lzPWMyZChz eXNbLFRdKSIpOwpAQCAtOTksNyArMTA1LDggQEAKIAogICBpZiAoIWlzX3NhbXBsZShUKSkKICAg ICBlcnJvcigic2FtcGxpbmcgcGVyaW9kIFQgbXVzdCBiZSBhIHBvc3RpdmUsIHJlYWwgc2NhbGFy Iik7Ci0gIGVsc2VpZiggISAoc3RyY21wKG9wdCwiZXgiKSB8IHN0cmNtcChvcHQsImJpIikgKSAp CisgIGVsc2VpZiggISAoc3RyY21wKG9wdCwiZXgiKSB8IHN0cmNtcChvcHQsImJpIikgfCBcCisJ ICAgICBzdHJjbXAob3B0LCJtYXRjaGVkIikgKSApCiAgICAgZXJyb3IoWyJpbnZhbGlkIG9wdGlv biBwYXNzZWQ6ICIsb3B0XSkKICAgZW5kaWYKIApAQCAtMTY3LDYgKzE3NCwyNCBAQAogICAgICAg RCA9IGQgKyAoYyppYWIpOwogICAgICAgc3RuYW1lZCA9IHN0cmFwcGVuZChzdG5hbWUsIl9kIik7 CiAgICAgICBkc3lzID0gc3Myc3lzKEEsQixDLEQsVCwwLHJvd3MoQSksc3RuYW1lZCxpbm5hbWUs b3V0bmFtZSk7CisgICAgZW5kaWYKKyAgZWxzZWlmKHN0cmNtcChvcHQsIm1hdGNoZWQiKSkKKyAg ICBpZihpc19kaWdpdGFsKHN5cykpCisgICAgICBlcnJvcigiYzJkOiBzeXN0ZW0gaXMgYWxyZWFk eSBkaWdpdGFsIik7CisgICAgZWxzZWlmKChsZW5ndGgoc3lzLmlubmFtZSkgIT0gMSkgfHwgKGxl bmd0aChzeXMub3V0bmFtZSkgIT0gMSkpCisgICAgICBlcnJvcigiYzJkOiBzeXN0ZW0gaW4gbm90 IHNpbmdsZSBpbnB1dCwgc2luZ2xlIG91dHB1dCIpOworICAgIGVsc2UKKyAgICAgIHN5cyA9IHN5 c3VwZGF0ZShzeXMsInpwIik7CisgICAgICBwID0gZXhwKHN5cy5wb2wqVCk7CisgICAgICB6ID0g ZXhwKHN5cy56ZXIqVCk7CisgICAgICBpbmZpbml0ZV96ZXJvcyA9IG1heChzaXplKHN5cy5wb2wp KS1tYXgoc2l6ZShzeXMuemVyKSktMTsKKyAgICAgIGZvciBpID0gMTppbmZpbml0ZV96ZXJvcwor CXogPSBbeiA7IC0xXTsKKyAgICAgIGVuZGZvcgorICAgICAgIyMgU2hvdWxkIHRoZSBmcmVhcXVl bmN5IHdlIGFkanVzdCBhcm91bmQgYWx3YXlzIGJlIDE/ICAgCisgICAgICBbY21hZyxjcGhhc2Us Y3ddID0gYm9kZShzeXMsMSk7CisgICAgICBbZG1hZyxkcGFoc2UsZHddID0gYm9kZSh6cDJzeXMo eixwLDEsVCksMSk7CisgICAgICBkc3lzID0genAyc3lzKHoscCxjbWFnL2RtYWcsVCk7CiAgICAg ZW5kaWYKICAgZWxzZQogICAgIGVycm9yKFsiQmFkIG9wdGlvbj0iLG9wdF0pCg== --Boundary-=_nWlrBbmQBhCDarzOwKkYHIDdqSCD-- ----------------------------------------------------------------------- Octave is freely available under the terms of the GNU GPL. Octave's home on the web: http://www.che.wisc.edu/octave/octave.html How to fund new projects: http://www.che.wisc.edu/octave/funding.html Subscription information: http://www.che.wisc.edu/octave/archive.html -----------------------------------------------------------------------