2018-08-18

第二章

连续时间系统的时域分析 (上)

  • 1、建立和求解线性微分方程的过程
    • 建立数学模型
      • KCL
      • KVL
      • 线性非时变系统的微分方程
        • \frac{d^n}{dt^n}r(t) + a_{n-1}\frac{d^{n-1}}{dt^{n-1}}r(t) + ...+a_1\frac{d}{dt}r(t)+a_0r(t)
        • = b_m\frac{d^m}{dt^m}e(t) +b_{m-1}\frac{d^{m-1}}{dt^{m-1}}e(t)+...+b_1\frac{d}{dt}e(t) +b_0e(t)
        • 初始条件:r^{n-1}(0),r^{n-2}(0),...,r(0)是几阶微分方程就需要几个初始条件
    • 求解
      • 1、时域法
        • 将响应分为通解和特解
        • e.g. r''(t) + 3r'(t) - 4r(t) = e(t)
        • e(t) = e^{-2t}
        • r(0) = 1,r'(0) = 1
        • \lambda^2 + 3\lambda - 4 = 0特征方程
        • \lambda_1 = 1,\lambda_2 = -4特征根
        • r(t) = C_1e^t + C_2e^{-4t}通解
        • r(t) = A e^{-2t}特解(猜测)
        • r(t) = C_1e^t + C_2e^{-4t}+ A e^{-2t}
        • C_1 + C_2 +A = 1
        • C_1 + -4C_2 - 2A = 1
        • (C_1e^t + C_2e^{-4t} + Ae^{-2t})''+3(C_1e^t + C_2e^{-4t} + Ae^{-2t})'-4(C_1e^t + C_2e^{-4t} + Ae^{-2t}) = e^{-2t}
      • 通解 - 自由响应
      • 特解 - 受迫响应
      • 2、卷积法(算子法)
        • 将响应分为两部分
          • 零输入响应:系统在没有输入激励的情况下,仅仅由系统的初始状态引起的响应r_{zi}(t)
          • 零状态响应:(没有初始储能)的条件下,仅仅由输入信号引起的响应r_{zs}(t)
          • 方程1:零状态响应:(经典发)
          • \frac{d^n}{dt^n}r_1(t) + a_{n-1}\frac{d^{n-1}}{dt^{n-1}}r_1(t) + ...+a_1\frac{d}{dt}r_1(t)+a_0r(t)
          • = b_m\frac{d^m}{dt^m}e(t) +b_{m-1}\frac{d^{m-1}}{dt^{m-1}}e(t)+...+b_1\frac{d}{dt}e(t) +b_0e(t)
          • 初始条件:r^{n-1}_1(0) = r^{n-2}_1(0) = 0,...,r_1(0) = 0
          • 方程2:零状态响应:(卷积积分法,激励信号是一个有始信号)
          • \frac{d^n}{dt^n}r_2(t) + a_{n-1}\frac{d^{n-1}}{dt^{n-1}}r_2(t) + ...+a_1\frac{d}{dt}r_2(t)+a_0r(t) = 0
          • 初始条件:r^{n-1}_2(0) =r^{n-1}(0) ,r^{n-2}_2(0) =r^{n-2}(0),...,r_2(0) = r(0)
          • r(t) = r_1(t) + r_2(t)方程的解
  • 2、微分算子
    • p = \frac{d}{dt},p^n = \frac{d^n}{dt^n},\frac{1}{p} = \int_{-\infty}^{t}()dt变限积分
    • 电感电容的伏安特性
      • U_L = L \frac{d}{dt}i_L(t) = L\cdot P\cdot i_L
      • U_C = \frac{1}{C}\int_{-\infty}^{t}i_C(t)dt = \frac{1}{C\cdot P}\cdot i_C
      • 即可以将电感和电容记成阻值为L\cdot P,\frac{1}{C\cdot P}的电阻
    • 利用算子可以将微分方程表示为:
    • p^nr(t) + a_{n-1}p^{n-1}r(t) + ...+a_1pr(t)+a_0r(t)
    • = b_mp^me(t) +b_{m-1}p^{m-1}e(t)+...+b_1pe(t) +b_0e(t)
    • 提取公因子,可以将上式化简为:
    • (p^n +a_{n-1}p^{n-1}+...+a_1p + a_0)r(t) =
    • (b_mp^m + b_{m-1}p^{m-1} + ...+b_1p+b_0)e(t)
    • r(t) = \frac{(b_mp^m + b_{m-1}p^{m-1} + ...+b_1p+b_0)}{(p^n +a_{n-1}p^{n-1}+...+a_1p + a_0)}e(t)
    • H(p) =\frac{(b_mp^m + b_{m-1}p^{m-1} + ...+b_1p+b_0)}{(p^n +a_{n-1}p^{n-1}+...+a_1p + a_0)} = \frac{N(p)}{D(p)}
    • r(t) = H(p) e(t)
    • 微分方程的一种简单记法,并不代表可以这样计算
    • 算子运算法则:
      • mp + np = (m+n)p,m,n为任意整数
      • p^mp^n = p^{m+n},m,n同为任意正整数(负整数)
      • p\frac{1}{p} = 1
      • 注:1、\frac{1}{p}p不一定等于1-微分和积分的次序不能交换
      • 2、如果px(t) = py(t),不一定能够推出x(t) = y(t),只能得到x(t) = y(t) + C
  • 3、零输入响应求解
    • D(p)r(t) =(p^n +a_{n-1}p^{n-1}+...+a_1p + a_0)r(t) = 0
      • 经典解法
        • 1、确定系统的自然频率,令D(p) = 0,将p看成一个代数量,解得其n个特征根\lambda_1,\lambda_2,...,\lambda_n
        • 2、确定零输入响应的形式解
          • 1)、如果没有重根
            • r_{zi}(t) = C_1e^{\lambda_1t} + C_2e^{\lambda_2t} +...+ C_ne^{\lambda_nt} = \sum_{i = 1}^n C_ie^{\lambda_it}
            • 其中C_1.C_2,...,C_n为待定常数
          • 2)、如果有重根,假设\lambda_1是一个k重根,即\lambda_1 = \lambda_2 = ... = \lambda_k则形式解为:
          • r_{zi}(t) = C_1e^{\lambda_1t} +C_2te^{\lambda_2t} +C_3t^2e^{\lambda_1t}+...+C_k t^{k-1} e^{\lambda_1t} + C_{k + 1}e^{k-1}t + ...+C_ne^{\lambda_nt}
          • = \sum_{i = 1}^{k}C_it^{i - 1}e^{\lambda_1t} + \sum_{i = k +1 }^{n}C_ie^{\lambda_it}
        • 3、根据初始条件,确定待定系数.一般的初始条件为已知零时刻的响应及各阶导数r(0),r'(0),r''(0),...,r^{n-1}(0),带入形式解中就可以确定待定系数.
          • r(0) = C_1 +C_2 +...+C_n
          • r'(0) = \lambda_1C_1 +\lambda_2C_2 +...+\lambda_nC_n
          • r''(0) = \lambda_1^{2}C_1 +\lambda_2^{2}C_2 +...+\lambda_n^{n}C_n
          • ......
          • r^{n-1}(0) =\lambda_1^{n-1}C_1 +\lambda_2^{n-1}C_2 +...+\lambda_{n}^{n-1}C_n
          • 或者记为矩阵的形式
          • \begin{bmatrix} r(0)\\ r'(0)\\r''(0)\\...\\r^{(n-1)}(0) \end{bmatrix} = \begin{bmatrix}1 & 1 & 1 &...& 1 \\ \lambda_1 & \lambda_2 & \lambda_3 &...& \lambda_n\\ \lambda_1^2 & \lambda_2^2 & \lambda_3^2 &...& \lambda_n^2\\ ... & ...&...&...&...\\ \lambda_n^{n-1} & \lambda_n^{n-1} & \lambda_n^{n-1} &...& \lambda_n^{n-1} \end{bmatrix} \begin{bmatrix} C_1\\ C_2\\C_3\\...\\C_n \end{bmatrix}
  • 4、零状态响应(利用线性系统的齐次和叠加原理)
    • 将信号分解为一系列"标准统一”的子信号之和(或积分)
    • 求线性系统对各个子信号的响应
    • 将各个子信号的响应相叠加,从而得到系统对激励信号的响应
  • 5、奇异函数
    • 找子信号
      • 完备性:任意函数都可以分解为该子信号的和
      • 简单性:容易求得子信号的响应
      • 相似性:不同子信号的响应具有内在联系
    • 奇异函数:阶跃函数,冲击函数
      • 阶跃函数\varepsilon(t) = \begin{cases}1,t\geq0 \\0,{\text{其他}} \end{cases}
        • 任意函数乘以\varepsilon(t),t<0部分等于零,成为有始函数,u(t)表示也可以
      • 冲激函数\delta(t)方波或者三角波的极限
        • 1)\delta(t) = \frac{d}{dt}\varepsilon(t)
          • 1、只在t = 0处非零值,其他各处都为零
          • 2、\delta(t),\varepsilon(t)互为微分和积分\varepsilon(t) = \int_{-\infty}^{t}\delta(t)dt
          • 3、\int_{-\infty}^{+\infty}\delta(t)dt = 1
          • 4、\delta(t) = \delta(-t)
          • 5、 \delta(t)f(t) = \delta(t)f(0)
          • 或者 \delta(t-t_0)f(t) = \delta(t-t_0)f(t_0)
          • 6、\int_{-\infty}^{+\infty}f(t)\delta(t)dt = f(0)
          • \int_{-\infty}^{+\infty}f(t-t_0)\delta(t)dt = f(t_0)冲击函数的取样特性
          • 7、\delta(at) = \frac{1}{|a|}\delta(t)
          • 冲击函数的图形表示方法:位置和强度
        • 2)定义分配函数
          • 通过函数对另外一个函数的作用来定义这个函数,利用冲激函数的抽样特性作为冲激函数的定义,即对于任意的函数f(t),都满足:
            \int_{-\infty}^{+\infty}f(t)\delta(t-t_0)dt = f(t_0)
  • 5、信号的时域分解(分割求近似,求和取极限
    • 1、将任意信号近似表示为多个阶跃函数之和,各个子信号为:
      • f_0(t) = f(0)\varepsilon(t)
      • f_1(t) = [f(\Delta t) - f(0)]\varepsilon(t-\Delta t)
        • = \frac{f(\Delta t)- f(0)}{\Delta t}\varepsilon(t - \Delta t)\Delta t
      • f_k(t) = [f(k\Delta t) - f((k - 1)\Delta t)]\varepsilon(t-k\Delta t)
        • \frac{f(k\Delta t) - f((k - 1)\Delta t)}{\Delta t}\varepsilon(t-k\Delta t)\Delta t
      • 总和为:
      • f(t)\approx f_0(t) + f_1(t) +... + f_k(t) +...
        • = f(0)\varepsilon(t) + \frac{f(\Delta t)- f(0)}{\Delta t}\varepsilon(t - \Delta t)\Delta t + \frac{f(2\Delta t)- f(\Delta t)}{\Delta t}\varepsilon(t - 2\Delta t)\Delta t + ... + \frac{f(k\Delta t)- f((k-1)\Delta t}{\Delta t}\varepsilon(t - k\Delta t)\Delta t + ...
        • \Delta t\to d\tau(或0)则
        • \frac{f(k\Delta t)- f((k-1)\Delta t)}{\Delta t}\varepsilon(t - k\Delta t)\Delta t = f'(k\Delta t)\varepsilon (t - k\Delta t)
      • f(t) = f(0)\varepsilon (t) + \int_0^tf'(\tau)\varepsilon(t - \tau)d\tau
      • 如果f(t)t = 0处连续可导,则:
        • f(t) = \int_{0}^{t}f'(\tau)\varepsilon(t-\tau)d\tau
    • 2、将任意信号近似表示为多个冲激函数之和,各个子信号为:
      • 将任意函数近似表示为一系列矩形脉冲函数之和:
        • u_{\tau}(t) = \begin{cases}1,0<t<\tau \\0,{\text{其他}} \end{cases}
          • f_0(t) = f(0)u_{\Delta t}(t)
          • f_1(t)= f(\Delta t)u_{\Delta t}(t - \Delta t)
          • f_2(t) = f(2\Delta t)u_{\Delta t}(t - 2\Delta t)
          • f_k(t) =f(k\Delta t)u_{\Delta t}(t - k\Delta t)
        • f(t) = f_0(t) + f_1(t)+...+f_k(t) + ...
          • = f(0)u_{\Delta t}(t)+ f(\Delta t)u_{\Delta t}(t - \Delta t)+...+f(k\Delta t)u_{\Delta t}(t - k\Delta t)
          • = f(0)\frac{u_{\Delta t}(t)}{\Delta t}\Delta t + f(\Delta t)\frac{u_{\Delta t}(t-\Delta t)}{\Delta t}\Delta t +...+f(k\Delta t)\frac{u_{\Delta t}(t-k\Delta t)}{\Delta t}\Delta t
          • = \sum_{k =0}^\infty f(k\Delta t)\frac{u_{\Delta t}(t-k\Delta t)}{\Delta t}\Delta t
            • \Delta t\to d\tau(或0)
            • f(t) = \int_{0}^{\infty}f(\tau)\delta(t - \tau)d\tau
            • \frac{u_{\Delta t}(t-k\Delta t)}{\Delta t}\to \delta(t - k\Delta t)方波脉冲的定义解释(面积为1的方波脉冲)
        • f(t) = \int_{-\infty}^{+\infty}f(\tau)\delta(t - \tau)d\tau更具有普遍性
        • f(t) = \int_{0}^{t}f(\tau)\delta(t - \tau)d\tau
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 218,284评论 6 506
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,115评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,614评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,671评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,699评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,562评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,309评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,223评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,668评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,859评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,981评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,705评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,310评论 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,904评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,023评论 1 270
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,146评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,933评论 2 355

推荐阅读更多精彩内容