9
实验四：线性方程组的求解
一、实验目的
1．掌握恰定方程组的解法。
2．了解欠定方程组、超定方程组的解法。 3．掌握左除法求解线性方程组。 4．学会测试程序的运行时间。
二、实验内容
1．用两种方法求下列方程组的解，并比较两种方法执行的时间。
?7x1 ?14x2?9x3?2x4?5x5?3x ?15x?13x?6x?4x12345????11x1 ?9x2?2x3?5x4?7x5?5x ?7x?14x?16x?2x2345?1???2x1 ?5x2?12x3?11x4?4x5?100?200?300?400?500
2．判定下列方程是恰定方程组、欠定方程组还是超定方程组，并求其解。
?6x1?9x2?14x3?11x4?5x5?68?x?14x?7x?15x?6x?294?12345???2x1?x2?7x3?12x4?x5??441??6x1?11x2?11x3?9x4?13x5?103
3．用网孔电流法求如下电路的各支路电流。
《Matlab语言及其在电子信息科学中的应用》实验指导书
10
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
《Matlab语言及其在电子信息科学中的应用》实验指导书
11
第二题：
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]; X=A\\B
结果： X =
0 -26.3759 -5.4126 -38.3960 -8.2392
第三题：
A=[1,1,-1,1;0,20,40,0;60,0,40,0;0,0,40,40]; B=[0;10;50;-40]; X=A\\B X =
0.7857 0.3571 0.0714 -1.0714
《Matlab语言及其在电子信息科学中的应用》实验指导书
12
实验五：函数编写与程序设计
一、实验目的
1．掌握函数的编写规则。 2．掌握函数的调用。
3．会用Matlab程序设计实现一些工程算法问题。
二、实验内容
1．编写一个[y,y1,y2]=mwave(f1,m1,f2,m2)函数，实现以下功能，并绘出y1、y2、y在t∈[0，2π]区间500个样点的图形。（其中调用参数2 ≤ f1、f2 ≤ 20 HZ；0.5 ≤ m1、m2 ≤ 2）
y1?m1sin(2?f1t);y2?m2sin(2?f2t); y?y1?y22．程序设计：相传古代印度国王要褒奖他的聪明能干的宰相达依尔（国际象棋发明者），问他要什么？达依尔回答：“陛下只要在国际象棋棋盘的第一个格子上放一粒麦子，第二个格子上放二粒麦子，以后每个格子的麦子数都按前一格的两倍计算。如果陛下按此法给我64格的麦子，就感激不尽，其他什么也不要了。”国王想：“这还不容易！”让人扛了一袋麦子，但很快用光了，再扛出一袋还不够，请你为国王算一下共要给达依尔多少小麦？（1袋小麦约1.4×108粒）。 3．程序设计：公元前五世纪我国古代数学家张丘建在《算经》一书中提出了“百鸡问题”：鸡翁一值钱五，鸡母一值钱三，鸡雏三值钱一。百钱买百鸡，问鸡翁、母、雏各几何？
三、设计提示
1．函数名和函数文件名应相同；对调用参数的取值范围要检验是否符合要求，如不符合要求，应给出出错信息（用if和error函数实现）。
2．程序设计——“百鸡问题”答案不唯一。提示：设x：鸡翁数，则x的范围：0～19；y：鸡母数，则y的范围：0～33；z：鸡雏数，则z的范围：0～100。
四、实验报告要求
1．编写实现第二节实验内容（1）所使用的函数文件。 2．程序设计用M文件编写，并记录执行结果。“百鸡问题”答案不唯一，要给出所有答案。
3．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。
五、参考答案
《Matlab语言及其在电子信息科学中的应用》实验指导书