利用bode图求传递函数例题
利用bode图求传递函数例题
第一篇：利用bode图求传递函数例题
例题：已知最小相位系统开环对数频率特性曲线如图所示。试写出开环传递函数Gk(s)。
解：
1)ω<ω1的低频段斜率为[-20]，故低频段为K/ｓ。
ω增至ω1，斜率由[-20]转为[-40]，增加[-20]，所以ω1应为惯性环节的转折频率，该环节为11。
11s 1ω增至ω2，斜率由[–40]转为[–20]，增加[+20]，所以ω2应为一阶微分环节的转折频率，该环节为 2s 1。
11ω增到ω3，斜率由[-20]转为[-40]，该环节为,ω>ω3，斜率保持不变。
31s 1故系统开环传递函数应由上述各典型环节串联组成，即
K(Gk(s) s(2)确定开环增益K 当ω=ωc时，A(ωc)=1。
2s 1)1
1 1s 1)( 3s 1)K(所以 A( c) 1 2 c)2 11 K1 21 c 1
c(1 1 c)2 ( 3 c)2 1 c 1 c故 K  2 c 所以，Gk(s) 11 1s(s 1)(s 1) 1 3
2 c1(s 1) 1 2
练习：
最小相位系统的对数幅频特性如下图所示，试分别确定各系统的传递函数。
(a)
（c）
a：G(s) 10s(s 1)
b：G(s) 100(10s 1)(s 1)
rank函数的用法cG(s) 100(0.5s 1)(0.2s 1)
(b)
第二篇：利用MATLAB求线性方程组
《MATLAB语言》课成论文
利用MATLAB求线性方程组
姓名：***
学号：***********
专业：通信工程
班级：2010级通信工程一班
指导老师：***
学院：物电学院
完成日期：2011年12月17日
利用MATLAB求解线性方程组
（郭亚兰 12010245331 2010 级通信一班）
【摘要】在高等数学及线性代数中涉及许多的数值问题，未知数的求解，微积分，不定积分，线性方程组的求解等对其手工求解都是比较复杂，而MATLAB语言正是处理线性方程组的求解的很好工具。线性代数是数学的一个分支，它的研究对象是向量，向量空间（或称线性空间），线性变换和有限维的线性方程组。因而，线性代数被广泛地应用于抽象代数和泛函分析中；由于科学研究中的非线性模型通常可以被近似为线性模型，使得线性代数被广泛地应用于自然科学和社会科学中。线性代数是数学的一个分支，它的研究对象是向量，向量空间（或称线性空间），线性变换和有限维的线性方程组。因而，线性代数被广泛地应用于抽象代数和泛函分析中；由于科学研究中的非线性模型通常可以被近似为线性模型，使得线性代数被广泛地应用于自然科学和社会科学中。线性代数是讨论矩阵理论、与矩阵结合的有限维向量空间及其线性变换理论的一门学科。
【关键字】线性代数
MATLAB语言
秩
矩阵
解
一、基本概念
1、N级行列式A：A等于所有取自不同性不同列的n个元素的积的代数和。
2、矩阵B：矩阵的概念是很直观的，可以说是一张表。
3、线性无关：一向量组（a1,a2,„,an）不线性相关，既没有不全为零的数k1,k2,„„„kn使得：k1*a1+k2*a2+„„„+kn*an=0
4、秩：向量组的极在线性无关组所含向量的个数成为这个向量组的秩。
5、矩阵B的秩：行秩，指矩阵的行向量组的秩；列秩类似。记：R(B)
a11*x1 a12*x2  a1n xn b1 a*x a*x  a*x b 2222nn26、一般线性方程组是指形式： 211
      as1*x1 as2*x2  asn*xn bs
