数学建模matlab例题参考及练习

合集下载

数学建模Matlab练习4

数学建模Matlab练习4

MATLAB练习四1. 绘制两个一元函数的曲线图(例如y=sin(x)、y=cos(x)),加上标注说明,包括坐标、抬头、曲线说明等。

%lianxi41.mx=0:pi/30:2*pi;y=sin(x);plot(x,y,'-ob','markersize',10);hold ony1=cos(x);plot(x,y1,'-.dm','markersize',10);legend('sinx','cosx');title('sin和cos函数的曲线图')xlabel('x');ylabel('y');2. 寻找一组数据拟合多项式,写出M文件,包括拟合的次数,数值的计算,数学表达式的呈现,图形的呈现等。

%lianxi42x=0:pi/20:pi;y=sin(x);n=input('n=');p=polyfit(x,y,n);x1=0:pi/20:pi*2;y1=sin(x1);y2=polyval(p,x1)plot(x1,y1,'-ob',x1,y2,':dm');legend('原曲线','拟合曲线');xlabel('x');ylabel('y');title('sin的拟合曲线和原曲线');axis([-0.2,6.5,-1.2,1.2])3. 求解一个代数方程,写出M文件,应该包括函数式M文件和命令式M文件。

%fc3.mfunction f=fc3(x)f='3*x-1=0';%lianxi43clc;clear;close;sprintf('原方程和其解析解的结果为:')s=solve(fc3)sprintf('保留10位有效数字的结果为:')digits(10)sprintf('数值解的结果为:')x=vpa(s)4. 求解一个代数方程组,写出M文件,应该包括函数式M文件和命令式M文件。

数学建模题目2道 解答 matlab

数学建模题目2道 解答 matlab

1.某市为了创建园林型城市,计划在该市的望湖公园里栽花1800支,由于施工队加强了施工力量,实际每天栽花的数量是原计划的3倍,最终提前2天完成了栽花任务,
问题一:原计划每天栽花的数量是多少支?
问题二:用Matlab求出问题一中的值。

解:问题一:设原计划每天栽花的数量是x支,由题意可得等量关系为:
原计划所用天数-实际所用天数=2,
所以根据等量关系可以列出方程:
18001800
-=
2
3
x x
x =600
故原计划每天栽花的数量是600支。

问题二:
输入命令:
syms x
x=solve(1800/x-1800/(3*x)==2,x) %解方程
运行结果:
x =600
2.一牧场每天支出15元经费用于养一头牛,估计可使一头250公斤重的生牛每天增加5公斤。

目前生牛市场售价为每公斤70元,但是预测每天会降低0.5元,问该牧场应该什么时候出售这样的生牛,可以使获得的利润最大。

(要求用Matlab求解,并附有Matlab 截图)
解:设在第x天出售这样的生牛(初始重250公斤的牛)可以获得的利润为y元。

每头牛投入:15x元
产出:(70-0.5x)(250+5x)元
利润:y=(70-0.5x)(250+5x)-15x=-2.5x^2+210x+17500
使用Matlab软件计算,得
当天数x=42时,利润y取得最大值21910(元)。

故该牧场应该在第42天时候出售这样的生牛,可以使获得的利润最大。

附Matlab计算截图:
①利润y取得最大值的Matlab计算截图
②天数x=42天的Matlab计算截图。

MATLAB习题及参考答案经典.doc

MATLAB习题及参考答案经典.doc

习题:1, 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。

2, 对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。

3, 已知:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察其结果。

4, 角度[]604530=x ,求x 的正弦、余弦、正切和余切。

(应用sin,cos,tan.cot)5, 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵: (1)组合成一个4⨯3的矩阵,第一列为按列顺序排列的a 矩阵元素,第二列为按列顺序排列的b 矩阵元素,第三列为按列顺序排列的c 矩阵元素,即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574(2)按照a 、b 、c 的列顺序组合成一个行矢量,即 []2965318772546, 将(x -6)(x -3)(x -8)展开为系数多项式的形式。

(应用poly,polyvalm)7, 求解多项式x 3-7x 2+2x +40的根。

(应用roots)8, 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。

(应用poly,polyvalm)9, 计算多项式9514124234++--x x x x 的微分和积分。

(应用polyder,polyint ,poly2sym)10, 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。

(应用x=a\b)11, 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。

(应用pinv) 12, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ,计算a 的行列式和逆矩阵。

(应用det,inv)13, y =sin(x ),x 从0到2π,∆x =0.02π,求y 的最大值、最小值、均值和标准差。

matlab建模习题

matlab建模习题

