MATLB实验报告1.doc

实验一 MATLAB运算基础 一、实验目的 1、熟悉启动和退出MATLAB的方法。 2、熟悉MATLAB命令窗口的组成。 3、掌握建立矩阵的方法。 4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。 二、实验内容 1、用逻辑表达式求下列分段函数的值。 ，其中t=0:0.5:2.5。 编程：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)) 程序运行结果如下：t =0 0.5000 1.0000 1.5000 2.0000 2.5000 y =0 0.2500 0 1.2500 1.0000 2.2500 2、求[100,999]之间能被21整除的数的个数。 编程：p=rem([100:999],21)==0; sum(p) 程序运行结果如下：ans =43 3、建立一个字符串向量，删除其中的大写字母。 编程：ch=ABcdefGHd,k=find(ch=Ach=Z),ch(k)=[] 程序运行结果如下：ch =cdefd 4、输入矩阵A=，并找出A中大于或等于5的元素（用行列表示）。 编程：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 程序运行结果如下：m =7 5 8 6 9 5、求矩阵A=的行列式值、逆和特征根。 编程：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) 若在命令窗口输入：A = 1 2 3 4 则程序运行结果如下：DA = -2 IA =-2.0000 1.0000 1.5000 -0.5000 EA =-0.3723 5.3723 若在命令窗口输入：A = 11 12 13 14 则运行结果如下：DA =-2 IA = -7.0000 6.0000 6.5000 -5.5000 EA =-0.0797 25.0797 6、不采用循环的形式求出和式S=的数值解。 编程：sum(2.^[0:63]) 程序运行结果如下：ans =1.8447e+019 实验二 M函数与M文件的编写与应用 一、实验目的 1、熟悉MATLAB环境与工作空间。 2、熟悉变量与矩阵的输入、矩阵的运算。 3、熟悉M文件与M函数的编写及应用。 4、熟悉MATLAB控制语句与逻辑运算。 二、实验内容 1、1行100列的Fibonacc数组a,a(1)=a(2)=1,a(i)=a(i-1)+a(i-2),用for循环指令来寻找该数组中第一个大于10000的元素，并指出其位置i。 编程：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 程序运行结果如下： ans =10946 i =21 2、编写M脚本文件，定义下列分段函数，并分别求出当、和时的函数值。 编程：function[p]=ff(x1,x2) If x1+x21 p=0.5457*exp(-0.75*x2^2-3.75*x1^2-1.5*x1); elseif x1+x2=-1 p=0.5457*exp(-0.75*x2^2-3.75*x1^2+1.5*x1); 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.0029;若在命令窗口输入：x1=1,x2=0.5,ff(x1,x2)程序运行结果如下：ans=0.0024 3、编写M函数表达曲线以及它的包络线，并从命令窗口输入命令语句绘制曲线。t的取值范围是。