拉普拉斯元件实现声阻运算（pspice) - 基础理论室 - 声学楼论坛

声学楼论坛 → 技术交流区 → 基础理论室 → 拉普拉斯元件实现声阻运算（pspice)

共有26540人关注过本帖树形打印复制链接主题：拉普拉斯元件实现声阻运算（pspice)

liheng_rd

小大 1楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

微型扬声器
等级：职业侠客帖子：121 积分：1927 威望：0 精华：0 注册：2007-5-31 19:07:18

拉普拉斯元件实现声阻运算（pspice) 发帖心情

Post By：2013-5-7 0:02:17 [只看该作者]

如题，

拉普拉斯元件实现声阻运算，参数如何设置

nisa

小大 2楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

无名小卒
等级：职业侠客帖子：114 积分：2635 威望：0 精华：0 注册：2006-11-12 0:48:14

Post By：2013-5-7 11:29:59 [只看该作者]

Modeling a frequency dependent inductor for AC sweep analysis

PROBLEM:
How do I model a frequency dependent inductor for AC sweep analysis?

SOLUTION:
The impedance of an inductor is:

Z = L * s

where

s = j * 2 * PI * Frequency

In a frequency dependent inductor, the "L" in the above formula changes with frequency.

If you know L as a math function of frequency, then:

Use a GLAPLACE device.
Connect the positive and negative inputs to the positive and negative outputs, respectively.
Write the transform expression as 1/(<L>*s), where <L> is replaced by the desired math formula, with frequency replaced by abs(s)/6.283185.

For example, if L = 1 / sqrt(Frequency), then write the transform expression as:

1/(1 / sqrt(abs(s)/6.283185)*s) or sqrt(abs(s)/6.283185)/s

http://www.orcad.com/documents/community.faqs/pspice/020020.aspx

更详细的介绍查阅实例(要翻越和谐之墙): Google搜索"Modeling a frequency dependent inductor for AC sweep analysis"的第二个结果:

http://acousticians.blogspot.com/2012/10/modeling-frequency-dependent-inductor.html

[此贴子已经被作者于2013-05-07 20:03:13编辑过]

jma

小大 3楼 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

小马哥
等级：版主帖子：3126 积分：27130 威望：2 精华：8 注册：2006-1-19 20:29:17

Post By：2013-5-7 20:09:03 [只看该作者]

1) 将GLAPLACE的输入端正负极与输出端正负极分别对接起来;

2) 改写转换函数为1/<R>形式, 其中<R>为计算后得到的频率f的函数;

3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28 .

独学而无友，则孤陋寡闻！！！
E-Mail:jma168@163.com 广告勿扰！

liheng_rd

小大 4楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

微型扬声器
等级：职业侠客帖子：121 积分：1927 威望：0 精华：0 注册：2007-5-31 19:07:18

Post By：2013-5-7 22:30:39 [只看该作者]

非常感谢楼上二位的帮助，谢谢

liheng_rd

小大 5楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

微型扬声器
等级：职业侠客帖子：121 积分：1927 威望：0 精华：0 注册：2007-5-31 19:07:18

Post By：2013-5-8 0:30:01 [只看该作者]

举个实例吧：

单孔的声阻：R=sqrt(2ηρw)*(l/a+2)/(π*a^2)

假设:η=1/2;ρ=1/2;l=1；a=1;

则计算的R=[sqrt(π*frequency)*3]/π

那么XFORM的参数是1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)这个表达式吗？

jeno

小大 6楼 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游民帖子：10 积分：247 威望：0 精华：0 注册：2006-8-15 11:50:23

Post By：2013-5-9 9:37:09 [只看该作者]

学习啦

nisa

小大 7楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

无名小卒
等级：职业侠客帖子：114 积分：2635 威望：0 精华：0 注册：2006-11-12 0:48:14

Post By：2013-5-9 11:32:50 [只看该作者]

以下是引用liheng_rd在2013-05-08 00:30:01的发言：

