上机说明系统仿真.doc

上机说明：
按上机要求完成实验一和实验二的题目，其余时间自由安排练习书上的例题，题目后面标有“作业”的，做为本次课的作业。完成时需要文字说明的用WORD编辑，WORD文档以“学号-姓名-作业1”命名；不需要文字说明的，给出题目中要求的程序（要求程序能正常运行），在命令窗口要输入的语句，写在“学号-姓名-作业1”文档中（标明题号），WORD文档及M文件压缩成压缩文件，以“学号-姓名-作业1”命名。发送给邮箱：wmq_117@
实验一 Matlab语言基本操作
上机实验内容：
1、如何在命令窗口输入命令函数并执行，主要体会“；”和“，”的不同意义；
2、通过帮助文件help查单相关函数的使用方法，并完成相应题目；
3、练习文本文件和函数文件的建立与使用，体会两者的不同；
4、在命令窗口中敲入Demo，查看相关范例。
相应题目：
1、在命令窗口中输入如下语句，查看显示结果
（1）>> x1=1:0.1:1.2
>> y1=2*x1
（2）>> x2=1:0.1:1.2；
>> y2=2*x2
（3）>> x3=1:0.1:1.2；
>> y3=2*x3；
（4）如果想查看x2，y1的内容如何实现？
>> x2 回车
>> y1 回车
（5）>> who 回车
（6）>> whos 回车
（7）>> save mydata x1 y1 x2 y2

>> load mydata x1 y1 x2 y2
>> y1 回车
>> x1 回车
2、某二阶欠阻尼系统的单位阶跃响应为：

试绘制响应曲线。
以文本文件形式输入：（练习进入文本文件编辑器的几种方法）
存盘后，运行（练习两种运行方法）
可通过help plot查看plot绘图功能的使用。
回到命令窗口查看t和c两个变量内容“；”和回车键的显示不同。
在命令窗口中敲入>> mp=max(c) 回车
显示结果：mp =1.1623 ％显示的是曲线中的峰值
3、函数文件的例子（作业）
（1）在m－file编辑器中输入：
%sq1.m
function[c]=sq1(a,b)
c=sqrt((a.^2)+(b.^2));
存盘，回到命令窗口中键入：
>> a=4;
>> b=3;
>> [c]=sq1(a,b)
c =5
说明此函数的功能是什么？
（2）试编写一个函数文件实现：

在x1＝3，x2＝5，x3＝1时的值。

实验二 熟悉MATLAB语言（运算）
一、熟悉MATLAB语言中矩阵的运算。
1、，求
2、，求，，，
3、，，求，，
4、规律填数：试输入1到100以5为公差的等差数列。
5、已知：A=[1 4 7;10 13 16;19 22 25]，k＝5。试求：B＝A+k、C=A-k、D＝A－k；
6、已知：A=[1 4 7;10 13 16;19 22 25]，k＝6。试求：B＝k.*A、C=k*A、D＝A*k、E＝A./k、F=A.\k；
7、已知：A=[1 4 7;10 13 16;19 22 25]，B=[6 6 6;3 3 3;1 1 1]。试求：C＝A.*B、D=B.*A、E＝A*B、F=B*A、G＝A./B、H=A.\B；
8、试计算矩阵a=[6 8 0;5 3 2;1 4 -3]与b=[1 1 -1 3;-1 -1 2 2;2 5 2 4;1 2 3 2]行列式的值。
9、已知：a =[1 1 -1 3;-1 -1 2 22 5 2 4;1 2 3 2]，试求a矩阵的转置，对a矩阵进行上下、左右、旋转转换。
10、求矩阵a＝[1 3 5;2 4 6;7 9 13]的特征值、特征向量及矩阵的秩。
程序：
11、在命令窗口实现以下操作：（作业）
（1）生成一个三行四列的单位矩阵A；
（2）将上述A阵中加入一新行（第四行）元素为：[8 4 6 7]，将新矩阵赋给矩阵B；
（3）在B矩阵中的第二行第三列元素改为10，将新矩阵赋给C;
（4）给C矩阵加入一新列元素4，5，3，6，分为两种情况：
①插入到C矩阵的最后一列，新矩阵赋给矩阵D；
②插入到C矩阵的任何列中，新矩阵赋给矩阵E。
二、练习MATLAB语言中控制语句及绘图函数的使用。
1、饱和非线性方程为：
试编写程序并绘制图形。
2、绘制图形，要求用红色点划线，在横轴上标注“x axis”，在纵轴上标注“y axis”，在图形适当地方写上字符串“y=sinx+cosx”，写出程序（以M文件写程序）。
3、在同一坐标系下绘制两条曲线，，要求用绿色点划线绘制，用黄色星号线绘制，在坐标系正上方标注题目名称“Figure Example”，写出程序（以M文件写程序）。（作业）
4、验证教材相关例题