MATLAB基础及其上机实习题1
matlab上机习详解_试题答案
P 第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=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 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-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.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y为纵坐标的曲线，要求有图形标注。
clc clearx=-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 xclc cleara=[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上机实验练习题及答案
MATLAB上机实验练习题及答案09级MATLAB上机实验练习题1、给出一个系数矩阵A[2 3 4;5 4 1;1 3 2]，U=[1 2 3]，求出线性方程组的一个精确解。
2、给出两组数据x=[0 0.3 0.8 1.1 1.6 2.3]’y=[0.82 0.72 0.63 0.60 0.55 0.50]’，我们可以简单的认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。
3、解线性方程4、通过测量得到一组数据：5、已知一组测量值6、从某一个过程中通过测量得到：分别采用多项式和指数函数进行曲线拟合。
7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近插值法，双线性插值法，二重三次方插值法）8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。
分别使用不同的颜色，线形和标识符。
9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额占总销售额的比例，并分离第三种产品的切片。
X= 19.3 22.1 51.634.2 70.3 82.4 61.4 82.9 90.8 50.5 54.9 59.1 29.4 36.3 47.010、对应时间矢量t ，测得一组矢量y采用一个带有线性参数的指数函数进行拟合，y=a 0+a 1e -t +a 2te -t ，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。
11、请创建如图所示的结构数组（9分）12、创建如图所示的元胞数组。
（9分）13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。
要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。
14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。
MATLAB上机习题一
MATLAB上机习题一请按以下步骤完成上机实验：1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内；2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面；3）如果习题是问答题，请将答案写在题目后；4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后；5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”；6）上传该文件到FTP的相关目录。
1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。
注：不需截图，但要说明如何操作的2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。
3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。
4. 通过使用帮助确定内置变量ispc 的含义。
只需用简单语言说明。
6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。
只需用简单语言说明。
5. 创建double类型的变量，并进行计算1）a=87，b=190，计算a+b 、a-b 、a*b ；2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致，并说明为什么。
6. 计算如下表达式：1）()sin 602）3e3）3cos 4π⎛⎫ ⎪⎝⎭4）27562323336+⨯-⨯⨯7. 设u=2，v=3，计算：1）4log uv v2）()22u e v v u+- 38. 计算如下表达式：1）()()3542i i -+2）()sin 28i -9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：1）420<2）420≤3）420==4）4~20=5）''''b B <10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b >2）a b <3）&&a b b c >>4）a d ==5）|a b c >6）~~d11. 判断下列变量名哪些是合法的MATLAB 变量：1）fred7）fred! 2）book_18）book-1 3）2ndplace9）Second_Place 4）#110）No_1 5）vel_511）vel.5 6）tan12）while 提示，可以使用isvarname 对上述变量名进行检验。
matlab上机实验题
“MATLAB及其在通信中的应用”上机实验一1姓名学号班级一、MATLAB环境与帮助1．简述command window的作用。
（5分）答案:输入命令，输出结果。
具体有：（1）输入语句（2）输入变量名（3）输入系统命令/函数（4）输入.m 文件名（5）输出上述的运算结果，包括数值或图形2．（1）简述command history的作用；（2）列举出四个在command window中快速操作command history内容的快捷键。
（5分）答案：（1）存储运行过的命令，便于以后调用或修改（2）快捷键：方向键上、下、左、右。
3．简述workspace的作用，并说明workspace的菜单上的五个按钮的功能。
（5分）答案：（1）列出程序创建的所有变量（2）加入新变量，打开所有变量，装载数据文件，保存空间或变量，删除变量4．简述清除command window、command history、workspace中内容的方法（5分）答案：清除command window内容：在command window键入clc清除command history内容：command history窗口上点右键，选择clear command history清除工作空间的内容：在command window键入clear5．（1）简述current directory的作用；（2）在桌面新建一个文件夹，命名为自己的学号，如（1110920101），将current directory1完成后以附件形式发送到邮箱commu_*******************。
邮件主题为“班级”加“下划线”加“姓名”加“下划线”加“实验一”，如“通信1101_姓名_实验一”，word文件名类似，如“通信1101_姓名_实验一”。
改为新建文件夹的路径，简述操作方法,并对更改后的current directory截图，填入到下面空白处。
Matlab基础及应用第一章上机作业
1、新建一个文件夹（自己的名字命名）
答：
2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。
保存，关闭对话框。
使用path命令查看MATLAB搜索路径。
3、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye。
4、使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
5、编写一段简短的脚本文件，保存并显示运行结果。
6、创建两个double型变量a=32、b=5，并计算a+b、a-b、a*b、a/b、a\b；
创建两个int8型变量a、b，取同样数值，并计算a+b、a-b、a*b、a/b、a\b，对于计算结果与前次计算结果不同的情形请给出解释。
7、查看int16数据类型的取值范围(intmin,intmax)；查看单精度数据类型的取值范围和精度(realmin,realmax,eps)。
8、求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量：
(1)sin(60o) (2) e3 (3) cos(3/4π)
9、用两种定义复数的方法计算(直接定义，complex(a,b))
(1)(3-5i)(4+2i) (2) sin(2-8i)。
Matlab上机实验答案
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 >> n=100:999; >> l=find(rem(n,21)==0); >> length(l)
ans =
43
>> ch='aegbBOIEG0je23RGnc';
150 335 520 705 890
77 237 397 557 717
>> D=C(3:5,2:3)
D=
520 705 890 >> whos Name
397 557 717
Size
Bytes Class
Attributes
A B C D
5x5 5x3 5x3 3x2
200 double 120 double 120 double 48 double
(2) >> A*B ans =
68 309 154 >> A.*B ans =
44 -72 -5
62 596 241
(3) >> A^3 ans =
12 102 68 0 9 -130
4 261 49
37226 247370 78688
233824 149188 454142
48604 600766 118820
>> A.^3 ans = 1728 39304 27 39304 343 274625 -64 658503 343
(4) >> A/B ans =
16.4000 -13.6000 35.8000 -76.2000 67.0000 -134.0000 >> B\A ans =
matlab上机练习答案
实验一MA TLAB的基本命令与基本函数１已知矩阵a =11 12 13 1421 22 23 2431 32 33 3441 42 43 44求(1) A(:,1) (2) A(2,:)(3) A(:,2:3) (4) A(2:3,2:3)(5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:)(9) ones(2,2) (10) eye(2)(11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1)(14) diag(A,-1) (15) diag(A,2)2(1)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(6) 求Msin(L)M-13已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0=1.785×10-3，a=0.03368,b=0.000221,求水在0，20，40，80℃时的黏度。
程序如下：miu0=1.785e-3;a=0.03368;b=0.000221;t=0:20:80miu=miu0./ (1+a*t+b*t.^2)（2）一个长管，其内表面半径为a，温度为Ta ；外表面半径为Tb；则其径向和切向应力可分别表示为：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛+----=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛----=r b a b r b a b a a b v T T E r b a b r b a b a a b v T T E b a t b a r ln ln 11)/ln()1(2)(ln ln 1)/ln()1(2)(2222222222ασασ式中r 为管子的径向坐标，E 为管子材料的弹性模量，ɑ为热膨胀系数。
MATLAB上机练习一参考解答
上机练习一参考解答一、实验目的1、 熟悉Matlab 编程2、 体会数学上恒等，算法上不一定恒等二、实验内容1． Using the Taylor polynomial of degree nine and three-digit rounding arithmetic to find an approximationto 5-e by each of the following methods.(A) ∑=--≈905!)5(n n n e , (B) ∑=-≈=9055!5/11n nn e e An approximate value of 5-e correct to three digits is 31074.6-⨯. Which formula, (A) or (B), gives the most accuracy, and why?1） 算法基础利用x e 的Taylor 公式00!!n nk x n n x x e n n ∞===≈∑∑，x -∞<<+∞ （1）及001/1/1/!!n nk x x n n x x e e n n ∞-====≈∑∑，x -∞<<+∞， （2）其中k 是根据精度要求给定的一个参数。
在本题中将k 取为9， x 取为－5或5即可由公式（1）或（2）得到5-e 的近似计算方法（A ）或（B ）。
2） 程序下述程序用公式（A ）及（B ）分别在Matlab 许可精度下及限定在字长为3的算术运算情况下给出5-e的近似计算结果，其中results_1, results_2为用方法（A ）在上述两种情况下的计算结果，err_1, err_2为相应的绝对误差；类似的，results_3, results_4为用方法（B ）在上述两种情况下的计算结果，err_3, err_4为相应的绝对误差；具体程序如下：% Numerical Experiment 1.1 % by Xu Minghua, May 17, 2008 clc; %Initialize the data x=-5; k=9; m=3; %three-digit rounding arithmetic %------------------------------------ % Compute exp(x) by using Method (A) % with the computer precision results_1=1; power_x=1; for i=1:k factor_x=x/i; power_x=power_x*factor_x; results_1=results_1+power_x; end results_1 err_1=abs(exp(x)-results_1)%------------------------------------% Compute exp(x) by using Method (A) % with the 3-digits precisionresults_2=1;power_x=1;for i=1:kfactor_x=digit(x/i,m);power_x=digit(power_x*factor_x,m); results_2=digit(results_2+power_x,m); endresults_2err_2=abs(exp(x)-results_2)%------------------------------------% Compute exp(x) by using Method (B) % with the computer precisiont=-x;results_3=1;power_x=1;for i=1:kfactor_x=t/i;power_x=power_x*factor_x;results_3=results_3+power_x; endresults_3=1/results_3err_3=abs(exp(x)-results_3)%------------------------------------% Compute exp(x) by using Method (B) % with the 3-digits precisiont=-x; results_4=1;power_x=1;for i=1:kfactor_x=digit(t/i,m);power_x=digit(power_x*factor_x,m);results_4=digit(results_4+power_x,m); endresults_4=digit(1/results_4,m)err_4=abs(exp(x)-results_4)%------------------------------------上述主程序用到一个子程序digit.m, digit(x,m)的作用是将x四舍五入成m位数。
MATLAB上机实验1答案
实验1 Matlab 初步一、问题已知矩阵A 、B 、b 如下：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡-------------=031948118763812654286174116470561091143A ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡------=503642237253619129113281510551201187851697236421B []1187531=b应用Matlab 软件进行矩阵输入及各种基本运算。
二、实验目的学会使用Matlab 软件构作已知矩阵对应的行（列）向量组、子矩阵及扩展矩阵，实施矩阵的初等变换及线性无关向量组的正交规范化，确定线性相关相关向量组的一个极大线性无关向量组，且将其余向量用极大线性无关向量组线性表示，并能编辑M 文件来完成所有的实验目的。
三、预备知识1、 线性代数中的矩阵及其初等变换、向量组的线性相关性等知识。
2、 Matlab 软件的相关命令提示如下；（1） 选择A 的第i 行做一个行向量：ai=A(i,:)；（2） 选择A 的第j 行做一个列向量：ai=A(j,:)；（3） 选择A 的某几行、某几列上的交叉元素做A 的子矩阵：A([行号],[列号])；（4） n 阶单位阵：eye(n)；n 阶零矩阵：zeros(n)；（5） 做一个n 维以0或1为元素的索引向量L ，然后取A(:,L)，L 中值为1的对应的列将被取到。
（6） 将非奇异矩阵A 正交规范化，orth(A) ；验证矩阵A 是否为正交阵，只需做A*A'看是否得到单位阵E 。
（7） 两个行向量a1和a2的内积：a1*a2'。
（8） 让A 的第i 行与第j 列互换可用赋值语句：A([i,j],:)=A([j,i],:)；（9）让K乘以A的第i行可用赋值语句：A(i,:)=K*A(i,:)；（10）让A的第i行加上第j行的K倍可用赋值语句：A(i,:)=A(i,:)+K*A(j,:)；（11）求列向量组的A的一个极大线性无关向量组可用命令：rref(A)将A化成阶梯形行的最简形式，其中单位向量对应的列向量即为极大线性无关向量组所含的向量，其它列向量的坐标即为其对应向量用极大线性无关组线性表示的系数。
MATLAB上机操作题
实验一 MATLAB 基本操作一、目的和要求（1） 熟练掌握MA TLAB 变量的使用。
（2） 熟练掌握矩阵的创建。
（3） 熟练掌握MA TLAB 的矩阵和数组的运算。
二、内容和步骤1． 创建矩阵矩阵是包括n m ⨯个元素的矩阵结构，矩阵中的元素可以是实数或复数，单个元素构成的标量以及多个元素构成的行向量和列向量都是矩阵的特殊形式。
下面用多种方式创建矩阵。
（1） 直接输入。
例：>> a=[1 2 3;4 5 6;7 8 9] a =1 2 34 5 6 7 8 9练习：从键盘输入一个44⨯的矩阵，矩阵元素自己定义，变量名定义为m1。
（2） 用from:step:to 方式。
例：>> a=[1:3 ; 4:6 ; 7:9]a =1 2 34 5 6 7 8 9练习：从键盘输入一个向量，用from:step:to 方式，要求从零开始，步长为1.5，到20结束，变量名定义为m2。
（3） 用linspace 函数。
例：>> a=[linspace(1,3,3);linspace(4,6,3);linspace(7,9,3)] a =1 2 3 4 5 6 7 8 9 练习：使用logspace 函数创建π4~0的行向量，有20个元素，查看其元素分布情况，变量名定义为m3。
（4） 使用matlab 内部函数创建矩阵。
例：>> a=ones(3)a =1 1 11 1 11 1 1练习：建立一个44⨯的正态分布矩阵，要求其均值为2，方差为2，变量名定义为m4。
2．矩阵的操作（1）矩阵元素的提取。
例：>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9提取元素8，可输入a(3,2)或a(6)。
练习：建立一个44⨯的对角阵，提取它的第三列元素，并将矩阵的第四列元素改为全1。
（2）矩阵的裁剪和拼接。
例：>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9>> b=[2 3 4;5 6 7;8 9 10]b =2 3 45 6 78 9 10>> c=cat(1,a,b)c =1 2 34 5 67 8 92 3 45 6 78 9 10>> c=cat(2,a,b)c =1 2 3 2 3 44 5 6 5 6 77 8 9 8 9 10>> c(:,3:5)=[]c =1 2 44 5 7 7 8 10练习：建立两个任意的44⨯矩阵，将它们沿着第二维拼接成一个新矩阵，并将新矩阵的第3列和第5列删除。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
要求:每个周二MATLAB基础理论课时上交打印版作业，要求写清楚题目，程序以及运行结果。将得到的结果以截图或者复制粘贴的形式保存到word文档里。
1.建立一个以学号命名的目录(文件夹)，将这个目录添加(扩展)成MATLAB的搜索路径，截图表示。
2.在命令窗口输入：x=rand(2,3)
3.在命令窗口输入：y=rand(2,3)
4.在命令窗口输入：A=sin(x)
5.在命令窗口输入：B=sin(2*y)
6.在命令窗口输入：who
7.在命令窗口输入:whos
8.在命令窗口输入：pi
9.求
10.将所有变量存入matlab1.mat这个文件中
11.读取matlab1.mat中的所有变量
12.清除所有变量