1.用dsolve求常微分方程的初值问题dy/dx+3y=8,y|x=0=2的解>>r=dsolve('Dy+3*y=8','y(0)=2','x')r=8/3-2/3*exp(-3*x)2.用dsolve求常微分方程的初值问题(1+x^2)y’’=2xy’,y|x=0=1,y’|x=0=3的解>>r=dsolve('D2y*(1+x^2)=2*x*Dy','y(0)=1,Dy(0)=3','x')r=1+x^3+3*x3.用dsolve求微分方程y’’’’-2y'''+y’’=0的解>>r=dsolve('D4y-2*D3y+D2y=0','x')r=C1*exp(x)+C2*exp(x)*x+C3+C4*x4.用dsolve求微分方程组的特解{2*dx/dt+4*x+dy/dt-y=e^t,x|t=0=3/2;{dx/dt+3x+y=0,y|t=0=0;>>[X,Y]=dsolve('2*Dx+4*x+Dy-y=exp(t),Dx+3*x+y=0','x(0)=3/2','y(0)=0')X=-3/2*exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)-3/2*exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)+1/6*exp(t)+1/2 *7^(1/2)*exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)-1/2*7^(1/2)*exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)Y=exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)+exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)-2/3*exp(t)1.选择适当的ode函数,求常微分方程的初值问题dy/dx+3y=8,y|x=0=2的解dy/dx+3y=8即是y’=8-3*y;y|x=0=2即是y(0)=2;故建立m-file文件function test1ode45(@fun,[0,1],2)%---%function f=fun(x,y)f=8-3*y;run之后结果2.选择适当的ode函数,求常微分方程的初值问题(1+x^2)y’’=2xy’,y|x=0=1,y’|x=0=3的解。

matlab考试题目及答案

matlab考试题目及答案

matlab考试题目及答案1. 题目:编写一个MATLAB函数,实现计算并返回一个向量中所有元素的平方和。

答案:以下是一个简单的MATLAB函数,用于计算并返回一个向量中所有元素的平方和。

```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目:给定一个3x3的矩阵A,使用MATLAB编写代码,求出矩阵A 的转置。

答案:可以通过简单的转置操作来求得矩阵A的转置。

```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目:编写一个MATLAB脚本,实现对一个二维数组进行排序,并返回排序后的数组。

答案:以下是一个MATLAB脚本,用于对一个二维数组进行排序,并返回排序后的数组。

```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目:给定一个向量x,使用MATLAB编写代码,计算并返回向量x的元素个数。

答案:可以通过内置函数`numel`来计算向量x的元素个数。

```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目:编写一个MATLAB函数,实现计算并返回两个向量元素的点积。

答案:以下是一个简单的MATLAB函数,用于计算两个向量的点积。

```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目:给定一个矩阵B,使用MATLAB编写代码,求出矩阵B的行列式。

matlab试题及答案

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的矩阵。

数学建模基础练习一及参考答案

数学建模基础练习一及参考答案

练习1matlab练习一、矩阵及数组操作:1.利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵([-1,1]之间)、正态分布矩阵(均值为1,方差为4),然后将正态分布矩阵中大于1的元素变为1,将小于1的元素变为0。

2.利用fix及rand函数生成[0,10]上的均匀分布的10×10的整数随机矩阵a,然后统计a中大于等于5的元素个数。

3.在给定的矩阵中删除含有整行内容全为0的行,删除整列内容全为0的列。

4.随机生成10阶的矩阵,要求元素值介于0~1000之间,并统计元素中奇数的个数、素数的个数。

二、绘图:5.在同一图形窗口画出下列两条曲线图像,要求改变线型和标记:y1=2x+5;y2=x^2-3x+1,并且用legend标注。

6.画出下列函数的曲面及等高线:z=sinxcosyexp(-sqrt(x^2+y^2)).7.在同一个图形中绘制一行三列的子图,分别画出向量x=[158101253]的三维饼图、柱状图、条形图。

三、程序设计:8.编写程序计算(x 在[-8,8],间隔0.5)先新建的,在那上输好,保存,在命令窗口代数;9.用两种方法求数列:前15项的和。

10.编写程序产生20个两位随机整数,输出其中小于平均数的偶数。

11.试找出100以内的所有素数。

12.当)1(433221)(+⨯++⨯+⨯+⨯=n n n f 时,?)20()30()40(=+f f f 四、数据处理与拟合初步:13.随机产生由10个两位随机数的行向量A,将A 中元素按降序排列为B,再将B 重排为A 。

14.通过测量得到一组数据:t 12345678910y 4.8424.3623.7543.3683.1693.0383.0343.0163.0123.005分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合,并画出散点及两条拟合曲线对比拟合效果。

matlab模拟试题及答案

matlab模拟试题及答案

matlab模拟试题及答案MATLAB模拟试题及答案1. 编写一个MATLAB函数,计算并返回一个向量中所有元素的和。

```matlabfunction sumVector = sumVectorElements(vector)sumVector = sum(vector);end```2. 给定一个矩阵A,编写一个MATLAB脚本,找出矩阵中的最大元素及其位置。

```matlabA = [1, 2, 3; 4, 5, 6; 7, 8, 9];[maxValue, maxIndex] = max(A(:));[maxRow, maxCol] = ind2sub(size(A), maxIndex);disp(['最大元素为: ', num2str(maxValue), ',位于第',num2str(maxRow), '行,第', num2str(maxCol), '列。

']);```3. 编写一个MATLAB函数,实现一个向量元素的排序。

```matlabfunction sortedVector = sortVector(vector)sortedVector = sort(vector);end4. 给定两个向量X和Y,编写一个MATLAB脚本,计算这两个向量的点积。

```matlabX = [1, 2, 3];Y = [4, 5, 6];dotProduct = dot(X, Y);disp(['X和Y的点积为: ', num2str(dotProduct)]);```5. 编写一个MATLAB函数,计算一个矩阵的行列式。

```matlabfunction determinant = calculateDeterminant(matrix)determinant = det(matrix);end```6. 给定一个向量,编写一个MATLAB脚本,找出向量中的最大值和最小值。

MATLAB数学建模习题-推荐下载

MATLAB数学建模习题-推荐下载

e=
x=[x1,x2,x3,x4];
A=[e;x;x.^2;x.^3];
D=det(A);。

2.反正切函数的 Tylor 级数展开式
级数 4
1 1 1 1
算。 function z=adds(errors) n=1; an=1; f=1;

arctan x (1)n x2n1
对全部高中资料试卷电气设备,在安装过程中以及安装结束后进行高中资料试卷调整试验;通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术通关,1系电过,力管根保线据护敷生高设产中技工资术艺料0不高试仅中卷可资配以料置解试技决卷术吊要是顶求指层,机配对组置电在不气进规设行范备继高进电中行保资空护料载高试与中卷带资问负料题荷试2下卷2,高总而中体且资配可料置保试时障卷,各调需类控要管试在路验最习;大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试,况卷要下安加与全强过,看度并22工且22作尽22下可22都能22可地护以缩1关正小于常故管工障路作高高;中中对资资于料料继试试电卷卷保破连护坏接进范管行围口整,处核或理对者高定对中值某资,些料审异试核常卷与高弯校中扁对资度图料固纸试定,卷盒编工位写况置复进.杂行保设自护备动层与处防装理腐置,跨高尤接中其地资要线料避弯试免曲卷错半调误径试高标方中高案资等,料,编试要5写、卷求重电保技要气护术设设装交备备置底4高调、动。中试电作管资高气,线料中课并敷3试资件且、设卷料中拒管技试试调绝路术验卷试动敷中方技作设包案术,技含以来术线及避槽系免、统不管启必架动要等方高多案中项;资方对料式整试,套卷为启突解动然决过停高程机中中。语高因文中此电资,气料电课试力件卷高中电中管气资壁设料薄备试、进卷接行保口调护不试装严工置等作调问并试题且技,进术合行,理过要利关求用运电管行力线高保敷中护设资装技料置术试做。卷到线技准缆术确敷指灵设导活原。。则对对:于于在调差分试动线过保盒程护处中装,高置当中高不资中同料资电试料压卷试回技卷路术调交问试叉题技时,术,作是应为指采调发用试电金人机属员一隔,变板需压进要器行在组隔事在开前发处掌生理握内;图部同纸故一资障线料时槽、,内设需,备要强制进电造行回厂外路家部须出电同具源时高高切中中断资资习料料题试试电卷卷源试切,验除线报从缆告而敷与采设相用完关高毕技中,术资要资料进料试行,卷检并主查且要和了保检解护测现装处场置理设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况,然后根据规范与规程规定,制定设备调试高中资料试卷方案。

数学建模matlab例题参考及练习

数学建模matlab例题参考及练习

数学建模matlab例题参考及练习数学实验与数学建模实验报告学院:专业班级:姓名:学号:完成时间:年⽉⽇承诺书本⼈承诺所呈交的数学实验与数学建模作业都是本⼈通过学习⾃⾏进⾏编程独⽴完成,所有结果都通过上机验证,⽆转载或抄袭他⼈,也未经他⼈转载或抄袭。

若承诺不实,本⼈愿意承担⼀切责任。

承诺⼈:年⽉⽇数学实验学习体会(每个⼈必须要写字数1200字以上,占总成绩的20%)练习1 ⼀元函数的图形 1.画出x y arcsin =的图象.2.画出x y sec =在],0[π之间的图象. 3.在同⼀坐标系中画出x y =,2x y =,3x y =,3x y =,x y =的图象.4.画出3232)1()1()(x x x f ++-=的图象,并根据图象特点指出函数)(x f 的奇偶性.5.画出)2ln(1++=x y 及其反函数的图象.6.画出321+=x y 及其反函数的图象.练习2 函数极限1.计算下列函数的极限.(1)xxx1lim4-+π→.程序:sym x;f=(1+sin(2*x))/(1-cos(4*x)); limit(f,x,pi/4)运⾏结果:lx21ans =1(2).程序:sym x;f=(1+cos(x))^(3*sec(x)); limit(f,x,pi/2)运⾏结果:lx22ans =exp(3)(3)22)2xx-ππ→.程序:sym x;f=log(sin(x))/(pi-2*x)^2; limit(f,x,pi/2)运⾏结果:lx23ans =-1/8(4)212lim xxex→.程序:x xx sec32)sym x ;f=x^2*exp(1/x); limit(f,x,0) limit(f,x,0,'right') limit(f,x,0,'left')运⾏结果:lx24ans = NaNans = Infans = 0%左极限为零,存在,右极限为⽆穷⼤,在x 趋近于零时函数没有极限(5))215(lim 122x x x x +-∞→.程序:sym x ;f=5*x^2/(1-x^2)+2^(1/x); limit(f,x,inf)运⾏结果:>> lx25ans = -4(6)x x x x x -+-→32112lim .程序:sym x ;f=(x^2-2*x+1)/(x^3-x); limit(f,x,1)运⾏结果:>> lx26ans = 0(7)x x x 11lim 20-+→.程序:sym x ;f=(sqrt(1+x^2)-1)/x; limit(f,x,0))3sin(cos 21lim 3π--π→x x x . 程序:sym x ;f=(1-2*cos(x))/sin(x-pi/3); limit(f,x,pi/3)运⾏结果:>> lx28ans = 3^(1/2)(9)tgxx x )1(lim 0+→.程序:sym x ;f=(1/x)^tan(x); limit(f,x,0,'right')运⾏结果:>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序:sym x ;f=(2/pi*atan(x))^x; limit(f,x,inf,'left')运⾏结果:>> lx210ans =Inf2.解⽅程012=-?x x . 程序:sym x ;X=solve(x*2^x-1)运⾏结果:>> lx202 X =lambertw(0, log(2))/log(2)%⽅程有两个解X=solve(3*sin(x)+1-x)运⾏结果:>> lx203 X =-0.53847936154.解⽅程03=++q px x .(p 、q 为实数) 程序:X=solve('x^3+p*x+q=0','x')运⾏结果: X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) -((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)))/2 p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数.(1))11)(1(-+=x x y程序:sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1); diff(f)运⾏结果:>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) + 1)/(2*x^(3/2))(2)x x x y ln sin =程序:sym x ;f=x*sin(x)*log(x); diff(f)运⾏结果:>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数⽅程所确定的函数的导数.(1)??==t y t x 44程序:ans =1/t^3(2)??-=+=arctgt t y t x )1ln(2程序:sym t ;f1=log(1+t^2);f2=t-atan(t); diff(f2)/diff(f1)运⾏结果:>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t) 3.求下列隐函数的导数.(1)22ln y x xyarctg+=程序:syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运⾏结果;>> lx331 yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 + 1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2)) (2)x y y x=程序:syms x y ; f=x^y-y^xyx=-diff(f,x)/diff(f,y)运⾏结果:>> lx332 f =x^y - y^x yx =f=exp(x)*sin(x); diff(f,x,4)运⾏结果:>> lx34 ans =(-4)*exp(x)*sin(x)5.验证x e y xsin =满⾜关系式:022=+'-''y y y程序:sym x ;f=exp(x)*sin(x); y2=diff(f,x,2); y1=diff(f,x,1); y=f;y2-y1*2+2*y=='0' 运⾏结果:>> lx35ans =1%运⾏结果为1表⽰y2-y1*2+2*y=='0'成⽴6.设)ln(y x x u +=,求22x u ??,22y u,y x u 2. 程序:syms x y ; f=x*log(x+y); uxx=diff(f,x,2) uyy=diff(f,y,2) f1=diff(f,x); uxy=diff(f1,y)运⾏结果: >> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ,.(1)1cos cos cos 222=++z y x程序:syms x y z ;-(cos(x)*sin(x))/(cos(z)*sin(z)) zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z= 程序:syms x y z ; f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z) zy=-diff(f,y)/diff(f,z)运⾏结果:>> lx372 f =exp(z) - x*y*z zx =(y*z)/(exp(z) - x*y) zy =(x*z)/(exp(z) - x*y) 8.证明函数22)()(lnb y a x u -+-=(b a ,为常数)满⾜拉普拉斯⽅程:02222=??+??y u x u (提⽰:对结果⽤simplify 化简)练习4 积分计算1.计算下列不定积分.(1)?+dxx x 12 (2)+x xdx 2sin 12sin2.计算下列定积分.(1)?exdxx 1ln (2)ππ342sin dxxx3.求?+tdx x x x4.求摆线)cos 1(),sin (t a y t t a x -=-=的⼀拱(π≤≤20t )与x 轴所围成的图形的⾯积.5.计算⼆重积分 (1)??≤++122)(y x dxdyy x (2)??≤++xy x dxdyy x 22)(226.计算?+Ldsy x 22 L 为圆周)0(22>=+a ax y x7.计算?++-L dy y x dx y x )()(2222,其中L 为抛物线2x y =上从点(0,0)到点(2,4)的⼀段弧.练习5 matlab ⾃定义函数与导数应⽤1.建⽴函数x x a a x f 3sin 31sin ),(+=,当a 为何值时,该函数在3π=x 处取得极值,它是极⼤值还是极⼩值,并求此极值.2.确定下列函数的单调区间.(1)7186223---=x x x y (2))0(82>+=x xx y3.求下列函数的最⼤值、最⼩值.(1)2332x x y -=41≤≤-x(2)312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1.计算(1)???--521111204321+???21(2)??-01301213?03010*******????? ??-205101(3)52422??- 2.设????? ??-=243121013A ,??-=112111201B ,求满⾜关系B X A =-23的X .练习7 矩阵与线性⽅程组1.求下列矩阵的秩.(1)???-321110021 (2)4820322513454947513253947543173125 2.求下列矩阵的⾏列式,如可逆,试⽤不同的⽅法求其逆矩阵.(1)??--285421122 (2)??---6201111121324321 3.设X ????? ?-111012111==--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分⽅程与级数求1-6题微分⽅程的通解1.1222+='y y y x 2.x y x y dx dy -+= 3.x xx y y +='cos 4.1)2sin cos (='+y y y x 5.x e y y y x2cos 3=-'+'' 6.x x y y sin 14++=+'' 求7、8题初值问题的解7.==-++-+=10)2(212222x y dx dy x xy y y xy x8.===++==0000222,02V dt dx x x x a dt dxn dtx d t t9.给出函数x x e x f xx cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式.10.判别下列级数的敛散性,若收敛求其和.(1)+++311(2)∑∞=+112n nntgπ11.求幂级数∑∞=--22)1(nnnnnx的和函数.12.求函数项级数∑∞=-1nnnn xπ的和函数.。

matlab试题及答案开卷

matlab试题及答案开卷

matlab试题及答案开卷1. MATLAB基础操作在MATLAB中,如何创建一个名为“myMatrix”的3x3矩阵,其元素从1到9按顺序排列?答案:在MATLAB中,可以使用以下命令创建名为“myMatrix”的3x3矩阵:```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用内置函数`magic`:```matlabmyMatrix = magic(3);```2. 矩阵运算给定两个矩阵A和B,其中A = [2 3; 4 5],B = [6 7; 8 9],计算A+B的结果。

答案:在MATLAB中,可以使用加号`+`来计算两个矩阵的和:```matlabA = [2 3; 4 5];B = [6 7; 8 9];C = A + B;```计算结果C为:```matlabC = [8 10; 12 14];```3. 向量操作创建一个从0到1的等差数列,步长为0.1。

答案:在MATLAB中,可以使用`linspace`函数或`:`操作符来创建等差数列:```matlabv = 0:0.1:1;```或者使用`linspace`函数:```matlabv = linspace(0, 1, 11);```两种方法都将得到一个包含11个元素的向量,从0开始,到1结束,步长为0.1。

4. 条件语句编写一个MATLAB脚本,判断一个给定的数n是否为素数。

答案:可以通过以下MATLAB脚本来判断一个数n是否为素数:```matlabfunction isPrime = isPrimeNumber(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```调用此函数时,传入一个数值n,函数将返回一个布尔值,表示n 是否为素数。

matlab试题及答案

matlab试题及答案

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循环两种循环结构。

Matlab与数学建模综合练习3

Matlab与数学建模综合练习3

1、画sin(6x)的图,x ∈[0,68×Pi];>> x=0:1:68*pi; >> plot(x,sin(6*x))2、画],0[],,0[),sin()sin(2Pi p Pi t p t ∈∈⨯的三维图; >> t=0:0.1:pi;p=t;>> [T,P]=meshgrid(t,p);>> Z=sin(T).*sin(P.^2);>> surf(T,P,Z)>> mesh(T,P,Z)3、画图]6,6[],6,6[,222-∈-∈+y x ey x 和]6,6[],6,6[),sin(22-∈-∈+y x y x ,并合并这两个图。

>> x=-6:0.3:6;y=x;>> [X,Y]=meshgrid(x,y);>> Z1=2.*exp(sqrt(X.^2+Y.^2));>> surf(X,Y,Z1)>> hold on>> Z2=sin(sqrt(X.^2+Y.^2));>> surf(X,Y,Z2)>> x=-6:0.3:6;y=x;>> [X,Y]=meshgrid(x,y);>> Z1=2.*exp(sqrt(X.^2+Y.^2)); >> mesh(X,Y,Z1)>> hold on>> Z2=sin(sqrt(X.^2+Y.^2)); >> mesh(X,Y,Z2)4、 设⎪⎩⎪⎨⎧=+=+32/)7(11x x x x n n n ,数列}{n x 是否收敛?若收敛,其值为多少?精确到6位有效数字。

5、设 ,131211p p p n nx ++++= }{n x 是否收敛?若收敛,其值为多少?精确到17位有效数字。

MATLAB数学建模习题

MATLAB数学建模习题

MATLAB 数学建模习题1一、单项选择题(将选择答案写在答题纸上,每小题2分共20分)1.在MA TLAB 命令窗口中键入命令,Vname=prod(7:9)/prod(1:3),可计算组合数!6!3!939⨯=C ,如果省略了变量名Vname ,MA TLAB 表现计算结果将用下面的哪一变量名做缺省变量名A )ans ;B )pi ;C )NaN ;D )eps2.宝石切割问题中,石料左右长度、前后长度、上下高度分别为a 1、a 2、a 3,即a 1×a 2×a 3(cm 3),而精品尺寸为b 1×b 2×b 3(cm 3)。

操作时,同向切割连续两次再旋转刀具。

某一切割方案的切割面积依次为:2a 1a 2→ 2a 1b 3 → 2b 2b 3,则这一切割方案为A )左右→前后→上下;B )上下→前后→左右;C )前后→上下→左右;D )前后→ 左右→上下3.机场指挥塔位置:北纬30度35.343分,东经104度2.441分,在MA TLAB 中用变量B=[30 35.343]表达纬度,L=[104 2.441]表达经度。

将数据转化为以度为单位的实数,下面正确的语句是A ) P=B(1)+B(2)/60,Q=L(1)+L(2); B) P = 60*B(1) + B(2),Q=60*L(1)+L(2)C ) P = B(1) + B(2)/60,Q=L(1)+L(2)/60; D) P=B(1)+B(2),Q=L(1)+ L(2);。

4.用MA TLAB 随机产生60个1到365之间的正整数,应该使用下面的哪一条命令A ) fix(365*rand(1,60));B )1+fix(366*rand(1,60));C )1+fix(364*rand(1,60));D )1+fix(365*rand(1,60))5.用A 、B 、C 表示三角形的三条边,用MA TLAB 表示条件“任意两条边之和大于第三条边”的逻辑表达式应该用下面哪一行语句A ) A+B>C | A+C>B | B+C>A ; B ) A+B>=C | A+C>=B | B+C>=A ;C ) A+B>=C&A+C>=B&B+C>=A ;D ) A+B>C & A+C>B & B+C>A ;6.在MATLAB 命令窗口中,键入命令syms x ; y=int(6*x^4)。

数学建模matlab例题参考及练习讲课稿

数学建模matlab例题参考及练习讲课稿

数学建模m a t l a b例题参考及练习数学实验与数学建模实验报告学院:专业班级:姓名:学号:完成时间:年月日承 诺 书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成,所有结果都通过上机验证,无转载或抄袭他人,也未经他人转载或抄袭。

若承诺不实,本人愿意承担一切责任。

承诺人:年 月 日数学实验学习体会(每个人必须要写字数1200字以上,占总成绩的20%)练习1 一元函数的图形1. 画出x y arcsin =的图象.2. 画出x y sec =在],0[π之间的图象.3. 在同一坐标系中画出x y =,2x y =,3x y =,3x y =,x y =的图象.4. 画出3232)1()1()(x x x f ++-=的图象,并根据图象特点指出函数)(x f 的奇偶性.5. 画出)2ln(1++=x y 及其反函数的图象.6. 画出321+=x y 及其反函数的图象.练习2 函数极限1. 计算下列函数的极限. (1)x xx 4cos 12sin 1lim 4-+π→.程序:sym x ;f=(1+sin(2*x))/(1-cos(4*x));limit(f,x,pi/4)运行结果:lx21ans =1(2).程序:sym x ;f=(1+cos(x))^(3*sec(x));limit(f,x,pi/2)运行结果:lx22ans =exp(3)(3)22)2(sin ln lim x xx -ππ→.程序:sym x ;f=log(sin(x))/(pi-2*x)^2;limit(f,x,pi/2)运行结果:lx23ans =x x x sec 3 2 ) cos 1 ( lim + π →-1/8(4)212lim xxex→.程序:sym x;f=x^2*exp(1/x);limit(f,x,0)limit(f,x,0,'right')limit(f,x,0,'left')运行结果:lx24ans =NaNans =Infans =%左极限为零,存在,右极限为无穷大,在x趋近于零时函数没有极限(5))215(lim122xx xx+-∞→.程序:sym x;f=5*x^2/(1-x^2)+2^(1/x);limit(f,x,inf)运行结果:>> lx25ans =(6)x x x x x -+-→32112lim .程序:sym x ;f=(x^2-2*x+1)/(x^3-x);limit(f,x,1)运行结果:>> lx26ans =0 (7)x x x 11lim 20-+→.程序:sym x ;f=(sqrt(1+x^2)-1)/x;limit(f,x,0)运行结果:>> lx27ans =0 (8))3sin(cos 21lim 3π--π→x xx . 程序:sym x ;f=(1-2*cos(x))/sin(x-pi/3);limit(f,x,pi/3)运行结果:>> lx28ans =3^(1/2) (9)tgxx x )1(lim 0+→.程序:sym x ;f=(1/x)^tan(x);limit(f,x,0,'right')运行结果:>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序:sym x ;f=(2/pi*atan(x))^x;limit(f,x,inf,'left')运行结果:>> lx210ans =Inf2. 解方程012=-⋅x x .程序:sym x ;X=solve(x*2^x-1)运行结果:>> lx202X =lambertw(0, log(2))/log(2)%方程有两个解3. 解方程1sin 3+=x x .程序:sym x ;X=solve(3*sin(x)+1-x)运行结果:>> lx203X =-0.538470451711254993610615326557454. 解方程03=++q px x .(p 、q 为实数)程序:X=solve('x^3+p*x+q=0','x')运行结果:X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 +q^2/4)^(1/2) - q/2)^(1/3))p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数. (1))11)(1(-+=x x y程序:sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1);diff(f)运行结果:>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) +1)/(2*x^(3/2))(2)x x x y ln sin =程序:sym x ;f=x*sin(x)*log(x);diff(f)运行结果:>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数方程所确定的函数的导数.(1)⎩⎨⎧==t y t x 44程序:sym t ;f1=t^4;f2=4*t;diff(f2)/diff(f1)运行结果:>> lx3211/t^3(2)⎩⎨⎧-=+=arctgt t y t x )1ln(2程序:sym t ;f1=log(1+t^2);f2=t-atan(t);diff(f2)/diff(f1)运行结果:>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t)3.求下列隐函数的导数. (1)22ln y x x y arctg +=程序:syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运行结果;>> lx331yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 +1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2))(2)x y y x =程序:syms x y ;f=x^y-y^xyx=-diff(f,x)/diff(f,y)运行结果:>> lx332f =x^y - y^x(x^(y - 1)*y - y^x*log(y))/(x*y^(x - 1) - x^y*log(x))4.设x e y x cos =,求)4(y .程序:sym x ;f=exp(x)*sin(x);diff(f,x,4)运行结果:>> lx34ans =(-4)*exp(x)*sin(x)5.验证x e y x sin =满足关系式:022=+'-''y y y程序:sym x ; f=exp(x)*sin(x);y2=diff(f,x,2);y1=diff(f,x,1);y=f;y2-y1*2+2*y=='0'运行结果:>> lx35ans =1%运行结果为1表示y2-y1*2+2*y=='0'成立6.设)ln(y x x u +=,求22x u ∂∂,22y u ∂∂,y x u ∂∂∂2. 程序:syms x y ;f=x*log(x+y);uxx=diff(f,x,2)uyy=diff(f,y,2)f1=diff(f,x);uxy=diff(f1,y)运行结果:>> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ∂∂∂∂,.(1)1cos cos cos 222=++z y x程序:syms x y z ;f=(cos(x))^2+(cos(y))^2+(cos(z))^2-1;zx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx371zx =-(cos(x)*sin(x))/(cos(z)*sin(z))zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z =程序:syms x y z ;f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx372f =exp(z) - x*y*zzx =(y*z)/(exp(z) - x*y)zy =(x*z)/(exp(z) - x*y)8.证明函数22)()(ln b y a x u -+-=(b a ,为常数)满足拉普拉斯方程: 02222=∂∂+∂∂y u x u (提示:对结果用simplify 化简)练习4 积分计算1.计算下列不定积分. (1)⎰+dx x x 12 (2)⎰+x xdx 2sin 12sin2.计算下列定积分.(1)⎰e xdx x 1ln (2)⎰ππ342sin dxx x3.求⎰+t dx x x x 12)ln (ln 1并用diff 对结果求导.4.求摆线)cos 1(),sin (t a y t t a x -=-=的一拱(π≤≤20t )与x 轴所围成的图形的面积.5.计算二重积分(1)⎰⎰≤++122)(y x dxdy y x (2)⎰⎰≤++x y x dxdy y x 22)(22 6.计算⎰+L ds y x 22 L 为圆周)0(22>=+a ax y x7.计算⎰++-L dy y x dx y x )()(2222,其中L 为抛物线2x y =上从点(0,0)到点(2,4)的一段弧.练习5 matlab 自定义函数与导数应用1.建立函数x x a a x f 3sin 31sin ),(+=,当a 为何值时,该函数在3π=x 处取得极值,它是极大值还是极小值,并求此极值.2.确定下列函数的单调区间.(1)7186223---=x x x y (2))0(82>+=x x x y3.求下列函数的最大值、最小值.(1)2332xx y -=41≤≤-x (2)312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1. 计算(1)⎪⎪⎪⎭⎫ ⎝⎛--521111204321+⎪⎪⎪⎭⎫ ⎝⎛--232002101041221 (2)⎪⎪⎭⎫ ⎝⎛-01301213⎪⎪⎪⎪⎪⎭⎫ ⎝⎛030101020501⎪⎪⎪⎭⎫ ⎝⎛-205101 (3)52422⎪⎪⎭⎫ ⎝⎛- 2.设⎪⎪⎪⎭⎫ ⎝⎛-=243121013A ,⎪⎪⎪⎭⎫ ⎝⎛-=112111201B ,求满足关系B X A =-23的X .练习7 矩阵与线性方程组1.求下列矩阵的秩.(1)⎪⎪⎪⎭⎫ ⎝⎛-321110021 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛4820322513454947513253947543173125 2.求下列矩阵的行列式,如可逆,试用不同的方法求其逆矩阵.(1)⎪⎪⎪⎭⎫ ⎝⎛--285421122 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---6201111121324321 3.设X ⎪⎪⎪⎭⎫ ⎝⎛-111012111=⎪⎪⎪⎭⎫ ⎝⎛521234311求X .4.解下列线性方程组.(1)⎪⎪⎩⎪⎪⎨⎧=--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)⎪⎪⎩⎪⎪⎨⎧-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分方程与级数求1-6题微分方程的通解1.1222+='y y y x 2.x y x y dx dy -+= 3.x x x y y +='cos 4.1)2sin cos (='+y y y x 5.x e y y y x 2cos 3=-'+'' 6.x x y y sin 14++=+''求7、8题初值问题的解7.⎪⎩⎪⎨⎧==-++-+=10)2(212222x y dx dy x xy y y xy x8.⎪⎪⎩⎪⎪⎨⎧===++==0000222,02V dt dxx x x a dt dx n dt x d t t9.给出函数x x e x f x x cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式.10.判别下列级数的敛散性,若收敛求其和. (1)Λ++++7151311 (2)∑∞=+112n n n tg π11.求幂级数∑∞=--22)1(n nn n n x 的和函数. 12.求函数项级数∑∞=-1)2sin )1(n n n n x π的和函数.。

Matlab考试题库及答案(教师出卷参考专用)

Matlab考试题库及答案(教师出卷参考专用)

Matlab考试题库及答案(教师出卷参考专用)一、选择题1.以下哪个函数用于在Matlab中创建一个图形窗口?A. figureB. plotC. graphD. window答案:A2.在Matlab中,以下哪个选项可以用来定义一个矩阵?A. A = [1 2 3; 4 5 6]B. A = (1, 2, 3, 4, 5, 6)C. A = {1, 2, 3, 4, 5, 6}D. A = 1 2 3; 4 5 6答案:A3.以下哪个函数用于求解线性方程组Ax=b?A. solveB. linsolveC. solve(A, b)D. linsolve(A, b)答案:D4.在Matlab中,如何计算矩阵A和矩阵B的乘积?A. A BB. A \ BC. A . BD. A .\ B答案:A5.以下哪个函数用于在Matlab中绘制三维散点图?A. scatterB. scatter3C. plot3D. bar3答案:B二、填空题1.在Matlab中,要创建一个名为"myfig"的图形窗口,可以使用______函数。

答案:figure('Name', 'myfig')2.在Matlab中,要计算矩阵A的行列式,可以使用______函数。

答案:det(A)3.在Matlab中,若要计算变量x的平方,可以使用______运算符。

答案:.^24.在Matlab中,若要计算矩阵A的逆矩阵,可以使用______函数。

答案:inv(A)5.在Matlab中,要绘制一个正弦波形,可以使用______函数。

答案:plot(sin(x))三、判断题1.在Matlab中,矩阵的索引从1开始计数。

()答案:正确2.在Matlab中,可以使用逻辑运算符"&&"和"||"。

()答案:错误3.在Matlab中,矩阵乘法满足交换律。

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 + 1.0000i ans = 0 + 1.0000i ans =2.2204e-016 ans =Inf ans = NaN ans =3.1416 ans =1.7977e+308 ans =2.2251e-3083.试写一函数 regPolygon(n),其功能为画出一个圆心在 (0, 0)、半径为 1 的圆,并在圆内画出一个内接正 n 边形,其中一顶点位于 (0, 1)。

例如 regPolygon(8) 可以画出如下之正八边型:解:新建regPolygon.m文件如下:function y=regPolyfon(n)n=8;%要画的n边形R=1; %圆的半径t=0:0.01: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 的脚本 plotParam.m,画出此曲线在 XY 平面的轨迹。

解:新建plotParam.m :t = linspace(0, 4*pi);x = sin(t);y = 1-cos(t)+t/10;plot(x, y, '-o'); -1-0.8-0.6-0.4-0.200.20.40.60.8100.511.522.535. 当一个小圆轮沿着一条曲线行进时,轮缘任一点的轨迹就会产生变化丰富的摆线。

MATLAB考试试题及答案

MATLAB考试试题及答案

MATLAB考试试题及答案一、选择题(每题5分,共25分)1. 在MATLAB中,下列哪个命令用于创建一个行向量?A. v = [1; 2; 3]B. v = [1 2 3]C. v = [1, 2, 3]D. v = (1, 2, 3)答案:B2. 在MATLAB中,下列哪个命令用于计算矩阵A的行列式?A. det(A)B. det(A')C. det(inv(A))D. det(A^2)答案:A3. 在MATLAB中,下列哪个命令用于计算矩阵A的逆?A. inv(A)B. A^(-1)C. pinv(A)D. A\B答案:A4. 在MATLAB中,下列哪个命令用于求解线性方程组Ax= b?A. A\bB. A/BC. B/AD. A^-1b答案:A5. 在MATLAB中,下列哪个命令用于绘制二维图形?A. plot(x, y)B. scatter(x, y)C. bar(x, y)D. pie(x, y)答案:A二、填空题(每题5分,共25分)6. 在MATLAB中,可以使用______命令创建一个等差数列。

答案:linspace7. 在MATLAB中,可以使用______命令创建一个等比数列。

答案:logspace8. 在MATLAB中,可以使用______命令计算矩阵A的特征值。

答案:eig(A)9. 在MATLAB中,可以使用______命令计算矩阵A的特征向量。

答案:eigenvector(A)10. 在MATLAB中,可以使用______命令计算矩阵A的奇异值。

答案:svd(A)三、解答题(每题25分,共75分)11. 编写MATLAB程序,求解以下线性方程组:2x + 3y - z = 1x - y + 2z = 03x + 2y - 4z = -3答案:```A = [2 3 -1; 1 -1 2; 3 2 -4];b = [1; 0; -3];x = A\b;disp('解为:');disp(x);```12. 编写MATLAB程序,绘制以下函数的图形:y = sin(x) + cos(x),x ∈ [0, 2π]答案:```x = linspace(0, 2pi, 100);y = sin(x) + cos(x);plot(x, y);title('y = sin(x) + cos(x)');xlabel('x');ylabel('y');grid on;```13. 编写MATLAB程序,计算以下矩阵的特征值和特征向量:A = [1 2 3; 4 5 6; 7 8 9]答案:```A = [1 2 3; 4 5 6; 7 8 9];[V, D] = eig(A);disp('特征值:');disp(diag(D));disp('特征向量:');disp(V);```14. 编写MATLAB程序,使用牛顿迭代法求解方程f(x) = x^3 - 4x + 2 = 0在x = 1附近的根。

Matlab与数学建模综合练习测验题

Matlab与数学建模综合练习测验题

1. 求032=-x e x 地所有根.(先画图后求解)>>syms x;>>f=exp(x)-3*x^2;>>ezplot(f)>>solve('exp(x)-3*x^2')ans =-2*lambertw(-1/6*3^(1/2))-2*lambertw(-1,-1/6*3^(1/2))-2*lambertw(1/6*3^(1/2))2. 求下列方程地根.1)0155=++x x>>solve('x^5+5*x+1')ans =1.1044655068824455162575638841973+1.0598296691525201166749456468980*i -1.0044974557968355184823910746206+1.0609465064060406435760940804509*i -.19993610217121999555034561915339-1.0044974557968355184823910746206-1.0609465064060406435760940804509*i 1.1044655068824455162575638841973-1.0598296691525201166749456468980*i 2) 至少三个根)(021s i n =-x x >>solve('x*sin(x)-1/2')ans =-.74084095509549062101093540994313>>f=x*sin(x)-1/2;>>ezplot(f)3) 所有根0c o s s i n2=-x x x>>f=sin(x)*cos(x)-x^2;>>solve('sin(x)*cos(x)-x^2')ans =0.>>ezplot(f)3. 求解下列各题:1)30sin lim x x x x ->- >>syms x;>>limit((x-sin(x))/x^3,x,0)ans =1/62) )10(,cos y x e y x 求=>>syms x;>>diff(exp(x)*cos(x),x,10)ans =-32*exp(x)*sin(x)3))位有效数字精确到⎰2/1017(2dx e x>>syms x;>>f=int(f,x,0,1/2);>>vpa(f,17)ans =.544987104183622204)⎰+dx xx 24425 >>syms x;>> int(x^4/(25+4*x^2),x)ans =1/12*x^3-25/16*x+125/32*atan(2/5*x)5) )(最高次幂为展开在将801=+x x>>syms x;>>f=sqrt(x+1);>>taylor(f,9)ans =1+1/2*x-1/8*x^2+1/16*x^3-5/128*x^4+7/256*x^5-21/1024*x^6+33/2048*x^7-429/32768*x^86) )2()3(1s i n y e y x 求=>>f=diff(exp(sin(1/x)),x,3)f =cos(1/x)/x^6*exp(sin(1/x))+6*sin(1/x)/x^5*exp(sin(1/x))+3*sin(1/x)/x^6*cos(1/x)*exp(sin(1/x))-6*cos(1/x)/x^4*exp(sin(1/x))-6*cos(1/x)^2/x^5*exp(sin(1/x))-cos(1/x )^3/x^6*exp(sin(1/x))>>limit(f,x,2)ans =-23/64*cos(1/2)*exp(sin(1/2))+3/64*exp(sin(1/2))*sin(1/2)*cos(1/2)-3/16*exp(sin(1/2))*cos(1/2)^2+3/16*exp(sin(1/2))*sin(1/2)-1/64*exp(sin(1/2))*cos(1/2)^34. 1)求矩阵⎪⎪⎪⎭⎫ ⎝⎛--=314020112A 地逆矩阵1-A 及特征值和特征向量. >> a=[-2 1 1;0 2 0;-4 1 3];>> inv(a)ans =-1.5000 0.5000 0.50000 0.5000 0-2.0000 0.5000 1.00002)求点(1,1,4)到直线l: (x-3)/-1 =y/0=(z+1)/2地距离.5. 已知,21)(222)(σμσπ--=x e x f 分别在下列条件下画出)(x f 地图形:);(在同一坐标系上作图,,=时=、);(在同一坐标系上作图,-,=时、421,0)2(110,1)1(σμμσ=、 6. 画 下列函数地图形:(1)202004cos sin ≤≤≤≤⎪⎪⎩⎪⎪⎨⎧===u t t z tu y t u x (1) >>t=[0:0.01:20];>>u=[0:0.01:2];>>[T,U]=meshgrid(t,u);>>x=U.*sin(T);>>y=U.*cos(T);>>mesh(x,y,z)(2) 30,30)sin(≤≤≤≤=y x xy z(2)>>x=[0:0.01:3];>>y=[0:0.01:3];>>[X,Y]=meshgrid(x,y);>>z=sin(X.*Y);>>mesh(x,y,z)(3)sin (3cos )02cos (3cos )02sin x t u t y t u u z u ππ=+⎧≤≤⎪=+⎨≤≤⎪=⎩.(3)>>t=[0:0.1:2*pi];>>u=[0:0.1:2*pi];>>[T,U]=meshgrid(t,u);>>x=sin(T*(3+cos(U)));>>y=cos(T*(3+cos(U)));>>z=sin(U);>>mesh(x,y,z) 7、 设⎪⎩⎪⎨⎧=+=+32/)7(11x x x x n n n ,数列}{n x 是否收敛?若收敛,其值为多少?精确到6位有效数字.8、设 ,131211p p p n nx ++++= }{n x 是否收敛?若收敛,其值为多少?精确到17位有效数字.注:学号为单号地取7=p ,学号为双号地取.8=p9、编程找出 5,1000+=≤b c c 地所有勾股数,并问:能否利用通项表示 },,{c b a ? 同时讨论c-b=2地情况.10、编程找出不定方程 )35000(122<-=-y Dy x 地所有正整数解.(学号为单号地取D=2, 学号为双号地取D=5)11、设 ⎩⎨⎧==+=--1,12121a a ma a a n n n , 编程计算.100a (学号为单号地取m=2, 学号为双号地取m=1)12、求下列极限:(1)1lim(39)x x xx →∞+; >>syms x;>>limit((3^x+9^x)^(1/x),x,inf)ans =9(2)2325(2)(3)lim (5)x x x x x x x +++→∞+++; >>syms x;>>limit(((x+2)^(x+2)*(x+3)^(x+3))/((x+5)^(2*x+5)),x,inf)ans =exp(-5)(3)00x y →→. 13、求下列函数地导数:(1)()y x =;>> syms x;>> diff(sqrt(x*sin(x)*sqrt(1-exp(x))))ans =1/2/(x*sin(x)*(1-exp(x))^(1/2))^(1/2)*(sin(x)*(1-exp(x))^(1/2)+x*cos(x)*(1-exp(x))^(1/2)-1/2*x*sin(x)/(1-exp(x))^(1/2)*exp(x))(2)()y t =14、已知参数方程ln cos cos sin x t y t t t=⎧⎨=-⎩,试求出dy dx 和223|t d y dx π=. 15、求积分:(1)0∞⎰; >> syms x;>> int(cos(x)/sqrt(x),x,0,inf)ans =1/2*2^(1/2)*pi^(1/2)(2)214011x dx x ++⎰. 16、求下列函数在0x =处地5阶泰勒展式:(1)1ln()1x x -+;(2)0sin x t dt t⎰. 17、求级数地和:(1)11116611(54)(51)n n ++++⨯⨯-⨯+;>> syms n;>> f=1/((5*n-4)*(5*n+1));>> s=symsum(f,n,1,inf)s =1/5(2)22111111()()()232323n n +++++++. >> syms n;>> f=1/2^n+1/3^n;>> s=symsum(f,n,1,inf)s =3/218、解微分方程:(1)222(1)2(1)0d x dx t t n n x dt dt --++=;(2)22222()0d x dx t t t n x dt dt++-=; 19、423dx x xy dt dy xy y dt⎧=-⎪⎪⎨⎪=-⎪⎩,初值为:(0)2,(0)3x y ==.版权申明本文部分内容,包括文字、图片、以及设计等在网上搜集整理.版权为个人所有This article includes some parts, including text, pictures, and design. Copyright is personal ownership.dvzfv。

matlab考试题B卷及答案

matlab考试题B卷及答案

matlab考试题B卷及答案1. MATLAB基础操作题:请在MATLAB中创建一个名为`myVector`的向量,包含元素1, 2, 3, 4, 5,并计算其元素之和。

答案:首先,创建向量`myVector`可以使用以下代码:```matlabmyVector = [1 2 3 4 5];```然后,计算元素之和可以使用`sum`函数:```matlabsumOfElements = sum(myVector);```最终,`sumOfElements`的值将为15。

2. MATLAB矩阵运算题:给定一个3x3的矩阵A,其元素为A=[1 2 3;4 5 6; 7 8 9],请计算矩阵A的转置。

答案:计算矩阵A的转置可以使用`transpose`函数或者简单的`.'`操作符。

以下是使用`transpose`函数的代码:```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = transpose(A);```或者使用`.'`操作符:```matlabA_transpose = A.';```得到的转置矩阵`A_transpose`将是:```1 4 72 5 83 6 9```3. MATLAB编程题:编写一个MATLAB函数,该函数接受一个整数n作为输入,并返回一个n阶的魔方矩阵。

答案:魔方矩阵是一个方阵,其中每一行、每一列以及两条对角线上的元素之和都相等。

以下是一个生成n阶魔方矩阵的MATLAB函数:```matlabfunction magicMatrix = createMagicSquare(n)if mod(n, 2) == 0error('n must be an odd number.');endmagicMatrix = zeros(n);num = 1;for i = 1:nfor j = 1:nif i == 1j = j;elsej = (j + 2 * (i - 1) - n) % n + 1;endmagicMatrix(i, j) = num;num = num + 1;endendend```调用此函数并传入一个奇数n,例如`createMagicSquare(3)`,将返回一个3阶魔方矩阵。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学实验与数学建模实验报告学院:专业班级:姓名:学号:完成时间:年月日承 诺 书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成,所有结果都通过上机验证,无转载或抄袭他人,也未经他人转载或抄袭。

若承诺不实,本人愿意承担一切责任。

承诺人:年 月 日数学实验学习体会(每个人必须要写字数1200字以上,占总成绩的20%)练习1 一元函数的图形1. 画出x y arcsin =的图象.2. 画出x y sec =在],0[π之间的图象.3. 在同一坐标系中画出x y =,2x y =,3x y =,3x y =,x y =的图象. 4. 画出3232)1()1()(x x x f ++-=的图象,并根据图象特点指出函数)(x f 的奇偶性. 5. 画出)2ln(1++=x y 及其反函数的图象. 6. 画出321+=x y 及其反函数的图象.练习2 函数极限1.计算下列函数的极限.(1)xxx4cos12sin1lim4-+π→.程序:sym x;f=(1+sin(2*x))/(1-cos(4*x));limit(f,x,pi/4)运行结果:lx21ans =1(2).程序:sym x;f=(1+cos(x))^(3*sec(x));limit(f,x,pi/2)运行结果:lx22ans =exp(3)(3)22)2(sinlnlimxxx-ππ→.程序:sym x;f=log(sin(x))/(pi-2*x)^2;limit(f,x,pi/2)运行结果:lx23ans =-1/8(4)212lim xxex→.程序:x xx sec32)cos1(lim+π→sym x ;f=x^2*exp(1/x);limit(f,x,0)limit(f,x,0,'right')limit(f,x,0,'left')运行结果:lx24ans =NaNans =Infans =%左极限为零,存在,右极限为无穷大,在x 趋近于零时函数没有极限(5))215(lim 122x x x x +-∞→.程序:sym x ;f=5*x^2/(1-x^2)+2^(1/x);limit(f,x,inf)运行结果:>> lx25ans =-4(6)x x x x x -+-→32112lim .程序:sym x ;f=(x^2-2*x+1)/(x^3-x);limit(f,x,1)运行结果:>> lx26ans =0 (7)x x x 11lim 20-+→.程序:sym x ;f=(sqrt(1+x^2)-1)/x;limit(f,x,0)运行结果:>> lx27ans =0 (8))3sin(cos 21lim 3π--π→x xx .程序:sym x ;f=(1-2*cos(x))/sin(x-pi/3);limit(f,x,pi/3)运行结果:>> lx28ans =3^(1/2) (9)tgxx x )1(lim 0+→.程序:sym x ;f=(1/x)^tan(x);limit(f,x,0,'right')运行结果:>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序:sym x ;f=(2/pi*atan(x))^x;limit(f,x,inf,'left')运行结果:>> lx210ans =Inf2. 解方程012=-⋅x x .程序:sym x ;X=solve(x*2^x-1)运行结果:>> lx202X =lambertw(0, log(2))/log(2)%方程有两个解3. 解方程1sin 3+=x x .程序:sym x ;X=solve(3*sin(x)+1-x)运行结果:>> lx203X =-0.53847936154. 解方程03=++q px x .(p 、q 为实数)程序:X=solve('x^3+p*x+q=0','x')运行结果:X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3))p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数. (1))11)(1(-+=x x y程序:sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1);diff(f)运行结果:>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) + 1)/(2*x^(3/2))(2)x x x y ln sin =程序:sym x ;f=x*sin(x)*log(x);diff(f)运行结果:>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数方程所确定的函数的导数.(1)⎩⎨⎧==t y t x 44程序:sym t ;f1=t^4;f2=4*t;diff(f2)/diff(f1)运行结果:>> lx321ans =1/t^3(2)⎩⎨⎧-=+=arctgt t y t x )1ln(2程序:sym t ;f1=log(1+t^2);f2=t-atan(t);diff(f2)/diff(f1)运行结果:>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t)3.求下列隐函数的导数. (1)22ln y x x y arctg +=程序:syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运行结果;>> lx331yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 + 1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2))(2)xy y x =程序:syms x y ;f=x^y-y^xyx=-diff(f,x)/diff(f,y)运行结果:>> lx332f =x^y - y^xyx =(x^(y - 1)*y - y^x*log(y))/(x*y^(x - 1) - x^y*log(x))4.设x e y x cos =,求)4(y .程序:sym x ;f=exp(x)*sin(x);diff(f,x,4)运行结果:>> lx34ans =(-4)*exp(x)*sin(x)5.验证x e y x sin =满足关系式:22=+'-''y y y程序: sym x ;f=exp(x)*sin(x);y2=diff(f,x,2);y1=diff(f,x,1);y=f;y2-y1*2+2*y=='0'运行结果:>> lx35ans =1%运行结果为1表示y2-y1*2+2*y=='0'成立6.设)ln(y x x u +=,求22x u ∂∂,22y u∂∂,y x u ∂∂∂2. 程序:syms x y ;f=x*log(x+y);uxx=diff(f,x,2)uyy=diff(f,y,2)f1=diff(f,x);uxy=diff(f1,y)运行结果:>> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ∂∂∂∂,.(1)1cos cos cos 222=++z y x程序:syms x y z ;f=(cos(x))^2+(cos(y))^2+(cos(z))^2-1;zx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx371zx =-(cos(x)*sin(x))/(cos(z)*sin(z))zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z =程序:syms x y z ;f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx372f =exp(z) - x*y*zzx =(y*z)/(exp(z) - x*y)zy =(x*z)/(exp(z) - x*y)8.证明函数22)()(ln b y a x u -+-=(b a ,为常数)满足拉普拉斯方程: 02222=∂∂+∂∂y u x u(提示:对结果用simplify 化简)练习4 积分计算1.计算下列不定积分. (1)⎰+dx x x 12 (2)⎰+x xdx 2sin 12sin2.计算下列定积分.(1)⎰e xdx x 1ln (2)⎰ππ342sin dx x x 3.求⎰+t dx x x x 12)ln (ln 1并用diff 对结果求导.4.求摆线)cos 1(),sin (t a y t t a x -=-=的一拱(π≤≤20t )与x 轴所围成的图形的面积.5.计算二重积分(1)⎰⎰≤++122)(y x dxdy y x (2)⎰⎰≤++x y x dxdy y x 22)(22 6.计算⎰+L ds y x 22 L 为圆周)0(22>=+a ax y x7.计算⎰++-L dy y x dx y x )()(2222,其中L 为抛物线2x y =上从点(0,0)到点(2,4)的一段弧.练习5 matlab 自定义函数与导数应用1.建立函数x x a a x f 3sin 31sin ),(+=,当a 为何值时,该函数在3π=x 处取得极值,它是极大值还是极小值,并求此极值.2.确定下列函数的单调区间.(1)7186223---=x x x y (2))0(82>+=x x x y3.求下列函数的最大值、最小值.(1)2332xx y -=41≤≤-x (2)312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1. 计算(1)⎪⎪⎪⎭⎫ ⎝⎛--521111204321+⎪⎪⎪⎭⎫ ⎝⎛--232002101041221 (2)⎪⎪⎭⎫ ⎝⎛-01301213⎪⎪⎪⎪⎪⎭⎫ ⎝⎛030101020501⎪⎪⎪⎭⎫ ⎝⎛-205101 (3)52422⎪⎪⎭⎫ ⎝⎛- 2.设⎪⎪⎪⎭⎫ ⎝⎛-=243121013A ,⎪⎪⎪⎭⎫ ⎝⎛-=112111201B ,求满足关系B X A =-23的X .练习7 矩阵与线性方程组1.求下列矩阵的秩.(1)⎪⎪⎪⎭⎫ ⎝⎛-321110021 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛4820322513454947513253947543173125 2.求下列矩阵的行列式,如可逆,试用不同的方法求其逆矩阵.(1)⎪⎪⎪⎭⎫ ⎝⎛--285421122 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---6201111121324321 3.设X ⎪⎪⎪⎭⎫ ⎝⎛-111012111=⎪⎪⎪⎭⎫ ⎝⎛521234311求X . 4.解下列线性方程组.(1)⎪⎪⎩⎪⎪⎨⎧=--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)⎪⎪⎩⎪⎪⎨⎧-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分方程与级数 求1-6题微分方程的通解1.1222+='y y y x 2.x y x y dx dy -+= 3.x x x y y +='cos 4.1)2sin cos (='+y y y x 5.x e y y y x 2cos 3=-'+'' 6.x x y y sin 14++=+'' 求7、8题初值问题的解7.⎪⎩⎪⎨⎧==-++-+=10)2(212222x y dx dy x xy y y xy x8.⎪⎪⎩⎪⎪⎨⎧===++==0000222,02V dt dxx x x a dt dx n dt x d t t9.给出函数x x e x f x x cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式.10.判别下列级数的敛散性,若收敛求其和.(1)++++7151311(2)∑∞=+112n nntgπ11.求幂级数∑∞=--22)1(nnnnnx的和函数.12.求函数项级数∑∞=-1)2sin)1(nnnn xπ的和函数.。

相关文档
最新文档