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

文件历史

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

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

以下2个页面使用本文件:

全域文件用途

以下其他wiki使用此文件: