线性代数实验报告

线性代数实验报告

姓名班级学号得分

2013年12月24日

数学实验报告题目

一、 实验目的

1．熟悉MATLAB 的矩阵初等运算；

2．掌握求矩阵的秩、逆、化最简阶梯形的命令； 3．会用MABLAB 求解线性方程组

二、 实验问题

1． 已知??

??

??????--=351503224A ，??????????---=112302431

B ，在MATLAB 命令窗口中建立A 、B 矩阵并对其进行以下操作：

(1) 计算矩阵A 的行列式的值det()?A =

(2) 分别计算下列各式：B A -2 、 B A *和B A *.、 1-AB 、 B A 1-、 2A 、 T A

2． 在MATLAB 中分别利用矩阵的初等变换及函数rank 、函数inv 求下列矩阵的秩和逆：(1) ??????????----=4211104532361A 求 Rank(A)=? (2)3501120010201202B ??

????=??????

求?1=-B 3. 在MATLAB 中判断下列向量组是否线性相关，并找出向量组中的一个最大线

性无关组：()11,1,3,2α'=,()2

1,1,1,3α'=--,()35,2,8,9α'=-,()41,3,1,7α'=-

4、在MATLAB 中判断下列方程组解的情况，若有多个解，写出通解：

（1）???????=+--=-++=+--=-+-0

61230

2730

20

24432143214

3214321x x x x x x x x x x x x x x x x (2) ??????

?-=+-=-+-=+-=++69413283542432321321321321x x x x x x x x x x x x 5、化方阵222254245A -??

?=- ? ?--??

为对角阵.

6、求一个正交变换，将二次型222

123121323553266f x x x x x x x x x =++-+-化为标准型。

7、判定三元二次方程的空间图形：C x x x x x x =--+32212

221442（分C=0,>0,<0三种情况讨论）。

三、 实验过程及结果分析

1．已知

?

?

?

?

?

?

?

?

?

?

-

-

=

3

5

1

5

3

2

2

4

A，

?

?

?

?

?

?

?

?

?

?

-

-

-

=

1

1

2

3

2

4

3

1

B，在MATLAB命令窗口中建立A、

B矩阵并对其进行以下操作：

(1)计算矩阵A的行列式的值det()?

A=

【程序设计】：

【结果分析】：

用det(A)算出矩阵A的行列式的值：

(2) 分别计算下列各式：B A -2 、 B A *和B A *.、 1-AB 、 B A 1-、 2A 、

T A

【程序设计】：

【结果分析】：

A ’表示矩阵A 的转置；

A^n 表示方阵A 的n 次方幂；

A/B 在矩阵B 可逆的情况下，表示1-AB ； A\B 在矩阵A 可逆的情况下，表示B A 1-；

2． 在MATLAB 中分别利用矩阵的初等变换及函数rank 、函数inv 求下列矩阵的

秩和逆：(1) ????

??????----=4211104532361A 求 Rank(A)=? (2)3

5011

20010201

20

2B ?????

?=??????

求?1=-B 【程序设计】

：

【结果分析】：

用rank(A)算出矩阵A 的秩； 用inv(B)算出矩阵B 的逆；

3.在MATLAB 中判断下列向量组是否线性相关，并找出向量组中的一个最大线性无关组：()11,1,3,2α'=,()2

1,1,1,3α'=--,()35,2,8,9α'=-,()41,3,1,7α'=-

【程序设计】：

【结果分析】：

观察得知由4,3,2,1αααα组成的矩阵A 化成的标准阶梯型的秩为3，3<4，所以它们线性

相关；

又因为r=3，所以3

,2,1ααα组成的向量组是最大的线性无关组。

4.在MATLAB 中判断下列方程组解的情况，若有多个解，写出通解：

（1）??????

?=+--=-++=+--=-+-0612302730

2024432143214

3214321x x x x x x x x x x x x x x x x (2) ??????

?-=+-=-+-=+-=++6

94132835424

323213213

21321x x x x x x x x x x x x

【程序设计】

：

【结果分析】：

根据下面的结果：

(1)由A 的标准阶梯型可知，A 为满秩矩阵，,04321====x x x x 是原方程组的唯一解；

(2)秩为2，2<3,所以原方程组由无数多组解，通解为?

???? ??-+????? ??-=????? ??021112321k x x x ，k 为任意常

数。

5、化方阵

222

254

245

A

-

??

?

=-

?

?

--

??

为对角阵.

