(完整版)matlab上机练习题答案
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上机答案
X
1
2
3
4
5
6
7
8
9
10
Y
16
32
70
142 260 436 682 1010 1432 1960
>> x=1:10;y=[16 32 70 142 260 436 682 1010 1432 1960];
>> p1=polyfit(x,y,1)
>>
p1 =
p2=polyfit(x,y,2),y2=polyval(p2,9.5)
0.01
0.005
0
-0.005
-0.01
-0.015
0
1
2
3
4
5
6
7
（3）大气压强 p 随高度 x 变化的理论公式为
，为验证这一公式，
测得某地大气压强随高度变化的一组数据如表所示。试用插值法和拟合法进行计算并绘图，
看那种方法较为合理，且总误差最小。
高度/m
0
300
600
1000
1500
2000
压强/Pa
equally spaced points between X1 and X2. 以 X1 为首元素，X2 为末元素平均生成 100 个元素的行向量。
LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2.
ans =
pi =
ans =
5
0
3.1416
答：3 次执行的结果不一样。exist（）函数是返回变量搜索顺序的一个函数。在第一次
MATLAB期末上机试题带答案
MATLAB 期末上机考试试题带答案版姓名： 学号： 成绩：1.请实现下图：xyy=sin(x)x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图：x=linspace(0,2*pi,100); y1=sin(x);subplot(2,2,1)plot(x,y1,'k--')grid onxlabel('x')ylabel('y')title('sin(x)')legend('y=sin(x)')y2=cos(x);subplot(2,2,2)plot(x,y2,'r--')grid onxlabel('x')ylabel('y')title('cos(x)')legend('y=cos(x)')y3=tan(x);subplot(2,2,3)plot(x,y3,'k-')grid onxlabel('x')ylabel('y')title('tan(x)')legend('y=tan(x)')y4=cot(x);subplot(2,2,4)plot(x,y4)grid onxlabel('x')ylabel('y')title('cot(x)')legend('y=cot(x)')3.解方程组：a=[3 2 1;1 -1 3;2 4 -4]；b=[7;6;-2] ；x=a\b4.请实现下图：yxx=linspace(0,4*pi,1000);y1=sin(x);y2=sin(2*x);plot(x,y1,'--',x,y2,'b*')grid onxlabel('x');ylabel('y');title('耿蒙蒙')legend('sin(x)','sin(2*x)')5.请在x，y在（-2，2）内的z=xexp (-x2-y2) 绘制网格图[x,y]=meshgrid(-2:0.1:2);z=x.*exp (-x.^2-y.^2);mesh(x,y,z)6.请实现peaks函数：-55xPeaksy[x,y]=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z) surf(x,y,z) shading flataxis([-3 3 -3 3 -8 8])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。
matlab试题及答案
matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。
答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。
矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。
矩阵的转置使用单引号`'`，例如`A'`。
矩阵的求逆使用`inv`函数，例如`inv(A)`。
2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。
`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。
A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。
答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。
MATLAB上机答案
熟悉Matlab工作环境1、熟悉Matlab的5个基本窗口思考题：(1)变量如何声明，变量名须遵守什么规则、是否区分大小写。
答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。
变量名要遵守以下几条规则：变量名必须以字母开头，只能由字母、数字或下划线组成。
变量名区分大小写。
变量名不能超过63个字符。
关键字不能作为变量名。
最好不要用特殊常量作为变量名。
(2)试说明分号、逗号、冒号的用法。
分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。
逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。
冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素(3) linspace()称为“线性等分”函数，说明它的用法。
LINSPACE Linearly spaced vector. 线性等分函数LINSPACE(X1, X2) generates a row vector of 100 linearlyequally spaced points between X1 and X2.以X1为首元素，X2为末元素平均生成100个元素的行向量。
LINSPACE(X1, X2, N) generates N points between X1 and X2.For N < 2, LINSPACE returns X2.以X1为首元素，X2为末元素平均生成n个元素的行向量。
如果n<2,返回X2。
Class support for inputs X1,X2:float: double, single数据类型：单精度、双精度浮点型。
(4)说明函数ones ()、zeros ()、eye ()的用法。
ones ()生成全1矩阵。
答：3次执行的结果不一样。
exist ()函数是返回变量搜索顺序的一个函数。
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上机练习一答案
Matlab上机练习⼀答案Matlab 上机练习⼀班级学号姓名按要求完成题⽬，并写下指令和运⾏结果。
（不需要画图）1、求下列联⽴⽅程的解=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59942、设 ??++=)1(sin 35.0cos 2x x x y ，把x=0-2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')3、产⽣8×6阶的正态分布随机数矩阵R1, 求其各列的平均值、和、乘积。
并求该矩阵全体数的平均值。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)或者u=reshape(a,1,48);p1=mean(u)p2=var(u)4、绘制)(222y x e x z +-=在定义域x=[-2,2],y=[-2,2]内的曲⾯。
（利⽤meshgrid ）x=-2:2;y=x;[X,Y]= meshgrid(x,y);Z=X^2*exp(-(X^2+Y^2));mesh(X,Y,Z)5、求代数⽅程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。
（利⽤roots 函数）p=[3 4 7 2 9 12];roots(p)6、把1开五次⽅，并求其全部五个根。
MATLAB上机答案
一熟悉Matlab工作环境1、熟悉Matlab的5个基本窗口思考题：（1）变量如何声明，变量名须遵守什么规则、是否区分大小写。
答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。
变量名要遵守以下几条规则：➢变量名必须以字母开头，只能由字母、数字或下划线组成。
➢变量名区分大小写。
➢变量名不能超过63个字符。
➢关键字不能作为变量名。
➢最好不要用特殊常量作为变量名。
（2）试说明分号、逗号、冒号的用法。
分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。
逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。
冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。
（3）linspace（）称为“线性等分”函数，说明它的用法。
LINSPACE Linearly spaced vector. 线性等分函数LINSPACE(X1, X2) generates a row vector of 100 linearlyequally spaced points between X1 and X2.以X1为首元素，X2为末元素平均生成100个元素的行向量。
LINSPACE(X1, X2, N) generates N points between X1 and X2.For N < 2, LINSPACE returns X2.以X1为首元素，X2为末元素平均生成n个元素的行向量。
如果n＜2，返回X2。
Class support for inputs X1,X2:float: double, single数据类型：单精度、双精度浮点型。
（4）说明函数ones（）、zeros（）、eye（）的用法。
ones（）生成全1矩阵。
zeros（）生成全0矩阵。
eye（）生成单位矩阵。
2、Matlab的数值显示格式思考题：（1）3次执行exist（’pi’）的结果一样吗？如果不一样，试解释为什么？>> pians =3.1416 >> sin(pi); >> exist('pi') ans =5 >> pi=0;>> exist('pi')ans =1>> pipi =>> clear>> exist('pi')ans =5>> pians =3.1416答：3次执行的结果不一样。
MATLAB上机习题一 解答
9. 判断下面语句的运算结果，并与 MATLAB 计算结果做比较： 1） 4 < 20 2） 4 ≤ 20 3） 4 == 20 4） 4 ~ = 20 5） ' b ' < ' B ' 解答： 1）true 2）true 3）false 4）true 5）false
注意，字符比较实际上比较的是字符的 ASCII 码的大小。
提示，可以使用 isvarname 对上述变量名进行检验。
解答：1、2、5、6、9、10 是合法的，其它不能作为变量名。 注意，第 6 问的函写脚本程序，命名为 barbell.m，完成如下计算： 如图 1 所示的杠铃，
每个球的半径是 10cm，两个求直接链接杆的长度是 15cm，杆的直径为 1cm，计 算杠铃的体积和表面积；使用 disp 函数显示计算结果。 （提示：disp 函数的用法 参考教材的例 2-31， 面积计算时相交部分球冠可近似为圆， 体积计算时可忽略重 叠的球冠体积） 解答： barbell.m（拷贝文件内容到如下方框中） ： %注意，将圆柱与球相交部分面积近似为圆形 clear all clc R=10; %球半径,单位cm l=15; %杆长度,单位cm d=1; %杆直径,单位为cm
n2a + p nRT (V − nb ) = V2
该方程新增了 a 和 b 两个变量，用来表示纯净气体的属性。 1）编写脚本程序，命名为 state_equation.m，利用如下数据： 压强，p 摩尔数，n 体积，V a b 理想气体常数，R 220 mbar 2 mol 1L 5.536 L2bar/mol2 0.03049 L/mol 8.314 J/(mol.K)
1. 运行 MATLAB 软件，拖放、关闭界面上的子窗口，并恢复到原始视图。 解答：略。注意采用下图方式恢复视图
matlab上机实验答案
matlab上机实验答案三、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数⾏提取出来，赋给B 矩阵，⽤magic(8)A =命令⽣成A 矩阵，⽤上述命令检验⼀下结果是不是正确。
>> A=magic(8) A =64 2 3 61 60 6 7 57 9 55 54 12 13 51 50 16 17 47 46 20 21 43 42 24 40 26 27 37 36 30 31 33 32 34 35 29 28 38 39 25 41 2322 44 45 19 18 48 49 15 14 52 53 11 10 56 8 58 59 5 4 62 63 1 >> B=A(2:2:end,:)B =9 55 54 12 13 51 50 16 40 26 27 37 36 30 31 33 41 23 22 44 45 19 18 48 8 58 59 5 4 62 63 1五、选择合适的步距绘制出下⾯的图形。
（1）)/1sin(t ，其中)1,1(-∈t ；（2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 。
1．>> t=[-1:0.0001:1];y=sin(1./t);plot(t,y) Warning: Divide by zero. >>2．>> t=[-pi:0.001:pi];y=sin(tan(t))-tan(sin(t));plot(t,y) >>七、试求出如下极限。
（1）x xx x 1)93(lim +∞→；（2）11lim00-+→→xy xy y x ；（3）22)()cos(1lim222200yx y x ey x y x +→→++-。
（1）>> syms x;f=(3^x+9^x)^(1/x);limit(f,x,inf)ans =9（2）>> syms x y;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0) ans =2（3）>> syms x y;f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2);limit(limit(f,x,0),y,0) ans =0九、假设?-=xytt ey x f 0d ),(2，试求222222yf yx f xf y x ??+-??。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(完整版)matlab 上机练习题答案1.计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积>〉 a=［6 9 3;2 7 5]; 〉〉 b=[2 4 1;4 6 8］; 〉〉 a 。
*b ans =12 36 3 8 42 402。
对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
〉〉 A=［4 9 2；7 6 4；3 5 7]; 〉> B=[37 26 28］’； >〉 X=A\B X = -0.5118 4.0427 1.33183。
⎥⎦⎤⎢⎣⎡-=463521a ，⎥⎦⎤⎢⎣⎡-=263478b ,观察a 与b 之间的六种关系运算的结果 >〉 a=[1 2 3;4 5 6]; 〉> b=［8 –7 4；3 6 2]； >〉 a 〉b ans =0 1 0 1 0 1 〉> a 〉=b ans =0 1 0 1 0 1 >> a 〈b ans =1 0 1 0 1 0 〉> a<=b ans =1 0 1 0 1 0 >〉 a==b ans =0 0 0 0 0 0 〉〉 a~=b ans =1 1 1 1 1 14计算多项式乘法（x 2+2x +2)（x 2+5x +4)>> c=conv （[1 2 2]，[1 5 4］) c =1 7 16 18 8 5计算多项式除法（3x 3+13x 2+6x +8）/(x +4） 〉〉 d=deconv （［3 13 6 8],［1 4］) d =3 1 26求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解〉〉 a=［2 4 7 4;9 3 5 6］； 〉> b=[8 5］’; 〉> x=pinv （a ）*b x =—0.2151 0.4459 0。
7949 0。
27077用符号函数法求解方程a t 2+b*t +c=0>〉 r=solve （’a＊t^2+b ＊t+c=0’,’t'） r =[ 1/2/a ＊（—b+（b^2—4＊a ＊c ）^(1/2))] [ 1/2/a ＊(—b-(b^2-4＊a*c)^(1/2))］ 8求矩阵⎥⎦⎤⎢⎣⎡=22211211a a a a A 的行列式值、逆和特征根 〉〉 syms a11 a12 a21 a22； >〉 A=［a11，a12;a21，a22]〉〉 AD=det(A) ％ 行列式 〉〉 AI=inv(A ） % 逆 〉〉 AE=eig （A) % 特征值 A =[ a11, a12］ [ a21, a22］ AD =a11*a22-a12＊a21 AI =［ —a22/(—a11＊a22+a12＊a21)， a12/（-a11*a22+a12*a21)］［ a21/（-a11＊a22+a12*a21)， -a11/(—a11*a22+a12＊a21)] AE =［ 1/2*a11+1/2*a22+1/2＊（a11^2—2＊a11＊a22+a22^2+4＊a12＊a21）^(1/2）][ 1/2*a11+1/2＊a22-1/2＊(a11^2—2＊a11＊a22+a22^2+4*a12*a21)^(1/2)］ 9因式分解:6555234-++-x x x x〉〉 syms x;(完整版)matlab 上机练习题答案>〉 f=x^4-5*x^3+5*x^2+5＊x-6; >〉 factor(f) ans =（x-1）＊(x —2)＊（x —3）＊(x+1)10⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x e x x af ax ，用符号微分求df/dx 。
〉> syms a x;>〉 f=[a, x^2, 1/x ； exp （a*x ）， log （x), sin(x)］; 〉〉 df=diff （f ） df =[ 0, 2＊x, -1/x^2] ［ a ＊exp(a*x)， 1/x, cos(x)］ 11求.arctan 2⎰xdx x 和.)(102⎰-dx x x程序如下: 〉〉 syms x;>> int （x^2＊atan （x ）,'x’） ans =1/3＊x^3＊atan(x ）-1/6*x^2+1/6＊log （x^2+1) 〉> simple(ans) 结果如下: ans =1/3＊x^3*atan(x)-1/6＊x^2+1/6＊log （x^2+1） -—-—-—-- 程序如下: >> syms x ；〉〉 int （x —x^2，’x’，0，1） 结果如下： ans = 1/612微分方程x e y y y x 2cos 52=+'-''的通解. 程序如下： 〉〉 syms x y ；〉〉 dsolve （'D2y —2＊Dy+5＊y=exp （x ）＊cos(2*x)’） ans =exp （t)*sin （2*t)*C2+exp （t ）＊cos （2*t)＊C1+1/5＊exp （x)＊cos （2＊x ）13求代数方程组⎪⎩⎪⎨⎧=+=++002y x c by ax 关于x ，y 的解〉〉 S=solve('a*x^2+b*y+c=0'，'b*x+c=0'，’x’,’y');〉> disp('S 。
x=') ， disp （S.x) 〉〉 disp('S.y=') ， disp(S 。
y ） S 。
x= —c/b S 。
y=-c ＊(a ＊c+b^2)/b^314求方程组⎪⎪⎩⎪⎪⎨⎧=---=++=+--=--+.0532,0375,023,02432143243214321x x x x x x x x x x x x x x x 和⎪⎪⎩⎪⎪⎨⎧=+-=+-=++-=++-53323221242143143214321x x x x x x x x x x x x x x程序如下：>> a1=[1 1 —2 —1]； >〉 a2=[3 —1 —1 2］; 〉> a3=［0 5 7 3］； 〉> a4=[2 —3 -5 -1]；>〉 linsolve([a1;a2；a3;a4］,[0 0 0 0］’）ans =0 0 0 0该方程组无解． --—-—-- 程序如下：〉〉 a1=[1 —1 2 1]; >〉 a2=［2 —1 1 2］; 〉> a3=[1 0 —1 1]; >〉 a4=[3 —1 0 3］；〉〉 null(［a1;a2;a3；a4］,'r’) ans = 1 -1 3 0 1 0 0 1所以该方程组的通解为： （其中k1 k2为任意常数） 15符号函数绘图法绘制函数x=sin(3t ）cos （t)，y=sin(3t ）sin(t)的图形， t 的变化范围为［0,2］〉> syms t〉〉 ezplot （sin （3*t ）*cos(t),sin(3*t)*sin(t ），[0,2*pi ］）(完整版)matlab 上机练习题答案16有一组测量数据满足-at e =y ，t 的变化范围为0～10，用不同的线型和标记点画出a=0.1、a=0。
2和a=0.5三种情况下的曲线，并加入标题和图列框（用代码形式生成)>〉 t=0:0.5：10；〉> y1=exp(—0。
1＊t)； 〉> y2=exp(—0.2＊t ）； 〉〉 y3=exp(—0。
5＊t); >> plot(t ，y1，’—ob'，t ，y2,’：*r'，t ，y3,'—。
^g ’）〉〉 title('\ity\rm=e^{-\itat ｝’,'FontSize ’,12）>〉 legend （'a=0。
1’，'a=0.2'，'a=0。
5’）17 x= [66 49 71 56 38]，绘制饼图并将第五个切块分离 >> x=［66 49 71 56 38］;〉〉 L=[0 0 0 0 1］； >〉 pie （x,L) 18 22y xxe z --=，当x 和y 的取值范围均为—2到2时,用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图〉> [x,y]=meshgrid([-2:。
2:2］);>〉 z=x 。
*exp （-x 。
^2-y.^2）; 〉> mesh （x ，y ，z) 〉> subplot （2,2，1）， plot3（x,y,z ） >> title(’plot3 （x ，y ，z ）’)〉> subplot （2，2，2）, mesh （x,y,z) >> title(’mesh (x,y ，z)'）>> subplot （2，2,3）, surf （x ，y,z) >> title('surf (x,y ，z)’）>〉 subplot(2，2，4), surf(x ，y,z), shading interp>> title(’surf （x,y,z ）， shading interp ’）19 在区间]1,1[-画出函数xy 1sin=的图形 程序如下：〉> fplot(’sin(1/x）',[-pi/12，pi/12]) 〉> grid 〉〉 title('graph of sin(1/x)'） 结果如下：(完整版)matlab 上机练习题答案graph of sin(1/x)20 分别画出坐标为)10,,2,1(),4,(),,(3222 =+i i i i i i 的散点图, 并画出折线图程序如下：>> for i=1：10plot(i,i 。
^2，’。
'）; hold onplot(i 。
^2,4*i.^2+i.^3，’.'); end>> x=1:10； 〉〉 y=x 。
^2; >> plot(x ，y);〉> plot(x.^2,4＊x 。
^2+x 。
^3）； 〉> axis （[0,105,0，1450］)21在区间]4,4[-上作出函数xx xx x f --=39)(的图形, 并计算)(lim x f x ∞→ 和 ).(lim 1x f x →程序如下: 〉〉 syms x;〉> f=(x^3—9*x)/（x^3-x ）； >〉 limit （f ，x,inf)ans = 1〉> limit （f ，x,1）ans = NaN绘制()f x 的图形程序如下:〉〉 f=inline(’(x.^3—9*x ）。