Matlab数学实验一2015(标准答案版)

合集下载

实验一 MATLAB软件应用复习含具体答案

，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证
。 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该9种商品的总收入和总利润。表1.1
货号
1
2
3. 用两种方法在同一个坐标下作出这四条曲线的图形，并要求用两种方法在图上加各种标注。 4．用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题， 1）概率曲线； 2）四叶玫瑰线； 3）叶形线
4）曳物线
。 5．作出下列曲面的3维图形， 1）； 2）环面
（2）新建一个M文件输入命令如下： v=0:0.01:2*pi; u=v; [U,V]=meshgrid(u,v); X=(1+cos(U)).*cos(V); Y=(1+cos(U)).*sin(V); Z=sin(U); mesh(X,Y,Z) 保存并运行，运行结果如下：
T6 新建一个M文件输入命令如下： function shuixianhuashu for m=100:999; a=fix(m/100); b=rem(fix(m/10),10); c=rem(m,10); if m==a.^3+b.^3+c.^3 disp(m) end end 保存，再在命令窗口输入shuixianhuashu结果如下： >> shuixianhuashu 153
探究实验自由发挥：自己提出问题，实验探索，广泛联想，发现规律，大胆猜想。比如函数cos(1/x)在 x=0附近的振荡现象，有无规律可寻？
T1
新建一个M文件输入命令如下： E=eye(3,3); R=rand(3,2); O=zeros(2,3); S=diag([3,7]); A=[E,R;O,S]; AA=A*A; B=[E,R+R*S;O,S*S]; if B==AA disp('命题成立'); else disp('命题不成立'); end 保存并运行，结果如下：命题成立由运行结果可知，得以验证。

数学实验(MATLAB)课后习题答案

数学实验练习2.1画出下列常见曲线的图形。

