a =

1 2 3 4 5 6 7 8 9 10 b =

5 3 1 -1 -3 -5 1. 代码:

linspace(1,10,10) linspace(5,-5,6) ans =

1 2 3 4 5 6 7 8 9 ans =

5 3 1 -1 -3 -5 2. 代码：

logspace(2,4,10) ans =

1.0e+004 *

0.0100 0.0167 0.0278 0.0464 0.0774 0.5995 1.0000 3. 代码:

D=10*rand(3,5) 结果: D =

0.1527 9.3181 8.4622 6.7214 6.8128 7.4679 4.6599 5.2515 8.3812 3.7948 4.4510 4.1865 2.0265 0.1964 8.3180 5.代码:

A=magic(5); e=diag(A); E=diag(e) 结果； E =

17 0 0 0 0 0 5 0 0 0 0 0 13 0 0 0 0 0 21 0 0 0 0 0 9 6.代码:

AA=magic(3); BB=rot90(AA,2); CC=AA+j*BB; DD=CC.'; EE=conj(CC); A=abs(CC) B=angle(CC) C=abs(EE)

10 0.1292 0.2154 0.3594

D=angle(EE) 结果: A =

8.2462 9.0554 7.2111 7.6158 7.0711 7.6158 7.2111 9.0554 8.2462 B =

0.2450 1.4601 0.5880 1.1659 0.7854 0.9828 0.1107 C =

8.2462 9.0554 7.6158 7.0711 7.2111 9.0554 D =

-0.2450 -1.4601 -1.1659 -0.7854 -0.9828 -0.1107 7.代码:

f=20*0.5.^(0:9); g=1:3:28; s=dot(f,g) 结果: s =

158.6719 8.代码： A=magic(9); M=A(4:6,4:6); sum(A,1) sum(A,2) 结果： ans =

369 369 369 ans = 369 369 369 369 369 369 369 369 369 9．代码：

0.4049 1.3258 7.2111 7.6158 8.2462 -0.5880 -0.4049 -1.3258 369 369 369 369 369 369 T=[1,2,3,4;2,3,4,5;3,4,5,6;4,5,6,7] A=rot90(T); B=tril(A);

T1=rot90(B,-1) 结果： T =

1 2 3 4 2 3 4 5 3 4 4 5 T1 =

1 2 2 3 3 4 4 0 5 6 6 7 3 4 4 0 0 0 0 0

实验四：线性方程组的求解

一、实验目的

1．掌握恰定方程组的解法。

2．了解欠定方程组、超定方程组的解法。 3．掌握左除法求解线性方程组。 4．学会测试程序的运行时间。

二、实验内容

1．用两种方法求下列方程组的解，并比较两种方法执行的时间。 2．判定下列方程是恰定方程组、欠定方程组还是超定方程组，并求其解。 3．用网孔电流法求如下电路的各支路电流。 4．用结点电压法求如下电路的结点电压un1、un2。

三、设计提示

1．在计算程序的执行时间之前，应注意用clear命令将内存变量清空。 2．求得线性方程组的解之后，代入原方程验证是否正确。

四、实验报告要求

1．编写实现第二节实验内容中所使用的函数命令，并记录相应的生成结果。 2．对于电路的求解，应列出相应的网孔方程和结点方程，并注意方向。

3．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。

五、参考答案 第一题：

A=[7,14,-9,-2,5;3,-15,-13,-6,-4;-11,-9,-2,5,7;5,7,14,16,-2;-2,5,12,-11,-4]; B=[100;200;300;400;500]; X=A\\B 结果: X =

370.9455 -224.7276 238.3333 -138.1891 503.6378

第二题：

A=[6,9,14,-11,5;1,14,-7,-15,-6;-2,1,-7,12,-1;6,11,11,-9,-13]; B=[68;294;-441;103];