举个实例吧：

单孔的声阻：R=sqrt(2ηρw)*(l/a+2)/(π*a^2)

假设:η=1/2;ρ=1/2;l=1；a=1;

则计算的R=[sqrt(π*frequency)*3]/π

那么XFORM的参数是1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)这个表达式吗？

要把分母最后的s去掉, 即XFORM= 1/Z =1/R= 1/([sqrt(pi*abs(s)/6.28)*3]/pi)

假设单孔的声质量表达式子L=sqrt(π*frequency)*3]/π ,

那么XFORM的参数才是XFORM= 1/Z =1/(L*s)=1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)

1/Z中的Z表示阻抗

lhyjs_63

小大 8楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游侠帖子：20 积分：1152 威望：0 精华：0 注册：2006-9-24 15:53:19

[讨论]

Post By：2013-6-28 15:33:50 [只看该作者]

小马哥，nisa您们好：

第3点中，3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28

有不明白请指教，s = jω = j·2πf，那么 f = s / j·2π, 又因为 j = √-1 所以 f = s / (√-1 ·2π) = -s·√-1 / 2π, 如何得到 f = abs(s)/6.28 的呢，谢谢。

nisa

小大 9楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

无名小卒
等级：职业侠客帖子：114 积分：2635 威望：0 精华：0 注册：2006-11-12 0:48:14

Post By：2013-6-29 7:54:44 [只看该作者]

以下是引用lhyjs_63在2013-06-28 15:33:50的发言：
小马哥，nisa您们好：

第3点中，3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28

有不明白请指教，s = jω = j·2πf，那么 f = s / j·2π, 又因为 j = √-1 所以 f = s / (√-1 ·2π) = -s·√-1 / 2π, 如何得到 f = abs(s)/6.28 的呢，谢谢。

你想多了; abs（）是求模值，即abs(s)=abs(jω)=ω, ω=2πf → f=ω/2π;

mustang

小大 10楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游民帖子：15 积分：231 威望：0 精华：0 注册：2012-11-16 19:57:09

Post By：2013-7-5 19:18:37 [只看该作者]

学习

jqr01

小大 11楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：侠之大者帖子：194 积分：1483 威望：0 精华：0 注册：2008-10-17 22:51:47

Post By：2013-7-8 2:41:17 [只看该作者]

阻抗Z= 单孔的声质量 L *S 这里的S是代表什么呢,为什么要把分母最后的s去掉 , 即XFORM= 1/Z =1/R= 1/([sqrt(pi*abs(s)/6.28)*3]/pi)

jqr01

小大 12楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：侠之大者帖子：194 积分：1483 威望：0 精华：0 注册：2008-10-17 22:51:47

Post By：2013-7-8 2:42:29 [只看该作者]

为什么这里又会出现S:那么XFORM的参数才是XFORM= 1/Z =1/(L*s)=1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)

1/Z中的Z表示阻抗

nisa

小大 13楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

无名小卒
等级：职业侠客帖子：114 积分：2635 威望：0 精华：0 注册：2006-11-12 0:48:14

Post By：2013-7-8 20:29:24 [只看该作者]

二、三楼的内容说得够详细了啊；自己多动脑子思考才好理解透。

katlin

小大 14楼 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游民帖子：10 积分：149 威望：0 精华：0 注册：2014-11-10 10:38:06

Post By：2015-5-6 10:07:44 [只看该作者]

为什么我将参数XFORM设置成 1/([sqrt(pi*abs(s)/6.28)*3]/pi)，提示出错呢，说是无效的传输函数？？

**** INCLUDING ceshi1-SCHEMATIC1.net ****

* source CESHI1

C_Cmd N16341 N11893 4.5m

E_pSd N03344 0 V2 N06338 3.8969e-5

L_Ma1f OUT 0 163

X_vSd N11893 N03344 N06338 V2 SCHEMATIC1_vSd

R_Rmd V1 N16144 4.0m

C_C2 N06242 0 0.094p