（其中a=1,b=2,c=3）1、立方抛物线3xy=解：x=-5:0.1:0;y=(-x).^(1/3);y=-y;x=0:0.1:5;y=[y,x.^(1/3)];x=[-5:0.1:0,0:0.1:5];plot(x,y)2、高斯曲线2x e=y-解：fplot('exp(-x.^2)',[-5,5])3、笛卡儿曲线)3(13,1333222axy y x t at y t at x =++=+=解：ezplot('x.^3+y.^3-3*x*y',[-5,5])xyx.3+y.3-3 x y = 0或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)4、蔓叶线)(1,1322322xa x y t at y t at x -=+=+=解：ezplot('y.^2-x.^3/(1-x)',[-5,5])xyy.2-x.3/(1-x) = 0或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi;x=t-sin(t); y=2*(1-cos(t)); plot(x,y)6、星形线)(sin ,cos 32323233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;plot(x,y)或ezplot('x.^(2/3)+y.^(2/3)-1',[-1,1])xyx.2/3+y.2/3-1 = 07、螺旋线ct z t b y t a x ===,sin ,cos 解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z) grid on8、阿基米德螺线θa r = 解：x =0:0.1:2*pi; r=x; polar(x,r)902701809、对数螺线θa e r = 解：x =0:0.1:2*pi; r=exp(x); polar(x,r)90270180010、双纽线))()((2cos 22222222y x a y x a r -=+=θ 解：x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-1,1]) grid onxy(x.2+y.2).2-(x.2-y.2) = 011、双纽线)2)((2sin 222222xy a y x a r =+=θ 解：x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-2*x*y',[-1,1]) grid onxy(x.2+y.2).2-2 x y = 012、心形线)cos 1(θ+=a r 解：x =0:0.1:2*pi; r=1+cos(x); polar(x,r)90270练习2.21、求出下列极限值。

MATLAB全部实验及答案

MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤1、命令窗口的简单使用（1）简单矩阵的输入（2）求[12+2×(7-4)]÷32的算术运算结果2、有关向量、矩阵或数组的一些运算（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,clca=[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);endc（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 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9（7）有一段指令如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量（8）使用三元组方法，创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一：clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a（9）写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B3、已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。

MATLAB数学实验答案（全）

MATLAB数学实验答案（全）第⼀次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作⼆维、三维⼏何图形，能够⽤Matlab 软件解决微积分、线性代数与解析⼏何中的计算问题。

补充命令vpa(x,n) 显⽰x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形在下⾯的题⽬中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//⾼阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项（要求将结果⽤向量的形式给出）。

（完整版）MATLAB）课后实验答案[1]

（完整版）MATLAB）课后实验答案[1]实验⼀ MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。

(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。

解：(1) 结果：(2). 建⽴⼀个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验⼆ MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证2 2E R RS A O S +??=。

解: M ⽂件如下；5. 下⾯是⼀个线性⽅程组：1231112340.951110.673450.52111456x x x ??=???(1) 求⽅程的解。

(2) 将⽅程右边向量元素b 3改为0.53再求解，并⽐较b 3的变化和解的相对变化。

(3) 计算系数矩阵A 的条件数并分析结论。

解： M ⽂件如下：123d4e56g9实验三选择结构程序设计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 ⽂件如下：2. 输⼊⼀个百分制成绩，要求输出成绩等级A、B、C、D、E。

MATLAB实验及答案详解

《MATLAB原理及应用》实验报告实验一MATLAB环境熟悉及基本操作一、实验目的熟悉MATLAB软件的工作环境和练习MATLAB命令窗口的基本操作。

二、实验内容1.命令窗口（Command Window）运行如入门【实验1-1】求23)]47(212[÷-⨯+的算术运算结果。

用键盘在MATLAB命令窗口中命令提示符“>>”后输入一下内容：(12+2*(7-4))/(3^2)在上述表达式输入完成后，按【Enter】键，该命令被执行。

在命令执行后，MATLAB命令窗口中将显示结果。

ans =2【实验1-2】简单矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=15524974111415412869811105132316A的输入步骤。

(1)用键盘在MATLAB命令窗口中命令提示符“>>”后输入一下内容：A=[16 3 2 13 -1 -4;5 10 11 8 7 9;9 6 8 12 -4 2;4 15 14 1 -5 15](2)按【Enter】键，该命令被执行。

(3)在命令执行后，MATLAB命令窗口中将显示结果.(4)采用分行输入，在命令窗口输入：A=[16 3 2 13 -1 -45 10 11 8 7 99 6 8 12 -4 24 15 14 1 -5 15]【实验1-3】用简短命令计算并绘制在0≤x≤6范围内的sin(Ax)、sinAx2、sin2(Ax)。

( A代表你的学号）在命令窗口输入：x=0:0.01:6;y=[sin(1*x);sin(1*x.^2);sin(1*x).^2];plot(x,y)(2) 按【Enter】键，该命令执行结果如图1图1 函数曲线图(3)在命令窗口输入：x=0:0.01:6;y1=sin(1*x);y2=sin(1*x.^2);y3=sin(1*x).^2; subplot(2,2,1),plot(x,y1),title('sin(1*x)') subplot(2,2,2),plot(x,y2),title('sin(1*x.^2)') subplot(2,2,3),plot(x,y3),title('sin(1*x) .^2') (4) 按【Enter】键，该命令执行结果如图2图2子图形式函数曲线图2.命令窗口（Command Window）常用的基本操作【实验1-4】命令行操作过程示例。

Matlab实验及答案

实验一、MATLAB基本操作一、实验目的2．学习使用图形函数计算器命令funtool及其环境。

3. 学习使用help命令进行帮助4. 掌握向量与矩阵的创建以及矩阵的基本操作5. 掌握数组与矩阵的概念二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；掌握数组与矩阵的概念；学会使用help命令进行帮助；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool；1.命令窗口的简单使用（1）简单矩阵的输入（自由创建）x=[1 3 5;2 4 6]x =1 3 52 4 6（2）求[12+2×(7-4)]÷32的算术运算结果，总结算术运算符先级[12+2*(7-4)]/3^2ans =22.有关向量、矩阵或数组的一些运算（1）设A=15；B=20；求C=A+B与c=a+b的区别A=15;B=20;>> C=A+BC =35>> c=a+bUndefined function or variable 'a'.（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是对应元素相乘)A=[1 2 3;4 5 6;7 8 9];B=[9 8 7;6 5 4;3 2 1];>> A*Bans =30 24 1884 69 54138 114 90>> A.*Bans =9 16 2124 25 2421 16 9（3）设a=10，b=20；求i=a/b与j=a\ba=10;>> b=20;>> i=a/bi =0.5000>> j=a\bj =2（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素的线性索引以及行列索引（sub2ind/ind2sub）。

重庆大学数学实验一 matlab的基本应用参考答案

《数学实验》第一次上机实验1．设有分块矩阵⎥⎦⎤⎢⎣⎡=⨯⨯⨯⨯22322333S O R E A ，其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证⎥⎦⎤⎢⎣⎡+=22S 0RS R E A 。

程序及结果：E=eye(3); %创建单位矩阵E% R=rand(3,2); %创建随机矩阵R% O=zeros(2,3); %创建0矩阵% S=diag(1:2); %创建对角矩阵% A=[E,R;O,S]; %创建A 矩阵%B=[E,(R+R*S);zeros(2,3),S^2] %计算等号右边的值%A^2 %计算等号左边的值%运行结果：B =1.00 0 0 1.632.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 ans =1.00 0 0 1.632.740 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.002．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。

表1.11)程序：a=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30]; b=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50]; c=[568 1205 753 580 395 2104 1538 810 694];s=sum((b-a).*c)i=b.*cmax((b-a).*c)min((b-a).*c)[m,n]=sort(b.*c)2）运行结果：s =4.6052e+004i =1.0e+004 *0.6305 1.8075 0.4518 0.9425 0.3911 3.8398 3.1990 1.95621.0757ans =1.3087e+004ans =1.2719e+003m =1.0e+004 *0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 1.9562 3.1990 3.8398n =5 3 1 4 9 2 8 7 63. 近景图将x的取值范围局限于较小的区间内可以画出函数的近景图,用于显示函数的局部特性。

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

