内容发布更新时间 : 2024/11/15 2:41:36星期一 下面是文章的全部内容请认真阅读。
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语言及其在电子信息科学中的应用》实验指导书