MATLAB实验指导教案资料Word格式文档下载.docx

1、3）结果如下：2、程序如下： p=rem（100:999,21）=0;sum（p）运行结果如下：3、程序代码如下：ch=ABcdefGHd,k=find（ch=A&ch=5）,for n=1:length（i）m（n）=A（i（n）,j（n）endM运行结果：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）4心得体会 通过本次的学习，使我对MATLAB的使用有了基本的了解。

2、熟悉了MATLAB的一系列运用和命令窗口的基本组成。掌握如何去建立矩阵及各种表达式的书写规则和常规函数的使用。对行列式的运用有了进一本的了解。实验二 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。2、编写M脚本文件，定义下列分段函数，并分别求出当（x1=1，x2=0.5）、（x1=-1，x

3、2=0）和（x1=0，x2=-0.5）时的函数值。3、编写M函数表示曲线sint以及其它的包络线，并从命令窗口输入命令语句绘制曲线。t的取值范围是0-4。4、设编写一个M文件，使的调用函数f（x）时，x可用矩阵代入，得出的f（x）为同价矩阵。5、实验程序及结果心得体会1、程序代码如下：n=100;a=ones（1,n）;for i=3:na（i）=a（i-1）+a（i-2）;if a（i）10000a（i）,break;end;end,i程序运行结果如下：2、程序代码如下：functionp=ff（x1,x2） if x1+x21p=0.5457*exp（-0.75*x22-3.75*12-1

4、.5*1）;elseif x1+x2=-1else p=0.7575*exp（-x22-6.*x12）;End在命令窗口输入x1=0.x2=-0.5ff（x1,x2）ans=0.5899输入 X1=-1,x2=0 Ff（x1,x2） 程序结果如下： Ans= 0.00229 输入 X1=1,x2=0.5程序结果如下Ans= 0.0224function y=ff（t ）y1=exp（-t/3）y2=exp（-t/3）.*sin（3*t）;y=y1;y2单击M文件编辑器的SAVE图标，使出现标准的文件保存对话框在文件对话框中选定目录填写文件名单击SAVE按钮把M文件保存到制定的文件目录中输入如下

5、命令 T=0:pi/100:4*pi; Plot（t,ff（tt）function f=fx（x） f=1./（x-2）.2+0.1）+1./（x-3）.4+0.01）4 实验心得 通过对M函数与M文件的编写与应用使我对MATLAB的环境与工作空进有了更近一步的熟悉，可以在文件保存对话框中选定目录，填写文件名并设置M文件保存在制定的目录中，熟悉了M文件与M函数的编写与应用和MATLAB的控制语句与逻辑运算。实验三 选择与循环结构程序的设计1、掌握if语句、switch语句、try语句的使用。2、掌握利用for语句、while语句实现循环结构的方法。1、根据求（1）y=3 my=y-f（n）%y

6、3时的最大y值 mn=n-1%y0,e=sum（p1）%正数的个数p2=fgrade=90 case 8B case 7C case 6D case num2cell（0:5）%数值向量必须转换为元胞数组才能作为判决条件E otherwise if per=100%grade=100 else grade=error%grade程序运行，在命令窗口输入成绩后可显示成绩等级。function y=ex3_4（x）for i=1:length（x） if （x（i）0）&（x（i）=-3）%如果x为一向量，x（x（i）5）&（x（i）=2）&（x（i）=3） y（i）=x（i）2-5*x（i）+6

7、else y（i）=x（i）2-x（i）-1y 程序运行如下：在命令窗口输入下列语句：X=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0,Y=ex3_4（x）有如下结果：y= 14.0000 11.0000 2.0000 1.0000 -0.2500 5.0000 19.0000 本实验是对循环与选择结构程序设计，在本实验中有对IF语句SWITCH语句TRY语句的使用。其中有当条件成立时，则执行语句组，执行完之后执行IF语句的后续语句，若条件不成立，则直接执行IF语句的后续语句。而SWITCH语句是更具表达式的取值不同，分别执行不同的语句。最后学习了循环结构有FOR WHILE循环

8、。，FOR循环变量可以为数据，WHILE循环以不定次数求一组命令的值，而BREAK语句与CONTINUE语句他们一般与IF语句配合使用。 实验四 绘图函数的应用1、熟悉二维和三维绘图函数的命令。2、熟悉图形修饰和控制方法。3、了解特殊坐标图形的绘制。1、在同一图形窗口绘制sinx，cosx曲线，要求用不同颜色和线型、数据点标记字符，X轴范围为0-2,y轴范围为-2-2，并加注解说明sinx、cosx以区分，图例注解放置在图形右下角。2、定义函数myfun：y（1,2）=200sinx/x，x2，绘制该函数在X=-20-20区间内的图形。3、试将图形窗口分割成三个区域，分别格绘制y=lgx在0-

9、100区间内的对数坐标，x半对数坐标及y半对数坐标，并加上标题添加图。4、t=-3:0.125:3，x=cos2t,y=cos2t，z=x2+2*y2，请绘制带有等高线的基于x,y,z的三维网格曲面图，并填充颜色。坐标轴范围为-1,1,-1,1,0,2.5、绘制出饱和非线性特性方程的曲线。1、输入的程序代码如下：0.05:2*pi;plot（t,sin（t）,r-.o,t,cos（t）,m-slegend（sinx,cosx,4）axis（0 2*pi -2 2）function Y=myfun（x）Y（:,1）=200*sin（x（:）./x（:）;,2）=x（:）.2;x=0:0.1:100;y=log10（x）;subplot（311）,loglog（x,y）gridtitle（loglog graphsubplot（312）,semilogx（x,y）semilogx graphsubplot（313）,semilogy（x,y）semilogy grapht=-3:3;x=sin（2.*t）;y=cos（2.*t）;x,y=meshgrid（x,y）;z=x.2+2*y.2;mesh（x,y,z）axis（-1