matlab语言实验指导书

matlab语言实验指导书
实验一matlab集成环境及数值运算
一、实验目的
1.熟悉启动和退出matlab的方法；
2.熟悉matlab的集成环境；
3.掌握建立矩阵的方法；
4.掌握MATLAB中各种表达式的编写规则和常用函数的使用。

二、实验内容
1.首先找到以下表达式的值，然后显示MATLAB workspace的用法并保存所有变量。

2sin85?（1）z1?21?e（2）z2?1?2i??21ln(x?1?x2)，其中
x250.45e0.3a?e?0.3a0.3?asin(a?0.3)?ln,a??3.0,?2.9,?,2.8,2.9,3.0（3）
z3?22提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。

?t20?t?1?2（4）z4??t?11?t?2，其中t=0:0.5:2.5
t2？2t？12? T3.提示：用逻辑表达式求分段函数的值。

2.已知：
12344131，b??20a??3478733?27???3657??求如下表达式的值：
（1） A+6*B和A-B+I（其中I是单位矩阵）（2）A*B和A*B
（3）a^3和a.^3（4）a/b和b\\a（5）[a,b]和[a([1,3],:);b^2]3.设有矩阵a和b
1234516 30? 67??? 891017-69 1112131415?, B023-4?
A.1617181920970 1311?? 2122232425??? 4.
（1）求他们的乘积c（2）将矩阵c的右下角3×2子矩阵赋给d（3）查看matlab工
作空间的使用情况。

4.完成下列操作：
（1）查找[100999]之间可被21除的数字数。

提示：先利用冒号表达式，再利用find和length函数。

（2）建立一个字符串向量，删除其中的大写字母。

提示：利用find函数和空矩阵。

实验二 matlab程序设计
一、实验目的
1.掌握建立和执行m文件的方法；
2.掌握利用if和switch语句实现选择结构的方法；
3.掌握利用for和while语句实现循环结构的方法；4,熟悉利用向量运算来代替循环操作的方法。

二、实验内容
1.求分段函数的值。

x2？十、6x？0&x？？3.Yx2？5倍？60? 十、5&x？2&x？三
x2x1other用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。

2.输入百分制分数，输出分数a、B、C、D、E，其中90~100分为a，80~89分为B，70~79分为C，60~69分为D，60分以下为E。

要求：
（1）分别用if语句和switch语句实现。

（2）输入百分制分数后，判断分数的合理性，对不合理的分数输出错误信息。

3.硅谷公司员工薪酬计算方法如下：
（1）工作时数超过120小时者，超过部分加发15%。

（2）工作时数低于60小时者，扣除700元。

（3）其余按每小时84元计发。

尝试编程，根据输入的工作编号和员工的工作时数计算应付工资。

4.根据
26111，求π的近似值。

当n分别去100,1000,10000时，结果是多少？22212n要求：分别用循环结构和向量运算来实现。

1115.根据y?1，求：
352n？1（1）y<3时的最大n值。

Y（5712）的值对应于991.44的值。

6、考虑下面的迭代公式：
xn?1?ab?xn其中a、b为正的常数。

（1）编写一个程序来查找迭代的结果。

迭代的终止条件是xn？1.xn？10? 5.迭代初始值x0=1.0，迭代次数不超过500次。

bb24a（2）如果迭代过程收敛于r，那么r的准确值是，当(a,b)的值取(1,1)、(8,3)、
2（10,0.1），分别将迭代结果与精确值进行比较。

实验三文件操作
一、实验目的
1.理解函数文件的概念；
2.掌握matlab函数的定义和调用方法；
3.掌握函数的递归调用。

二、实验内容
1.定义函数文件，找到给定复数的索引、对数、正弦和余弦，并在命令文件中调用函
数文件。

2.一物理系统可用下列方程组来表示：
m1cos？？msin？？1.0 m10m20？罪余弦？？罪余弦？0 a1？？0 a？？mg？0 2 1.0 n1？？0 1?? n2？？m2g？输入M1、M2，然后从键盘θ中找到A1、A2、N1和N2的值。

其中G为9.8，输入θ的角度。

要求：定义一个求解线性方程组ax=b的函数文件，然后在命令文件中调用该函数文件。

3.一个自然数是素数，且它的数字位置经过任意对换后仍为素数，则称是绝对素数。

例如13是绝对素数。

试求所有两位绝对素数。

要求：定义一个判断素数的函数文件。

