home *** CD-ROM | disk | FTP | other *** search
/ Piper's Pit BBS/FTP: ibm 0020 - 0029 / ibm0020-0029 / ibm0028.tar / ibm0028 / MCAD254.ZIP / BODE.MCD < prev    next >
Encoding:
Text File  |  1990-02-01  |  2.2 KB  |  91 lines
  1. .MCD 25000 0
  2. .CMD SURFACEFORMAT rot=10 tilt=35 vScale=20 size=15,30
  3. .CMD SKETCHFORMAT mag=1.000000,1.000000 center=0.500000,0.500000 size=15,30 box=y
  4. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  5. .CMD FORMAT  rd=d ct=56 im=i et=3 zt=15 pr=3 mass length time charge
  6. .CMD SET ORIGIN 0
  7. .CMD SET TOL 0.001000
  8. .CMD MARGIN 0
  9. .CMD LINELENGTH 78
  10. .CMD SET PRNCOLWIDTH 8
  11. .CMD SET PRNPRECISION 4
  12. .TXT 1 0 1 13 
  13. a1,12,10,11
  14. BODE PLOTS
  15. .TXT 2 10 4 64 
  16. a4,63,61,214
  17. This document examines the dynamic response of a second-order 
  18. system, and presents it in the form of a Bode plot. The 
  19. results demonstrate how to use the Laplace domain and 
  20. semi-log and log-log plots in MathCAD.
  21. .TXT 5 -10 1 10 
  22. a1,10,8,9
  23. Example:
  24. .TXT 2 10 2 52 
  25. a2,51,54,66
  26. First, define resonant frequency, ▀, and damping 
  27. coefficient, α:
  28. .EQN 3 0 4 9 
  29. α:\2/2
  30. .EQN 1 20 1 15 
  31. ▀:10*hertz
  32. .TXT 4 -20 1 41 
  33. a1,40,38,39
  34. The second order system is defined as:
  35. .EQN 2 0 5 27 
  36. f(s):▀^2/(s^2+2*α*▀*s+▀^2)
  37. .EQN 7 0 1 9 
  38. n:50
  39. .EQN 0 20 1 12 
  40. i:0;n
  41. .TXT 2 -20 3 61 
  42. a3,60,58,158
  43. We will be working with a system with resonant frequency   
  44. 10 rad/sec, so the minimum and maxumum frequency of 
  45. interest for plotting are 1 and 100 rad/sec.
  46. .EQN 4 0 1 10 
  47. min:1
  48. .EQN 0 20 1 12 
  49. max:100
  50. .TXT 2 -20 3 58 
  51. a3,58,56,124
  52. r is a log ratio of maximum and minimum frequencies. The 
  53. values of s computed below will be equally spaced on a 
  54. log plot.
  55. .EQN 4 0 3 14 
  56. r:ln(max/min)
  57. .EQN 4 0 5 22 
  58. s[i:min*e^(i*r/(n))*hertz
  59. .TXT 6 -10 1 20 
  60. a1,20,18,19
  61. Global definitions
  62. .EQN 2 10 1 10 
  63. sec~1t
  64. .EQN 0 20 3 13 
  65. hertz~(2*π)/sec
  66. .EQN 0 20 3 11 
  67. deg~π/180
  68. .EQN 4 -40 1 37 
  69. Φ(Θ)~angle(Re(Θ),Im(Θ))-360*deg
  70. .TXT 3 -10 1 1 
  71. x1,1,0,0
  72. .TXT 2 30 1 35 
  73. a1,34,32,33
  74. Bode Plot of Second Order System
  75. .EQN 2 -25 17 73 
  76. 10&10^-2&|f(s[i*1j){67,66,15,60,O}@10^2*hertz&1*hertz&s[i
  77. .TXT 18 20 1 36 
  78. a1,36,34,35
  79. Phase angle of second order system
  80. .EQN 2 -20 15 74 
  81. 0*deg&-180*deg&Φ(f(s[i*1j)){2,194,13,60,l}@10^2*hertz&10^0*hertz&s[i
  82. .EQN 17 -5 2 19 
  83. Θ[i:Φ(f(s[i*1j))
  84. .TXT 0 23 1 34 
  85. a1,33,31,32
  86. Nyquist Plot of system function
  87. .EQN 3 -23 2 20 
  88. Mag[i:|f(s[i*1j)
  89. .EQN 0 28 14 51 
  90. 1&-1&-sin(Θ[i)*Mag[i{2,2,12,35,l}@1&-1&cos(Θ[i)*Mag[i
  91.