【程序设计】：

【结果分析】：

通过将矩阵A化成标准阶梯型而化成对角阵。

6、求一个正交变换，将二次型222

123121323553266f x x x x x x x x x =++-+-化为标准型。

【程序设计】：

【结果分析】：

由下面算出的矩阵得知22

2194y y f +=

7、判定三元二次方程的空间图形：C x x x x x x =--+32212221

442（分C=0,>0,<0三种情况讨论）。

【程序设计】：

【结果分析】：

由D 可以得知方程对应矩阵的特征值为-2、1、4； 所以标准型为C y y y =++-23

22

21

42；

从而分如下三种情况讨论：

(ⅰ)C=0时,此三元二次方程的空间图形为开口沿1y 方向的椭圆锥面；

(ⅱ)C>0时,此三元二次方程的空间图形为开口沿1y 方向的单叶双曲面； (ⅲ)C<0时,此三元二次方程的空间图形为开口沿1

y 方向的双叶双曲面。

四、 实验总结与体会

在平时的线性代数运算中，时常会遇到繁琐的计算，费时费力，而MATLAB 提供了方便快捷的运算，大大地减少了题目的运算量，使我受益匪浅。

通过本次试验，我学习到多种MATLAB 有关线性代数运算的指令，主要学习运用MATLAB 解决矩阵除法，线性方程组的通解，矩阵相似对角化问题，二次型化为标准型，计算矩阵特征值等等。熟悉了MATLAB 的矩阵初等运算、掌握求矩阵的秩、逆、化最简阶梯形的命令，会用MATLAB 求解线性方程组，并综合运用多种指令解决应用题，十分方便准确快捷。在此次实验学习实践的过程中，加深了对线性代数和MATLAB 的理解，也产生了对本学科更深的兴趣。相信在以后更多的实践中能够更加熟练地运用MATLAB 解决实际问题，并继续深入学习。

线性代数(李建平)习题答案详解__复旦大学出版社

线性代数课后习题答案 习题一 1.2.3（答案略） 4. (1) ∵ (127435689)415τ=+= (奇数) ∴ (127485639)τ为偶数 故所求为127485639 (2) ∵(397281564)25119τ=+++= （奇数） ∴所求为397281564 5.(1)∵(532416)421106τ=++++= (偶数) ∴项前的符号位()6 11-=+ （正号） （2）∵325326114465112632445365a a a a a a a a a a a a = (162435)415τ=+= ∴ 项前的符号位5(1)1-=- （负号） 6. (1) (2341)(1)12n n τ-?L L 原式=(1)(1)!n n -=- （2）()((1)(2)21) 1(1)(2)21n n n n n n τ--??---??L L 原式=(1)(2) 2 (1) !n n n --=- （3）原式=((1)21) 12(1)1(1) n n n n n a a a τ-?--L L (1) 2 12(1)1(1)n n n n n a a a --=-L 7.8（答案略） 9. ∵162019(42)0D x =?-?+?--?= ∴7x = 10. （1）从第2列开始，以后各列加到第一列的对应元素之上，得 []11(1)1110 01(1)1110 (1)1 1 (1)1 1 1 x x n x x x n x x x n x x n x x +-+--==+-+--L L L L L L L L L L L L L L L L L L L L L []1(1)(1)n x n x -=+-- （2）按第一列展开： 11100000 (1)(1)0 0n n n n n y x y D x x y x y x y -++=?+-=+-L L L L L L L L

用MATLAB解决线性代数问题实验报告

实验三使用MATLAB解决线性代数问题学院：数计学院班级：1003班姓名：黄晓丹学号：1051020144 实验目的： 学习MATLAB有关线性代数运算的指令，主要学习运用MATLAB解决矩阵除法，线性方程组的通解，矩阵相似 对角化问题，以及解决投入产出分析等应用问题。 实验内容： 矩阵转置：A=[1 2;3 4];B=[4 3;2 1]; >> A',B' ans = 1 3 2 4 ans = 4 3 3 1 矩阵加减：A-B ans= -3 -1 1 3 矩阵乘法：A*B,A.*B(数组乘法)||比较矩阵乘法与数组乘法的区别ans= 8 5 20 13 ans= 4 6 6 4 矩阵除法：A\B,B./A ans=

