************************ MATLAB上机指导书
************************
昆明理工大学机电学院
刘孝保
2013年3月
实验三符号计算
一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。
1.findsym：帮助我们获取系统定义的自变量
f= sym('sin(a*x+b*y)'); findsym(f)
2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系
数获得符号表达式）
[n,d]=numden(sym('x*x+y')+sym('y^2'))
p=sym('2*x^3+3*x^2+4'); sym2poly(p)
x=[2,3,0,4]; poly2sym(x)
3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式
分解；simplify: 对表达式化简
syms x y; collect(x^2*y+y*x-x^2-2*x)
collect((x+y)*(x^2+y^2+1), y)
syms x y; expand((x-2)*(x-4))
syms x;horner(x^3-6*x^2+11*x-6)
syms x;factor(x^3-6*x^2+11*x-6)
syms x;simplify((x^2+5*x+6)/(x+2))
4. finverse ：求得符号函数的反函数。
syms x y; finverse(1/tan(x))
f= x^2+y; finverse(f,y)
finverse(f)
https://www.doczj.com/doc/f53348457.html,pose 求符号函数的复合函数
syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g)
6. subs :表达式替换。
syms a b;subs(a+b,a,4)
subs(a+b,4)
subs(cos(a)+sin(b),{a,b},{sym('alpha'),2}) subs('x^2+2*y',{'x','y'},{3,4})
7.极限：limit
syms x t;limit(sin(x)/x)
limit((x-2)/(x^2-4),2)
limit((1+2*t/x)^(3*x),x,inf)
limit(1/x,x,0,'right')
8. 微分：diff
syms a b c x;f=sym('a*x^2+b*x+c')
diff(f)
diff(f,2)
diff(f,a)
diff(f,a,2)
9. 积分：int
syms a b c x
f=sym('a*x^2+b*x+c')
int(f)
int(f,x,0,2)
int(f,a)
int(int(f,a),x)
10级数：symsum,taylor
syms k
symsum(1/k,k,1,inf)
symsum(1/(k*(k+1)),k,1,inf)
syms x
taylor(sin(x),x,10)
subs(ans,x,pi/2)
11.solve 解符号方程式
f=sym('a*x^2+b*x+c');solve(f)
f1=sym('x+y+z-10=0');
f2=sym('x-y+z=0');
f3=sym('2*x-y-z+4=0');
solve(f1,f2,f3);
[x,y,z]=solve(f1,f2,f3)
12解微分符号方程式：dsolve,Dy代表dy/dt, D2y代表d2y/dt2 。dsolve('Dy=5')
dsolve('Dy=x','x')
dsolve('D2y=1+Dy')
dsolve('D2y=1+Dy','y(0)=1','Dy(0)=0')
[x,y]=dsolve('Dx=x+y, Dy=2*x','x(0)=0','y(0)=1')
二、使用部分：利用matlab求解，并记录运行结果
1、计算符号表达式
()cos()sin()
f x x x x
=+-在1
xπ
=-处的值
2、设x为符号变量，
42
()21
f x x x
=++，32
()635
g x x x x
=+++，试进行如下运算。
（1）
()() f x g x
+，
（2）
()() f x g x
?，
（3）对
()
f x进行因式分解
3、指出下面各条语句的输出结果>> f=sym('2*u');
>> subs(f,'u',2)
>> f2=subs(f,'u','u+2')
>> a=3;
>> subs(f2,'u',a+2)
>> subs(f2,'u','a+2')
>> syms x y
>> f3=subs(f,'u',x+y)
>> subs(f3,[x,y],[1,2])
>> subs(f3,[x,y],[x+y,x+y])
>> f=sym('2*u');
4、化简
5、 因式分解：
6555234-++-x x x x
6、求极限 ：
7、已知f(x)=ax2 +bx+c,求f(x)的微分
8、已知f(x)=ax2 +bx+c,求f(x)的积分。
9、
三、思考题：
圆形凸轮机构如下图所示，其中凸轮半径为100mm ，滚子半径20mm ，旋转偏心距为
50mm ，顶针长度为120mm ，以旋转中心作为坐标原点，请用符号函数表示出顶针上端顶点的位移、速度和加速度表达式，求出当旋转角度=0:5:360度时对应的位移、速度、加速度值，并试图画出位移图速度图、加速度图。
实验四 矩阵运算
一、实验目的：掌握基本的矩阵运算及常用的函数。 二、实验内容：
??????=654321a ??????-=531142b ??
???
?????-=201c ????
??????=063258741d
1、下列运算是否合法，为什么？如合法，结果是多少？ (1) result1 = a' (2) result2 = a * b
(3) result3 = a + b
()f x =
lim()x
x x a x a
→∞+-11111111............23k 22x33x4k(k+1)++++++++++和的和。
(4) result4 = b * d
(5) result5 = [b ; c' ] * d (6) result6 = a . * b (7) result7 = a . / b (8) result8 = a . * c (9) result9 = a . \ b (10) result10 = a . ^2 (11) result11 = a ^2 (12) result11 = 2 . ^ a
2、用MATLAB 求下面的的方程组。
(1)????
?
???????-=??????????????????????
???----01741323151122231592127
4321x x x x (2)???????=-++=--=-++=++5
6533332821
w z y x w y x w z y x z y x
3、已知????
?
????
???----=1323151122231592127
A (1)求矩阵A 的秩(rank)
(2)求矩阵A 的行列式(determinant) (3)求矩阵A 的逆(inverse)
(4)求矩阵A 的特征值及特征向量(eigenvalue and eigenvector) 4、关系运算与逻辑运算 已知a=20,b=-2,c=0,d=1
(1) r1 = a > b
(2) r2 = a > b & c > d (3) r3 = a == b* (-10) (4) r4 = ~b | c
三、思考题
10991010
10
22222++++==
---=∑ n n y ,求y=？（用format long 查看y 的值）
预习 MATLAB 图形绘制基础
一、实习目的
1、进一步熟悉M文件调试过程
2、熟练掌握MATLAB二维曲线的绘制
3、掌握图形的修饰
二、实习内容
1、基本命令：plot（数据作图、函数作图）, ezplot（简易作图），fplot(简易作图)
（1）plot(x,y,’s’) s用来表示色彩、线型、点型
例3.1 用plot绘制y=sinx/x的图形
x=-2*pi:pi/10:2*pi;
y=sin(x)./x; %准备数据
plot(x,y) %用plot作图
注意x=0处的图形
例3.2 x=-2*pi:pi/10:2*pi;
x=x+(x==0)*eps;
y=sin(x)./x; %（准备数据,处理分母零点))
plot(x,y) %（用plot作图）：
例3.3 画出()π,0包罗线。
t=(0:pi/100:pi);
y1=sin(t);y=-sin(t);
y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;
y3=sin(t3).*sin(9*t3);
plot(t,y1,'r:',t,y,'r:',t,y2,'b',t3,y3,'bo'),axis([0,pi,-1,1])
gtext(‘该点数值x’)
运行结果见图3-1。
（2）ezplot用作图：
ezplot(‘sin(x)/x’,[-4*pi,4*pi])
运行结果如图3-2
（3）fplot用作图：
fplot(‘sin(x)/x’,[-2*pi,2*pi])
图
3-1
图3-2
2、基本绘图控制参数：设置线型、线色、数据点形 （1）、曲线点型、线色允许设置值
（2（3（4）、图形标注：坐标轴名(label)、图形标题（title ）、图例
(legend)、图形说明（text,gtext ）,
格式如下：
title(‘s’) 书写图形标题 xlabel (‘s’) 横坐标轴名 ylabel(‘s’) 纵坐标轴名
legend(‘s1’,’s2’,…) 绘制曲线所用线型、色彩、或数据点型图例 text(x1,y1,’s’) 在图面（x1,y1）坐标处书写字符注释 gtext(‘s’) 用鼠标选择书写字符注释的位置 （5）、常用坐标控制命令
（hold on 使当前轴及图形保持而不被刷新，准备接受此后将绘制的新曲线。 hold off 使当前轴及图形不再具备不被刷的性质。 hold 是否具备不被刷的性质的双向开关。
subplot(m,n,k) 图形窗中有m×n幅子图，k是子图的编号。
3、三维曲线绘图
（1）基本命令：
plot3(x,y,z,’s’) % s为线色、线形、点形
plot3(x1,y1,z1,’s1’,…. xn,yn,zn,’sn’)
例3.5
t=-2*pi:0.01:2*pi;
x=2*cos(t);y=2*sin(t),z=3*t; %(准备数据)
plot3(x,y,z,’r:’); %（用plot3作图）
（2）三维曲面作图
1）三维网格命令（mesh）
例3.6 x=-1:0.05:1;y=x;
[x,y]=meshgrid(x,y);
z=x.^2+2*y.^2;
mesh(x,y,z)
运行结果如图3-4。
图3-3
2）三维表面命令surf
例3.7
x=-1:0.05:1;y=x;
[x,y]=meshgrid(x,y);
z=x.^2+2*y.^2;
surf(x,y,z)
运行结果如图3-5。
图3-4
3）基本三维命令的几个改进命令
meshc(带等高线的网线图)
surfc(带等高线的曲面图)
4）、三维视图的可视效果控制
view, rotate3d, rotate
5)、图形的透视（hidden）
图形颜色的遮掩shading。
hidden on 使图形消隐（透视）。
hidden off 关闭消隐（此为默认状态）。
shading 用于曲面颜色不匀时的均衡处理，使用格式为。
shading interp 对网眼内采用匀色处理，使色彩自然连贯。
shading faceted 对网眼颜色不做特殊处理，但加深网线黑色，这种格式有较强的表现力，是默认格式。
4、球面（sphere）及柱面（cylinder）的表达
（1）sphere(n)绘制分格线条数为n的单位球面。
例3.8 sphere(50) %绘制分格线条数为50的单位球面，运行结果如图3-6。
图3-5
（2）cylinder(r,n)绘制半径为r分格线条数为n的旋转柱面（曲面），纵坐标为0到1。R为母线的描述。
例3.9 cylinder(5,30)，运行结果如图3-7。
（3）[x,y,z]= sphere(n), [x,y,z]= cylinder(r,n)，分别为读取球面及柱面的数据。
图3-6
5、特殊图形：直方图hist、饼图pie、条形图bar、矢量图feather、等高线contour。
三、实习步骤
（1）plot指令的使用
在commmand窗口键入以下指令：
X=1:10
Y=[1 2 3 4 5 6 7 8]
Z=1:-1:-10
Plot(x)
Plot(y)
Plot(x,z)
观察指令窗及图形窗中的结果。
（2）subplot指令的使用；采样点数对绘图的影响的理解；title、axis指令的使用键入所以下指令：
t1=(0:11)/11*pi;
y1=sin(t1).*sin(9*t1);
t2=(0:100)/100*pi;
y2=sin(t2).*sin(9*t2);
subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图(1)')
subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图(2)')
subplot(2,2,3),plot(t1,y1,t1,y1,'r.')
axis([0,pi,-1,1]),title('子图(3)')
subplot(2,2,4),plot(t2,y2)
axis([0,pi,-1,1]),title('子图(4)')
观察指令窗及图形窗中的结果。
（3）图形交互指令的使用
①、在（2）的四个子图上用ginput指令分别获取其峰值；
②、在（2）的四个子图中的子图1上用gtext指令将“不太完整的离散图形”字符串
放置到图形上。在子图2上用gtext指令将“完整的离散图形”字符串放置到图形上。在子图3上用gtext指令将“失真的连续图形”字符串放置到图形上。在子图4上用gtext
指令将“正确的连续图形”字符串放置到图形上。 （4）其它绘图指令的使用及精细指令的使用
键入所以下指令：
clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t); tt=t(find(abs(y-1)>0.05));ts=max(tt); plot(t,y,'r-','linewidth',3) axis([-inf,6*pi,0.6,inf])
set(gca,'xtick',[2*pi,4*pi,6*pi],'ytick',[0.95,1,1.05,max(y)]) grid on
title('\ity=1-e^{-\alphat}cos{\omegat}') text(13.5,1.2,'\fontsize{12}{\alpha}=0.3') text(13.5,1.1,'\fontsize{12}{\omega}=0.7') hold on;plot(ts,0.95,'bo','markersize',10);hold off cell_string{1}='\fontsize{12}\uparrow';
cell_string{2}='\fontsize{16}\fontname{隶书}镇定时间'; cell_string{3}='\fontsize{6}';
cell_string{4}=['\fontsize{14}\rmt_{s}=' num2str(ts)]; text(ts,0.85,cell_string)
xlabel('\fontsize{14}\bft\rightarrow') ylabel('\fontsize{14}\bfy\rightarrow') 观察指令窗及图形窗中的结果。
（5）提高部分
本实习要求事先编写好解决下面问题的程序，然后上机输入程序并调试运行程序。 A 、绘制曲线s a (t)=sin(x)/x,其中]10,10[ππ-∈x （注意：x=0处需要特别处理） B 、画出t e
y t
cos 5.0-=的曲线（t 的取值范围是0到2π），在x 轴上标上0.5，1.2和
4.5三个刻度。在y 轴上标0.2，1.0，3.5三个刻度。用xlabel 在横轴下标上t;用ylabel
在y 轴左侧标y 在图上用“小红圈”标出第一次使y=0.6的那点位置，并要求在该点旁边用12号字，斜体标出其横坐标。
四、实习要求：
1、进一步熟悉和掌握MA TLAB 的编程及调试。
2、掌握二维图形的绘制。
3、掌握图形交互指令的使用。
实验一 Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句； 3、熟悉Matlab绘图命令及基本绘图控制； 4、熟悉Matlab程序设计的基本方法。 二、实验内容 1、帮助命令 2、矩阵运算 （1）矩阵的乘法和乘方 已知A=[1 2;3 4]:B=[5 5;7 8]:求A^2*B ( 2 )矩阵除法 已知A=[1 2 3:4 5 6:7 8 9]:B=[1 0 0:0 2 0:0 0 3],求矩阵左除A\B,右除A/B。 ( 3 )矩阵的转置及共轭转置 已知A=[5+i,2-i,1;6*i,4,9-i],求该复数矩阵的转置A',共轭转置A.' ( 4 )使用冒号选出指定元素 已知：A=[1 2 3:4 5 6:7 8 9];求A中第3列前2个元素;A中第2、3行元素。 ( 5 )方括号[] 用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式 （1）求多项式P(x)=x3-2x-4的根 ( 2 )已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4].,构造多项式，并计算多项式值为20的解。 4、基本绘图命令 ( 1 ) 绘制余弦曲线 ] 2,0[ ), cos(π ∈ =t t y 。 ( 2 ) 在同一坐标系中绘制曲线 ] 2,0[ ), 5.0 sin( ), 25 .0 cos( yπ ∈ - = - =t t y t 5、基本绘图控制 绘制 ] 4,0[π 区间上的y=10sint曲线，并要求： （1）线形为点划线，颜色为红色，数据点标记为加号； （2）坐标轴控制：显示范围，刻度线，比例，网络线； （3）标注控制：坐标轴名称，标题，相应文本。 6、基本程序设计 （1）编写命令文件：计算1+2+....+n<2000时的最大n值； (2)编写函数文件：分别用n和which循环结构编写程序，求2的0到n次幂的和. n=input('输入正数n:') ji=1: for i=1:n; ji=ji+2^i: end ji input('输入正数n:') ji-1:i-1: While i<=n ji=ji+2^i; i=i+1; end ji (3)如果想对一个变量x赋值，当从键盘输入y或Y时，x自动赋为1；当从键盘输入n或N时，x自 动赋为0；输入其他字符时终止程序。 k=input('shuruX:'.'s'): if k=='y' k=='Y' x=1 ; else k=='n' k=='N' x=0; else ruturn end >> n=input('输入正数n:') 输入正数n:20 n =20
信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：
实验一MATLAB 工作环境熟悉及简单命令的执行 一、 实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。 二、 实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单 （Menu ）、工具条（T oolbar ）、命令窗口(Command Window)、历史命令窗口、工作 空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、启动MATLAB ，熟悉MATLAB 的桌面。 2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。 （1）（365-52?2-70）÷3； （2）area=pi*2.5^2； （3）已知x=3，y=4，在MATLAB 中求z ： ()23 2y x y x z -= ； （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。 m1=???? ? ???? ???11514412679810115133216 执行以下命令 >>m1( 2 , 3 ) >>m1( 11 ) >>m1( : , 3 ) >>m1( 2 : 3 , 1 : 3 ) >>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) （5）执行命令>>help abs 查看函数abs 的用法及用途，计算abs( 3 + 4i )； （6）执行命令 >>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y) （7）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。 四、思考题 1、以下变量名是否合法？为什么？ （1）x2； （2）3col ； （3）_row ； （4）for ；
实验二 MATLAB矩阵及其运算 学号：3121003104 姓名：刘艳琳专业：电子信息工程1班日期：2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. （1）新建一个.m文件，验证书本第15页例2-1； （2）用命令方式查看和保存代码中的所有变量；
（3）用命令方式删除所有变量； （4）用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式
长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 （1）w=sqrt(2)*(1+0.34245*10^(-6)) （2）x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; （3）y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; （4）z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];
4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20]，对其进行如下操作：（1）输出A在[ 7, 10]范围内的全部元素； （2）取出A的第2，4行和第1，3，5列； （3）对矩阵A变换成向量B，B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]； （4）删除A的第2，3，4行元素； (1) (2)
实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i
1.0000 1.4142 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B 解：A=[1 2;3 4 ]; B=[5 5;7 8 ]; A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解：A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B （3）矩阵的转置及共轭转置
Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容：
实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量： 以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 的近似面积，并与它的精确值41288km2比较。
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因 c b a 、、的不同取值而定） ，这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 的值，调用该函数后，
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名： 学号： 指导教师：
一、实验名称 实验二 线性方程组求解和函数的数值逼近 二、实验目的 通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。 实验涉及的核心知识点：病态方程求解、矩阵分解和方程组求解、Lagrange 插值。 实验重点与难点：算法设计和MATLAB 编程 三、实验内容 1. 对高阶多项式 ()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 2. 对2 20n =，生成对应的Hilbert 矩阵，计算矩阵的条件数；通过先确定解获得常向量 b 的方法，确定方程组 ()n H x b = 最后，用矩阵分解方法求解方程组，并分析计算结果。 3. 对函数 ()2 1 125f x x = + []1,1x ∈- 的Chebyshev 点 ()()21cos 21k k x n π ?? -= ? ?+? ? ，1,2,,1k n =+ 编程进行Lagrange 插值，并分析插值结果。 四、实验数据及结果分析 1. 对高阶多项式
()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 p=[1,-1]; for i=2:20 n=[1,-i]; p=conv(p,n); % 求多项式乘积 end m=zeros(1,21); % m 的最高次幂为20，有21项 hold on x=1:20; d=[-1,0,0.1,0.5,1]; for i=1:5 delt=d(i); m(2)=delt; y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); end title('方程p(x)=0的解与扰动量delt 的关系') legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1') 2468101214161820 010 20 30 40 50 60 方程p(x)=0的解与扰动量delt 的关系 delt=-1delt=0delt=0.1delt=0.5delt=1
************************ MATLAB上机指导书 ************************ 昆明理工大学机电学院 彭用新 2015年3月
实验三符号计算 一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。 1.findsym：帮助我们获取系统定义的自变量 f= sym('sin(a*x+b*y)'); findsym(f) 2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系 数获得符号表达式） [n,d]=numden(sym('x*x+y')+sym('y^2')) p=sym('2*x^3+3*x^2+4'); sym2poly(p) x=[2,3,0,4]; poly2sym(x) 3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式 分解；simplify: 对表达式化简 syms x y; collect(x^2*y+y*x-x^2-2*x) collect((x+y)*(x^2+y^2+1), y) syms x y; expand((x-2)*(x-4)) syms x;horner(x^3-6*x^2+11*x-6) syms x;factor(x^3-6*x^2+11*x-6) syms x;simplify((x^2+5*x+6)/(x+2)) 4. finverse ：求得符号函数的反函数。 syms x y; finverse(1/tan(x)) f= x^2+y; finverse(f,y) finverse(f) https://www.doczj.com/doc/f53348457.html,pose 求符号函数的复合函数 syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g) 6. subs :表达式替换。 syms a b;subs(a+b,a,4)
数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)
cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);
昆明理工大学信息工程与自动化学院学生实验报告 （ 2011——2012 学年 第一学期 ） 课程名称：控制系统计算机辅助设计 开课实验室：信自楼 234 2011年10月28日 年级、专业、班 学号 姓名 成绩 实验项目名称 实验二 控制系统分析 指导教师 胡蓉 教师评 语 该同学是否熟悉实验内容： A.熟悉□ B.比较熟悉□ C.不熟悉□ 该同学的实验能力： A.强 □ B.中等 □ C.差 □ 该同学的实验是否达到要求 ： A.达到□ B.基本达到□ C.未达到□ 实验报告是否规范： A.规范□ B.基本规范□ C.不规范□ 实验过程是否详细记录： A.详细□ B.一般 □ C.没有 □ 注：5个A 为优，5个B 为中，介于二者间为良，5个C 为不及格，3个Ｂ以上为及格。 教师签名： 年 月 日 实验二 控制系统分析 一、 实验目的 1. 掌握如何使用Matlab 进行系统的时域分析。 2. 掌握如何使用Matlab 进行系统的频域分析。 3. 掌握如何使用Matlab 进行系统的根轨迹分析。 4. 掌握如何使用Matlab 进行系统的稳定性分析。 二、 实验内容 1.时域分析 （1）典型二阶系统传递函数为：
当ζ=0.7, ω 取2、4、6、8、10、12的单位阶跃响应。 n 程序为：>> num1=4;den1=[1,2.8,4];sys1=tf(num1,den1); >> num2=16;den2=[1,5.6,16];sys2=tf(num2,den2); >> num3=36;den3=[1,8.4,36];sys3=tf(num3,den3); >> num4=64;den4=[1,11.2,64];sys4=tf(num4,den4); >> num5=100;den5=[1,14,100];sys5=tf(num5,den5); >> num6=144;den6=[1,16.8,144];sys6=tf(num6,den6); >> step(sys1,sys2,sys3,sys4,sys5,sys6); 运行单位阶跃响应结果图为： （2）典型二阶系统传递函数为： 当ω =6, ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。 n
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).
西安邮电学院 《Matlab》 实验报告 （二） 2011- 2012 学年第 1 学期 自动化 专业： 自动0903 班级： 学号： 姓名： 2011 年10月15日
实验二 MATLAB 的基本计算 一、实验目的 1．掌握建立矩阵的方法。 2．掌握MATLAB 各种表达式的书写规则以及常用函数的使用。 3．能用MATLAB 进行基本的数组、矩阵及符号运算。 4．掌握矩阵分析的方法以及能用矩阵求逆法解线性方程组。 二、实验设备及条件 计算机一台（带有MATLAB6.5或以上版本的软件环境）。 三、实验内容 1．利用diag 等函数产生下列矩阵。 ??????????-=032570800a ???? ? ?????=8040 507 2 b 2．利用reshape 函数将1题中的a 和b 变换成行向量。
3．产生一个均匀分布在（-5，5）之间的随即矩阵（10×2），要求精确到小数点后一位。 4．已知： ???????? ??-=765 3877 34434 12A ???? ? ?? ???--=7312 033 2 1 B 求下列表达式的值：
(1) B A K *611+=和I B A K +-=12（其中I 为单位矩阵） (2) B A K *21=和B A K *.22= (3) 331^A K =和3.32^A K = (4) B A K /41=和A B K \42= (5) ],[51B A K =和]2:);],3,1([[52^B A K =
5．下面是一个线性方程组： ???? ? ?????=????????????????????52.067.095.03216/15 /14 /15/14/13/14/13/12/1x x x (1) 求方程的解 (2) 将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相对变 6．利用randn 函数产生均值为0，方差为1的6×6正态分布随机矩阵C ，然后统计C 中大于-0.3，小于0.3的元素个数t
《数学实验》报告实验名称 Matlab 二维曲线绘图 2011年 5月
一、【实验目的】 学习Matlab 绘图的运用，学会制作二维曲线，三维图形的绘画。 二、【实验任务】 P79 第3，5，9题。 1，在同一图形窗口画三个子图…… 2，绘制圆锥螺线的图像并加各种标注…… 3，画三维曲面z=5-x^2-y^2与平面z=3的交线。 三、【实验程序】 1. >> clear >> x=-pi:pi/50:4*pi; y1=x.*cos(x); y2=x.*tan(1./x).*sin(x.^3); y3=exp(1./x).*sin(x); subplot(3,1,1) plot(x,y1,'r*'),grid on title('y1=xcosx') xlabel('x轴'),ylabel('y轴') axis([-pi pi -pi pi]) gtext('y1=xcosx'),legend('y1=xcosx') subplot(3,1,2),plot(x,y2,'b'),grid on title('y=xtan(1/x)sin(x^3)') gtext('y=xtan(1/x)sin(x^3)') legend('y=xtan(1/x)sin(x^3)') axis([pi 4*pi -2 2]) subplot(3,1,3),plot(x,y3,'y'),grid on title('y=exp(1/x)sinx') xlabel('x轴'),ylabel('y轴') gtext('y=exp(1/x)sinx') legend('y=exp(1/x)sinx') axis([1 8 -3 3]) 2. >> clear >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2.*t; plot3(x,y,z) title('圆锥螺线') xlabel('x轴'),ylabel('y轴'),zlabel('z轴') >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t);
数学软件与数学实验2013-2014学年度秋季学期期末试卷 专业：统计学 班级：11级2班 学号：20110723 姓名：晏静
一、按要求计算出下列表达式的值 (1)318, 3 162 53 ?? + ? ?? , 21 eπ+, 2.5 tg, 2 log15； (2)给出π的9位和e的10位近似值； (3)求658和4102的最大公约数及35和25的最小公倍数； (4)产生10个0与10之间随机数的一个表; (5)求虚数1453 i i i i +- -的实部,虚部,模,共轭,辐角。 (6)自己运用Table建立两个表，并进行表运算，如连接、并集、交、排序等操作。
二、因式分解 22212321332112322 1 22(1)()()()4;(2)21;x x x x x x x x x x x x x x x +++++---- 解： 三、解方程(组) 1234234124234-2+344-+-3(1)+31-73+3 x x x x x x x x x x x x x -=??=? ? +=??+=-? 65432(2)5232002000.x x x x x x -+--++= 四、求极限 () 20 (1)1sin ;(2);(3)56! ctg x n x n n n Lim x Lim n n →→∞ →∞++
(1) (2) (3) 五、求导数 32 22(1)()=ln(x+1+);(2)()=cos 2,; (3)=log (),Z . x f x x f f x e y x y Z xy x y y ???求的导数已知求求关于的二阶导 (1) (2) (3) 六、求下列定积分与不定积分： ()()()12201+sin ln 1+(1);(2);(3)sin (1+cos ) +1(1+)(2+-) x x dx dx x x x x x x ? ? ?2 2-(4)=0,=1,==.y D D x y y x I x e d σ??设是由直线围成的区域，计算的值 (1) (2)
实验一 MATLAB语言平台及基本运算 一、实验目的 1、熟悉MATLAB基本的操作界面。 2、掌握MATLAB中变量、数组、向量等对象的生成方法。 3、掌握矩阵的生成和基本运算方法。 4、掌握MATLAB中的常用绘图命令使用方法 二、实验设备 计算机，MATLAB语言环境 三、实验指导原理 1、常见数学函数 如：输入 x=[-4.85 -2.3 -0.2 1.3 4.56 6.75]，则 ceil(x)= -4 -2 0 2 5 7 fix(x) = -4 -2 0 1 4 6 floor(x) = -5 -3 -1 1 4 6 round(x) = -5 -2 0 1 5 7 系统的在线帮助 （1） help 命令： 1）.当不知系统有何帮助内容时，可直接输入 help以寻求帮助: >> help（回车） 2）.当想了解某一主题的内容时，如输入： >> help syntax （了解 Matlab 的语法规定） 3）.当想了解某一具体的函数或命令的帮助信息时，如输入： >> help sqrt （了解函数 sqrt 的相关信息）
（2） lookfor 命令 现需要完成某一具体操作，不知有何命令或函数可以完成，如输入： >> lookfor line （查找与直线、线性问题有关的函数） 3、常量与变量 系统的变量命名规则：变量名区分字母大小写；变量名必须以字母打头，其后 可以是任意字母，数字，或下划线的组合。此外，系统内部预先定义了几个有特殊意义和用途的变量，见下表： （1）数值型向量（矩阵）的输入 1）．任何矩阵（向量），可以直接按行方式输入每个元素：同一行中的元素用逗号（，）或者用空格符来分隔；行与行之间用分号（；）分隔。所有元素处于一方括号（[ ]内；例1： >> Time = [11 12 1 2 3 4 5 6 7 8 9 10] >> X_Data = [2.32 3.43；4.37 5.98] 2）．系统中提供了多个命令用于输入特殊的矩阵： 上面函数的具体用法，可以用帮助命令 help 得到。如：help zeros ZEROS Zeros array. ZEROS(N) is an N-by-N matrix of zeros. ZEROS(M,N) or ZEROS([M,N]) is an M-by-N matrix of zeros.等等 4、数组（矩阵）的点运算 运算符：+（加）、-（减）、./（右除）、.\（左除）、.^（乘方）， 例4： >> g = [1 2 3 4]；h = [4 3 2 1]； >> s1 = g + h, s2 = g.*h, s3 = g.^h, s4 = g.^2, s5 = 2.^h 5、矩阵的运算
实验二 M A T L A B 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器 （Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 21 3105168421 63105168421 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 编写一个函数，计算下面函数的值，给出标量x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x) 选择一些数据测试你编写的函数。 5. 编写一个函数求向量x 中元素的平均值、最大值、最小值、均方根值。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0
实验报告说明： matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。 第一次实验内容： 实验一 MATLAB 运算基础 一、实验目的 1.熟悉启动和退出MA TLAB 的方法。 2.熟悉MA TLAB 命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握MA TLAB 各种表达式的书写规则以及常用函数的使用。 二、实验内容 1.先求下列表达式的值，然后显示MA TLAB 工作空间的使用情况并保存全部变量。 （1）2 2sin 8511z e ?= + （2 ）12ln(2 z x =+ ，其中2120.45 5i +? ? =? ?-?? （3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02 a a e e z a a --= +=--- 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。 （4）2 2 2 01 41 1221 23 t t z t t t t t ?≤
12344347873657A -????=??????，131203327B -???? =????-?? 求下列表达式的值： （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 12345678910111213141516171819202122232425A ????????=????????， 30 161769 23497041311B ?? ?? -?? ??=-?? ????? ? （1） 求它们的乘积C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D （3） 查看MA TLAB 工作空间使用情况。 4.完成下列操作： （1）求[100，999]之间能被21整除的数的个数。 提示：先利用冒号表达式，再利用find 和length 函数。 （2）建立一个字符串向量，删除其中的大写字母。 提示：利用find 函数和空矩阵。 第二次实验内容： 实验三 选择结构程序设计 一、实验目的 1. 掌握建立和执行M 文件的方法。 2. 掌握利用if 语句实现选择结构的方法。 3. 掌握利用switch 语句实现多分支选择结构的方法。 4. 掌握try 语句的使用。 二 、实验内容
Matｌab数学实验一——matlab初体验 一、实验目的及意义 ［1] 熟悉MAＴＬＡB软件的用户环境； ［2] 了解MAＴLAB软件的一般目的命令； ［３] 掌握MATLＡB数组操作与运算函数； 通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。 二、实验内容 １.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务 根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。 （1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝
matlab实验指导书 实验一 MATLAB运算基础及矩阵分析与处理 一实验目的 1．熟悉启动和退出MATLAB的方法。 2．熟悉MATLAB 命令窗口的组成。 3．掌握建立矩阵的方法。 4．掌握MATLAB各种表达式的书写规则以及常用函数的使用。 5. 掌握生成特殊矩阵的方法。 6. 掌握矩阵分析的方法。 7. 用矩阵求逆解线性方程组。二实验仪器和设备装有以上计算机一台三实验原理 MATLAB 是以复杂矩阵作为基本编程单元的一种程序设计语言。它提供了各种矩阵的运算与操作，并有较强的绘图功能。基本规则 一般MATLAB命令格式为 [输出参数1，输出参数2，……]= 输出参数用方括号，输入参数用圆括号如果输出参数只有一个可不使用括号。 %后面的任意内容都将被忽略，而不作为命令执行，一般用于为代码加注 释。 可用↑、↓键来重现已输入的数据或命令。用←、→键来移动光标进行修改。所有MATLAB命令都用小写字母。
大写字母和小写字母分别表示不同的 变量。 常用预定义变量，如pi 、Inf、NaN、ans 矩阵的输入要一行一行的进行，每行各元素用空格或“，”分开，每行用 “；”分开。如 123 A456789 MATLAB书写格式为A=[1 2 3 ；4 5 6 ；7 8 9] 在MATLABZ中运行如下程序可得到A 矩阵 a=[1 2 3;4 5 6;7 8 9] a = 1 2 3 4 5 6 7 8 9 需要显示命令的计算结果时，则语句后面不加“；”号，否则要加“；”号。运行下面两种格式可以看出它 们的区别： a=[1 2 3;4 5 6;7 8 9] a=[1 2 3;4 5 6;7 8 9]; a= 1 2 3 不显示结果 4 5 6 7 8 9 当输入语句过长需要换行时，应加上“…”后再回车，则可续行输入。文件管理常用命令帮助命令 MATLAB有很多命令，因此很不容易记忆。使用HELP命 令可以得到有关命令的屏幕帮助信息。 如在MATLAB环境下直接运行HELP命令就会在屏幕上