二、基本理论
三种基本变换：1，用一非零的数乘某一方程；2，把一个方程的倍数加到另一方程；3，互换两个方程的位置。以上称出等变换。消元法
首先用初等变换化线性方程组为阶梯形方程组：1，如果剩下的方程当中最后的一个等式等于一非零数，那么方程组无解；否则有解；2，如果阶梯形方程组中方程的个数r等于未知量的个数，那么方程组有唯一的解；3，如果阶梯形方程组中方程的个数r小于未知量的个数，那么方程组就有无穷个解。定理1：线性方程组有解的充要条件为：R(A)=R(A,b)线性方程组解的结构：
1：对齐次线性方程组，a:两个解的和还是方程组的解；b:一个解的倍数还是方程组的解。定义：齐次线性方程组的一组解u1,u2,„ui称为齐次线性方程组的一个基础解系，如果：齐次线性方程组的任一解都表成u1,u2,„ui的线性组合，且u1,u2,„ui线性无关。2：对非齐次线
性方程组
（1）方程组的两个解的差是（2）的解。
（2）方程组的一个解与（2）的一个解之和还是（1）的解。
定理2 如果R0是方程组（1）的一个特解，那么方程组（1）的任一个解R都可以表成;R=R0+V„.其中V是（2）的一个解，因此，对方称（1）的任一特解R0，当v取遍它的全部解时，（3）就给出了(1)全部解。
三：基本思路
线性方程的求解分为两类：一类是方程求唯一解或求特解；一类是方程组
求无穷解即通解。
（1）判断方程组解的情况。1:当R(A)=R(B)时，有解（R(A)=R(A，b)）>=n 唯一解，R(A)=R(A,b)(n,有无穷解)；2：当R(B)+1=R(A，b)时无解。
（2）求特解；
（3）求通解（无穷解），线性方程组的无穷解=对应齐次方程组的通解+非齐 次方程组的一个特解；
注：以上针对非齐次线性方程组，对齐次线性方程组，主要使用到（1）,(2)步！
四、基本方法
基本思路将在解题的过程中得到体现。
1、（求线性方程组的唯一解或特解），这类问题的求法分为两类：一类主要用于解低阶稠密矩阵——直接法；一类是解大型稀疏矩阵——迭代法。
2、利用矩阵除法求线性方程组的特解（或一个解）
方程：AX=b,解法：X=Ab，（注意此处’’不是’/’）
(1)求方程组
2x1 x2 x3 x4 2 x x 2x x 4 1234  4x 6x 2x 2x 4234 1  3x1 6x2 9x3 7x4 9 命令如下：
A=[2,-1,-1,1;1,1,-2,1;4,-6,2,-2;3,6,-9,7];%产生4x4阶系数矩阵 b=[2,4,4,9]';%对矩阵进行转置 x=Ab %进行左初运算 x = NaN Inf Inf-3.0000 曾介绍过利用矩阵求逆来解线性方程组，即其结果于使用左除是相同的。
2、利用矩阵的分解求线性方程组
矩阵分解是指根据一定的原理用某种运算将一个矩阵分解成若干个矩阵的乘积。常见矩阵分解如，LU,QR和Cholesky分解求方程组的解，这三种分解，再求大型方程组是很有用。其优点是运算速度快，可以节省磁盘空间，节省内存。（1）LU分解又称Gauss 消去分解，可把任意方阵分解为下三角矩阵的基本变
换形式（行变换）和上三角矩阵的乘积。即A=LU,L为下三角阵,U为上三角
阵。
则：A*X=b 变成L*U*X=b 所以X=U(Lb)这样可以大大提高运算速度。
命令 [L,U]=lu(A)在matlab中可以编如下通用m文件； 在MATLAB建立M文件如下 % exp1.m
A;b;[L,U]=lu(A);%产生一个三角矩阵A和一个变换形式的下三叫矩阵L（交 换行），使之满足A=LU X=U(Lb)%L右乘b的结果再右乘U得到x的值
x1 x2 x3 x4 1 (2)求方程组  2x1 5x2 3x3 2x4 3
7x 7x 3x x 6234 1命令如下：
A=[1,1,-1,-1;2,-5,3,2;7,-7,3,1];%产生3x4阶系数矩阵 b=[1,3,6]';%对矩阵进行转置
[L,U]=lu(A);%产生一个三角矩阵A和一个变换形式的下三叫矩阵L（交换行），使之满足A=LU x=U(Lb)%L右乘b的结果再右乘U得到x的值 x = 0.4286-0.4286 0 0 rank = 2, tol = 6.7642e-015.采用第二种格式分解，在MATLAB建立M文件如下 %exp1.m A;b;[L,U,P]=lu(A);X=U(LP*b) x1 x2 x3 x4 1 (3)求方程组  2x1 5x2 3x3 2x4 3
7x 7x 3x x 7234 1命令如下：
A=[1,1,-1,-1;2,-5,3,2;7,-7,3,1];%产生3x4阶矩阵 b=[1,3,7]';%对矩阵进行转置
[L,U,P]=lu(A);%产生一个三角矩阵A和一个下三角阵L以
及一个置换矩阵P,使之满足PA=LU x=U(LP*b)%x的值 x = 0.6667-0.3333 0 0
rank = 2, tol = 6.7642e-015.(II)Cholesky分解
若A为对成正定矩阵，则Cholesky分解可将矩阵A分解成上三角矩阵和其 转置的乘积，即:A=R’*R 其中R为上三角矩阵。
方程 A*X=b 变成 R’*R*X=b 所以 X=R(R’b)在MATLAB中建立M文件如下 %exp2.m A;b;[R’,R]=chol(A);%产生一个上三角矩阵R，使R’R=A X=R(R’b)%x的值  x1 x2 x3 x4 0 (4)求方程组  x1 x2 x3 3x4 1

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。