4.编写一个函数文件，查找小于任何自然数n的fibnacci序列项。

fibnacci序列定
义如下：
f11f21fffn1n2nn1n2
N25.首先用函数的递归调用定义一个函数文件
1++。

kkkk?1k?1k?121005010?ii?1nm，然后，调用该函数文件求：
实验四 matlab绘图
一、实验目的
1.掌握绘制二维和三维图形的常用功能；
2.掌握绘图的辅助操作；
3.掌握图形对象属
性的基本操作；
4.掌握利用图形对象进行绘图操作的方法。

二、实验内容
3sinx??1.设y??0.5?cosx，在x=0~2π区间取101点，绘制函数的曲线。

2?1?x??2.
已知y1?x2,y2?cos(2x),y3?y1?y2,完成下列操作：（1）在同一坐标下用不同的颜色和线
型绘制三条曲线；（2）以子图形式绘制三条曲线；
（3）分别用条形图、梯形图、条形图和填充图绘制三条曲线。

3.已知
xe2y1ln(x1x22x0x0
以-5为间隔绘制函数曲线≤ 十、≤ 5.
4.建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在
按下鼠标器的左键之后显示出leftbuttonpressed字样。

5.先用默认属性绘制曲线y？X2e2x，然后通过图形句柄操作更改曲线的颜色、线型和线宽，并使用文本对象向曲线添加文本注释。

6.利用曲面对象绘制曲面v(x,t)?10e?0.01xsin(2000?t?0.2x??)，要求与上题相同。

实验五 matlab符号运算
一、实验目的
1.掌握定义符号对象的方法；
2.掌握符号表达式的运算法则以及符号矩阵运算；
3.掌握求符号函数极限及导数的方法；
4.掌握求符号函数定积分和不定积分的方法。

二、实验内容
1.给定x=6和y=5，使用符号表达式来求Z？十、13? 十、Y
提示：定义符号常数x=sym(‘6’),y=sym(‘5’)。

2．分解因式。

（1） x4？y4（2）51353。

简化表达式。

（1）sin?1cos?2?cos?1sin?24.已知
010 100 A.P010 A.dp1？？1002 001??? 101??? gbehc？F我
4x2?8x?3（2）
2倍？1.完成以下操作：
(1)b?p1?p2?a（2）b的逆矩阵并验证结果（3）包括b矩阵主对角线的下三角阵（4）
b的行列式值5.用符号方法求下列极限或导数。

x（esinx？1）？2（etanx？1）？？arccosx（1）lim（2）lim3？十、0x？？1sinxx？
1.axt3？1.cos（2x）d2ad2a（3）y分别找到y'，y“（4）已知的a和2。

？，
xdxdtdt?tcosxlnx?（5）已知f(x,y)?(x?2x)e6.用符号方法求下列积分。

（1）?2?x2?y2?xy?y?2f，求、
十、十、yx？0，y？一
dxdx（2）?(arcsinx)21?x21?x4?x8??（3）
0x2？1dx（4）？ex（1？ex）2dx4x？10ln2实验6 Simulink模拟
一、实验目的
1.熟悉Simulink的运行环境；
2.掌握绘制系统模型的方法；
3.对简单系统所给出的数学模型能转化为系统仿真模型并进行仿真分析。

二、实验内容
1.利用simulink仿真下列曲线，去ω=2π。

1111x（？t）？罪Tsin3？Tsin5？Tsin7？Tsin9？T
35792.建立图1示的系统模型并进行仿真。

改变gain模块的增益，观察scope显示
波形的变化。

图1
3.将图1所示的scope模块换成output模块。

在configurationparameters对话框
中把时间和输出作为返回变量，分别设置变量名t和y。

运行仿真并用绘图命令画出曲线
t-y。

4.建立如图2所示的系统模型并进行仿真。

改变slidergain模型的增益，观察X-Y
波形的变化。

用两个输出模块替换图2中所示的xygraph模块。

在“配置参数”对话框中，将时间和输出作为返回变量，并分别设置变量名T和[Y1，Y2]。

运行模拟并使用绘图命令
绘制曲线t-y1、t-y2和y1-y2。

图2
5.建立如图3所示的系统模型并进行仿真。

图3
6.系统的微分方程为：
x’(t)=-4x(t)+2u(t)
U（T）是一个振幅为1、角频率为1rad/s的方波输入信号。

建立了系统模型并进行了仿真。