Matｌab数学实验一——matlab初体验
一、实验目的及意义
［1] 熟悉MAＴＬＡB软件的用户环境；
［2] 了解MAＴLAB软件的一般目的命令；
［３] 掌握MATLＡB数组操作与运算函数；
通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。

二、实验内容
１.认识ｍatｌab的界面和基本操作
2.了解mａtlab的数据输出方式（ｆormat)
3. MATLAＢ软件的数组(矩阵）操作及运算练习;
三、实验任务
根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会)
完成如下题目,并按照实验报告格式和要求填写实验报告
1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。

iｊeps inf nａn pi realｍaｘrealｍin
2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。

3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。

4. 写出在命令窗口中的计算步骤和运行结果。

（1)计算
1.22
10
(ln log)
81
e
ππ
+-
；
>＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns =
０．０3４8
(2)
>> x=2;y＝４；
>> z=x^２+exｐ(x+y)-y*lｏg(x)-3
z ＝
40１.65６2
（3)输入变量
13
5.3,
25
a b
⎡⎤
==⎢⎥
⎣⎦
，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0
１.ｍat”文件。

测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3
ａ＝
５.３０00
＞> b=［1 3; 2 5]
b =
1 ３
２５
>> who
Your ｖａr ｉables a ｒe:
a b
>> ｗhos
Nam ｅＳi ｚe By ｔe ｓ Clas ｓ
ａ 1x1 8 dou ｂle ar ｒa ｙ
ｂ２x2 32 ｄo ｕｂle array
Ｇｒan ｄ t ｏｔal is 5 ｅlements usi ｎg ４０ bytes
＞＞ save D:＼exe01
＞> c ｌear 清除内存中在全部变量
>> ｌo ａｄ D:＼exe0１
5. 对矩阵，求其行列式(det ）、逆矩阵(inv ）、矩阵的特征值和特征向量(eig)、矩阵的秩（r ａn ｋ)、矩阵的行最简形（rref ）、以该矩阵为系数矩阵的线性方程组Ａx=0的通解(ｎｕｌl);
①已知422134305,203153211A B -⎛⎫⎛⎫ ⎪ ⎪=-=-- ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭
，在ＭA ＴＬA Ｂ命令窗口中建立A 、B 矩阵并对其进行以下操作：
(1) 计算矩阵A 的行列式的值det()A
>> Ａ=[4,-2，2;-3,0,5；1,5,3］;
>> d ｅt(A ）
ans =
-15８
(2) 分别计算下列各式:1122,*,.*,,,,T A B A B A B AB A B A A ---
>> A=［４,-2,２;－3,０,5;1,5,3];Ｂ=[1，3，4;-2,0，－3;2，-１,1];
＞> 2*Ａ-Ｂ
ans ＝
７ -7 ０
-４ 0 13
0 11 5
>> A*B。