Matlab上机_测试题c答案讲课教案
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,采用不同插值方法进行插值,并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
数学建模2021c题解析用matlab
《数学建模2021C题解析用Matlab》一、引言数学建模是一门研究怎样应用数学知识和方法来解决实际问题的学科。
而在数学建模的实际应用中,Matlab是一个常用的数学建模工具。
本文将以2021年C题为例,介绍用Matlab进行数学建模的方法和步骤。
二、题目分析2021年C题的题目是关于某体育场馆的冷却系统优化问题。
通过分析题目,我们可以了解到需要解决以下几个问题:1. 如何建立冷却系统的数学模型?2. 如何优化冷却系统的参数以提高效率?3. 如何利用Matlab进行模拟实验和数据分析?三、建立数学模型在建立数学模型时,我们需要考虑以下因素:1. 建立冷却系统的热传导方程和流体力学方程;2. 考虑不同参数对于冷却系统的影响;3. 建立合适的边界条件和初始条件。
在Matlab中,我们可以通过编写相应的程序来建立数学模型,并进行模拟实验。
我们可以利用Matlab来解决热传导方程和流体力学方程,得到冷却系统的温度分布和流速分布。
我们可以通过改变不同参数,比如冷却系统中的换热器面积、流体的流速等,来观察参数变化对系统性能的影响。
四、优化冷却系统在优化冷却系统时,我们可以利用Matlab来进行参数优化。
通过设置合适的优化目标和约束条件,可以通过Matlab内置的优化函数来优化冷却系统的参数。
我们可以通过最小化能耗或最大化换热效率来优化冷却系统的参数。
在优化过程中,我们还可以利用Matlab来进行灵敏度分析,以了解不同参数对于系统性能的影响程度。
这将有助于我们更好地理解冷却系统的特性,并为优化提供更多的参考信息。
五、个人观点和理解通过上述分析和讨论,我认为Matlab作为数学建模的工具,具有很高的灵活性和可扩展性。
它不仅可以帮助我们建立复杂的数学模型,还可以进行模拟实验、数据分析和参数优化。
我相信在数学建模的实际应用中,Matlab将会发挥越来越重要的作用。
六、总结通过以上分析,我们可以清晰地了解了如何利用Matlab进行数学建模,尤其是在解决冷却系统优化问题时的具体方法和步骤。
大学matlab考试题及答案
大学matlab考试题及答案大学MATLAB考试题及答案一、选择题1. MATLAB的全称是什么?A. Matrix LaboratoryB. Microprocessor Application Tool SetC. Microsoft Advanced Tool SetD. Microprocessor Application Test System答案:A2. 在MATLAB中,以下哪个命令用于绘制三维图形?A. plotB. scatterC. surfD. bar答案:C3. MATLAB中用于求解线性方程组的函数是?A. solveB. linsolveC. linprogD. fsolve答案:A二、简答题1. 简述MATLAB的基本数据类型有哪些,并给出至少两个每种类型的示例。
答案:MATLAB的基本数据类型包括数值数组、字符数组和单元数组。
数值数组可以是向量、矩阵或多维数组。
例如,向量 `v = [1 2 3]`,矩阵 `M = [1 2; 3 4]`。
字符数组是由单引号或双引号括起来的字符序列,如 `C = 'Hello'`。
单元数组可以包含不同类型的数据,如`{1, 'text', [1; 2; 3]}`。
2. 描述如何在MATLAB中实现矩阵的转置和翻转。
答案:矩阵的转置可以通过 `'T'` 来实现,例如 `A'` 表示矩阵A 的转置。
矩阵的翻转可以通过 `flip` 函数实现,例如 `flip(A)` 可以翻转矩阵A的所有行和列,`flipud(A)` 仅翻转矩阵A的行,而`fliplr(A)` 仅翻转矩阵A的列。
三、编程题1. 编写一个MATLAB函数,该函数接受一个向量作为输入,并返回向量中所有元素的和以及平均值。
```matlabfunction [sumVal, avgVal] = calculateSumAndAverage(V)sumVal = sum(V);avgVal = mean(V);end```2. 设计一个MATLAB脚本来解决以下问题:给定一个3x3的矩阵,找出其中最大的元素,并显示其位置。
MATLAB 上机 习题及答案
15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5,要求先求得P(x)=P1(x)+P2(x),然后计算xi=0.2*i各点上的P(xi)(i=0,1,2,…,5)值。
p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵:A=[1 2]1234B⎡⎤=⎢⎥⎣⎦,求出矩阵A和B的乘积,并将结果赋给变量C。
>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息,了解inv命令的调用格式,并作简要说明。
help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法,然后利用plot命令绘制函数y=sin(x)的图形,其中0xπ≤≤。
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上机练习(附答案)
1.以下两种说法对吗?(1)MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
(2)MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同?DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上?“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗?4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组?用直接键入法?用数组编辑器?用M文件编辑器?5.运用数组算术运算符去掉下面程序里的for/end循环:x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符,重写下面的程序代码:x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7.某公司销售电脑打印机的价格方案如下:()如果顾客只买一台打印机,则一台的基本价格为$150。
MATLAB期末上机试题带答案
MATLAB 期末上机考试试题带答案版姓名:学号:成绩:1.请实现下图:50100150200250x y 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=[321;1-13;24-4];b=[7;6;-2];x=a\b4.请实现下图:x y x=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 (-x 2-y 2)绘制网格图[x,y]=meshgrid(-2:0.1:2);z=x.*exp (-x.^2-y.^2);mesh(x,y,z)6.请实现peaks 函数:-33x Peaksy [x,y]=meshgrid(-3:1/8:3);z=peaks(x,y);mesh(x,y,z)surf(x,y,z)shading flataxis([-33-33-88])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。
Matlab上机实验题及参考解答
Matlab上机实验题及参考解答目录实验一Matlab初步实验 (2)一matlab基本功能介绍 (2)二Matlab扩展功能 (2)三练习 (2)四练习题参考解答 (3)实验二概率模型实验 (5)一复习 (5)二事件的响应 (5)三Matlab中随机数字的生成与处理 (5)四练习 (5)五练习题参考解答 (5)实验三插值与拟合 (7)实验四线性规划与非线性规划 (8)4.1 实验目的 (8)4.2 实验内容 (9)4.3 综合练习 (10)4.4 课外作业 (11)实验五数值计算 (12)5.1 实验目的 (12)5.2 实验内容 (12)4.3 综合练习 (15)4.4 课外作业 (15)实验六计算机图像处理 (16)6.1 实验目的 (16)6.2 实验内容 (16)6.3 综合练习 (17)6.4 课外作业 (19)实验七综合练习 (19)7.1 实验目的 (19)7.2 实验内容 (19)7.3 综合练习 (20)7.4 课外作业 (21)实验一 Matlab 初步实验 一 matlab 基本功能介绍1 编程环境2语法规范:for … end; if …else if …end; 3 矩阵运算 4 图形绘制二 Matlab 扩展功能1 编程练习:(1) 绘出序列kk x x r r 0(1),0.2083=+=;(2) 绘出曲线rtx t x e t 0(),0=>2 扩展功能(1) 矩阵中全部数据、部分数据的截取、更改; (2) 矩阵的初始化与赋值如:A=zeros(5,5); A(2:2:)=[1,2 3 4 5] 3 微积分基础(见实验4) 符号计算三 练习(课上编程完成下列练习,课后上机验证) 1 求和S=1+2+3+…+100; 2 求和e 1111!2!10!1...=++++3求和S 1112310!1...=++++4设A 234576138⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦, 求A 的逆、特征值和特征向量;验证Ax=λx 5 画函数图()011mrtm x x t x e x -=⎛⎫+- ⎪⎝⎭6 展开 (x-1)(x-2)…(x-100)7 因式分解 x 8—y 8; 因数分解200520068 求极限312lim +∞→⎪⎭⎫⎝⎛++n n n n9 )](sin[cos 22x x y += 求dxdy10 求积分x xdx 10ln ⎰11 求积分3⎰并且画出所求的平面区域12 设x+2y=1, 2x+3y=6, y=2x 2, 画出各个方程图形,求出曲线交点.四 练习题参考解答%MatlabTrain1.m clear all % 2nd e=1; temp=1; for I=1:1:10temp=temp*I; e=e+1/temp; end e%%%%%%%%%%% clear all % 3nd S=0; temp=1;for I=1:1:100temp=temp*I; endfor J=1:1:temp S=S+1/J; end S%%%%%%%%%%%%%% clear all % 11ndx=linspace(0,4); y=1./sqrt(x.^5+1); plot(x,y) for t=1:0.1:3yt=1./sqrt(t.^5+1);hold online([t,t],[0,yt]);end%fill(t,yt,'b') %%%%%%%%%%%%% clear all% 12ndx=linspace(-2,2);y=[0.5-0.5*x; 2-2/3.*x; 2*x.^2]; plot(x,y)grid实验二概率模型实验一复习1 小结上次编程练习中存在的问题,讲述部分习题答案2 画图命令介绍:line二事件的响应(1) 获取鼠标的位置%MatlabTrain2.mclear all% 鼠标响应p=ginput(3)plot(p(:,1),p(:,2),'r*')(2) 键盘输入相应t=input('How many apples? t=');m=t+3三Matlab中随机数字的生成与处理1 随机数的生成2 产生随机数字3 产生某区间的整数4 生日模拟问题的Montecaro法设计技术、思路学生尝试编程四练习(1) 编程验证人数在不同年龄段的生日的概率计算(2) 编程实现游戏”聪明伶俐100分”(3) 编程实现两家电影院的座位数问题(4) 编程实现某图形面积的计算五练习题参考解答(1) 生日问题程序示例:%birthPro.mn=0;nStudents=30;for I=1:1000 %how many times testy=0;x=1+floor(365*rand(1,nStudents));%get nStudents random numbersfor J=1:nStudents-1for K=J+1:nStudentsif x(J)==x(K)y=1;break;endendendn=n+y;%count, n times of that there are two people's dirthday in the same dayendfreq=n/I % caculating the frequently(2) 编程实现游戏”聪明伶俐100分”参考答案%MatlabTrain2.mclear all% 鼠标响应x=floor(10*rand(1,4))t=input('填入四个数字[n1 n2 n3 n4]=');flag=0;A=0;B=0;for I=1:1:8flag=flag+1;A=0;B=0;if t==xswitch flagcase 1disp('聪明绝顶!');case 2disp('聪明!');case 3disp('有点聪明!');case 4disp('还可以!');case 5disp('聪明伶俐100分!');case 6disp('聪明伶俐90分!');case 7disp('聪明伶俐85分!');case 8disp('聪明伶俐80分!');otherwisedisp('赫赫!');endbreak;endfor J=1:1:4for K=1:1:4if x(J)==t(K) & J==KA=A+1;else if x(J)==t(K) & J~=KB=B+1;endendendends='AABB';s(1)=INT2STR(A);s(3)=INT2STR(B);disp(s);t=input('不重复填入四个数字[n1 n2 n3 n4]=');endif flag>0disp('太烂了! 正确答案是:');xend实验三插值与拟合一复习讲述聪明伶俐100分的编程中的问题二插值三拟合课堂练习2 某之股票价格from 2003 09 01 to 2004 01 02,试进行插值、拟合%TimerS.m%from 2003 09 01 to 2003 01 02clear all;dataST=[15.09 14.7514.95 14.722.88 21.8619.82 19.09];plot(dataST)四课外练习112)进行多项式拟合,求出拟合多项式,并求出多项式在t=4, 5处的值.实验四线性规划与非线性规划4.1 实验目的1 用Matlab求解线性规划2 用Matlab求解非线性规划4.2 实验内容4.2.1 线性规划求解实用格式:x=lp(c, A, b, xLB,xUB,x0,nEq)可以求解下列线性规划模型:min f=c’xs.t. Ax=<=b(其中前nEq个约束为等式约束,即等式约束的个数,其余是不等式约束<=) xLB<=x<=xUB函数中x0参数是算法迭代的初始点,任意取值例1 求解下列线性规划1)123123123123min2..360210200,1,2,3jz x x xs t x x xx x xx x xx j=--+⎧⎪++≤⎪⎪-+≤⎨⎪+-≤⎪≥=⎪⎩,2)1235635623416367min..3621060,1,,7jz x x x x xs t x x xx x xx xx x xx j=-++-⎧⎪++=⎪⎪+-=⎪⎨-+=⎪⎪++=⎪≥=⎪⎩例1求解示例c=[-2 -1 1]';%book page 72 Number 16-1A=[3 1 1;1 -1 2;1 1 -1];b=[60 10 20]';xlb=[0 0 0]';xub=[inf inf inf]';x0=[0 0 0]'; x=lp(c,A,b,xlb,xub,x0,0)% x=(15 5 0)'例2 求解示例c2=[1 -1 1 0 1 -1 0]';%book page 72 Number 16-3A2=[0 0 3 0 1 1 0;...0 1 2 -1 0 0 0;...-1 0 0 0 0 1 0;...0 0 1 0 0 1 1];b2=[6 10 0 6]';xlb2=[0 0 0 0 0 0 0]';xub2=[inf inf inf inf inf inf inf]';x02=[0 0 0 0 0 0 0]';x2=lp(c2,A2,b2,xlb2,xub2,x02,4)% unbounded4.2.2 非线性规划1)命令格式1:[X, OPTIONS]=constr(‘FUN’, X, OPTIONS,VLB,VUB)2)命令格式2:X=FMINCON(FUN,X0,A,B,Aeq,Beq)% minimizes FUN subject to the linear equalities% Aeq*X = Beq as well as A*X <= B. (Set A=[] and B=[] if no inequalities exist.)例2 求解非线性规划y x x x x s t x3211221min22 ..1=++-≤-求解示例%unconop.mfunction y=unconop(x)y=x(1).^3+2*x(1).*x(2)+2*x(2).^2;%book page 148 ex.7-1 后建立调用函数xx=fmincon('unconop',[0 0]',[-1 0],-1,[],[])%book page 148 ex.7-1 4.3 综合练习学生独立编写程序,求解一个含有2个变量的线性规划问题,要求:1)编写程序,把可行域画上阴影;2)求出最优解,在可行域上标出最优解;3)求出基本解,并在上图中表示出来;4)求出基本可行解,观察单纯形方法迭代时,顶点的变化.可行域画图与表出阴影示例:syms x y[u(1),v(1)]=solve('y=x+2','y=2*x');%求出交点坐标[u(2),v(2)]=solve('y=-x+2','y=2*x');[u(3),v(3)]=solve('y=x+2','y=-x+2');x=linspace(0,3,5); %直线作图y=[2*x;-x+2;x+2];line(x,y); gridpatch(double(u),double(v),'b'); 运行结果:4.4 课外作业1 求解线性规划131223min ..250.530,1,2,3i x x s t x x x x x i +⎧⎪+≤⎪⎨+=⎪⎪≥=⎩ (1) 求解线性规划;x *=()(2) 目标函数中c 1由1变为(-1.25)时求最优解;(3) 目标函数中c 1由1变为(-1.25),c 3由1变为2时求最优解;(4) 约束条件中53b ⎛⎫= ⎪⎝⎭变为21b -⎛⎫'= ⎪⎝⎭时,求解;(5) 约束条件中53b ⎛⎫= ⎪⎝⎭变为23b ⎛⎫'= ⎪⎝⎭时,求解[刁在筠,运筹学(第二版),高等教育出版社,2004,01 p74第20题]2 求解非线性规划y x x x x x x x 3221122233min 2223=++++ 注:无约束非线性规划问题, 命令:fminunc子函数% unconop.mfunction y=unconop(x)y=x(1).^2+2*x(1).*x(2)+2*x(2).^2+2*x(2).*x(3)+3*x(3).^2;%book page 148 ex.7-1 主函数:xx=fminunc('unconop',[0.1 0.1 1]')思考:绘出两个变量的线性规划问题的可行域、标出可行的整数解和求出可行解;演示单纯形方法的迭代过程,如j z x x s t x x x x x j 121212min 2..360200,1,2=--⎧⎪+≤⎪⎪+≤⎨⎪⎪≥=⎪⎩实验五 数值计算5.1 实验目的1 掌握代数数值计算2 掌握常微分方程数值计算5.2 实验内容5.2.1 关于多项式设多项式1110()n n n n p x a x a x a x a --=++++表示为110[,,,,]n n p a a a a -=1)求多项式的根 roots(p) %求出p(x)=0的解。
matlab上机试验[教学]
Matlab 上机实验一、 实验目的1、 掌握绘制MATLAB 二维、三维和特殊图形的常用函数;2、 熟悉并掌握图像输入、输出及其常用处理的函数。
二、 实验内容1 绘制函数的网格图和等高线图。
422cos cos y x yex z +-=其中x 的21个值均匀分布在[-5,5]范围,y 的31个值均匀分布在[0,10],要求将产生的网格图和等高线图画在同一个图形窗口上。
2 绘制三维曲面图,使用纯铜色调色图阵进行着色,并进行插值着色处理。
⎪⎩⎪⎨⎧===s z t s y ts x sin sin cos cos cos230,20ππ≤≤≤≤t s3 已知⎪⎪⎩⎪⎪⎨⎧>++≤+=0),1ln(210,22x x x x e x y π在-5<=x<=5区间绘制函数曲线。
4 已知y1=x2,y2=cos(2x),y3=y1*y2,其中x 为取值-2π~2π的等差数列(每次增加0.02π),完成下列操作:a) 在同一坐标系下用不同的颜色和线型绘制三条曲线,给三条曲线添加图例;b) 以子图形式,分别用条形图、阶梯图、杆图绘制三条曲线,并分别给三个图形添加标题“y1=x^2”,“y2=cos(2x)”和“y3=y1*y2”。
5 在xy 平面内选择区域[][],,-⨯-8888,绘制函数z =的三种三维曲面图。
6 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x 轴,y 轴,标题为“正弦余弦函数图象”.7 分别用线框图和曲面图表现函数z=cos(x)sin(y)/y ,其中x 的取值为[-1.5pi,1.5pi],y=x ,要求:要有标题、坐标轴标签8 有一组测量数据满足-ate =y ,t 的变化范围为0~10,用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,并加入标题和图列框(用代码形式生成)9 22y x xez --=,当x 和y 的取值范围均为-2到2时,用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图10 x= [66 49 71 56 38],绘制饼图,并将第五个切块分离出来。
Matlab上机题库及详细答案_Tonyxie
ax 2 + bx + c 0.5 ≤ x < 1.5 = y a sin c b + x 1.5 ≤ x < 3.5 c ln b + 3.5 ≤ x < 5.5 x
clc;clear; a=input('Please input a= '); b=input('Please input b= '); c=input('Please input c= '); x=input('Please input x= '); disp('#if#'); if x>=0.5&x<1.5 y=a*x^2+b*x+c elseif x<3.5&x>=1.5 y=a*sin(b)^c+x elseif x>=3.5&x<5.5 y=log(abs(b+(c./x))) else y='ERROR!' end disp('#switch#') switch(x*10/5) case {1 2} y=a*x^2+b*x+c case {3 4 5 6} y=a*sin(b)^c+x case {7 8 9 10} y=log(abs(b+(c./x))) otherwise y='ERROR!' end
பைடு நூலகம்
8、数值与符号计算
ex (1)求极限 lim x →+∞ a + be x
(3)已知线性方程组 Ax=b,其中
(2)求不定积分
∫ xe
ax
Matlab上机题库及详细答案
解:c=input('请输入一个字符','s');
if c>='A' & c<='Z'
disp(setstr(abs(c)+abs('a')-abs('A')));
elseif c>='a'& c<='z'
disp(setstr(abs(c)- abs('a')+abs('A')));
(2)均值为0.6方差为0.1的5阶正态分布随机矩阵
解:>>x=20+(50-20)*rand(5);
>>y=0.6+sqrt(0.1)*randn(5)
例13:将101~125等25个数填入一个5行5列的表格中,使其每行每列及对角线的和均为565。
解:M=100+magic(5)
M =117 124 101 108 115
-0.4606
0.3848
例2、用简短命令计算并绘制在0x6范围内的sin(2x)、sinx2、sin2x。
解:x=linspace(0,6)
y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2;
plot(x,y1,x, y2,x, y3)
例3:画出指数衰减曲线y1=exp(-t/3)*sin(3*t)和它的包络y2=exp(-t/3),t的取值范围是(0,4pi)。
matlab上机指导书
matlab上机指导书⽬录实习1. MATLAB软件环境实习2. 矩阵的建⽴和运算实习3. 编程操作练习实习4. 图形操作练习实习5. 符号运算练习实习6. 数值运算练习实习7. GUI设计⼊门实习8. SIMULINK建模练习实验1 矩阵的建⽴与运算1.上机⽬的1.1掌握MATLAB 软件的安装与运⾏1.2通过上机启动MATLAB 软件,熟悉软件的界⾯和菜单 1.3掌握变量名的命名规则及常⽤熟悉函数 1.4熟练数据操作2.上机内容例1:数学函数 A=cos(3.14); B=log(2.71928); C=exp(3.14); D=fix(3.14); E=ceil(3.14);例2:数据操作A=input('what's your name?\n','s') Emily %通过键盘输⼊ ans=Emily例3:数据处理求解⽅程02=++c bx ax 的根,其中,a=1,b=2,c=3。
a=1;b=2;c=3;d=sqrt(b^2-4*a*c); x1=(-b+d)/(2*a); x2=(-b-d)/(2*a); disp(x1) disp(x2)3.上机练习题3.1 使⽤gcd 函数求15和35的最⼤公约数。
(提⽰:先⽤help 命令查看函数⽤法);3.2 已知a=4.96,b=8.11,计算)ln(b a e ba +-的值;3.3 定义三个变量a=1,b=2,c=3,全部存⼊⼀个⽂件中,再把a ,b 存⼊另⼀个⽂件中,清空⼯作空间后,检查⼯作空间,再调⼊变量a ,再检查⼯作空间。
实验2 矩阵的建⽴与运算1.上机⽬的1.1熟练矩阵的建⽴1.2掌握矩阵的关系运算和掌握常⽤的矩阵运算指令1.3熟悉⽂件的操作2.上机内容例1:矩阵输⼊Time = [11 12 1 2 3 4 5 6 7 8 9 10]X_Data = [2.32 3.43;4.37 5.98]例2:矩阵运算g = [1 2 3 4];h = [4 3 2 1];s1 = g + h, s2 = g.*h, s3 = g.^h, s4 = g.^2, s5 = 2.^h例3:矩阵运算A=[2 0 –1;1 3 2]; B=[1 7 –1;4 2 3;2 0 1];M = A*B % 矩阵A与B按矩阵运算相乘X = A/B % A/B = A*B-1,即XB=A,求XY = B\A % B\A = B-1*A,即BY=A,求Y3.上机练习题3.1 输⼊A=[7 1 5;2 5 6;3 1 5],B=[1 1 1; 2 2 2; 3 3 3],在命令窗⼝中执⾏下列表达式,掌握其含义: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;3.2 查找已创建变量的信息,删除⽆⽤的变量;3.3 利⽤函数产⽣3X4阶单位矩阵和全部元素都是4.5的4X4阶常数矩阵;3.4 利⽤画图软件画⼀副画,存盘后,读⼊MATLAB⼯具空间,将它改名字改格式另存。
MatLab考试题题库(带答案)
MatLab 考试题题库(必做题)(带答案)一,1.请登陆美国MathWorks公司的网站,查看看现在大概有多少本MATLAB-based books (以MATLAB为基本软件,来说明各个专业领域的教科书或工具书)。
哪一个领域的MATLAB-based books最多中文书共有几本答: 1612 本,数学方面的最多,中文书共有37 本。
2.请在 MATLAB 中直接输入下列常数,看它们的值是多少:a.ib.jc.epsd.infe.nanf.pig.realmaxh.realmin依次解为: ans = 0 +ans = 0 +ans =ans =Inf ans = NaN ans =ans = +308ans =3.试写一函数regPolygon(n) ,其功能为画出一个圆心在(0, 0)、半径为 1 的圆,并在圆内画出一个内接正n 边形,其中一顶点位于(0, 1)。
例如regPolygon(8)可以画出如下之正八边型:解:新建文件如下:function y=regPolyfon(n)n=8;%要画的 n 边形R=1; %圆的半径t=0::2*pi;x=R*cos(t);y=R*sin(t);m=linspace(pi/ 2,5/2*pi,n+1);xz=R*cos(m);yz=R*sin(m);hold onplot(x,y,xz,yz);axis 'equal';4.一条参数式的曲线可由下列方程式表示:x= sin(t), y = 1 - cos(t) + t/10当 t 由 0 变化到4*pi时,请写一个MATLAB 的脚本,画出此曲线在XY 平面的轨迹。
解:新建:t = linspace(0, 4*pi);x = sin(t);y = 1-cos(t)+t/10;plot(x, y, '-o');32.521.510.5-1-0.8-0.6-0.4-0.200.20.40.60.815.当一个小圆轮沿着一条曲线行进时,轮缘任一点的轨迹就会产生变化丰富的摆线。
matlab上机教学
a b c d 1 2 3 4 A 12 15 22 17 e f g h d e a b 1 2 3 4 B 12 15 22 17 5 7 1 6
1 2 b f c
习题 11页 第3题的(1)-(3) 24页第1题 38页第1题(5)-(7) 49页第1题(1)(2) 49页第2题(1)-(3) 126页第1题(1)-(2)
Warning: Could not get change notification handle for local C:\Program Files\MATLAB\R2007b\work. Performance degradation may occur due to on-disk directory change checking. >> >> >> a=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4]; >> b=rref(a) %不打;则计算机将显示b的值 b= 1.0000 0 0.3333 0 1.7778 0 1.0000 0.6667 0 -0.1111 0 0 0 1.0000 -0.3333 0 0 0 0 0 >>
上机界面
MATLAB中基本代数运算符
运算 加法,a+b 减法,a-b 乘法,a×b 除法,a÷b 乘幂,ab 符号 + * / or \ ^ 举例 5+3 5-3 5*3 48/4=4\48=12 5^2=25
变量及数组输入
2:矩阵方式输入(请记下红色部分!) a=[1,2,3;4,5,6;7,8,0] %矩阵输入 (a为3阶方阵) b=[366;804;351] %列矩阵输入 c=[366;804;351] ‘ %行矩阵(转置)输入 det(a) %方阵行列式 inv(a) %方阵的逆 x=a\b %ax=b方程组的解 y=inv(a)*b %与x相同 r=rank(A), %矩阵的秩 disp([a,b,x]) %显示矩阵
完整版本matlab上机练习试题含答案
1.计算 a6 9 3 2 4 1 的数组乘积1 71618 82 7 与 b6 8 5 计算多项式除法 (3x 3+13 x 2+6x+8)/( x+4)54>> a=[6 9 3;2 7 5]; >> d=deconv([3 13 6 8],[1 4])>> b=[2 4 1;4 6 8]; d =>> a.*b 31 2ans =6 求欠定方程组2 4 7 4x 8 的最小范数解12 3639 35 6 584240>> a=[2 4 7 4;9 3 5 6];4 9 2 37>> b=[8 5]';>> x=pinv(a)*b2.对于 AXB ,假如 A7 6 4 ,B 26,求解 X 。
3 5 728x =-0.2151>> A=[4 9 2;7 6 4;3 5 7]; 0.4459 >> B=[37 26 28] ’;0.7949 >> X=A\B 0.2707X =7 用符号函数法求解方程 at 2+b* t+c=0-0.5118 >> r=solve('a*t^2+b*t+c=0','t') 4.0427 r =1.3318[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] 1 2 5 , b8 7 4 [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]3. a64 362 ,察看 a 与 b 之间的3a 11 a 12六种关系运算的结果8求矩阵 A的队列式值、逆和特点根a 21 a 22 >> a=[1 2 3;4 5 6];>> b=[8 –7 4;3 6 2]; >> syms a11 a12 a21 a22; >> a>b >> A=[a11,a12;a21,a22]ans =>> AD=det(A) % 队列式 0 1 0 >> AI=inv(A) % 逆 11>> AE=eig(A) % 特点值>> a>=b A =ans =[ a11, a12] 0 1 0 [ a21, a22] 11AD =>> a<b a11*a22-a12*a21 ans =AI =1 0 1 [ -a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21)]1[ a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21)] >> a<=b AE =ans =1 0 1 [ 1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]1>> a==b [ 1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]ans =9 因式分解: x 45x 3 5 x 2 5x 60 0 0 >> syms x;>> f=x^4-5*x^3+5*x^2+5*x-6; >> a~=b >> factor(f) ans =ans =1 1 1 (x-1)*(x-2)*(x-3)*(x+1) 111ax 214 计算多项式乘法 (x 2+2x+2)(x 2+5x+4),用符号微分求 df/dx 。
上机习题1教案
1.简述MATLAB的主要功能。
主要功能有:数值计算,符号计算,绘图,编程以及应用工具箱。
2.在命令窗口输入“w=3+2”,然后依次使用clear和clc命令,分别观察命令窗口、工作空间窗口和历史命令窗口的变化。
答:使用clear命令:命令窗口以及历史命令窗口没有变化,工作空间窗口清零。
使用clc命令:命令窗口全部清除3.将硬盘上一已有目录,加入到搜索路径,并将其设置为当前工作目录。
答:选择MATLAB主窗口中的【File】|【Set Path】菜单,进入【设置搜索路径】对话框。
点击folder添加文件夹即可设置为当前工作目录。
或者在命令窗口中执行“addpath 目录- end”,将新的目录加到整个搜索路径的末尾。
4.通过命令窗口,查询函数sin()的用法。
答:在命令窗口中输入help sin即可查询到sin()函数及其用法。
5.通过联机帮助系统,查询函数inv()和eig()的用法。
答:直接单击MATLAB主窗口菜单的?键,在此菜单中选择search再输入inv()和eig()进行查询,在打开的右窗口中即可查询其用法。
6.通过联机演示系统,查询并运行“Control Systems Toolboxes”下的“Case Studies”中的“Yaw Damper for a 747 Aircraft”演示程序。
答:直接在帮助页面上选择Demos,在Demos页面中打开contents,在其中找Control Systems Toolboxes”下的“Case Studies”中的“Yaw Damper for a 747 Aircraft”,单击右边窗口中的运行这个程序即可。
7.在命令窗口依次执行“w=5;”、“p=2*w”和“q=p+w”。
答:在命令窗口中输入w=5;回车后无显示。
再输入p=2*w按下回车,得到结果10.输入q=p+w按下回车得到结果15.或输入代码w=5;p=2*w;q=p+w再回车8.在命令窗口同时执行下述代码:w=5;p=2*wq=p+w答:要同时执行则在输入下一条命令时,按住shift和enter键进行换行。
Matlab考试题库+答案(教师出卷参考专用)
Matlab考试题库+答案(教师出卷参考专用)[1 2;3 4][1 0;0 4]21. A=[1 2 3;6 5 4]; reshape(A,3,2)=______________.[1 5;6 3;2 4]22. A=[1,2;3,1]; B=[1,0;0,1]; A~=B=___________.[0 1;1 0]23. A=[1,2,3;1,4,9;1,7,27]; triu(A,-1)=_ ____.提取矩阵上三角矩阵的函数[1 2 3;1 4 9;0 7 27] 24. A=[4,15,-45,10,6;56,0,17,-45,0]; find(A>=10 & A<=20)=___ __________.3 6 725.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。
A &B 的运行结果_______;A ./B 的运行结果_______。
ans=[0 1 1;1 1 0]ans= [0 2.25 2;1 0.6 NaN]26. 产生4阶全0方阵的命令为_______;产生3阶全1方阵的命令为。
zeros(4)ones(3)27. A=rand(2,5);b=size(A);c=length(A);则b和c的值分别为_______ 。
2 5和 528. 用if判断语句判断80≥ x >60,在MATLAB中if语言后的判断应写为_______。
if x<=80 & x>6029.为了使两个pl ot的图形在同一个坐标显示,可以使用_______命令进行图形保持;可以使用_______命令为图形添加网格。
hold ongrid on30. MATLAB的工作空间中有三个变量v1, v2, v3,写出把它们保存到文件my_data.mat中的指令_______;save my_date31. MATLAB提供了两种运行方式,是按_______和把_______文件copy到matlab里,回车。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
<<<<<<精品资料》》》》》
Matlab 上机测试题答案
班级 姓名 学号
要求: (1)请将程序和运行结果粘贴到题目下面,写清楚哪部分是程序,哪部分是运行结果,运行结果包括图形,图形请适当缩小,然后将word 文档传到教育在线,注意保存格式是2003的doc ,文件名”自111_110110_张三”
(2)在最开始前请将输入Matlab 命令窗口中输入”diary ”命令,随后在workspace 窗口中可以发现一个名为“diary ”的文件。
在老师要求结束练习后,在命令窗口输入”diary off ”命令,然后将文件“diary ”上传到教育在线。
(3)注意题号不要弄乱。
---------------------------------------------------------------------
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、已知矩阵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
--------------------------------------------------------------------- 3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵,并计算两者的乘积。
解:命令为: A=eye(3,4) B=8*ones(4) C=A*B
---------------------------------------------------------------------
4、创建矩阵a=⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601
,取出其前两列构成的矩阵b ,取出前两行构成矩阵c ,转置矩阵b
构成矩阵d ,计算a*b 、c<d ,c&d, c|d ,~c|~d 解:命令为:
a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7] b=a(:,[1 2]) c=a([1 2],:) d=b ’ e=a*b f=c<d g=c&d h=c|d i=~c|~d
--------------------------------------------------------------------- 5、求
!20
1
∑=n n
解:命令文件为
sum=0; s=1;
for n=1:20 s=n*s;
sum=sum+s; end sum
---------------------------------------------------------------------
<<<<<<精品资料》》》》》
6、求a aa aaa aa a S n ++++=得值,其中a 是一个数字,由键盘输入,表达式中位数最多项a 的个数,也由键盘输入。
例如:输入a=3,n=4,
解:命令文件为
sum=0; s=0;
a=input('请输入a 的值:') n=input('请输入n 的值:') for i=1:n
s=s+a*10^(i-1); sum=sum+s;
end sum
--------------------------------------------------------------------- 7、计算下列分段函数的值:
⎪
⎩⎪⎨⎧>≤<-≤=1021054350)(x x
x x x x f
解:命令文件为
x=input('输入x 的值:') if x<=5 f=0;
elseif x>5&x<=10 f=(x-3)/4.0; else
f=2*x; end f
--------------------------------------------------------------------- 8、有一群鸡和兔子,加在一起头的数量是36,脚的数量是100,编程序解答鸡和兔子数量各是多少? 解:命令文件为 for x=0:36
for y=0:36-x
if (x+y==36)&(2*x+4*y==100) disp('鸡有:') disp(x)
disp('兔子有') disp(y) end
end end
--------------------------------------------------------------------- 9、(1) 化简表达式2
cos 22sin f x x =+; (2)对表达式62
4--x x 进行因式分解。
解:命令文件为 syms x
f1=cos(2*x)+(sin(x))^2 f2=x^4-x^2-6 t1=simple(f1) t2=factor(f2)
--------------------------------------------------------------------- 10、求下列极限值(1)0sin 2lim sin 5x x
x
->
(2
)x →∞
解:命令文件为
syms x
f1=sin(2*x)/sin(5*x)
f2=sqrt(x^2+x)-sqrt(x^2-x) t1=limit(f1,x,0) t2=limit(f2,x,inf)
--------------------------------------------------------------------- 11、求下列级数之和
23(1)123(2)21491610000
n s x x x nx s =+++=++++
解:命令文件为: syms x n m s1=n*x^n s2=n^2
t1=symsum(s1,n,1,n) t2=symsum(s2,n ,1,100)
--------------------------------------------------------------------- 12、在区间1<x<4,同一窗口的两个子窗口内分别绘制函数x x +-12和)35sin(2
+x 曲线,第一个子窗口的采用红色点线,并且要求坐标轴采用等刻度,第二个子窗口绿色点化线,并且要求坐标轴为正方形。
请将程序与结果包括图形均粘贴到下面.
<<<<<<精品资料》》》》》
解: figure;
subplot(121); axis equal;
fplot('sqrt(x^2-1)+x',[1 4],'r:') subplot(122);
fplot('sin(5*x^2+3)',[1 4],'-.g') axis square;
--------------------------------------------------------------------- 13、x []2,2-∈,y ]2,2[-∈,自己设定步长,产生三维网格X ,Y ,并在一个图形窗口的四个子窗口分别绘制2
2y x e
x z --*=的三维网格图,带底座的三维网格
图,三维曲面图和带等高线的三维曲面图。
解:
figure;
[X,Y]=meshgrid(-2:0.1:2);Z=X.*exp(-X.^2-Y .^2);
subplot(2,2,1)
mesh(X,Y ,Z) %mesh 绘制三维网格图,meshz 绘制带底座的三维网格图title ‘三维网格图’ subplot(2,2,2) meshz(X,Y,Z)
title ‘带底座的三维网格图’
subplot(2,2,3) surf(X,Y ,Z)
title ‘三维曲面’
subplot(2,2,4)
surfc(X,Y ,Z)
title ‘带等高线的三维曲面’。