-6 -5 5 4 ans= 4 1.5 0.6667 0.25 特殊矩阵生成：zeros（m，n）||生成m行n列的矩阵 ones（m，n）||生成m行n列的元素全为一的矩阵 eye（n）||生成n阶单位矩阵 rand（m，n）||生成m行n列[0 ，1]上均匀分布随 机数矩阵 zeros(2,3) ans = 0 0 0 0 0 0 >> ones(3,3) ans = 1 1 1 1 1 1 1 1 1 >> eye(3)

ans = 1 0 0 0 1 0 0 0 1 >> rand(2,4) ans = Columns 1 through 3 0.9501 0.6068 0.8913 0.2311 0.4860 0.7621 Column 4 0.4565 0.0185 矩阵处理：trace（A）||返回矩阵的迹 diag（A）||返回矩阵对角线元素构成的向量 tril（A）||提取矩阵的下三角部分 triu（A）||提取矩阵的上三角部分 flipud（A）||矩阵上下翻转 fliplr（A）||矩阵左右翻转 reshape(A,m,n)||将矩阵的元素重排成m行n列矩阵A=[1 2 3;4 5 6;7 8 9]; >> t=trace(A),d=diag(A),u=triu(A)

信号与系统实验报告总结

信号与系统实验 实验一常用信号的观察 方波： 正弦波： 三角波： 在观测中，虚拟示波器完全充当实际示波器的作用，在工作台上连接AD1为示波器的输入，输入方波、正弦波、三角波信号时，可在电脑上利用软件观测到相应的波形，其纵轴为幅值可通过设置实现幅值自动调节以观测到最佳大小的波形，其横轴为时间，宜可通过设置实现时间自动调节以观测到最佳宽度的波形。

实验四非正弦周期信号的分解与合成 方波DC信号： DC信号几乎没有，与理论相符合，原信号没有添加偏移。 方波基波信号： 基波信号为与原方波50Hz信号相对应的频率为50Hz的正弦波信号，是方波分解的一次谐波信号。 方波二次谐波信号: 二次谐波信号频率为100Hz为原方波信号频率的两倍，幅值较一次谐波较为减少。

方波三次谐波信号： 三次谐波信号频率为150Hz为原方波信号的三倍。幅值较一二次谐波大为减少。方波四次谐波信号: 四次谐波信号的频率为200Hz为原方波信号的四倍。幅值较三次谐波再次减小。方波五次谐波信号： 五次谐波频率为250Hz为原方波信号的五倍。幅值减少到0.3以内，几乎可以忽略。 综上可知：50Hz方波可以分解为DC信号、基波信号、二次、三次、四次、五次谐波信号…，无偏移时即无DC信号，DC信号幅值为0。分解出来的基波信号即一次谐波信号频率与原方波信号频率相同，幅值接近方波信号的幅值。二次谐波、三次谐波、四次谐波、五次谐波依次频率分别为原方波信号的二、三、四、五倍，且幅值依次衰减，直至五次谐波信号时几乎可以忽略。可知，方波信号可分解为多个谐波。

方波基波加三次谐波信号： 基波叠加上三次谐波信号时，幅值与方波信号接近，形状还有一定差异，但已基本可以看出叠加后逼近了方波信号。 方波基波加三次谐波信号加五次谐波信号： 基波信号、三次谐波信号、五次谐波信号叠加以后，比基波信号、三次谐波信号叠加后的波形更加接近方波信号。 综上所述：方波分解出来的各次谐波以及DC信号，叠加起来以后会逼近方波信号，且叠加的信号越多，越是接近方波信号。说明，方波信号可有多个谐波合成。

线性代数知识点总结汇总

线性代数知识点总结 1 行列式 （一）行列式概念和性质 1、逆序数：所有的逆序的总数 2、行列式定义：不同行不同列元素乘积代数和 3、行列式性质：（用于化简行列式） （1）行列互换（转置），行列式的值不变 （2）两行（列）互换，行列式变号 （3）提公因式：行列式的某一行（列）的所有元素都乘以同一数k，等于用数k 乘此行列式 （4）拆列分配：行列式中如果某一行（列）的元素都是两组数之和，那么这个行列式就等于两个行列式之和。 （5）一行（列）乘k加到另一行（列），行列式的值不变。 （6）两行成比例，行列式的值为0。 （二）重要行列式 4、上（下）三角（主对角线）行列式的值等于主对角线元素的乘积 5、副对角线行列式的值等于副对角线元素的乘积乘 6、Laplace展开式：（A是m阶矩阵，B是n阶矩阵），则 7、n阶（n≥2）范德蒙德行列式

