MATLAB实验指导课件.docVIP

实验一 MATLAB 运算基础 1.实验目的 1）熟悉启动和退出 MA TLAB 的方法。 2）熟悉 MATLAB 命令窗口的组成。 3）掌握建立矩阵的方法。 4）掌握 MATLAB 各种表达式的书写规则以及常用函数的使用。 2.实验内容 1、用逻辑表达式求下列分段函数的值。 2 t , 0 t 1 2 t t 1, 1 t 2, 其中 t 0 : 0.5 : 2. 5. 2 2 t 3 t 2t 1, 2、求[100-999]之间能被 21 整除的数的个数。 3、建立一个字符串向量，删除其中的大写字母。 1 2 3 4、输入矩阵 A 4 5 6 ，并找出 A 中大于或等于 5 的元素（用行列表示） 。 7 8 9 5、求矩阵 a a 11 12 A 的行列式值、逆和特征根。 a a 21 22 6、不采用循环的形式求出和式 63 i S 2 的数值解。 i 0 3.实验程序及结果 1、程序如下： t=0:0.5;2.5... y=t.^2*((t=0)(t1))+(t.^2-1).*((t=1)(t2))+(t.^2-2*t+1).*((t=2)(t3)) 结果如下： 2、程序如下： p=rem([100:999],21)==0; sum(p) 运行结果如下： 3、程序代码如下： ch=ABcdefGHd,k=find(ch=Ach=Z),ch(k)=[] 4、程序代码如下： A=[1 2 3;4 5 6;7 8 9],[i,j]=find(A=5), for n=1:length(i) m(n)=A(i(n),j(n)) end M 运行结果： 5、程序代码： a11=input(a11=),a12=input(a12=), a21=input(a21=),a22=input(a22=), A=[a11,a12,a21,a22], DA=det(A),IA=inv(A),EA=eig(A) 运行结果如下： 6、程序代码如下： sum(2.^[0:63]) 结果如下： 四．心得体会 通过本次的学习，使我对 MATLAB 的使用有了基本的了解。熟悉了 MATLAB 的一系列运用和命令窗口的基本组成。掌握如何去建立矩阵及各种表 达式的书写规则和常规函数的使用。对行列式的运用有了进一本的了解。 实验二 M 函数与 M 文件的编写与应用 1.实验目的 1、熟悉 MATLAB 的环境与工作空间。 2、熟悉变量与矩阵的输入、矩阵的运算。 3、熟悉 M 文件与 M 函数的编写及应用。 4、熟悉 MATLAB 控制语句与逻辑运算。 2.实验内容 1、1 行 100 列的 Fibonacc 数组 a，a（1）=a(2)=1,a(i)=a(i-1)+a(i-2), 用 for 循环指令来寻求该 数组中第一个大于 10000 的元素，并指出其位置 i。 2、编写 M 脚本文件，定义下列分段函数，并分别求出当（ x1=1，x2=0.5 ）、（x1=-1 ，x2=0） 和（x1=0，x2=-0.5 ）时的函数值。 0.5457 e 0 .75 2 2 x 3 .75 2 x 1 1.5 x 1 , x 1 x 2 1 P( x , 1 x ) 0.7575 e 2 2 x 2 2 6x 1 , 1 x 1 x 2 1 x 1 x 2 1 0.5457 e 0 .75 2 2 x 3 .75 2 x 1 1.5 x 1 , 3、编写 M 函数表示曲线 t/ 3 y e sint 以及其它的包络线 t /3 y e ，并从命令窗口输入命令 1 语句绘制曲线。 t 的取值范围是 [0-4 ]。 1 1 4、设 f ( x) ,编写一个 M 文件，使的调用函数 f（x）时，x 2 x 4 (x 2) 0.1 ( 3) 0.01 可用矩阵代入，得出的 f（x）为同价矩阵。 \\ 5、实验程序及结果心得体会 3.实验程序及结果 1、程序代码如下： n=100;a=ones(1,n); for i=3:n a(i)=a(i-1)+a(i-2); if a(i)10000 a(i), break; end; end,i 程序运行结果如下： 2、程序代码如下： function[p]=ff(x1,x2) if x1+x21 p=0.5457*exp(-0.75*x2^2-3.75*1^2-1.5*1); elseif x1+x2=-1 p=0.5457*exp(-0.75*x2^2-3.75*1^2-1.5*1); else p=0.7575*exp(-x2^2-6.*x1^2); End 在命令窗口输入 x1=0.x2=-0.5 ff(x1,x2) 运行结果： ans=0.5899 输入 X1=-1,x2=0 Ff(x1,x2) 程序结果如下： Ans= 0.002