MATLAB 测试题
- 格式:doc
- 大小:130.50 KB
- 文档页数:6
第一章测试1.MATLAB、FORTRAN、MAPLE并称为三大数学软件,以研究和解决各种具体的数学问题和工程问题。
A:对B:错答案:B2.下列属于MATLAB独特优势的是()。
A:强大的科学计算及数据处理能力B:简单易用的程序语言C:友好的工作平台和编程环境D:出色的图形处理功能答案:ABCD3.MATLAB界面有()个窗口。
A:4B:3C:5D:6答案:D4.MATLAB、FORTRAN、Mathematica并称为三大数学软件,以研究和解决各种具体的数学问题和工程问题。
A:错B:对答案:A5.在MATLAB中,可以在命令行窗口中通过帮助命令来查询帮助信息,最常用的帮助命令是()。
A:helpwinB:helpC:whichD:what答案:B第二章测试1.下列可作为MATLAB合法变量名的是()。
A:@hB:xyz_2aC:合计D:123答案:B2.设A=[1,2;3,4],B=[5,6;7,8],则A*B=()。
A:[19,22; 21,32]B:[19,22;43,50]C:[5,12;43,50]D:[5,12;21,32]答案:B3.输⼊字符串时,要⼊()将字符括起来。
A:{ }B:C:[ ]D:” ”答案:B4.使⼊语句t=0:7⼊成的是7个元素的向量。
()A:对B:错答案:B5.下列语句中正确的是()。
A:x==y==3B:y=3,x=yC:x=y=3D:x=y==3答案:ABD第三章测试1.sym函数可以建立一个符号量。
()A:错B:对答案:B2.下列指令可用于对符号运算表达式进行化简的是()。
A:expand(S)B:factor(S)C:collect(S)D:horner(S)答案:ABCD3.MATLAB用于符号常微分方程求解的函数是()。
A:solverB:solveC:dsolverD:dsolve答案:D4.ezplot函数可用于绘制符号函数的三维曲线。
()A:错B:对答案:A5.MATLAB将函数展开为幂级数,所使用的函数是()。
PART I 教材复习要点第一章 MATLAB 入门测试1.11.MATLAB 的命令窗口的作用是什么?编辑/调试窗口?图象窗口? MATLAB 命令窗口是MATLAB 起动时第一个看到的窗口,用户可以在命令窗口提示符">>"后面输入命令,这些命令会被立即执行。
命令窗口也可以用来执行M 文件。
编辑/调试窗口是用来新建,修改或调试M 文件的。
图像窗口用来显示MATLAB 的图形输出。
2.列出几种不同的得到MATLAB 帮助的方法。
(1)在命令窗口中输入help <command_name>,本命令将会在命令窗口中显示关于些命令的有信息。
(2)在命令窗口中输入lookfor <keyword>,本命令将会在命令窗口中显示所有的在第一注释行中包含该关键字的命令和函数。
(3)通过在命令窗口输入helpwin 或helpdesk 启动帮助浏览器,或者是在启动板中选择"Help"。
帮助浏览器包含了基于超文本的MATLAB 所有特性的描述,HTML 或PDF 格式的在线手册,这是MATLAB 最全面的帮助资源。
3.什么是工作区?在同一工作区内,你如何决定它里面存储了什么?工作区是命令、M 文件或函数执行时被MATLAB 使用的变量或数组的收集器,所有命令都在命令窗口(所有的脚本文件也是从命令窗口执行)共享公共工作区,因此它们也共享所有变量,工作区的内容可以通过whos 命令来查看,或者通过工作区浏览器来图形化地查看。
4.你怎样清空MATLAB 工作区内的内容?要清除工作区的内容,只需在命令窗口中输入clear 或clear variables 即可。
第二章MATLAB 基础测试2.11.数组,矩阵,向量有什么区别?数组是在内存中被组织成行和列的数据集合,只有一个名称,数据要通过在数组名后面圆括号里加上表示数据所在行和列的数字来访问。
matlab考试题及答案【篇一:matlab期末考试试卷及参考答案】lass=txt>1.设有程序a=[2,0,1,0;1,0,1,5;0, 6 ,2, 3];b=[1,2,-1;1,3,1];c=a(1:3,[1,2]);d=size(c)- length(b(2,:));f= c.*b’将下列命令的运行结果填在横线上0*、02. a=[2,3,2,1;1,2,2,1; 1,2,3,-2];a1=sum(aa(2,1));a(2,:)=[],a2=a.^2-4a1=; a2=;3.p=[1,2,5,4;2,2,1,3];max(p;mean(p3??x?sint(1).?(0?t?2?);3??y?cost二、(本题满分10分)试编写计算程序画出曲线与曲面图形(2).z?(x2?1)e?x?y22?xy,?2?x?2,?2?y?2.2%1 %2t=0:0.1: 2*pi; [x,y]=meshgrid(-2:0.1:2);x= sin(t).^3;z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; y= cos (t).^3;mesh(x,y,z) plot(x,y)三、(本题满分12分)编写程序解决下列问题:(1)试产生100?4阶矩阵使其元素在0—100之间的随机整数,可将此矩阵理解为100名同学的四门课程(课程1, 课程2, 课程3, 课程4)的成绩. (2)计算出每门课程和每位同学的平均成绩; (3)统计出总的优秀人次(90分及以上为优秀);(4)统计出四门课程成绩至少有两门课程不及格的人数. 解:(1) a=fix((100-0+1)*rand(100,4));(2) mean(a)mean(a,2) (3)sum(sum(a=90,2)==4)(4)sum(sum(a60,2)=2)四、(本题满10分)n设f(x)??k?1(?1)kk?1sinkx,返回该函数并保存y?fun1(x,n));(1)试编写函数式文件程序(程序的第一句为:function(2)用子图画出该函数在[0,2?]内图形(分别取n?20,50).(1). function y=fun1(x,n)(2). x=0:0.1:2*pi; y=0; y1= fun1 (x,20); for k=1:n y2= fun1 (x,50);y=y+ ((-1)^(k+1)*sin(k*x))/k; subplot(2,1,1),plot(x,y1)endsubplot(2,1,2),plot(x,y2) y五、(本题满分10分)试利用微分方程的数值解法写出下列微分方程的求解步骤和matlab程序. x???3tx??x?t,x(0)?2,x?(0)?0.(t的变化区间为[0,10]). 第一步:先将高阶微分方程转化为一阶微分方程组:选择状态变量x1?x?,x2?x,则原方程化为???3tx1?x2?t,?x1???x1.?x2x1(0)?0,x2(0)?2dxdt = [-3*t*x(1)+ x(2)+t; x(1)];第三步:求解微分方程,命令如下:[t,x]=ode45(@fun2,[0,10],[0;2])六、运用符号运算功能计算1.求方程ex下列各题,:?x2?2的根;?xsinxa???xe??lnx?da21?,并计算,?a(x)dx,a;dxs?x??;t?x?2e,x(0)?x?(0)?x??(0)?x???(0)?1的解2.生成符号矩阵3.计算极限limcosx?x(esinxcosx?x)x?04.用dsolve函数求解微分方程x六、参考程序: 1. solve(exp(x)=x^2+2) 2.syms x sa=[x*sin(x),log(x);exp(-x),1/(s+x)] diff(a) int(a) a*a 3.syms x limit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1))4.dsolve(d4x+x=2*exp(t),x(0)=1,dx(0)=1,d2x(0)=1,d3x(0)=1)七.用fsolve函数以及符号运算功能22??4x1?x2?2x1x2?x2?2?0?22??2x1?x2?3x1x2?3?0求解下面的非线性方程组:解一:function y=fc(x)y(1)=4*x(1).^2+x(2).^2+2*x(1).*x(2)-x(2)-2y(2)=2*x(1).^2+x(2).^2+3*x(1).*x(2)-3 y=[y(1) y(2)];x0=[1 1];fsolve(@fc,x0)解二: syms x1 x2[x1 x2]=solve(‘4*x1^2+ x2^2+2* x1* x2- x2-2=0’,’ 2*x1^2+x2^2+3* x1* x2-3=0’,’ x1’,’ x2’)《matlab程序设计》参考程序一、 (30分)1. d =0 -1 f =2 0 2 0 0 6 2. a1 =1 3 3 0 a2 =0 5 0 -3-3 0 5 0 3. max(p,[ ],1)=2 2 5 4mean(p,2)=32 4. x=15 5.0 1 4 -1 0 1 -4 -1 0 6.z =2010 7.285二、(10分)%1t=0:0.1: 2*pi; x= sin(t).^3; y= cos (t).^3; plot(x,y) %2[x,y]=meshgrid(-2:0.1:2);z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; mesh(x,y,z)三、(16分)syms xlimit((x-1)*log(x)/(exp(x-1)-1),x,1) %2solve(exp(x)-3*cos(x)+1=0) %3 syms xm=[ exp(-x), x;sin(x),exp(-0.1*x)* cos(3*x)] diff(m,2) int(m,x,0,pi) %4 syms t xdsolve(dx-6*x=exp(-t) ,x(0)=1)四、(12分)x=0.5:9.5;y=[ 0.16, 0.40, 0.66, 0.75, 0.88, 1.0, 1.1, 1.1, 1.2, 1.3];x1=1./x;y1=1./y; p=polyfit(x1,y1,1); a=p(2); b= p(1);ye=x./(a*x+b); plot(x,y,ro, x, ye, b-)grid onxlabel(自变量x ) ylabel(因变量y )title(因变量 y与自变量 x的拟合曲线图)五、(12分)%参考程序cj=round(normrnd(70,10,100,4)); for i=1:100 for j=1:4if cj(i,j)100 cj(i,j)=100;elseif cj(i,j)0 cj(i,j)=0;end end endkmean=round(mean(cj))disp(课程1课程2课程3课程4)disp(sprintf( %4d%7d%8d%8d ,kmean(1),kmean(2),kmean(3),k mean(4))) studentmean=round(mean(cj,2));studentnumber=[1:100];[studentnumber, studentmean] c=sum(sum(cj=90))【篇二:matlab习题及答案】241?b?1,计算a??与?468?的数组乘积。
1、求以下变量的值,并在MATLAB 中验证。
( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ;a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 82、下列运算是否合法,为什么?如合法,结果是多少? >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 25 8 11 >> result4=b*d result4 =31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13-5 -8 7>> result6=a.*b result6 =2 8 -34 15 30 >> result7=a./b result7 =0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 =2 4 816 32 64 3、用MATLAB 求解下面的的方程组。
(1)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----017413231511222315921274321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332821w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6]>> B2=[1;8;3;5] >> x2=inv(A1)*B24、已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=1323151122231592127A (1)求矩阵A 的秩(rank) (2)求矩阵A 的行列式(determinant) (3)求矩阵A 的逆(inverse) (4)求矩阵A 的特征值及特征向量(eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 5、109910101022222++++==---=∑n n y ,求y=?(运行format long g 命令后,查看y 的值) m1=0; for m=-10:10 m1=m1+2^m; end m1 m1 = 2047.9990234375 6、求分段函数的值。
第二题求解系统的模态的matlab程序:m=1200j=1100k1=3000k2=1500l1=2l2=1K=[k1+k2 -k1*l1+k2*l2;k1*l1+k2*l2 -k1*l1^2+k2*l2^2] M=[m 0;0 j][eigve,eigva]=eig(inv(M)*K)运行结果为:eigve =0.8492 0.32360.5281 0.9462eigva =1.4178 00 -7.2133龙格库塔法求解程序:function dy=kjt001(t,y)w1=1;w2=1F=1;T=1a1=1200;b1=1100;a2=4500;b2=7500;a3=-4500;b3=-10500;dy = zeros(4,1);dy(1)=y(2);dy(2)=(F*sin(w1*t)-a2*y(2)-a3*y(3))/a1;dy(3)=y(4);dy(4)=(T*sin(w2*t)-b2*y(2)-b3*y(3))/b1;%[t,y] = ode45(@kjt001,[0 30],[0 0 0 0]);%plot(t,y(:,2),'-',t,y(:,4),'-.')求解结果为:第三题中心差分法求解系统方程,matlab运行程序为:clear%输入参数%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=linspace(0,10,290);P=[2*sin(2*t); -2*cos(3*t); sin(3.5*t)];k=50;m=2;time=10; %积分时间dim=3; %维数C=zeros(3,3);M=diag([m m m],0); %质量K=k*[2 -1 0;-1 2 -1;0 -1 1]; %刚度%求解特征值问题%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [eigve,eigva]=eig(inv(M)*K);%%%%%%%%%%%%%%CENTRAL INTEGRATION%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% calculate the time stepnum=1; %时间步长取临界步长的 1/num倍或者num倍check=1; %求解模式:0-时间步长小于临界时间步长;(给模式求解了num=1,2)%1-时间步长大于临界时间步长 (该模式求解了num=2,5,7)if check==0dt=0.1elsedt=num/((eigva(3,3))^0.5*pi);endnn=floor(time/dt+1)-1;%INITIAL CONDITION AND STARTING PROCEDURE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%a=zeros(3,nn);a(:,1)=0;da0=[0 0 0]';dv0=inv(M)*(P(:,1)-K*a(:,1)-C*a(:,1));a_dt=a(:,1)-da0*dt+0.5*dv0*dt^2;%%%%%%%%%% THE CONST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%c0=1/(dt)^2;c1=1/(2*dt);c2=2*c0;c3=1./c2;EM=c0*M+c1*C; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%for i=2:nnif i==2Q=P(:,i)-(K-c2*M)*a(:,i-1)-(c0*M-c1*C)*a_dt;elseQ=P(:,i)-(K-c2*M)*a(:,i-1)-(c0*M-c1*C)*a(:,i-2);enda(:,i)=inv(EM)*Q;endfiguretim=[0:dt:(nn-1)*dt];plot(tim,a(1,:),tim,a(2,:),tim,a(3,:))xlabel('t/s (时间步长为临界步长的num倍)')ylabel('distance/m')title('位移时程响应-central difference')legend('mass 1','mass 2','mass 3')运行结果为:。
《Matlab 基础》测试题2学号: 姓名:要求:将解答命令及结果写在答题纸上,试卷当场回收. 同时将详细的过程及运行结果保存在word 文档中,在本周五12:00之前发送至邮箱gaoxiawang@,邮件主题为: 考试,附件文件名为:学号+姓名,要求试卷解答要与电子文本一致.另外请用自己的电子邮箱发送.一.建立如下矩阵(用命令方式或编程方式,而非直接输入)?1010101010101010101010⨯--⎛⎫ ⎪-- ⎪ ⎪ ⎪ ⎪--⎝⎭二. 设有对称矩阵 ⎪⎪⎪⎭⎫ ⎝⎛=1849424942A , 求其特征根a 及对应的特征向量v . 说明特征根的和仍为原矩阵的迹, 特征根的积仍为原矩阵的行列式, 并验证A*v 1-a 1*v 1=0.三、数字图像处理中使用的Butterworth 低通滤波器的数学模型为2/),(11),(D v u D v u H n +=,其中,2020)()(),(v v u u v u D -+-=,0D 为给定的区域半径,0u 和0v 为区域中心. 假设2,2000==n D ,试绘制该滤波器图形.四. 编写函数M 文件,描述如下分段函数:2sin ,0()1,0x x x f x x x e x ⎧⎪>=⎨⎪+≤⎩要求编写的函数文件对向量的输入能够产生向量的输出;利用plot 命令作出函数)(x f 在区间]8,8[-上的图形;利用quad 命令求定积分⎰-32)(dx x f .五. 假如你有一组实测数据,例如:x=[75 86 95 108 112 116 135 151 155 160 163167 171 178 185];y=[10 12 15 17 20 22 35 41 48 50 5154 59 66 75];1. 对上述数据做二次多项式拟合.并写出拟合得到的函数表达式.2. 做出数据点与拟合曲线的图形.3. 计算绝对误差.4. 求出在数据点 80 与100 处的估计值.。
2016 /2017第 1 页共8 页第 2 页 共 8 页13. 已知A 是矩阵,求A 的对角矩阵B 的命令是_______,求A 的下三角矩阵C 的命令是_______。
14. 生成正态分布4⨯6随机矩阵A 的命令是 ,删除矩阵A 的第一列和第三列命令是 __________________15. 已知A=[1 5 3;4 2 6; 11 1 0]; A(7)=__________。
16. A=[1 3 2 4; 5 4 6 7]; reshape(A,4,2); MATLAB 命令窗口执行上述命令的结果是: ___________ ___.17.已知A=[1,2;3,1];B=[1,0;0,1];,则A~=B=___________.18. A=[4,-15,15, 6, 10;56, 17, 0 ,-45,0]; x=find(A>=10&A<=20),则x=_____________。
19 已知A=[9 6 0;1 3 0];B=[1 4 3;1 5 0]; 写出下列各指令运行的结果。
A~=B 的运行结果_______;20 已知A=[0 9 6;1 0 3];B=[1 4 3;1 0 0]; 写出下列各指令运行的结果。
A./B 的运行结果_______;21. 产生5阶全0方阵的命令为_______;产生5阶单位方阵的命令为 。
22. A=rand(4,3); b=size(A); c=length(A);则b= 和c=_______。
23. 用if 判断语句判断80≥x>60,在MATLAB 中if 语言后的判断应写为_______。
24. 为了使两个plot 的图形在同一个坐标显示,可以使用_______命令进行图形保持;可以使用_______命令为图形添加网格。
25. MATLAB 的工作空间中有三个变量v1,v2,v3,写出把它们保存到文件my_data.mat 中的指令_______;26. MATLAB 中clf 用于 、clc 用于 、 clear 用于 。
matlab试题及答案一、单项选择题(每题2分,共10分)1. MATLAB中用于创建向量的命令是:A. vectorB. arrayC. linspaceD. colon答案:D2. 在MATLAB中,以下哪个函数用于计算矩阵的行列式?A. detB. rankC. invD. eig答案:A3. MATLAB中用于进行矩阵转置的运算符是:A. 'B. .C. ^D. !答案:A4. 若A是一个3x3的矩阵,执行命令A(2,:)=0;后,矩阵A的第二行将变为:A. [0 0 0]B. [1 0 0]C. [0 1 0]D. [0 0 1]答案:A5. MATLAB中,以下哪个函数用于绘制三维曲面图?A. plotB. surfC. barD. hist答案:B二、填空题(每题3分,共15分)1. MATLAB中,使用________命令可以清除所有变量。
答案:clear2. 若要在MATLAB中创建一个从0到1的100个元素的向量,可以使用命令________。
答案:linspace(0,1,100)3. MATLAB中,使用________函数可以计算矩阵的特征值。
答案:eig4. 在MATLAB中,________函数用于计算两个矩阵的点乘。
答案:dot5. 若要在MATLAB中绘制一个圆,可以使用________函数。
答案:plot三、简答题(每题5分,共20分)1. 请解释MATLAB中矩阵索引的概念。
答案:在MATLAB中,矩阵索引指的是通过行号和列号来访问矩阵中特定元素的过程。
例如,A(2,3)表示访问矩阵A的第二行第三列的元素。
2. MATLAB中如何实现矩阵的元素乘法?答案:在MATLAB中,矩阵的元素乘法可以通过使用点乘运算符(.*)来实现。
例如,C = A .* B,其中A和B是相同大小的矩阵。
3. 请说明MATLAB中如何使用循环结构。
答案:MATLAB中可以使用for循环和while循环两种循环结构。
2014年春季MATLAB 基础测试作业题
学号:2012101204
姓名:黄辉
说明:要求用A4纸打印,第13周周1晚19:00交到L2408教室; 提前交,可以交到L1602.
1、 与其它高级语言相比,如下四点是否都是有MATLAB 的特点?
你觉得它还有哪些显著特点? (10分)
① 可靠的数值计算和符号计算功能 ② 强大的绘图功能 ③ 简单易学的语言体系
④ 为数众多的MATLAB 应用工具箱
答:(不少于100字)
以上四点均为 MATLAB 的特点。
其他显著特点有 :
1、MATLAB 支持与 C++等多种语言的混合编程;
2、MATLAB 确切说不是一个程序设计语言,而是一个数学运算工具。
3、MATLAB 里数据存储的基本单元是矩阵,即使一个数,内部也是看做一个矩阵。
数组也是矩阵。
所以对矩阵的运算是最快的。
4、MATLAB 的特点就是它将矩阵作为基本存储单元。
5、MATLAB 的工具箱也很丰富,在图像处理、信号处理、小波、人工智能、经济数学等方面的工具箱里的工具、示例非常多,功能非常强大。
2、 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥
⎦
⎤
⎢⎢⎢⎣⎡638976242之和。
(5分)
命令如下,请运行并给出结果
>> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6]; >> a+b 结果:ans = 7
7 7
9 14 13
15 12 14 3、 计算⎥⎦⎤⎢
⎣⎡=572396a 与⎥⎦
⎤
⎢⎣⎡=864142b 的数组乘积。
(5分)
命令如下,请运行并给出结果
>> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; >> a.*b 结果:ans =12
36 3
8 42 40
4、对于B AX =,如果⎥⎥
⎥⎦
⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤
⎢⎢⎢⎣⎡=282637B ,求解X 。
(5分)
命令如下,请运行并给出结果
>> A=[4 9 2;7 6 4;3 5 7]; >> B=[37 26 28]’; >> X=A\B 结果:X =
-0.5118483412322 4.0426540284360 1.3317535545023
5、已知:⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察
其结果。
(5分)
命令如下,请运行并给出结果
>> a=[1 2 3;4 5 6;7 8 9]; >> a.^2 >> a^2 结果:ans=
1 4 9
16 25 36 49 64 81
ans =
30 36 42
66 81 96 102 126 150
分析:a.^2 只是把矩阵里面的每个元素进行平方处理 ,得到结果为第一个, 而 a^2 为 矩阵相乘,结果为第二个。
6、绘制饼图。
把图截下来,并说明都表示些什么?(10分)
程序为:
x=[190 33 45 42 45]; explode=[0 1 0 0 0]; figure(1)
subplot(2,1,1) colormap hsv pie(x,explode) gtext('生活费') gtext('资料费') gtext('电话费') gtext('衣服') gtext('其它') title('二维饼图')
subplot(2,1,2) colormap hsv pie3(x,explode) title('三维饼图') 图形:
说明:二位饼图通过每个区域的扇形面积来表 示该部分 占整体的比例,而三
维饼图则是通过 每个区域的体积比例来表 示该部分 占整体的比例 。
7、如下函数文件,可求小于任意自然数n 的Fibnacci 数列各项。
Fibnacci 数列定义如下:
⎪⎩⎪
⎨⎧>+=====--2
,2,
11,
12121n f f f n f n f n n n
请将如下函数文件保存到你的“work ”,然后取一适当的n 进行调用
运算,并记录你的结果。
(10分)
function f=fib(n) f(1)=1; if n>=2 f(2)=1; end
for i=3:n
f(i)=f(i-1)+f(i-2);
end
结果:user_entr y=input (' Please enter the number you want to know of
the Fibonacci :’);
Please enter the number you want to know of the Fibonacci: 5 >> if user_entry == l || user_entry==2
Result =l else
a= l;
b=l ;
for i=3:user_entry
Result =a+b; a=b ;
b=Result ; end
Result =Result end
Result = 5
8、求方程1
sin 02
x x -
=的根。
(至少3个,先画图后求解,要有程序与截图) (10分)
答:命令:ezplot(‘sin(x)’, ‘0.5/x ’,[-10,2*pi]) 图形:
结果:-0.1425 0.129 -0.3048
9、根据1252
2
22=-+a
y a x 绘制平面曲线,并分析参数a 对其形状的影响。
(请截图,并附分析文字)(10分)
程序如下:
a=1:4;
t=0:0.01:2*pi; for i=1:4
x=a(i)*cos(t);
y=sqrt(25-a(i)^2)*sin(t); subplot(2,2,i);plot(x,y); end
图形:
-1-0.500.51
-2-1012
-50
5
-4-2
024
分析:当a 小于a0√(12.5),椭圆随着a 的增大而渐渐地变圆 : 当a 等于a0√
( 12.5),椭圆变为圆 ; 当a 大于a0√( 12.5),椭圆随着a 的增大而渐渐地变扁.
10、求解下列各题(m 取成自己学号的后三位)。
给出命令语句与答案。
(30分。
每题6分)
1)3
0sin lim x mx mx
x ->- 命令: syms x y;limit((204*x-sin(204*x))/(x^3),x,0) 我的学号为 :2012101204 m=204 答案 :ans= 1414944
2) (10)
cos ,x y e x y =求
命令: syms x ;diff(exp(x)*cos(x),x,10) 答案:-32*exp(x)*sin(x)
3
)2
1/2
mx e dx
⎰ (精确到17位小数)
命令: clear; x=0:0.0001:0.5;
>>y=exp(204*x.^2); z=trapz(x,y)
答案:
4) 在0x =展开(最高次幂为8)
命令: syms x;taylor((204/1000+x)^(1/2),0,8) (M 取值为 204)
答案:
1/100*2370^(1/2)+5/237*2370^(1/2)*x-1250/56169*2370^(1/2)*x^2+625000/13312053*2370^(1/2)*x^3-390625000/3154956561*2370^(1/2)*x^4+273437500000/747724704957*2370^(1/2)*x^5-68359375000000/5
9070251691603*2370^(1/2)*x^6+5371093750 00000/13999649650909911*2370^(1/2)*x^7
5)
1
s i n
x
y e
=求(3)()
y m(精确到17位有效数字)
命令:clear;x=237;
>> diff(exp(sin(1/x)),x,3) (M 取值为204) 答案:Empty array: l-by-l-by-0
6)求矩阵
-211
=020
-41
100
A
m
⎛⎫
⎪
⎪
⎪
⎪
⎪
⎝⎭
的逆矩阵1-
A及特征值和特征向量。
命令:>>A=[-21 1;020;-4 1237/100];
>> format rat
>> D=inv(A)
C=eig(A)
答案:D= -237174 137/148 50/37
0 1/2 0
-200/37 50/37 100/37 C= 804/1157
804/755
2。