数学归纳法证明 ★8、对角线的元素为a，其余元素为b的行列式的值： （三）按行（列）展开 9、按行展开定理： （1）任一行（列）的各元素与其对应的代数余子式乘积之和等于行列式的值（2）行列式中某一行（列）各个元素与另一行（列）对应元素的代数余子式乘积之和等于0 （四）行列式公式 10、行列式七大公式： （1）|kA|=k n|A| （2）|AB|=|A|·|B| （3）|A T|=|A| （4）|A-1|=|A|-1 （5）|A*|=|A|n-1 （6）若A的特征值λ1、λ2、……λn，则 （7）若A与B相似，则|A|=|B| （五）克莱姆法则 11、克莱姆法则： （1）非齐次线性方程组的系数行列式不为0，那么方程为唯一解

（2）如果非齐次线性方程组无解或有两个不同解，则它的系数行列式必为0 （3）若齐次线性方程组的系数行列式不为0，则齐次线性方程组只有0解；如果方程组有非零解，那么必有D=0。 2 矩阵 （一）矩阵的运算 1、矩阵乘法注意事项： （1）矩阵乘法要求前列后行一致； （2）矩阵乘法不满足交换律；（因式分解的公式对矩阵不适用，但若B=E,O,A-1，A*,f(A)时，可以用交换律） （3）AB=O不能推出A=O或B=O。 2、转置的性质（5条） （1）（A+B）T=A T+B T （2）（kA）T=kA T （3）（AB）T=B T A T （4）|A|T=|A| （5）（A T）T=A （二）矩阵的逆 3、逆的定义： AB=E或BA=E成立，称A可逆，B是A的逆矩阵，记为B=A-1 注：A可逆的充要条件是|A|≠0 4、逆的性质：（5条） （1）（kA）-1=1/k·A-1 (k≠0) （2）（AB）-1=B-1·A-1 （3）|A-1|=|A|-1 （4）（A T）-1=（A-1）T （5）（A-1）-1=A

北大版 线性代数第一章部分课后答案详解

习题1.2： 1 .写出四阶行列式中 11121314212223243132333441 42 43 44 a a a a a a a a a a a a a a a a 含有因子1123a a 的项 解：由行列式的定义可知，第三行只能从32a 、34a 中选，第四行只能从42a 、44a 中选，所以所有的组合只有() () 13241τ-11233244a a a a 或() () 13421τ-11233442a a a a ，即含有因子1123a a 的项 为11233244a a a a 和11233442a a a a 2. 用行列式的定义证明111213141521 22232425 31 3241425152 000000000 a a a a a a a a a a a a a a a a =0 证明：第五行只有取51a 、52a 整个因式才能有可能不为0，同理，第四行取41a 、42a ，第三行取31a 、32a ，由于每一列只能取一个，则在第三第四第五行中，必有一行只能取0.以第五行为参考，含有51a 的因式必含有0，同理，含有52a 的因式也必含有0。故所有因式都为0.原命题得证.。 3.求下列行列式的值： （1）01000020;0001000 n n -L L M M M O M L L （2）00100200100000 n n -L L M O M O M L L ； 解：（1）0100 0020 0001 000 n n -L L M M M O M L L =()()23411n τ-L 123n ????L =()1 1!n n --

数值分析实验报告1

实验一误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中ε（1.1）和（1.221,,,a a 的输出b ”和“poly ε。 （1（2 （3）写成展 关于α solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程： 程序： a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);

ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n

很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号：06450210 姓名：万轩 实验二插值法

信号与系统实验报告1

学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟

20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求

实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。

实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;

线性代数实验题04-交通网络的流量分析

数学实验报告 学号: , 姓名: , 得分: 实验内容:实验题：交通网络流量分析问题（线性方程组应用） 城市道路网中每条道路、每个交叉路口的车流量调查，是分析、评价及改善城市交通状况的基础。 问题：某城市有下图所示的交通图，每条道路都是单行线，需要调查每条道路每小时的车流量。图中的数字表示该条路段的车流数。如果每个交叉路口进入和离开的车数相等，整个图中进入和离开的车数相等。 求（1）建立确定每条道路流量的线性方程组； （2）分析哪些流量数据是多余的； （3）为了唯一确定未知流量，需要增添哪几条道路的流量统计。 解： （1）由题意得：x1+ x7=400 x1+ x9= x2+300 x2+100=300+ x11 x3+ x7=350+ x8 x4+ x10= x9+ x3 x11+500= x4+ x12 x8+ x5=310 x6+400= x10+ x5 x12+150= x6+290