R_Re N02561 N02595 31.2

C_C1 N03778 0 0.093p

X_BLi N16272 N03050 V1 N03236 SCHEMATIC1_BLi

R_Ra11 OUT N04192 10

X_BLv 0 N03236 N03050 0 SCHEMATIC1_BLv

L_Le N02595 N16272 80uH

R_Ra1f 0 OUT 190M

L_Mmd N16144 N16341 19uH

R_Ra1r N06242 N06138 246M

G_G1 N03778 N03852 LAPLACE { V(N03778, N03852) } {

+ 1/[sqrt(abs(s)/6.28)*0.09] }

---------------------------------------------------------$

ERROR -- Invalid transfer function

L_Maf V2 N03778 38.6

R_Ra31 N06138 N05374 10

L_Mar N06338 N06242 38.6

G_G2 0 N05072 LAPLACE { V(0, N05072) } { 1/[sqrt(abs(s)/6.28)*0.09] }

-----------------------------------------------$

ERROR -- Invalid transfer function

R_R20 N03778 0 100G

L_Ma1 N03852 N04192 210

R_R21 N06242 0 100G

L_Ma1r N06138 N06242 186

V_Eg N02561 0 DC 0Vdc AC 1.79

L_Ma3 N05072 N05374 455

Acoustics

小大 15楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

电声工程师
等级：版主帖子：3520 积分：26241 威望：2 精华：19 注册：2006-2-13 13:44:44

Post By：2015-5-6 15:25:57 [只看该作者]

* Example Schematics Netlist *

R_R2 $N_0002 $N_0001 0.001

V_V3 $N_0002 0 DC 0V AC 1v

R_R3 $N_0002 0 10meg

R_Rsim 0 $N_0003 10meg

G_Rvar2 $N_0001 0 LAPLACE { V($N_0001, 0) } {
+ 1/(k2*sqrt(S/(2*3.14159265))) }

R_R1 $N_0004 0 16000

G_Lvar $N_0003 0 LAPLACE { V($N_0003, 0) } {
+ 1/(S*1e-6*(k3-(k4*LOG(k5*(S/(2*3.14159265)))))) }

G_Rvar1 $N_0001 $N_0003 LAPLACE { V($N_0001, $N_0003) } {
+ 1/(k1*sqrt(S/(2*3.14159265))) }

C_C $N_0001 $N_0004 .64pF

To view the specific effective simulation results

Include the following Macros in your Probe trace analysis to see frequency vs. inductance, impedance, phase angle (in degrees), and Q factor:

PI = 3.14159265

L=(IMG(V(R3:1)/I(R2))/(2*pi*FREQUENCY)

Z = V(R3:1)/I(R2)

ANG = (180/PI)*ARCTAN((IMG(V(R3:1)/I(R2)))/(R(V(R3:1)/I(R2))))

QFACT = ABS((IMG(V(R3:1)/I(R2)))/(R(V(R3:1)/I(R2))))

由Wei Nisa于 24th October 2012发布

学无止境！

vickyzhang1981

小大 16楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：论坛游侠帖子：35 积分：247 威望：0 精华：0 注册：2016-12-5 11:01:42

Post By：2017-1-20 16:48:52 [只看该作者]

我的权限也不够，也来狂顶楼主一通图片点击可在新窗口打开查看

XYO104249

小大 17楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：业余侠客帖子：78 积分：486 威望：0 精华：0 注册：2012-10-11 9:52:28

Post By：2017-12-25 14:56:47 [只看该作者]

太深奥了~~收藏学习

eie05031

小大 18楼 | 信息 | 搜索 | 邮箱 | 主页 | UC

加好友

发短信

等级：侠之大者帖子：172 积分：938 威望：0 精华：0 注册：2018-9-28 14:50:33

Post By：2023-5-17 13:12:40 [只看该作者]

么的深奥啊，闷啦

下一主题 >> << 上一主题

主题：拉普拉斯元件实现声阻运算（pspice)