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

实验一 Matlab使用方法和程序设计

一、 实验目的

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

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

二、 实验内容：

1、帮助命令

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

SQRT Square root.

SQRT(X) is the square root of the elements of X. Complex results are produced if X is not positive.

See also sqrtm, realsqrt, hypot.

Overloaded methods: darray/sqrt sym/sqrt

Reference page in Help browser doc sqrt 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]; S=A^2*B S =

105 115 229 251 （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]; s1=A\\B s2=A/B s1 =

1.0e+016 *

-0.4504 1.8014 -1.3511 0.9007 -3.6029 2.7022 -0.4504 1.8014 -1.3511 s2 =

1.0000 1.0000 1.0000 4.0000 2.5000 2.0000 7.0000 4.0000 3.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.0000i 0 - 6.0000i 2.0000 + 1.0000i 4.0000 1.0000 9.0000 + 1.0000i s2 =

5.0000 + 1.0000i 0 + 6.0000i 2.0000 - 1.0000i 4.0000 1.0000 9.0000 - 1.0000i （4） 使用冒号表达式选出指定元素 已知： A=[1 2 3;4 5 6;7 8 9];

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

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

B1 = 3 6 B2 =

4 5 6 7 8 9 C =

16 2 3 5 11 10 9 7 6 4 14 15 3、多项式

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

2.0000 -1.0000 + 1.0000i -1.0000 - 1.0000i

（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ， 求矩阵A的特征多项式;

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

>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]; p=poly(A) c=polyvalm(p,A) p =

1.0000 -6.9000 -77.2600 -86.1300 604.5500 c =

1.0e-011 *

0.1364 -0.3350 0.0080 0.1421 -0.0455 -0.0114 -0.1627 -0.0909 -0.3084 0.1705 -0.0227 -0.2217 0.0171 -0.1265 -0.0327 0.0227 4、基本绘图命令

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

10.80.60.40.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),hold on plot(t,y2)

10.80.60.40.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') grid on xlabel('t'); ylabel('x1');

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

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

6、基本程序设计

（1）编写命令文件：计算1+2+…+n<2000 时的最大n值； >> s=0;m=0;

while (s<2000),m=m+1; s=s+m; end

s=s-m;m=m-1; [s,m] ans =

1953 62

（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。 For循环结构的函数：

function [s]=shiyan162a(k) s=0;

for i=1:k, s=s+2^i; end

While循环结构的函数：

function [m,s]=shiyan162b(k) s=0;m=0; while (m

（3）如果想对一个变量x自动赋值。当从键盘输入y或Y时（表示是），x自动赋为1；

当从键盘输入n或N时（表示否），x自动赋为0；输入其他字符时终止程序。

>> str=input('输入字符：','s'); if str=='y'|str=='Y' x=1

elseif str=='n'|str=='N' x=0

end

本文来源：https://www.bwwdw.com/article/1zx8.html