整理得： x 1+ x 7=400 x 1- x 2+ x 9=300 x 2+ x 11=200 x 3+ x 7- x 8=350 -x 3+x 4+ x 10- x 9=0 -x 4+x 11- x 12=-500 x 5 +x 8=310 - x 5+x 6- x 10=-400 -x 6+ x 12= 140 将方程组写成矩阵向量形式为AX = b 1 0 0 0 0 0 1 0 0 0 0 0 400 x 1 1 -1 0 0 0 0 0 0 1 0 0 0 300 x 2 0 1 0 0 0 0 0 0 0 0 1 0 200 x 3 A= 0 0 1 0 0 0 1 -1 0 0 0 0 b= 350 X= x 4 0 0 -1 1 0 0 0 0 -1 1 0 0 0 x 5 0 0 0 -1 0 0 0 0 0 0 1 -1 -500 x 6 0 0 0 0 1 0 0 1 0 0 0 0 310 x 7 0 0 0 0 -1 1 0 0 0 -1 0 0 -400 x 8 0 0 0 0 0 -1 0 0 0 0 0 1 140 x 9 x 10 x 11 x 12 在MATLAB 环境中，首先输入方程组的系数矩阵A 和方程组右端向量b A=[1,0,0,0,0,0,1,0,0,0,0,0;1,-1,0,0,0,0,0,0,1,0,0,0;0,1,0,0,0,0,0,0,0,0,1,0;0,0,1,0,0,0,1,-1,0,0,0,0;0,0,-1,1,0,0,0,0,-1,1,0,0;0,0,0,-1,0,0,0,0,0,0,1,-1;0,0,0,0,1,0,0,1,0,0,0,0;0,0,0,0,-1,1,0,0,0,-1,0,0;0,0,0,0,-1,0,0,0,0,0,1] b = [400;300;200;350;0;500;310;-400;140] 解得 x 1=- x 9+500 x 2=200 x 3=- x 9+ x 10- x 12

信号系统实验报告

电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现

一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数

数值分析实验报告

实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p Λ 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a Λ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a Λ 的全部根；而函数 poly(v)b = 的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots +

线性代数课后习题答案全)习题详解

线性代数课后习题答案全)习题详解 第一章 行列式 1.利用对角线法则计算下列三阶行列式： （1）381141102---； （2）b a c a c b c b a ; （3）222111c b a c b a ； （4）y x y x x y x y y x y x +++. 解 （1）=---3 811411 02811)1()1(03)4(2??+-?-?+?-?)1()4(18)1(2310-?-?-?-?-??- =416824-++-=4- （2）=b a c a c b c b a cc c aaa bbb cba bac acb ---++3333c b a abc ---= （3）=2 221 11c b a c b a 222222cb ba ac ab ca bc ---++))()((a c c b b a ---= （4）y x y x x y x y y x y x +++yx y x y x yx y y x x )()()(+++++=333)(x y x y -+-- 33322333)(3x y x x y y x y y x xy ------+= )(233y x +-=

2.按自然数从小到大为标准次序，求下列各排列的逆序数： （1）1 2 3 4； （2）4 1 3 2； （3）3 4 2 1； （4）2 4 1 3； （5）1 3 … )12(-n 2 4 … )2(n ； （6）1 3 … )12(-n )2(n )22(-n … 2. 解（1）逆序数为0 （2）逆序数为4：4 1，4 3，4 2，3 2 （3）逆序数为5：3 2，3 1，4 2，4 1,2 1 （4）逆序数为3：2 1，4 1，4 3 （5）逆序数为 2 ) 1(-n n ： 3 2 1个 5 2，5 4 2个 7 2，7 4，7 6 3个 ……………… … )12(-n 2，)12(-n 4，)12(-n 6，…，)12(-n )22(-n )1(-n 个 （6）逆序数为)1(-n n 3 2 1个 5 2，5 4 2个 ……………… … )12(-n 2，)12(-n 4，)12(-n 6，…，)12(-n )22(-n )1(-n 个 4 2 1个 6 2，6 4 2个 ……………… … )2(n 2，)2(n 4，)2(n 6，…，)2(n )22(-n )1(-n 个 3.写出四阶行列式中含有因子2311a a 的项.

数值线性代数第二版徐树方高立张平文上机习题第一章实验报告(供参考)

上机习题 1.先用你所熟悉的的计算机语言将不选主元和列主元Gauss 消去法编写成通用的子程序；然后用你编写的程序求解84阶方程组；最后将你的计算结果与方程的精确解进行比较，并就此谈谈你对Gauss 消去法的看法。 Sol ： （1）先用matlab 将不选主元和列主元Gauss 消去法编写成通用的子程序，得到P U L ,,： 不选主元Gauss 消去法：[])(,A GaussLA U L =得到U L ,满足LU A = 列主元Gauss 消去法：[])(,,A GaussCol P U L =得到P U L ,,满足LU PA = （2）用前代法解()Pb or b Ly =，得y 用回代法解y Ux =，得x 求解程序为()P U L b A Gauss x ,,,,=（P 可缺省，缺省时默认为单位矩阵） （3）计算脚本为ex1_1 代码 %算法（计算三角分解：Gauss 消去法） function [L,U]=GaussLA(A) n=length(A); for k=1:n-1 A(k+1:n,k)=A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n); end

U=triu(A); L=tril(A); L=L-diag(diag(L))+diag(ones(1,n)); end %算法计算列主元三角分解：列主元Gauss消去法） function [L,U,P]=GaussCol(A) n=length(A); for k=1:n-1 [s,t]=max(abs(A(k:n,k))); p=t+k-1; temp=A(k,1:n); A(k,1:n)=A(p,1:n); A(p,1:n)=temp; u(k)=p; if A(k,k)~=0 A(k+1:n,k)=A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n); else break; end end L=tril(A);U=triu(A);L=L-diag(diag(L))+diag(ones(1,n));

