MTALAB试验一 Matlab使用方法和程序设计VIP

本文格式为Word版，下载可任意编辑——MTALAB试验一Matlab使用方法和程序设计试验一Matlab使用方法和程序设计

一、试验目的

1、把握Matlab软件使用的基本方法；

2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法

二、试验内容：

1、帮助命令

使用help命令，查找sqrt（开方）函数的使用方法；>>helpsqrt

SQRTSquareroot.

SQRT(X)isthesquarerootoftheelementsofX.ComplexresultsareproducedifXisnotpositive.

Seealsosqrtm,realsqrt,hypot.

Overloadedmethods:darray/sqrtsym/sqrt

ReferencepageinHelpbrowserdocsqrt2、矩阵运算（1）矩阵的乘法

已知A=[12;34];B=[55;78];求A^2*B

>>A=[12;34];

B=[55;78];S=A^2*BS=

105115229251（2）矩阵除法

已知A=[123;456;789];B=[100;020;003];求A\\B,A/B

>>A=[123;456;789];B=[100;020;003];s1=A\\Bs2=A/Bs1=

1.0e+016*

-0.45041.8014-1.35110.9007-3.60292.7022-0.45041.8014-1.3511s2=

1.00001.00001.00004.00002.50002.00007.00004.00003.0000

（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.',A'

>>A=[5+i,2-i,1;6*i,4,9-i];s1=A's2=A.'s1=

5.0000-1.0000i0-6.0000i2.0000+1.0000i4.00001.00009.0000+1.0000is2=

5.0000+1.0000i0+6.0000i2.0000-1.0000i4.00001.00009.0000-1.0000i（4）使用冒号表达式选出指定元素已知：A=[123;456;789];

求A中第3列前2个元素；A中所有列第2，3行的元素；方括号[]

用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列>>A=[123;456;789];B1=A([1,2],[3])B2=A([2,3],:)C=magic(4);C(:,4)=[]

B1=36B2=

456789C=

162351110976414153、多项式

（1）求多项式p(x)?x3?2x?4的根>>P=[10-2-4];roots(P)ans=

2.0000-1.0000+1.0000i-1.0000-1.0000i

（2）已知A=[1.2350.9;51.756;3901;1234]，求矩阵A的特征多项式;

把矩阵A作为未知数代入到多项式中；

>>A=[1.2350.9;51.756;3901;1234];p=poly(A)c=polyvalm(p,A)p=

1.0000-6.9000-77.2600-86.1300604.5500c=

1.0e-011*

0.1364-0.33500.00800.1421-0.0455-0.0114-0.1627-0.0909-0.30840.1705-0.0227-0.22170.0171-0.1265-0.03270.02274、基本绘图命令

（1）绘制余弦曲线y=cos(t)，t∈[0，2π]>>t=[0:0.05:2*pi];y=cos(t);plot(t,y)

0.20-0.2-0.4-0.6-0.8-101234567

（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]

>>t=[0:0.05:2*pi];y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1),holdonplot(t,y2)

0.20-0.2-0.4-0.6-0.8-101234567

5、基本绘图控制

绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；>>t=[0:0.1:4*pi];x1=10*sin(t);plot(t,x1,'r-.x')gridonxlabel('t');ylabel('x1');

title('x1=10sin(t)函数图象')

x1=10sin(t)函数图象108642x10-2-4-6-8-100246t8101214

6、基本程序设计

（1）编写命令文件：计算1+2+…+n>s=0;m=0;