MATLAB 第一次实验课课堂作业
北科大matlab第一次作业
《数学实验》报告实验名称MATLAB基础知识______ 学院专业班级姓名学号2012年4月一、【实验目的】1.掌握MATLAB的具体操作与操作键；2.掌握MATLAB中的常用函数与变量、表达式的定义方法；3.掌握MATLAB常用的工作方式M文件的编程工作方式；4.掌握MATLAB语言中的程序结构。
二、【实验任务】完成课后第一章第4题，第二章第1、2、3题，第三章第1题。
三、【实验程序】n? n=20编写函数文件，计算时的表达式的值。
，并求出当1.!k1k?functionf=first(n)f=0;i=1;while i<=ns=prod(1:i);i=i+1;f=f+s;endf2.用不同的数据格式显示自然底数e的值，并分析各个数据格式之间有什么相同于不同之处。
e=exp(1)format short e,eformat long,eformat long e,eformat short g,eformat rat,eformat hex,eformat bank,e134826?????????6AB?6545 3. A*B,A.*B，，计算，并比较两者的区别。
矩阵5????????827392????>> A=[1 2 3;4 5 6;7 8 9]A =21.002.003.004.005.006.007.00 8.00 9.00>> B =4.00 6.00 8.005.00 5.006.003.00 2.00 2.00>> A*B>> A.*B5212????，做简单的关系运算并做逻辑运算4.已知矩阵，?A?B,B,A?,A?BA?B????9921????).?ABB)&((B??)&(A?B),A??(A>> A=[5 2;9 1];B=[1 2;9 2];>>A>B>>A==B>>A<B>> (A==B)&(A<B)>> (A==B)&(A>B)?111?6?为止。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/f6c08c5902d276a201292ece.png)
1 + e2 (2) z = 1 ln( x + 1 + x 2 ) ，其中 x = ⎡⎢ 2⎣-0.45 ⎦2 2 ⎪t 2 - 2t + 1 2 ≤ t <3 ⎨实验一MATLAB 运算基础1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1) z = 2sin 8501221 + 2i ⎤5 ⎥(3) z = e 0.3a - e -0.3asin(a + 0.3) + ln 0.3 + a ,a = -3.0, - 2.9, L , 2.9, 3.03⎧t 2 0 ≤ t < 1 (4) z = ⎪t 2 - 11 ≤ t <2 ，其中 t=0:0.5:2.5 4⎩解：M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4.完成下列操作：(1)求[100,999]之间能被21整除的数的个数。
(2)建立一个字符串向量，删除其中的大写字母。
解：(1)结果：m=100:999;n=find(mod(m,21)==0);length(n)ans=43(2).建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch=⎣O2⨯3⎥，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩S⎦阵和对角阵，试通过数值计算验证A=⎢⎥。
matlab上机实验答案-整理版
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y)subplot(2,1,1)plot(x,y)subplot(2,1,2)plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi;rho=a*cos(b+n*theta);polar(theta,rho)3. 列出求下列空间曲面交线的程序clcclearx=[-5:0.5:5];[X,Y]=meshgrid(x);z1=X.^2-2*Y.^2;z2=X.*2-Y.*3;xlabel('x')ylabel('y')zlabel('z')surf(X,Y,z1)hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5);x1=X(k)y1=Y(k)z3=x1.^2-2*y1.^2hold onplot3(x1,y1,z3,'*')⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2212y x z -=y x z 322-=4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clcclearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2));plot(x,y,'b*-');title('绘图');xlabel('x 坐标');ylabel('y 坐标');legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x clccleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8];c=b/a;x=c(1,1)y=c(1,2)z=c(1,3)w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
实验一MATLAB基本操作并且附有答案.
实验一、MA TLAB基本操作一、实验目的熟悉MA TLAB软件环境,掌握命令窗口的使用。
二、实验内容及步骤1、命令窗口的简单使用(1简单矩阵的输入1 2 3A = 4 5 67 8 9>> A=[1 2 3 ;4 5 6 ;7 8 9 ](2求[12+2×(7-4]÷32的算术运算结果>>[12+2*(7-4]/3^22、有关向量、矩阵或数组的一些运算(1设A=15;B=20;求C=A+B与c=a+b?(2设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与A.*B? (3设a=10,b=20;求i=a/b与j=a\b?(4设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全下标的形式,并将其单下标转换成全下标。
>> a=[1 -2 3;4 5 -4;5 -6 7];>> [i,j]=find(a<0>>b=find(a<0 >> a(b >> [i,j]=ind2sub(size(a,b (5在MA TLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?>> >>A=[1,2;3,4]+i*[5,6;7,8]A =1.0000 + 5.0000i2.0000 + 6.0000i3.0000 + 7.0000i4.0000 + 8.0000i当输入A=[1,2;3,4]+i[5,6;7,8]显示出错(6请写出完成下列计算的指令:a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?>>a=[1 2 3;3 4 2;5 2 3];a^2 代表两个矩阵相乘ans =22 16 1625 26 2326 24 28>> a.^2 代表A矩阵元素的平方ans =1 4 99 16 425 4 9(7有一段指令如下,请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : 访问向量是列优先。
同济matlab选修课第一次实验内容参考答案
z1=poly2str(z,'x'); w=polyval(z,x); plot(x,y,'o',x,w,'r')
实验体会与总结
-3-Βιβλιοθήκη 的值，8 3 2 8, b2 3, r2 15
1 10, b1 , r1 28
在三维空间上画出解轨线。 解：
equation.m equation_run.m
function dx=equation(t,x,a,b,c); dx=zeros(3,1); dx=[-a*x(1)+a*x(2);c*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-b*x(3)];
实 验 报 告
班级 实 验 名 称 学号 姓名 Matlab 第一次试验 1．熟悉 Matlab 的操作环境和基本操作； 实 验 目 的 2. 掌握随机矩阵的构建； 3．掌握微分方程数值求解； 4. 掌握曲线拟合 参加实验时间 实 验 地 点 年 月 日 7、8 节课 电信楼 313
实验内容及其结果（写出源程序及其输出结果）
>> x0=[0,1e-10,1]; a1=10;b1=8/3;c1=28; [t,x]=ode45(@equation,[0,100],x0,[],a1,b1,c1); a2=8;b2=3;c2=15; [t,y]=ode45(@equation,[0,100],x0,[],a2,b2,c2); subplot(1,2,1);plot3(x(:,1),x(:,2),x(:,3)); subplot(1,2,2);plot3(y(:,1),y(:,2),y(:,3))
%
[i,j]= find(abs(A)>3);
MATLAB)课后实验答案
实验一 MATLAB 运算基础1、 先求下列表达式得值，然后显示MATLAB 工作空间得使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0、5:2、5 解：4、 完成下列操作：(1) 求[100,999]之间能被21整除得数得个数。
(2) 建立一个字符串向量，删除其中得大写字母。
解：(1) 结果：(2)、 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果就是：实验二 MATLAB 矩阵分析与处理1、 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵与对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5、 下面就是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。
(2) 将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。
(3) 计算系数矩阵A 得条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1、 求分段函数得值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y 值。
实验一-熟悉MATLAB-环境---答案
实验一熟悉MATLAB环境本人自己做的，并不是标准答案，仅供参考。
［实验目的］1.熟悉MATLAB主界面，并学会简单的菜单操作。
2.学会简单的矩阵输入与信号输入。
3.掌握部分绘图函数。
［实验原理］MATLAB是以复杂矩阵作为基本编程单元的一种程序设计语言。
它提供了各种矩阵的运算与操作，并有较强的绘图功能。
用户第一次使用MATLAB时，建议首先在屏幕上键入demo命令，它将启动MATLAB的演试程序，用户可在此演示程序中领略MATLAB所提供的强大的运算与绘图功能。
也可以键入help进行进一步了解。
MATLAB启动界面如图所示：操作界面主要的介绍如下：口指令窗(CommandWindow),在该窗可键入各种送给MATLAB运作的指令、函数、表达式，并显示除图形外的所以运算结果。
口历史指令窗(CommandHistory)，该窗记录已经运行过的指令、函数、表达式；允许用户对它们进行选择复制、重运行，以及产生M文件。
口工作空间浏览器(WorkspaceBrowser)，该窗口罗列出MATLAB工作空间中所有的变量名、大小、字节数；并且在该窗中，可对变量进行观察、编辑、提取和保存其它还有当前目录浏览器（CurrentDirectoryBrowser）、M文件编辑/调试器（Editor/Debugger）以及帮助导航/浏览器（HelpNavigator/Browser）等，但通常不随操作界面的出现而启动。
利用File菜单可方便对文件或窗口进行管理。
其中File|New的各子菜单，M-file（M文件）、Figure（图形窗口）、或Model（SIMULINK编辑界面）分别可创建对应文件或模块。
Edit菜单允许用户和Windows的剪切板交互信息。
MATLAB语言最基本的赋值语句结构为：变量名列表=表达式。
表达式由操作符或其它字符，函数和变量名组成，表达式的结果为一个矩阵，显示在屏幕上，同时输送到一个变量中并存放于工作空间中以备调用。
大学数学matlab实验作业1
重庆大学学生实验报告实验课程名称数学实验开课实验室学院年级专业班学生姓名学号开课时间至学年第学期数理学院制开课学院、实验室：实验时间：年月日2）比较函数33)(,)(,)(xxhxxxgxxf=+==在x→0时函数的性态。
观察到什么现象？从观察到的现象，反映了什么结论。
程序：x=-1:0.0001:1;y1=x;y2=x.^3;y3=y1+y2;plot(x,y1,x,y2,x,y3)结果：结论：当x→0时，f(x)与g(x)很接近，而h(x)与前两个函数都不接近。
3）比较函数33)(,)(,)(xxhxxxgxxf=+==在x→∞时函数的性态。
程序如下所示：x=linspace(-100000,100000,30); y1=x;y2=x+x.^3;y3=x.^3;subplot(2,2,1),plot(x,y1),title('f(x)=x'),xlabel('x');ylabel('f(x)'); grid; subplot(2,2,2),plot(x,y2),title('g(x)=x+x^3'),xlabel('x');ylabel('g(x)'); grid;subplot(2,2,3),plot(x,y3),title('h(x)=x^3'),xlabel('x');ylabel('h(x)');grid ；结果：4）在日常生活中我们有这样的经验：与幂函数相比，指数函数是急脾气，对数函数是慢性子。
这就是说，当x→∞时，再小的指数函数也比幂函数变化快，再大的对数函数也比幂函数变化慢。
当x→∞时,比较10xy=与xy1.1=的大小.当x→∞时,比较001.0xy=与xy lg1000=的大小.程序如下所示：x=linspace(5000,8000,500);y1=x.^10;y2=1.1.^x;Subplot(1,2,1),plot(x,y1),xlabel('x');ylabel('y)'); grid; title('y=x^1^0');Subplot(1,2,2),plot(x,y2),xlabel('x');ylabel('y)'); grid; title('y=1.1^x');结果：从上图可以看出来指数函数变化快程序如下所示：x=linspace(5000,8000,500);y1=x.^0.001;y2=1000.*log(x);Subplot(1,2,1),plot(x,y1),xlabel('x');ylabel('y)'); grid; title('y=x^0.001');Subplot(1,2,2),plot(x,y2),xlabel('x');ylabel('y)'); grid; title('y=1000.*log(x)');结果：分析：由以上函数图形可知对数函数变化比幂函数慢。
【原创】MATLAB实验报告-第一次-MATLAB语言程序设计基础-北京交通大学
输入what<目录名>，表示将分类显示所选择目录下的所有相关文件；
输入ver singal，表示列出singal工具箱的版本信息。
2、在命令窗口中输入 命令，执行后用方向 或 键（历史命令）将该命令调出，在行内编辑，将其改为6阶魔术矩阵命令。
输入magic（4），显示4阶魔术矩阵，按上方向键调出magic（4）命令，并将4改为6，便显示出6阶魔术矩阵。
3、用Matlab语言绘制函数 在 内的曲线。
输入>> t=0:0.1:2*pi;y=t.*cos(t);plot(t,y)，显示出图像。
4、选择适当的坐标范围，用Matlab语言绘制函数： 的三维曲线。
输入
>> [x,y]=meshgrid(-3:0.1:3);
>> z=(x.^2+y).*exp(-x.^2-y.^2/2);
输入addpath E:\qinmin，将E:\qinmin目录添加到Matlab包含的路径下：
>> addpath E:\qinmin
（3）显示格式
输入format long，表示输出双精度数时小数点后保留14位有效数字，输出单精度数是小数点后保留7位有效数字；
输入format short e，表示输出数字用5位科学技术表示；
>> surf(x,y,z),colorbar
显示三维图像：
三、实验感想
本次Matlab上机实验是我第一次亲自操作Matlab软件并实现一些简单的功能。虽然只是实验了一些基本的命令，显示了魔术矩阵和绘制了简单的函数图像，但是这次实验也使我意识到了Matlab软件的强大。Matlab是著名的数学工具，在各科研领域都具有不可忽视的影响力。作为通信工程的本科生，我们也不可避免地需要使用Matlab来处理一些实验数据，所以熟练使用Matlab也是我们必备的技能之一。
MATLAB实验作业一
MATLAB的基本操作(一)实验内容及步骤1.进入MATLAB工作环境，熟悉各窗口的功能。
（1）双击桌面图标，或从“开始”菜单打开MATLAB.exe,启动MATLAB。
（2）查看MATLAB界面各窗口的布局、了解其功能，并完成各窗口之间的切换。
（3）设置当前工作目录。
在D盘创建mymatlab目录，并将其作为当前工作目录。
今后的实验过程中以此目录作为当前工作目录。
2.计算5.542+cos(35π)√17.89÷3.5−e2的值。
（1）在命令窗口（Command Window）中输入程序：a=5.54^2+cos(3/5*pi)*sqrt(17.89)/3.5-exp(2)（2）按回车键运行，如果出现“Error”（出错信息），则应找出原因并改正，再运行。
（3）运行结果22.9291过程： >> a=5.54^2+cos(3/5*pi)*sqrt(17.89)/3.5-exp(2)a =22.9291>>3．输入矩阵A=[428256319], B=[111222333], 在命令窗口中执行下列表达式，掌握其含义：A(2,3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A（1）在命令窗口键入A=[4 2 8;2 5 6;3 1 9] ,生成矩阵A;键入B=[1 1 1;2 22;3 3 3] ，生成矩阵B。
>>A=[4 2 8;2 5 6;3 1 9] A =4 2 82 5 63 1 9>> B=[1 1 1;2 2 2;3 3 3]B =1 1 12 2 23 3 3>> A(2,3)ans =6>> A(6)ans =1>> A(:,2)ans =251>> A(3,:)ans =3 1 9>> A(:,1:2:3)4 82 63 9>> A(:,3).*B(:,2)ans =81227>> A(:,3)*B(2,:)ans =16 16 16 12 12 12 18 18 18>> A*Bans =32 32 32 30 30 30 32 32 32>> A.*Bans =4 2 8 4 10 12 9 3 27>> A^2ans =44 26 116 36 35 10041 20 111>> A.^2ans =16 4 644 25 369 1 81>> B/Aans =0.5000 0.0769 -0.38461.0000 0.1538 -0.76921.5000 0.2308 -1.1538>> B./Aans =0.2500 0.5000 0.12501.0000 0.4000 0.33331.0 3.0000 0.3333>> B.\Aans =4.0000 2.0000 8.00001.00002.50003.00001.0000 0.3333 3.0000>>4. 产生一个5阶魔方矩阵，将矩阵的第3行4列元素赋值给变量a；将由矩阵第2,3,4行第2,5列构成的子矩阵赋值给变量b。
大连东软信息学院Matlab实验报告-第一次作业
MATLAB教程练习题1．利用Matlab中的函数生成3×3单位阵，2×3全0阵和4×3全1阵；2．利用Matlab中的函数生成10个在区间[14,99]上均匀分布的随机整数。
3. 计算向量（2，4，6，8）的最大值，最小值，平均值，中值【查一下这个函数】；4. 在行向量（2，4，6，8）和（1，3，5，7）之间实施加减乘除及幂运算；5. 在向量（1，3，5）与标量2之间实施加减乘除及幂运算。
6. 求 121211130⎛⎫ ⎪ ⎪ ⎪⎝⎭－的转置、行列式【查一下这个函数】、逆矩阵【查一下这个函数】。
7.要求利用diag 函数和其他函数生成数组A =。
8.将上面生成的数组A = ，转换成行向量B。
9. 产生一均匀分布在（-5，5）之间的5*3阶随机矩阵C，并要求所有的数值为整数。
10. 利用函数产生一个（0,1）间均匀分布的6*6随机矩阵D，然后统计D中大于等于0.6的元素个数。
11.产生一个4*5的随机矩阵E，要求为整数，且矩阵的元素范围为[0,100]，计算出矩阵的最大值所在的位置。
12. 产生100到1的等差数列，数据点为20个，赋值给F（1）取F的第10个单元的值，赋值给F1；（2）取F的前10个单元的值，赋值给F2；（3）取F的后5个单元的值，赋值给F3；（4）取F的偶数单元的值，赋值给F4；（5）将F的数据倒排，赋值给F5。
13. 定义矩阵G=[1 2 3;4 5 6;7 8 9]；（1）抽取A中的第1、2行，第2、3列元素，赋值给G1；（2）抽取A中的第1、4个元素作为新矩阵G2的第1行，抽取第2、3个元素作为新矩阵G2的第2行；resleep(1 4 2 3 , 2,2)（3）对矩阵G1和G2进行求和、求差、求乘和点乘。
matlab第一次作业
《信号与系统》课程研究性学习手册姓名吕航(第十一组)学号09213044同组成员马琳刘煜晨朱志杰王琛珏指导教师魏杰时间2010/10研讨题目：1．M2-3题目分析：题目考查连续信号的编写，和信号的尺度变换、翻转、时移（平移）。
仿真程序：x.mfunction yt=x(t)yt=t.*(t>=0&t<2)+2.*(t>=2&t<=3)-1.*(t>=3&t<=5);untile1.mt=-1:0.001:6;plot(t,x(t))title('x(t)')axis([-1,6,-2,3])untile2.mt=-1:0.001:11;subplot(2,2,1)plot(t,x(0.5*t))title('x(0.5*t)')axis([-1,11,-2,3])t=-7:0.001:5subplot(2,2,2)plot(t,x(2-0.5*t))title('x(2-0.5*t)')axis([-8,6,-2,3])仿真结果：结果分析：符合题意，结果正确。
自主学习内容：课本P63 2.6.3阅读文献：【1】陈后金信号与系统【M】高等教育出版社2007. 【2】魏杰ppt 2010.发现问题(专题研讨或相关知识点学习中发现的问题)：在函数定义中一开始写的是：x.mfunction yt=x(t)yt=t.*(t>=0&t<=2)+2.*(t>=2&t<=3)-1.*(t>=3&t<=5); untile1.mt=-1:0.001:6;subplot(2,1,1)plot(t,x(t))title('x(t)')axis([-1,6,-2,3])得到的是如图：在t=2时不是连续的。
后来发现t.*(t>=0&t<=2)中t<=2与后面的表达式在t=2点重复定义。
MATLAB实验一(课后练习答案)
求： ， ， ， ， ， ，
答：
>> a=[1 2;3 4;5 6],b=[5 6;3 4;1 2],c=[1 2 3;4 5 6]
a =
1 2
3 4
5 6
b =
5 6
3 4
1 2
c =
1 2 3
4 5 6
>> [mA,nA]=size(a),[mB,nB]=size(b)
mA =
3
nA =
-3 6 -5 -9
2 -7 -12 -8
B =
5 4 3 -2
6 -2 3 -8
-1 3 -9 7
>> C1=A*B',C2=A'*B,C3=A.*B
C1 =
19 -82 30
12 27 3
-38 54 29
C2 =
-15 16 -24 36
63 -17 93 -105
22 6 117 -60
19 46 84 -10
2
mB =
3
nB =
2
>>a+b
ans =
6 8
6 8
6 8
>>a-b
ans =
-4 -4
0 0
4 4
>>a*c
ans =
9 12 15
19 26 33
29 40 51
>>a*5
ans =
5 10
15 20
25 30
>>a/2
ans =
0.5000 1.0000
1.5000 2.0000
2.5000 3.0000
Matlab实验第一次实验答案
实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；解：sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = 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.RemarksSee sqrtm for the matrix square root.Examplessqrt((-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）求多项式 42)(3--=x x x p 的根解：A=[1 0 -2 -4];B=roots(A)（2）已知A=[ 3 5 ;5 5 6;3 9 0 1;1 2 3 4] ，求矩阵A 的特征多项式;解：A=[ 3 5 .9 ; 5 5 6 ;3 9 0 1 ;1 2 3 4];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)grid5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；解法1：s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2：s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
MATLAB - 第一次试验结果
>> Ex01_02
3． 画出衰减振荡曲线 y = e
−
t 3
sin 3t 及其包络线 y 0 = e 3 ，y0 = e 3 。t 值的范围是 [0,4π ] 。
%Ex01_03 function y=Ex01_03
−
t
t
t=0:0.1:4*pi; y=sin(3*t).*exp(-t/3); y00=-exp(-t/3); y01=exp(-t/3); plot(t,y,'-',t,y00,'--',t,y01,'-.'); title('衰减震荡曲线及其包络线'); xlabel('t'); ylabel('y');
实 验 报 告
班级 机械 2 班 学号 081704 姓名 刘鑫 实 验 名 称 Matlab 第一次试验 1．掌握 Matlab 绘制曲线； 实 验 目 的 2．掌握 Matlab 绘制曲面； 3、掌握构造分段函数的方法； 参加实验时间 实 验 地 点 实验结果 − 4 − 2 0 2 4 1．找出数组 A = 中所有绝对值大于 3 的元素。 − 3 − 1 1 3 5
%Ex01_02 function Ex01_02 %画出所表示的三维曲面，用暖色调。 x=[-8:0.1:8]; y=[-8:0.1:8]; [xx,yy]=meshgrid(x,y);%设置栅格坐标矩阵 zz=sin(sqrt(xx.^2+yy.^2))./sqrt(xx.^2+yy.^2); surf(xx,yy,zz);%画三维曲面图 shading interp; colormap(hot);%暖色
-3-
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/d27ebaa958fb770bf78a55b8.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下： 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
MATLAB-第一次实验课课堂作业
M A T L A B-第一次实验课课堂作业(总1页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--MATLAB 第一次实验课课堂作业（4学时）学号 2 姓名陈丹注：1）此课堂作业作为本课程结业成绩的重要依据，请同学们认真、独立完成，不得抄袭。
2）请在授课教师规定的时间内完成；3）完成作业后，请以word格式保存，文件名为：学号+姓名1. 请列出 d:\ 文件夹下所有的文件或子文件夹信息。
dir d:\$ RECYCLER$360Rec System Volume Information36dc61c8857c08269dda9343f78d03 zemax2003b9952393a57d6eerecycled2. 请写出文件 fft..m的路径。
C:\MATLAB7\toolbox\matlab\datafun\3. 请使用帮助命令查询的使用方法。
完成如下问题：1）简要说明fft的中文使用方FFT(X):对于变量x的离散傅里叶变换的一种快速算法。
Fft(x)的计算长度为向量x的长度；fft(x,n)完成指定长度为n的DFT变换。
其中x序列被自动补零或截短到长度n。
FFT(X)是 X的快速傅立叶变换，FFT(X,N)是X的n点快速傅立叶变换，FFT(X,[],DIM)是X的多维快速傅立叶变换，FFT(X,N,DIM)是X的n维度的快速傅立叶变换。
2）计算x=[1 0 2 3]的离散傅里叶变换的值X。
X(0)= X(1)= + X(2)= X（3）= -3）把X保存为数据文件 X..matsave X2。
Matlab实验作业
《MATLAB 语言及其应用》实验指导书目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；解：sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = 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.RemarksSee sqrtm for the matrix square root.Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.41422、矩阵运算（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）求多项式 42)(3--=x x x p 的根解：A=[1 0 -2 -4];B=roots(A)（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ， 求矩阵A 的特征多项式;解：A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4]; 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(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π] t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1，t,y2)grid5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；解法1：s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2：s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB 第一次实验课课堂作业（4学时）
学号100401202 姓名陈丹
注：1）此课堂作业作为本课程结业成绩的重要依据，请同学们认真、独立完成，不得抄袭。
2）请在授课教师规定的时间内完成；
3）完成作业后，请以word格式保存，文件名为：学号+姓名
1. 请列出d:\ 文件夹下所有的文件或子文件夹信息。
dir d:\
$RECYCLE.BIN.exe RECYCLER recycled.exe
$Recycle.BIN RECYCLER.exe sysanti.exe
360Rec System V olume Information syslive.exe
36dc61c8857c08269dda9343f78d03 autorun.inf zemax2003
Desktop.ini b9952393a57d6ee
Folder.htt recycled
2. 请写出文件fft..m的路径。
C:\MA TLAB7\toolbox\matlab\datafun\fft.m
3. 请使用帮助命令查询fft.m 的使用方法。
完成如下问题：
1）简要说明fft的中文使用方
FFT(X):对于变量x的离散傅里叶变换的一种快速算法。
Fft(x)的计算长度为向量x的长度；fft(x,n)完成指定长度为n的DFT变换。
其中x序列被自动补零或截短到长度n。
FFT(X)是X的快速傅立叶变换，FFT(X,N)是X的n点快速傅立叶变换，
FFT(X,[],DIM)是X的多维快速傅立叶变换，FFT(X,N,DIM)是X的n维度的快速傅立叶变换。
2）计算x=[1 0 2 3]的离散傅里叶变换的值X。
X(0)= 6.0000 X(1)= -1.0000 + 3.0000i X(2)= X（3）=-1.0000 - 3.0000i
3）把X保存为数据文件X..mat
save X。