还剩18页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
实验1熟悉软件环境和基本的操作
一、实验目的熟悉MATLAB运行环境和了解基本操作
二、实验内容MATLAB的启动、操作界面组成
1.熟悉MATLAB图形界面打开MATLAB,单击命令窗口菜单栏中的各个下拉菜单按钮,试使用各个按钮引出的选项;把光标移动到工具栏中各个图标上不要按下,查看它们与菜单选项的对应情况
2.熟悉MATLAB的基本命令在命令窗口中分别键入以下内容,以建立若干变量A=[l2;34;56]B=[7,8,9;10,ll,12]C=[567;l83];D=B+C问题1如何输入一个矩阵变量的行元素和列元素?问题2:观察每行命令后是否加,对显示执行结果有什么区别?键入以下命令或执行操作,查看效果,并体会命令功能1工作空间管理whowhosclear A2路径编辑试用菜单File/Set Path将D盘根目录及其下的所有子目录和文件夹包含进来,设为搜索路径问题3:当前路径是什么?问题4:搜索路径是什么意思?3联机帮助
三、思考题2=
1.0013821p=cos
3.52sin3p二---------------------e
1.取合适的0范围,在同一图形窗口绘制下4幅极坐标图4p=1-37^
2.某校共有1560名学生,其中计算机系有213名学生,外语系有387COS名学生,音乐系有220名学生,美术系有280名学生,中文系有280名学生,理科系180名学生,分别画出饼图、条形图示意学生分布
一、实验目的熟悉MATLAB的三维图形绘制函数
二、实验内容
1.标准球面绘制程序1subplot2,2,1,sphere3;titlen=3,axis equalsubplot2,2,2,sphere6;title,n=6,,axis equalsubplot2,2,3,sphere10titlen=10,axis equalsubplot2,2,4,sphere15;titlen=15,axis equal
2.标准球面绘制程序2ot=linspacepi/2,
3.5*pi,50;R=cost+2;subplot2,2,l;cylinderR,3,title*n=3!;subplot2,2,2cylinderR,6,title,n=6f;subplot2,2,3cylinderR,title,n=20,subplot2,2,4cylinderR,50,titlen=
503.绘制三维网格图绘制下面给出的二元函数的网格图sin Jx2+y2z二——.卜+y2[x,y]=meshgrid-8:
0.5:8,-10:
0.5:10;R=sqrtx.A2+y.A2+eps;z=sinR./R;meshx,y,z;
4.绘制三维网格图x=-pi:
0.15:pi;y=sinx;subplot2,1,1,H=bar x,y;xx=get H,xdata,;yy=get H,ydata,;subplot2,1,2,plot xx,yy;
三、思考题
1.试绘制出的三维曲面图和三视图
2.绘出函数的慧星效果图
一、实验目的熟悉和掌握MATLAB的m文件程序和程序设计
二、实验内容
1.for循环结构for i=l:10;xi=i;end;x
2.while循环结构Fibonaci数组的元素满足Fibonacci规则且现要求该数组中第一个大于10000的元素al=l;a2=l;i=2;owhile ai=10000ai+l=ai-l+ai;i=i+l;end;i,ai
3.if-else-end结构cost=10;number=12;if number8sums二number*.95*cost;end,sums
4.M文件的一般结构A编写一个画任意半径彩色线型的圆.B完整函数文件的基本结构.C函数文件各基本组成部分的作用.function sa=functiondemor,s%CIRCLEif nargin2error输入变量太多’;endif nargin==ls二b,;end;elf;t=0:pi/100:2*pi;x=r*expi*t;if nargin==2plot x,s;elsesa=pi*r*r;fill realx,imagx,sendaxissquare5
三、思考题
1.编写函数文件,分别用for和while循环语句编写程序,求出之和,并比较算法的运算时间
2.用MATLAB语言实现下面分段函数X由键盘输入
3.用for-end循环求出从100-200之间的所有素数
4.用while-end循环求1-100的所有数值和
一、实验目的熟悉和掌握MATLAB程序编写和调试的常用函数
二、实验内容
1.input:请求用户输入信息reply=input请输入一个字符;s;
2.menu:为用户输入选择菜单t=0:
0.1:2*pi;yl=sint;y2=cos t;y3=yl.*y2;k=menu菜单,菜单1,菜单2,菜单3;if k=lplot t,yl,r elseifk==2plot t,y2,g elseplott,y3,b end
3.pause:程序暂停执行pause5%程序暂停执5秒
4.global:定义全局变量假设工作空间中已有全局变量a,执行完下述函数后,观看工作空间中的变量有何变化.function z=example9_lglobal a;a=5;b=10;z=a+b;
三、思考题
1.输入一行字符,分别统计出其中英文字母、空格、数字和其它字符的个数
2.求之值,其中a是一个数字例如2+22+222+2222+22222(此时n=5),a和n由键盘输入
3.编写程序,判断某一年是否为闰年,闰年的条件是
(1)能被4整除,但不能被100整除的年份都是闰年,如1996年和2004年;
(2)能被100整除,又能被400整除的年份是闰年,如1600年和2000年不符合这两个条件的年份不是闰年(提示rem命令可以计算两数相除后的余数)实验9MATLAB应用
一、实验目的了解MATLAB在曲线拟合、插值和解线性方程方面的应用
二、实验内容
1.曲线拟合已知离散点上的数据集[x1,y1,x2,y2,•••,xn,yn],求得一解析函数y=fx,使fx在原离散点xi尽可能接近给定yi值,这一过程叫曲线拟合,最常用的曲线拟合是最小二乘拟合,拟合结果可使误差的平方和最小,即找出使最小的fxx=[
0.5,
1.0,
1.5,
2.0,
2.5,3,0];y=[
1.75,
2.45,
3.81,
4.80,
7.00860];p=polyfitx,y,2;xl=
0.5:
0.05:
3.0;yl=polyvalp,xl;plotx,y;*r,xl,yl,-b;
2.一维插值year=1900:10:2010product=[
75.995,
91.972,
105.711,
123.203,
131.669,...
150.697/
79.232,
203.212,
226.505,
249.633,
256.344,
267.893];p2005=interp1year,product,2005x=1900:l:2010y=interpl year,products,cubic*;plotyear,product;,x,y
3.解常微分方程二阶微分方程x+x2-lx+x=0可改写成两个耦合的一阶微分方程,玉=~°一一马、文-Xl2模拟这个系统的第一步是建立微分方程的M文件,随意取名作difeq.m,其内容是function xdot=difeqt,xxdot=zeros2,l;xdotl=xl.*l-x
2.A2-x2;xdot2=xl;为了在的范围内模拟该方程,可引用ode45函数求解并画出其图形tO=O;tf=20;xO=[O
0.25];[t,x]=ode23difeq,t0,tf,x0;plott,x
三、思考题f41,求定积分Q=[fxdx JOInx2当x V2式中fM=M16当x22+sinx+l
712.求解描述振荡器的经典的Ver derPol微分方程,,3,已知x=[
1.
21.
82.
12.
42.
63.
03.3],y=[
4.
855.
25.
66.
26.
57.
07.5],求对x和y进行6阶多项式拟合的系数附录MATLAB实验报告格式报告人年级学号实验日期报告完成日期
一、实验名称XXXXXX
二、实验目的XXXXXX
三、实验内容(结合实验指导书和课本,列出自己在本次实验中练习或验证了什么)
四、回答问题(回答实验指导书中提出的问题)
五、思考题(回答或编程求解每个实验后列出的思考题)
六、遇到的问题及解决(列出实验和求解思考题的过程中遇到的问题及解决方法)
七、体会(小结出从本次实验中学会了什么,难点是什么,有什么想法和体会)help pausehelpwin4窗口清理先画出正弦函数在0-2冗之间的图形,再用以下各种窗口清理命令,看每项命令都清除了什么figureplotsin0:
0.1:
6.28claelfclose注意figure为打开一幅图形图像窗口close为关闭当前图形图像窗口,而close all为关闭所有已打开的图形图像窗口5MATLAB基本矩阵操作演示playshow intro6MATLAB图形绘制演示playshow buckydem7MATLAB数学功能演示快速傅氏变换playshow fftdemo8MATLAB三维造型演示茶壶演playshow teapotdemo
3.打开MATLAB命令窗口,键入demos,观看演示程序
三、思考题
1.将pi分别用15位数字格式、分数格式、十六进制格式、5位数字的科学计数法显示
一、实验目的熟悉和掌握MATLAB中关于矩阵变换以及矩阵运算的各种命令
二、实验内容
1.数、数组、矩阵的输入1数的输入a=5b=2-5i2数组的输入c=[l,3,5,7,9,ll]%元素之间要用逗号用空格分开d=l:2:lle=linspacel,ll,6体会以上输入有什么区别和联系3矩阵的输入A=[2,3,5;l,3,5;6,9,4]%行之间要用分别隔开
2.矩阵大小的测试和定位A=[3,5,6;,2,5,8;3,5,9;3,7,9][n,m]=sizeAAl,
33.矩阵的块操作A2,:A[l,3],A2:3,l:2问题
2.1如何将A的2,3列互换?
4.矩阵的四则运算A=[3,5,8;-2,3,6;l,4,9]B=rand3,3C=A+B D=A-B E=A*B问题
2.2E为矩阵A.B的乘积运算结果,如果要求E的结果为A和B对应元素相乘的结果,应输入什么命令?F=A/B命令问题
2.3如果要求F为A,B对应元素作除法运算的结果,应输入什么
5.矩阵的点运算A=[12;34];B=[56;78];A*BA.*BAA2A.A
26.矩阵的逻辑运算A=[l2;34];B=[06;08];A|BABxorA,Ba=-5;b=-10;b〜=0a/b5b==0||a/b0
三、思考题
1.输入一个矩阵A,取出A的第2行第1列的元素;取出A的第1,3,4列的所有元素;让A的第1列和第3列互换;删除A的第二列1233l+4i436723574+2iA=13576+7i32391894_
2.用MATLAB可以识别的格式输入下面两个矩阵再求出它们的乘积矩阵C,并将C矩阵的右下角2X3子矩阵赋给D矩阵,赋值完成之后,调用相应的命令,查看MATLAB工作空间的情况
3.已知,B=,求i A*B-2A iiA*B iiiA.*B ivA*B-B*A
4.已知x=
[123],y=
[456],试计算z=x.*y、*./丫和*.\丫
5.解线性方程
6.解方程组:
一、实验目的熟悉和掌握MATLAB中各种矩阵生成函数
二、实验内容
1.随机矩阵的产生A=rand5,5%产生5X5均匀分布随机矩阵B=randn5,5%产生5X5正态分布随机矩阵
2.单位矩阵的产生A=eye3,3%产生3X3单位阵
3.全零、全1矩阵的产生A=zeros3,3Boones3,
34.矩阵的翻转操作A=rand5,4flipudA%进行上下翻转fliplrA%进行左右翻转rot90A%逆时针旋转90度
5.矩阵线性代数方面的运算A=[1,3,5;4,9,3;2,0,5]A%A转置det A%A的行列式,A必须是方阵rank A%A的秩inv A%A的逆问题
3.1矩阵A在什么时候不能求逆?eigA%A的特征矢量[X,D]=eigA%A的特征矢量及特征值A~2/矩阵A的幕运算
6.构造多维数组Al=[l,2,3;4,5,6;7,
8.9];A2=Al A3=A1-A2cat1,Al,A2,A3cat2,Al,A2,A3cat3,Al,A2,,A
37.字符串变量及其处理s=Hello s=[s;,Hello,]whos sn=input JHow manyapples%先显示单引号中的字符串,向用户提示,再把%用户键入的数字或表达式赋给n
三、思考题
1.产生3义4维全1矩阵,产生4X2维的随机矩阵,产生4X4的单位矩阵.
2.将A的第2行元素扩大2倍后作为A的第3行元素
3.输入任意矩阵A.B它们的元素个数相等,命令A:和人:尸B会产生什么结果?
4.A=[l,3,5;5,8,3;6,1,6],B=[3,6;9,3;4,7],C=[3,7,9;4;5,7],D=2:6,体会命令[A,B],[A;C],[A,B;D]所产生的结果,学习由小矩阵生成大矩阵的方法
5.计算表达式+233log5-tan21的值
26.已知矩阵A=[52;91],B=[12;92],做简单的关系运算AB,A==B,AB,并做逻辑运算A==BAB,A==BAB实验4图形绘制
(1)
一、实验目的熟悉和掌握MATLAB基本的二维图形绘制函数
二、实验内容
1.绘制简单的二维图形t=0:
0.l:2*piy=sint plott,y
2.一个坐标系绘制多幅图形t=0:
0.1:2*pi;yl=sint;y2=cost;y3=yl.*y2;plott,y1;-r;t,y2;:g;t,y3;x*plott,y l;-c\t,y2;-.y;t,y3;d,
3.图形标识和坐标控制t=0:
0.l:2*piy=sintplott,ygrid on,xlabelC时间’,ylabel幅值,title正弦曲线’axis[-l,8,-
1.2,
1.2]
4.交互式图形指令axis[0,5,0,5];hold on;box on;x=[];y=[];wh i1e1[xl,yl,button]=ginput1;ifbutton~=l break;endplot xl,yl,o;x=[x,xl];y=[y,yl];endlinex,y;hold off;gtext用左键取点,然后划线;
三、思考题
1.在同一坐标系绘制t3「t2,t2sint在[0,22内的曲线图
2.在同一图形窗口画三个子图,要求使用指定gtext axis、legend、title xlabel和ylabel:1y=XCOSX,X G一肛71132y=xtan—sinx,x,4xG3y=ex smx.x e[1,8]
3.绘制的图像,要求用蓝色的星号画图;并且在同一坐标轴中画出其包络线的图像,用红色的点划线画图实验5图形绘制2
一、实验目的熟悉和掌握MATLAB的多种二维图形绘制函数
二、实验内容
1.二维图形绘制函数X—2:
0.1:2;y=sinx;subplot221,stairsx,y,titlea stairssubplot222,compasscosx,y,title,b compass1yl=randn1,10000;subplot223,histy,20,title,c histsubplot224,[u,v]=meshgrid-2:
0.2:21:
0.15:1;z=u.*exp-u.A2-v.A2;[px,py]=gradientz,
0.2,
0.15;contouru,v,z,hold onquiveru,v,px,py,hold off,axis imagetitled quiver*
2.误差限图绘制函数x=-2:
0.2:2;y=sinx;L=rand1Jengthx/10;U=randl Jengthx/10;errorbarx,y,L U7:
93.复数图绘制函数z=[2+3i,2+2i,l-2i,4i,-3];x=[2,2,1,0,-3];y=[3,2,-2,4,0];subplot1,2,1,compass z,r subplot1,2,2,feather x,y,b
4.条形图与直方图绘制函数x=-pi:
0.15:pi;y=sinx;subplot2,1,1,H=bar x,y;xx=get H,xdata,;yy=get H,ydata;subplot2,1,2,plot xx,yy;。