Matlab与科学计算样题 (加主观题答案)
Matlab与科学计算考试样题（客观题）
1 下面的MATLAB语句中正确的有：
a) 2a＝pi;
b) record_1=3+4i
c) a=2.0,
d) c=1+6j
2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t为20，30，40度时的粘度分别是：
为0℃水的黏度，值为；a、b为常数，分别为0.03368、0.000221。
（a）0.0018 0.0010 0.0007
(b) 0.0010 0.0007 0.0005 (0.0010 0.0008 0.0007)
(c) 1.7850e-003 1.0131e-003 6.6092e-004
(d) 1.0131e-003 6.6092e-004 4.6772e-004
(1.0131e-003 8.0795e-004 6.6092e-004)
a=0.03368;b=0.000221;u0=1.785e-3;
t=[20 30 40];u=u0./(1+a*t+b*t.^2)
>>format short %format short e
>>u
3.请补充语句以画出如图所示的图形：
[x,y]=meshgrid(-2:0.1:2, -2:0.1:2);
Z=x.*exp(-x.^2-y.^2);
;
a) Plot3(x,y,Z)
b) plot3(x,y,Z)
c) mesh(x,y,Z)
d) plot3(x,y,z)
2
x 0.5 1.0 1.5 2.0 2.5 3.0
y 1.75 2.45 3.81 4.80 8.00 8.60
a) 0.4900 1.2501 0.8560
b) 0.8560 1.2501 0.4900
c) -0.6341 3.8189 -3.7749
d) 3.8189 -3.7749 2.8533
解释说明：
>> x=0.5:0.5:3.0;
>> y=[1.75,2.45,3.81,4.80,8.00,8.60]; >> a=polyfit(x,y,2) a =
0.4900 1.2501 0.8560
>> x1=[0.5:0.25:3.0];
>> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on
>> plot(x1,y1,'--r')
5. 求方程在x=0.5附近的根.
a) 0.6180
b) -1.1719e-25 c) －1 d) -1.6180
6. 用Newton-Cotes 方法计算如下积分
（a ）133.6625 (b) 23.8600 (c) 87.9027 (d) -1.6180
7. y=ln(1+x)，求x=1时y"的近似值。 a) -0.25 b) 0.5 c) -0.6137 d) -1.6137
8. 某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘米。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262
计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p ＝1.3431e-005，没有显著差异。 b) p ＝0.9688，没有显著差异。 c) p ＝0.4956，有显著差异。 function f=fun(x) f=x.*x.*sqrt(2*x+3)
quadl(‘fun ’,1,5,1e-10)
或quadl('x.*x.*sqrt(2*x+3)',1,5,1e-10)
或 fun=@(x)(x.*x.*sqrt(2*x+3)); quadl(fun,1,5,1e-10)
>> X=[0.236 0.238 0.248 0.245 0.243;0.257 0.253 0.255 0.254 0.261;0.258 0.264 0.259 0.267 0.262];
>> P=anova1(X')
fzero('x.^2+x-1',0.5)
syms x y=log(1+x) f=diff(y,2) subs(f,1)
d) p＝0.9688，有显著差异。
9. 求解如下非线性方程组在（x=-1，y=-1）附近的解
a) 0.5671 0.5671
b) 无解
c) 无穷解
d) 0 0
10. 采用ODE45求解如下多阶常微分方程，并求出当x＝1.8505时的函数值。
a) 31.6441
b) 74.6907
c) 118.7862
d) 63.2564
=
11. 求解下列方程组。
a) -0.1819 -1.6630 2.2172 -0.4467
b) -0.7841 -0.0037 2.1994 -0.4226
c) -0.4467 2.2172 -1.6630 -0.1819
d) -0.4226 2.1994 -0.0037 -0.7841
12. 求极限
a）a) -1/6
function F=myfun(x)
F=[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];
x=fsolve('myfun',[-1,1])
或者
fsolve('[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))]',[-1 1])
建立求解函数文件myfun03
function dy=myfun03(x,y)
dy=zeros(3,1) %初始化变量dy，改行可以没有
dy(1)=y(2); %dy(1)表示y的一阶导数,其等于y的第二列值
dy(2)=y(3); %dy(2)表示y的二阶导数
dy(3)=2*y(3)/x^3+3*y(2)/x^3+3*exp(2*x)/x^3 %dy(3)表示y的
三阶导数
求解过程：
[x45,y45]=ode45('myfun03',[1,10],[1 10 30]);
查workspace中的矩阵框
a=[0.4096 0.1234 0.3678 0.2943;0.2246 0.3872
0.4015 0.1129;0.3645 0.1920 0.3781
0.0643;0.1784 0.4002 0.2786 0.3927];
b=[0.4043 0.1550 0.4240 -0.2557]’
x=a\b 或x=inv(a)*b
syms x
f=(x^3+x^2+x+1)^(1/3)-sqrt(x^2+x+1)*log(exp(x)+x)/x
limit(f,x,inf)
b) Inf
c) –Inf d）-1
有关上机考试说明：
（1）样题中每一题对应一组相似的题，每个人考试的时候每一组题目只会出现一道题，同组题可能会有一些细节的变化，比如说某个参数变化了或者某个积分函数发生了变化，但是所用到的基本命令是一样的。
（2）考试的时候可以启动Matlab运行以得到所需要的结果。
（3）考试采用闭卷考试，但是可以使用联机帮助。
Matlab 与科学计算考试样题（主观题）
考试要求：
1、要求独立完成不得与他人共享，答卷雷同将做不及格处理。
2、答卷用Word 文件递交，文件名为学号+姓名.doc ，试卷写上姓名及学号。
3、答卷内容包括： （1）程序；
（2）运行结果及其分析； （3）图也要粘贴在文档中。
1. 求van der Pol 方程y ''?μ (1? y 2) y '+ y = 0的数值解（μ=1），并作出y （x ）和y '（x ）的图形。（15’）
2. 在金属材料塑性变形时的流变应力σ与应变ε的近似表达式为
，对于某金属材料
function dy=vdp1(t,y);
dy=[y(2);(1-y(1)^2)*y(2)-y(1)]; 调用Matlab 函数。对于初值y (0) = 2, y
'(0) = 0 ，解为 [T,Y]=ode45('vdp1',[0 20],[2;0]); 观察结果。利用图形输出解的结果：
plot(T,Y(:,1),'-',T,Y(:,2),'--');
title('Solution of van der Pol Equation,mu=1'); xlabel('time t');
ylabel('solution y'); legend('y1','y2')
或者plot(T,Y(:,1),'-');hold on;plot(T,Y(:,2),'--'); title('Solution of van der Pol Equation,mu=1'); xlabel('time t'); ylabel('solution y')
应力σ925 1125 1625 2125 2625 3125 3625
应变ε0.11 0.16 0.35 0.48 0.61 0.71 0.85 计算参数k和n，并分别画出实验测试数据点和拟合曲线（15’）。
近似表达式可以写成
y=[925,1125,1625,2125,2625,3125,3625];
x=[0.11,0.16,0.35,0.48,0.61,0.71,0.85];
x1=log(x);
y1=log(y);
p=polyfit(x1,y1,1);
n=p(1),k=exp(p(2))
xi=linspace(0.1,0.9,800);xi=0.1:0.001:0.9;
yi=exp(polyval(p,log(xi)));yi=k*xi.^n;
plot(x,y,'o',xi,yi)
xlabel('\epsilon')
ylabel('\sigma')
legend('experimental','Fitting')
3. 在4个子图中绘制不同的三角函数图（10’）。
函数范围：x=0:0.1*pi:2*pi;
函数为：sin(x); cos(x)；sin(x)+cos(x)；sin(x).*cos(x)
>> x=0:0.1*pi:2*pi;
>> subplot(2,2,1);
>> plot(x,sin(x),'-*');
>> title('sin(x)');
>> subplot(2,2,2);
>> plot(x,cos(x),'-*');
>> title('cos(x)');
>> subplot(2,2,3);
>> plot(x,sin(x)+cos(x),'-*'); >> title('sin(x)+cos(x)');
>> subplot(2,2,4);
>> plot(x,sin(x).*cos(x),'-*'); >> title('sin(x)*cos(x)');
Matlab上机实验答案
Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -
>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +
Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50
+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =
MATLAB试题解答
2．一单位反馈控制系统，若要求： （1）跟踪单位斜坡输入时系统的稳态误差为2。 （2）设该系统为三阶，其中一对复数闭环极点为-1±j1 求满足上述要求的开环传递函数。 利用MATALAB 完成以下操作： （1）求控制系统的单位阶跃响应、单位斜坡响应、单位加速度响应； （2）绘制控制系统的根轨迹，并判断系统的稳定性； （3）绘制控制系统的奈奎斯特图和Bode 图。 解：由于控制系统为单位反馈系统 则 H （s ）=1 由条件一可知，输入信号R(s)= 21s 时，稳态误差为： 2) (10)()(1)(==→+=s sG s s H s G s sR Ess 由条件二可知，系统为三阶且有。两个极点s1=-1+j1 ；s2=-1-j1 则可得，开环传递函数为 ) 222(1)(++=s s s s G （1）求控制系统的单位阶跃响应： G=tf([1],[1, 2, 2, 0]); Step(G) title('单位阶跃响应'); grid
(2)控制系统的单位斜坡响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1},0]; sys1=tf(num,den) step(sys1)
(3)控制系统的单位加速度响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1} 0 0]; sys2=tf(num,den) step(sys2)
MATLAB全部实验及答案
MATLAB全部实验及答案 实验一、MATLAB基本操作 实验内容及步骤 4、有关向量、矩阵或数组的一些运算 （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？ A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？ （4）设a=[1 -2 3;4 5 -4;5 -6 7] 请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全 下标的形式），并将其单下标转换成全下标。 clear,clc a=[1 -2 3;4 5 -4;5 -6 7]; [x,y]=find(a<0); c=[]; for i=1:length(x) c(i,1)=a(x(i),y(i)); c(i,2)=x(i); c(i,3)=y(i); c(i,4)=(y(i)-1)*size(a,2)+x(i); end c
（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行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^2= 22 16 16 25 26 23 26 24 28 a.^2= 1 4 9 9 16 4 25 4 9 （7）有一段指令如下，请思考并说明运行结果及其原因 clear X=[1 2;8 9;3 6]; X( : ) 转化为列向量 （8）使用三元组方法，创建下列稀疏矩阵 2 0 8 0 0 0 0 1 0 4 0 0 6 0 0 0 方法一： clear,clc
matlab练习题及答案
第1章MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 答：MATLAB语言有一下特点： 1、起点高。 2、人机界面适合科技人员。 3、强大而简易的作图功能。 4、智能化程度高。 5、功能丰富，可扩展性强。 1.2 MATLAB系统由那些部分组成？ 答：MATLAB系统主要由开发环境、MATLAB数学函数库、MMATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。
1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器。 （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到
MATLAB精通科学计算_偏微分方程求解
一、Maple V 系统 Maple V是由Waterloo大学开发的数学系统软件，它不但具有精确的数值处理功能，而且具有无以伦比的符号计算功能。Maple V的符号计算能力还是MathCAD和MATLAB等软件的符号处理的核心。Maple提供了2000余种数学函数，涉及范围包括：普通数学、高等数学、线性代数、数论、离散数学、图形学。它还提供了一套内置的编程语言，用户可以开发自己的应用程序，而且Maple自身的2000多种函数，基本上是用此语言开发的。 Maple采用字符行输入方式，输入时需要按照规定的格式输入，虽然与一般常见的数学格式不同，但灵活方便，也很容易理解。输出则可以选择字符方式和图形方式，产生的图形结果可以很方便地剪贴到Windows应用程序内。 二、MATLAB 系统 MATLAB原是矩阵实验室（Matrix Laboratory）在70年代用来提供Linpack和Eispac k软件包的接口程序，采用C语言编写。从80年代出现3.0的DOS版本，逐渐成为科技计算、视图交互系统和程序语言。MATLAB可以运行在十几个操作平台上，比较常见的有基于W indows 9X/NT、OS/2、Macintosh、Sun、Unix、Linux等平台的系统。 MATLAB程序主要由主程序和各种工具包组成，其中主程序包含数百个内部核心函数，工具包则包括复杂系统仿真、信号处理工具包、系统识别工具包、优化工具包、神经网络工具包、控制系统工具包、μ分析和综合工具包、样条工具包、符号数学工具包、图像处理工具包、统计工具包等。而且5.x版本还包含一套几十个的PDF文件，从MATLAB的使用入门到其他专题应用均有详细的介绍。 MATLAB是数值计算的先锋，它以矩阵作为基本数据单位，在应用线性代数、数理统计、自动控制、数字信号处理、动态系统仿真方面已经成为首选工具，同时也是科研工作人员和大学生、研究生进行科学研究的得力工具。MATLAB在输入方面也很方便，可以使用内部的E ditor或者其他任何字符处理器，同时它还可以与Word6.0/7.0结合在一起，在Word的页面里直接调用MATLAB的大部分功能，使Word具有特殊的计算能力。 三、MathCAD 系统 MathCAD是美国Mathsoft公司推出的一个交互式的数学系统软件。从早期的DOS下的1. 0和Windows下的4.0版本，到今日的8.0版本，功能也从简单的数值计算，直至引用Map le强大的符号计算能力，使得它发生了一个质的飞跃。 MathCAD是集文本编辑、数学计算、程序编辑和仿真于一体的软件。MathCAD7.0 Profe ssional（专业版）运行在Win9X/NT下，它的主要特点是输入格式与人们习惯的数学书写格式很近似，采用WYSWYG（所见所得）界面，特别适合一般无须进行复杂编程或要求比较特殊的计算。MathCAD 7.0 Professional 还带有一个程序编辑器，对于一般比较短小，或者要求计算速度比较低时，采用它也是可以的。这个程序编辑器的优点是语法特别简单。 MathCAD可以看作是一个功能强大的计算器，没有很复杂的规则；同时它也可以和Wor d、Lotus、WPS2000等字处理软件很好地配合使用，可以把它当作一个出色的全屏幕数学公式编辑器。 四、Mathematica 系统 Mathematica是由美国物理学家Stephen Wolfram领导的Wolfram Research开发的数学系统软件。它拥有强大的数值计算和符号计算能力，在这一方面与Maple类似，但它的符
MATLAB)课后实验答案
实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012 2sin 851z e =+ (2) 21ln(2 z x =，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022 a a e e a z a a --+= ++=--L (4) 22 42011 122123t t z t t t t t ?≤
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。解：(1) 结果： (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是：
实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22 E R RS A O S +?? =???? 。 解: M 文件如下； 5. 下面是一个线性方程组： 1 231 1 12340.951110.673450.5211145 6x x x ?? ??????????????=??? ?????????????????? ? ch = 123d4e56g9
(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 值。 解：M 文件如下：
matlab考试题及答案
%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (i
科学计算与MATLAB语言(第四课)
第四讲绘图功能
作为一个功能强大的工具软件，Matlab 具有很强的图形处理功能，提供了大量的二维、三维图形函数。由于系统采用面向对象的技术和丰富的矩阵运算，所以在图形处理方面即常方便又高效。
4.1 二维图形 一、plot函数 函数格式：plot(x,y)其中x和y为坐标向量函数功能：以向量x、y为轴，绘制曲线。【例1】在区间0≤X≤2 内，绘制正弦曲线Y=SIN（X），其程序为： x=0:pi/100:2*pi; y=sin(x); plot(x,y)
一、plot函数 【例2】同时绘制正、余弦两条曲线Y1=SIN（X）和Y2=COS（X），其程序为： x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,x,y2) plot函数还可以为plot(x,y1,x,y2，x,y3，…)形式，其功能是以公共向量x为X轴，分别以y1，y2，y3，…为Y轴，在同一幅图内绘制出多条曲线。
一、plot函数 （一）线型与颜色 格式：plot(x,y1,’cs’,...) 其中c表示颜色，s表示线型。 【例3】用不同线型和颜色重新绘制例4.2图形，其程序为：x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,'go',x,y2,'b-.') 其中参数'go'和'b-.'表示图形的颜色和线型。g表示绿色，o表示图形线型为圆圈；b表示蓝色，-.表示图形线型为点划线。
一、plot函数 （二）图形标记 在绘制图形的同时，可以对图形加上一些说明，如图形名称、图形某一部分的含义、坐标说明等，将这些操作称为添加图形标记。 title(‘加图形标题'); xlabel('加X轴标记'); ylabel('加Y轴标记'); text(X,Y,'添加文本');
MATLAB实验题答案
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 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 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 -3 4 1 5 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 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x 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) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w 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)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A
matlab所有实验及答案
实验二 习题 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。取出子矩阵?? ? ? ??21301，并求该矩阵的维数。 解：命令为： Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) 2、 建立一个数值范围为0—100内4*5的整数随机矩阵。 求出大于50的元素的位置。 解：命令为： G=int8(100*rand(4,5)) find(G>50) 3、 已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B 解：命令为： A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B 4、 利用函数产生3*4阶单位矩阵和全部元素都为8的 4*4阶矩阵，并计算两者的乘积。 解：命令为： A=eye(3,4) B=8*ones(4) C=A*B 5、 创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c 2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来， 赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的 形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。 （1）)/1sin(t ，其中)1,1(-∈t ； （2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三 视图 7. 试求出如下极限。 （1）x x x x 1)93(lim +∞ →； （2）1 1lim 0-+→→xy xy y x ； （3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。 （1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ； （2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。 （1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。 《Matlab与科学计算》作业 第一章MATLAB环境 1、MATLAB通用操作界面窗口包括哪些？命令窗口、历史命令窗口、当前目录窗口、工作空间窗口各有哪些功能？ 答：MATLAB通用操作界面窗口包括：命令窗口、历史命令窗口、当前目录浏览器窗口、工作空间窗口、变量编辑器窗口、M文件编辑/调试器窗口、程序性能剖析窗口、MATLAB帮助。 命令窗口是MATLAB命令操作的最主要窗口，可以把命令窗口当做高级的“草稿纸”。在命令窗口中可以输入各种MATLAB的命令、函数和表达式，并显示除图形外的所有运算结果。 历史命令窗口用来记录并显示已经运行过的命令、函数和表达式，并允许用户对它们进行选择、复制和重运行，用户可以方便地输入和修改命令，选择多行命令以产生M文件。 当前目录窗口用来设置当前目录，可以随时显示当前目录下的M、MKL等文件的信息，扬文件类型、文件名、最后个修改时间和文件的说明信息等，并可以复制、编辑和运行M文件及装载MAT数据文件。 工作空间窗口用来显示所有MATLAB工作空间中的变量名、数据结构、类型、大小和字节数。 2、熟悉课本中表格1.4、1.5、1.6、1.7、1.8的内容。 3、如何生成数据文件？如何把数据文件中的相关内容输入到工作空间中，用实例进行操作。 生成数据文件： 把数据文件中的相关内容输入到工作空间中： 结果： 4、在工作空间中可以通过哪些命令管理变量，写出每种语法的具体操作过程。答：（1）把工作空间中的数据存放到MAT数据文件。 语法：save filename 变量1 变量2 ……参数。 （2）从数据文件中取出变量存放到工作空间。 语法：load filename 变量1 变量2 ……。 实验四GUI 设计专业 学号姓名成绩电气工程及其自 动化201409140305 杨诚1.创建GUI 绘制方程c bx ax y ++=2图形，需要显示绘图结果的坐标系窗口，还能够输入 a ， b ， c 的值和x 取值范围（最大值和最小值）。 回调函数的编写： a=str2num(get(handles.edit1,'String'));b=str2num(get(handles.edit2,'String'));c=str2num(get(handles.edit3,'String'));xmin=str2num(get(handles.edit4,'String'));xmax=str2num(get(handles.edit5,'String'));x=xmin:0.1:xmax;y=a*x.^2+b*x+c;plot(x,y); 设计的运行界面截图： 2.设计一个GUI，完成画出y=sin(x)、y=cos(x)和y=x 3.的波形图。回调函数的编写： x=-5:0.1:5 plot(x,sin(x)) x=-5:0.1:5 plot(x,cos(x)) x=-5:0.1:5 y=x.^3 plot(x,y) 设计的运行界面截图： 3.创建一个GUI，含有下拉菜单，下拉菜单中有背景颜色选择。回调函数的编写： yanse=get(handles.popupmenu1,'value'); switch yanse case1 set(gcf,'color','r'); case2 set(gcf,'color','y'); case3 set(gcf,'color','g'); end 设计的运行界面截图： MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？ 进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？ 在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？ 在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。 7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。 相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？ 在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。 盐城师范学院《MATLAB与科学计算》期末论文 2016－2017学年度第一学期 用MATLAB解决解析几何的图形问题 学生姓名吴梦成 学院数学与统计学院 专业信息与计算科学 班级数15(5)信计 学号 15213542 用MATLAB 解决解析几何的图形问题 摘 要 将 MATLAB 的图形和动画功能都用于解析几何教学，可使教学形象生动。以图形问题为例，详细给出了实例的程序编写和动画实现过程 。在解析几何教学中有一定的应用价值。 【关键词】： MATLAB ； 解析几何 ；图形 ； 动 画；编程 1 引 言 在解析几何的教学中，使用传统的教学方法。许多曲线及曲面的形成过程与变换过程只通过传统的教师讲授静态图示就很难形象生动地表示出来 。在解析几何教学中使用MATLAB 软件辅助教学，不仅可以很容易绘制出复杂的立体图形。把曲线、曲面的形成和变化过程准确地模拟出来 ，而且还能够对它们进行翻转 、旋转 ，甚 至还能够轻而易举地实现图形的动画效果!这对提高教学效率和培养学生的空间想象能力可起到事半功倍的效果。下面结合实例从几个方面说明MATLAB 在解析几何画图方面的应用。 2 利用 MATLAB 绘制三维曲线 在空间解析几何中，各种曲线和曲面方程的建立都离不开图形 ，而空间曲线和曲面图形既难画又费时。借助MATLAB 的绘图功能 ，可以快捷 、 准确地绘出图形，使教学变得形象 、生动 。有利于学生观察三维空间图形的形状 ， 掌握图形的性质 。 一 般地 ，MATLAB 可用plot3，ezplot3，comet3等函数来各种三维曲线 。 例如画螺旋曲线的图形，其参数方程设为 ：t at cos x =，t b sin t y -=，ct =z 。使用 plot3语句画螺旋曲线图形的方法如下( 设a =2 ,b=4,c=3)： );*3),sin(*.*4),cos(*.*2(3;*10:50/:0t t t t t plot pi pi t -= MATLAB 用两条简单的语句就可以画出螺旋 曲线(图1)，但上述方法是静态的 ，为了体 Matlab 与科学计算考试样题（客观题） 1 下面的MATLAB 语句中正确的有： a) 2a ＝pi 。 b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 2 1at bt μμ=++0μ为0℃水的黏度，值为31.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2)。 Z=x.*exp(-x.^2-y.^2)。 。 a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明： >> x=0.5:0.5:3.0。 >> y=[1.75,2.45,3.81,4.80,8.00,8.60]。 >> a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]。 >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在 x=0.5附近的根. 21 x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes方法计算如下积分 1 5 x? （a）133.6625 (b)23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8.某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘M。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p＝1.3431e-005，没有显著差异。 实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0 matlab练习题和答案 控制系统仿真实验 Matlab 部分实验结果 目录 实验一MATLAB基本操作..................................................................... .......................1 实验二Matlab 编程..................................................................... ...............................5 实验三Matlab 底层图形控制..................................................................... .....................6 实验四控制系统古典分析...................................................................... .......................12 实验五控制系统现代分析..................................................................... . (15) 实验六PID 控制器的设计..................................................................... ......................19 实验七系统状态空间设计...................................................................... .......................23 实验九直流双闭环调速系统仿真..................................................................... . (25) 实验一MATLAB基本操作 实验一：Matlab操作环境熟悉 一、实验目的 1．初步了解Matlab操作环境。 2．学习使用图形函数计算器命令funtool及其环境。 二、实验内容 熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算： 1．单函数运算操作。 求下列函数的符号导数 (1) y=sin(x); (2) y=(1+x)^3*(2-x); 求下列函数的符号积分 (1) y=cos(x); (2) y=1/(1+x^2); (3) y=1/sqrt(1-x^2); (4) y=(x-1)/(x+1)/(x+2); 求反函数 (1) y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2)); 代数式的化简 (1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4); (2) sin(x)^2+cos(x)^2; (3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x); 2．函数与参数的运算操作。 从y=x^2通过参数的选择去观察下列函数的图形变化 (1) y1=(x+1)^2 (2) y2=(x+2)^2 (3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6) y6=x^2/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x)matlab习题及答案
《Matlab与科学计算》作业 2010010099
matlab实验四及其答案
matlab习题及答案2
《MATLAB与科学计算》期末论文
Matlab与科学计算样题(加主观题答案)
MATLAB实验题目及答案
matlab练习题和答案
Matlab实验指导书(含答案)汇总