信号与系统实验报告_1(常用信号的分类与观察)

实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：

图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)

线性代数课后习题答案

线性代数课后题详解 第一章 行列式 1.利用对角线法则计算下列三阶行列式： 相信自己加油 （1） 3811411 02 ---； （2）b a c a c b c b a （3） 2 2 2 111 c b a c b a ； （4） y x y x x y x y y x y x +++. 解 注意看过程解答（1）=---3 81141 1 2811)1()1(03)4(2??+-?-?+?-? )1()4(18)1(2310-?-?-?-?-??- =416824-++- =4- （2） =b a c a c b c b a cc c aaa bbb cba bac acb ---++ 3333c b a abc ---= （3） =2 2 2 1 11c b a c b a 222222cb ba ac ab ca bc ---++ ))()((a c c b b a ---= （4） y x y x x y x y y x y x +++ yx y x y x yx y y x x )()()(+++++=333)(x y x y -+-- 33322333)(3x y x x y y x y y x xy ------+= )(233y x +-= 2.按自然数从小到大为标准次序，求下列各排列的逆序数：耐心成就大业 （1）1 2 3 4； （2）4 1 3 2； （3）3 4 2 1； （4）2 4 1 3； （5）1 3 … )12(-n 2 4 … )2(n ； （6）1 3 … )12(-n )2(n )22(-n … 2. 解（1）逆序数为0

（2）逆序数为4：4 1，4 3，4 2，3 2 （3）逆序数为5：3 2，3 1，4 2，4 1,2 1 （4）逆序数为3：2 1，4 1，4 3 （5）逆序数为2 ) 1(-n n ： 3 2 1个 5 2，5 4 2个 7 2，7 4，7 6 3个 ……………… … )12(-n 2，)12(-n 4，)12(-n 6，…，)12(-n )22(-n )1(-n 个 （6）逆序数为)1(-n n 3 2 1个 5 2，5 4 2个 ……………… … )12(-n 2，)12(-n 4，)12(-n 6，…，)12(-n )22(-n )1(-n 个 4 2 1个 6 2，6 4 2个 ……………… … )2(n 2，)2(n 4，)2(n 6，…，)2(n )22(-n )1(-n 个 3.写出四阶行列式中含有因子 2311a a 的项. 解 由定义知，四阶行列式的一般项为 43214321)1(p p p p t a a a a -，其中t 为4321p p p p 的逆序数．由于3,121==p p 已固定， 4321p p p p 只能形如13□□，即1324或1342.对应的t 分别为 10100=+++或22000=+++ ∴44322311a a a a -和42342311a a a a 为所求. 4.计算下列各行列式： 多练习方能成大财 （1）?? ??????? ???711 00251020214214； （2）????? ? ??? ???-26 0523******** 12； （3）???? ??????---ef cf bf de cd bd ae ac ab ； （4）?? ??? ???????---d c b a 100 110011001 解 (1) 7110025102021421434327c c c c --0 1001423102 02110214--- =34)1(14 3102211014+-?---

信号实验报告

大连理工大学 本科实验报告 课程名称：信号与系统实验 学院(系）：电子信息与电气工程学部专业: 通信工程 班级: 1401班 学号：20148３091 学生姓名：李睿 20１6年 5 月2１日 ?实验项目列表

?大连理工大学实验预习报告 学院（系）：电信专业:通信工程班级:1401班 姓名:李睿学号：20148309１组:5 _＿_ 实验时间：2016、5、6 实验室：创新园大厦c0221 实验台： 5 指导教师签字：成绩： 信号得频谱图 一、实验目得与要求 1、掌握周期信号得傅里叶级数展开 ２、掌握周期信号得有限项傅里叶级数逼近 3、掌握周期信号得频谱分析 4、掌握连续非周期信号得傅立叶变换 5、掌握傅立叶变换得性质 二、实验用得mａｔlａb命令与例子

１、ａ:b：ｃ：产生一个从ａ到 c，间隔为ｂ得等间隔数列例：5：1:11,产生一个从 5 到11,间隔为 1 得等间隔数列 ２、ｑuａrｅ（t，dutｙ)：周期性矩形脉冲信号(duty 表示占空比）调用形式： y=squaｒe(ｔ，dｕtｙ）例：产生一个周期为2π，幅值为±１得周期性方波。y＝sｑuarｅ(2＊pi＊30*t,７5）； plot（ｔ,y）,grｉd on aｘis（[—０、１,０、1，—1、5，1、5]） 3、plot（）：ｍaｔlab 中二维线画图函数ｐloｔ（x，ｙ，’颜色与标识’）：若 y 与ｘ为同维向量，则以ｘ为横坐标，y 为纵坐标绘制连线图. 若x 就是向量，y 就是行数或列数与ｘ长度相等得矩阵,则绘制多条不同色彩得连线图，x 被作为这些曲线得共同横坐标.若 x 与 y 为同型矩阵，则以ｘ，ｙ对应元素分别绘制曲线,曲线条数等于矩阵列数. 例:在０≤x≤2π区间内，绘制曲线 y=2ｅ-0、5xcoｓ(４πx）。 ｘ=0：2＊pi； y＝２＊eｘp(－０、5*x）、＊cos（４＊ｐｉ＊x); plot（x，y) ‘’：y 黄ｍ紫 c 青 r 红 g 绿 b 蓝ｗ白 k 黑—实线、点 <小于号 :点线ｏ圆s 正方形 -、点划线x 叉号 d 菱形- -虚线 +加号h 六角星 *星号 p 五角星 v 向下三角形 ^向上三角形〉大于号 4、grｉd ｏｎ:有网格 gｒｉｄ off:关掉格网下面就是加上命令grｉd on后画得图，有网格. 5、 aｘis（[a b c d］）:表明图线得x轴范围为a～by轴范围为c～ｄ例：ｐlｏt(x,ｙ）ａｘis(［０ 1 ２３]) ｇｒｉｄ oｎ ６、 lengｔｈ（a)：表示矩阵ａ得最大得长度比如lｅngth([1 2 3；4 5 ６]) 等于3，因为2行与3列中最大就是3。当a就是向量时,即表示向量得元素个数，因为向量总就是１×n或ｎ×1得，而n一定大于或等于1、所以得到得结果一定就是n. ７、 1、/tan(pi、*x)：表示点乘。点乘就是值对值得运算上面得式子中 X 可能就是一个向量或矩阵，PＩ后面得点就是一个ＰI 与一个向量相乘，得到得也就是一个向量;1 后面乘得自然也就是个向量所以要加点,也就就是对应不同得Ｘ，有不同得 Y 值. 8.fｉｇurｅ就是建立图形得意思. 系统自动从 1，2，3,4、、、来建立图形，数字代表第几幅图形,figure（1)，figurｅ（２）就就是第一第二副图得意思,在建立图形

数值线性代数第二版徐树方高立张平文上机习题第三章实验报告

数值线性代数第二版徐树方高立张平文上机习题第三章实验报告

第三章上机习题 用 你所熟悉的的计算机语言编制利用QR 分解 求解线性方程组和线性最小二乘问题的通用子程序，并用你编制的子程序完成下面的计算任务： （1）求解第一章上机习题中的三个线性方程组，并将所得的计算结果与前面的结果相比较，说明各方法的优劣； （2）求一个二次多项式+bt+c y=at 2 ，使得在残向量 的2范数下最小的意义下拟合表3.2中的数据； （3）在房产估价的线性模型 11 1122110x a x a x a x y ++++= 中，11 2 1 ,,,a a a 分别表示税、浴室数目、占地面积、车库数目、房屋数目、居室数目、房龄、建筑类型、户型及壁炉数目，y 代表房屋价格。现根据表3.3和表3.4给出的28组数据，求出模型中参数的最小二乘结果。

（表3.3和表3.4见课本P99-100） 解 分析： （1）计算一个Householder 变换H ： 由于T T vv I ww I H β-=-=2，则计算一个Householder 变换H 等价于计算相应的 v 、β。其中 ) /(2,||||12v v e x x v T =-=β。 在实际计算中， 为避免出现两个相近的数出现的情形，当0 1 >x 时， 令 2 12 221||||)(-x x x x v n +++= ； 为便于储存，将v 规格化为1 /v v v =，相应的，β变为)/(22 1 v v v T =β 为防止溢出现象，用∞ ||||/x x 代替 （2）QR 分解： 利用Householder 变换逐步将n m A n m ≥?,转化为上三 角矩阵A H H H n n 11 -=Λ，则有

数值线性代数二版徐树方高立张平文上机习题第三章实验报告

- 1 - 第三章上机习题 用你所熟悉的的计算机语言编制利用QR 分解求解线性方程组和线性最小二乘问题的 通用子程序，并用你编制的子程序完成下面的计算任务： （1）求解第一章上机习题中的三个线性方程组，并将所得的计算结果与前面的结果相比较，说明各方法的优劣； （2）求一个二次多项式+bt+c y=at 2 ，使得在残向量的2范数下最小的意义下拟合表3.2中的数据； （3）在房产估价的线性模型 111122110x a x a x a x y ++++= 中，1121,,,a a a 分别表示税、浴室数目、占地面积、车库数目、房屋数目、居室数目、房龄、建筑类型、户型及壁炉数目，y 代表房屋价格。现根据表3.3和表3.4给出的28组数据，求出模型中参数的最小二乘结果。 （表3.3和表3.4见课本P99-100） 解 分析： （1）计算一个Householder 变换H ： 由于T T vv I ww I H β-=-=2，则计算一个Householder 变换H 等价于计算相应的v 、β。其中)/(2,||||12v v e x x v T =-=β。 在实际计算中， 为避免出现两个相近的数出现的情形，当01>x 时，令2 12221||||) (-x x x x v n +++= ； 为便于储存，将v 规格化为1/v v v =，相应的，β变为)/(221v v v T =β 为防止溢出现象，用∞||||/x x 代替 （2）QR 分解： 利用Householder 变换逐步将n m A n m ≥?,转化为上三角矩阵A H H H n n 11 -=Λ，则有

?? ? ???=0R Q A ，其中n H H H Q 21=，:),:1(n R Λ=。 在实际计算中，从n j :1=，若m j <，依次计算)),:((j m j A x =对应的)1()1()~ (+-?+-k m k m j H 即对应的j v ，j β，将)1:2(+-j m v j 储存到),:1(j m j A +，j β储存到)(j d ，迭代结束 后再次计算Q ，有??? ? ?? ??=-~001 j j j H I H ，n H H H Q 21=（m n =时1-21n H H H Q =） （3）求解线性方程组b Ax =或最小二乘问题的步骤为 i 计算A 的QR 分解； ii 计算b Q c T 11=，其中):1(:,1n Q Q = iii 利用回代法求解上三角方程组1c Rx = （4）对第一章第一个线性方程组，由于R 的结果最后一行为零，故使用前代法时不计最后一行，而用运行结果计算84x 。 运算matlab 程序为 1 计算Householder 变换 [v,belta]=house(x) function [v,belta]=house(x) n=length(x); x=x/norm(x,inf); sigma=x(2:n)'*x(2:n); v=zeros(n,1); v(2:n,1)=x(2:n); if sigma==0 belta=0; else alpha=sqrt(x(1)^2+sigma); if x(1)<=0 v(1)=x(1)-alpha; else v(1)=-sigma/(x(1)+alpha); end belta=2*v(1)^2/(sigma+v(1)^2); v=v/v(1,1); end end

信号与系统实验报告

实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10()0 t u t t >?=? 0); %定义函数体，即函数所执行指令