File:Second order transfer function.svg

原始文件 (SVG文件,尺寸为631 × 356像素,文件大小:80 KB)


摘要

描述
English: Step responses for a second order system defined by the transfer function:

where is the damping ratio and is the undamped natural frequency. The equations were obtained from here, plotted using maxima and edited in a text editor to insert the Greek alphabets in the plot. The equations are:

日期
来源 自己的作品
作者 Krishnavedala
Source code using w:python (programming language) with numpy and matplotlib toolboxes
from matplotlib.pyplot import *
from numpy import *

wt = linspace(0,15,100)
b = lambda z: sqrt(1. - z**2)
t = lambda z: arctan(b(z)/z)
h1 = lambda wt,z: 1. - exp(-z*wt)*sin(b(z)*wt+t(z))/b(z)
h2 = lambda wt: 1. - cos(wt)
h3 = lambda wt: 1. - exp(-wt)*(1.+wt)
s1 = lambda z: (z + sqrt(z**2-1.))
s2 = lambda z: (z - sqrt(z**2-1.))
h4 = lambda wt,z: 1. + ( (exp(-s1(z)*wt)/s1(z)) - \
	(exp(-s2(z)*wt)/s2(z)) ) / (2.*sqrt(z**2-1.))

fig = figure(figsize=(8,4))
ax = fig.add_subplot(111)
ax.grid(True)
ax.plot(wt,h2(wt),'g',label=r"undamped $(\zeta=0)$")
ax.plot(wt,h1(wt,.5),'b',label=r"under $(\zeta=0.5)$")
ax.plot(wt,h3(wt),'r',label=r"critical $(\zeta=1.0)$")
ax.plot(wt,h4(wt,1.5),'m',label=r"over $(\zeta=1.5)$")
ax.set_ylim(0,2)
ax.minorticks_on()

leg = ax.legend(frameon=False,handletextpad=.05)
setp(leg.get_texts(),fontsize=10)
ax.set_xlim(0,15)
ax.set_xlabel(r"$\omega t$",fontsize=15)
ax.set_ylabel("Step response",fontsize=12)
fig.savefig("Second_order_transfer_function.svg",bbox_inches="tight",\
	pad_inches=.15)
The maxima source code
beta(zeta) := sqrt(1-zeta^2);
theta(zeta) := atan(beta(zeta)/zeta);
h_under(wt) := 1 - beta(0.5)^-1*exp(-0.5*wt)*sin(wt*beta(.5)+theta(0.5));
h_un(wt) := 1 - cos(wt);
h_crit(wt) := 1 - exp(-wt) * (1+wt);
s1(zeta) := zeta+sqrt(zeta^2-1);
s2(zeta) := zeta-sqrt(zeta^2-1);
h_over(wt) := 1 + ((exp(-s1(1.5)*wt)/s1(1.5))-(exp(-s2(1.5)*wt)/s2(1.5)))/(2*sqrt(1.5^2-1));
load(draw);
draw2d(dimensions=[800,400],terminal=svg,
  user_preamble="set mxtics; set mytics;",
  grid=true, yrange=[0,2], xlabel="omega t",
  line_width=1.5, ylabel="Step response",
  key="under (zeta=0.5)",color=blue,explicit(h_under(wt),wt,0,15),
  key="critical (zeta=1)",color=red,explicit(h_crit(wt),wt,0,15),
  key="over (zeta=1.5)",color=magenta,explicit(h_over(wt),wt,0,15),
  key="undamped (zeta=0)",color=green,explicit(h_un(wt),wt,0,15)
 );

许可协议

我,本作品著作权人,特此采用以下许可协议发表本作品:
w:zh:知识共享
署名 相同方式共享
您可以自由地:
  • 共享 – 复制、发行并传播本作品
  • 修改 – 改编作品
惟须遵守下列条件:
  • 署名 – 您必须对作品进行署名,提供授权条款的链接,并说明是否对原始内容进行了更改。您可以用任何合理的方式来署名,但不得以任何方式表明许可人认可您或您的使用。
  • 相同方式共享 – 如果您再混合、转换或者基于本作品进行创作,您必须以与原先许可协议相同或相兼容的许可协议分发您贡献的作品。

说明

添加一行文字以描述该文件所表现的内容

此文件中描述的项目

描绘内容

创作作者 简体中文(已转写)

某些值没有维基数据项目

作者姓名字符串 简体中文(已转写):​Krishnavedala
维基媒体用户名 简体中文(已转写):​Krishnavedala

版权状态 简体中文(已转写)

版权所有 简体中文(已转写)

文件来源 简体中文(已转写)

上传者的原创作品 简体中文(已转写)

文件历史

点击某个日期/时间查看对应时刻的文件。

日期/时间缩⁠略⁠图大小用户备注
当前2011年6月7日 (二) 15:572011年6月7日 (二) 15:57版本的缩略图631 × 356(80 KB)Krishnavedalafigure correction after correcting the equations
2011年5月26日 (四) 19:252011年5月26日 (四) 19:25版本的缩略图800 × 400(36 KB)Krishnavedalaerroneous under-damped curve corrected.
2011年5月26日 (四) 18:282011年5月26日 (四) 18:28版本的缩略图800 × 400(36 KB)Krishnavedala

以下页面使用本文件:

全域文件用途

以下其他wiki使用此文件: