实验一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 +
0+
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）矩阵的转置及共轭转置
已知 A=[5+i,2-i,1;6*i,4,9-i];
求A.', A'
解： A=[5+1i,2-1i,1;6*1i,4,9-1i ];
A1=A.',A2=A'
（4）使用冒号表达式选出指定元素
已知： A=[1 2 3;4 5 6;7 8 9];
求A中第 3列前 2个元素； A中所有列第 2， 3行的元素；方括号 []
解： A=[1 2 3;4 5 6;7 8 9 ];
B=A([1,2],[3]),C=A(2:end, : )
用magic函数生成一个 4阶魔术矩阵，删除该矩阵的第四列3、多项式
（1）求多项式p( x) x32x 4的根
解： A=[1 0 -2 -4];
B=roots(A)
（2）已知A=[35;5 56;3901;1234]，
求矩阵 A的特征多项式 ;
解：
A=[35.9;5 56;3901;1234];
A=poly(A);
A=poly2sym(A)
把矩阵 A作为未知数代入到多项式中；
4、基本绘图命令
（1）绘制余弦曲线y=cos(t)，t∈[0，2π]
解： t=0:.1:2*pi;
y=cos(t);
plot(t,y),grid
（2）在同一坐标系中绘制余弦曲线y=cos和正弦曲线 y=sin ， t ∈[0 ，2π] t=0:.1:2*pi; y1=cos;
y2=sin;
plot(t,y1， t,y2)
grid
5、基本绘图控制
绘制 [0 ，4π] 区间上的 x1=10sint 曲线，并要求：
（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线
（3）标注控制：坐标轴名称、标题、相应文本；
t=0:.01:4*pi;
y=10*sin(t);
plot(t,y,'-.' ,t,y, 'r' )
grid
6、基本程序设计
（1）编写命令文件：计算 1+2+ +n<2000 时的最大 n值；解法 1：
s=0;i=0;
while (s<2000) i=i+1;s=s+i;end
s=s-i,i=i-1
解法 2： s=0;
for i=1:1000; s=s+i;if (s>2000) , break ;
end
end
s=s-i,i=i-1
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2的 0到 n次幂的和。
解法 1： function[ s i ] = mitifang1(n)
s=0;
for i=0:n
s=s+2^i; if (i>n),break ;
end
end
s,n
解法 2： function[ s i ] = mitifang2(n)
s=0;i=0;
while (i<=n) ,s=s+2^i;i=i+1;
end
s,n
（3）如果想对一个变量 x自动赋值。当从键盘输入y或 Y时（表示是）， x自动赋为 1；
当从键盘输入 n或N时（表示否）， x自动赋为 0；输入其他字符时终止程序。reply=input( '?? ê?è?×???: ' , 's' );
'N'
while reply== 'y' ||reply== 'Y' ||reply== 'n'
||reply==
if (reply=='y' ||reply== 'Y' )
x=1;
disp('x=' )
disp(x)
elseif (reply== 'n' ||reply== 'N' )
x=0;
disp('x=' )
disp(x)
else
disp( 'shu ru you wu' )
end
reply=input('?? ê?è?×???￡o ' , 's' );
end
disp( 'shu ru you wu ')
解：
三、实验报告要求：
编写实验内容中的相关程序在计算机中运行，程序、运行结果及相关图形一并写在报告上。
