MATLAB与数学实验艾冬梅版第三章习题答案
数学实验(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课后习题答案第三章
第3章数值数组及其运算习题3及解答1 要求在闭区间]2,0[ 上产生具有10个等距采样点的一维数组。
试用两种不同的指令实现。
〖目的〗●数值计算中产生自变量采样点的两个常用指令的异同。
〖解答〗%方法一t1=linspace(0,2*pi,10)%方法二t2=0:2*pi/9:2*pi %要注意采样间距的选择，如这里的2*pi/9.t1 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.2832t2 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.28322 由指令rng('default'),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。
〖目的〗●数组下标的不同描述：全下标和单下标。
●sub2ind, int2str, disp的使用。
●随机发生器的状态控制：保证随机数的可复现性。
〖解答〗rng('default')A=rand(3,5)[ri,cj]=find(A>0.5);id=sub2ind(size(A),ri,cj);ri=ri';cj=cj';disp(' ')disp('大于0.5的元素的全下标')disp(['行号 ',int2str(ri)])disp(['列号 ',int2str(cj)])disp(' ')disp('大于0.5的元素的单下标')disp(id')A =0.8147 0.9134 0.2785 0.9649 0.95720.9058 0.6324 0.5469 0.1576 0.48540.1270 0.0975 0.9575 0.9706 0.8003大于0.5的元素的全下标行号 1 2 1 2 2 3 1 3 1 3列号 1 1 2 2 3 3 4 4 5 5大于0.5的元素的单下标1 2 4 5 8 9 10 12 13 153 采用默认全局随机流，写出产生长度为1000的“等概率双位（即取-1，+1）取值的随机码”程序指令，并给出 -1码的数目。
matlab习题及答案3
>> ystd=std(y)
4.1 创建符号变量有几种方法？ MATLAB 提供了两种创建符号变量和表达式的函数：sym 和 syms。 sym 用于创建一个符号变量或表达式，用法如 x=sym(‘x’) 及 f=sym(‘x+y+z’)，syms 用于创建多个符号变量， 用法如 syms x y z。 f=sym(‘x+y+z’) 相当于 syms x y z f= x+y+z 4.2 下面三种表示方法有什么不同的含义？ （1）f=3*x^2+5*x+2 （2）f='3*x^2+5*x+2' （3）x=sym('x') f=3*x^2+5*x+2 （1）f=3*x^2+5*x+2 表示在给定 x 时，将 3*x^2+5*x+2 的数值运算结果赋值给变量 f，如果没有给定 x 则指示错误信息。 （2）f='3*x^2+5*x+2' 表示将字符串'3*x^2+5*x+2'赋值给字符变量 f，没有任何计算含义，因此也不对字符串中的内容做任何分析。 （3）x=sym('x') f=3*x^2+5*x+2 表示 x 是一个符号变量，因此算式 f=3*x^2+5*x+2 就具有了符号函数的意义，f 也自然成为符号变量了。 4.3 用符号函数法求解方程 at2+b*t+c=0。 >> r=solve('a*t^2+b*t+c=0','t')
5.1 绘制曲线 y x 3 x 1 ，x 的取值范围为[-5,5]。 >> x=-5:0.2:5; >> y=x.^3+x+1; >> plot(x,y)
matlab第三章课后部分答案
matlab第三章课后部分答案习题三3-2 从键盘输入一个三位整数，将它反向输出。
如输入639，输出为936程序如下：m=input('请输入一个三位整数：');m1=fix(m/100);%求m的百位整数m2=rem(fix(m/10),10);%求m的十位数字m3=rem(m,10);%求m的个位数字m=m3*100+m2*10+m1%反向输出m3-3 输入一个百分制成绩，要求输出成绩等级A,B,C,D,E。
其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。
要求：（1）分别用if语句和switch语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
程序如下：（1）if语句c=input('请输入成绩：');if c>=90&c<=100disp('A 成绩合理');elseif c>=80&c<=89disp('B 成绩合理');elseif c>=70&c<=79disp('C 成绩合理'); elseif c>=60&c<=69disp('D 成绩合理'); elseif c<60disp('E 成绩合理');elsedisp('成绩错误');end（2）switch语句c=input('请输入成绩：'); switch fix(c)case num2cell(90:100)disp('A 成绩合理'); case num2cell(80:89)disp('B 成绩合理'); case num2cell(70:79)disp('C 成绩合理'); case num2cell(60:69)disp('D 成绩合理'); case num2cell(0:59)disp('E 成绩合理');x=fix(rand(1,20)*89)+10;x1=fix(sum(x)/20);disp(['平均数是：',num2str(x1)])m=(rem(x,2)==0&x<x1);n=find(m);disp(['小于平均数的数是：',num2str(x(n))]); 3-6 输入20个数，求其中最大数和最小数。
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】命令行操作过程示例。
(完整版)第3章MATLAB矩阵处理_习题答案
3
B．y=cond(A) D．y=norm(A)
1．建立 3 阶幺矩阵的命令是
。ones(3)
2．设 A 为 2×3 矩阵，则用 zeros(size(A))建立的矩阵是
行
列的
矩阵。2，3，零
3．将 3 阶魔方矩阵主对角线元素加 10，命令是
。10*eye(3)+magic(3)
4．
可以用来描述矩阵的性能，它越接近于
D．B=diag(tril(A))
7．在 MATLAB 中定义 A=randn(5,4,3,2)，则下列关于 A 的操作中正确的是（
）。D
A．y=eig(A)
B．y=reshape(A,[4 3 6 7])
C．y=cond(A)
D．y=sin(A)
8．在命令行窗口中分别输入下列命令，对应输出结果正确的是（
7 8 9
10 11 12
，
B
10
7
4
1
11 8 5 2
12
3
9
，
C
6
6
9
3
12
2 5 8 11
1
4
7
10
第 1 题: (1)
a=1:12;
A=reshape(a,3,4)
(2)
B=rot90(A)
(3)
C=rot90(A,-1) 3．用矩阵求逆法求线性方程组的解。
4
x1
3x1
2x2 x3 x2 2x3
）。 C
A．命令 x=[-2:2]'，结果 x=[-2 -1 0 1 2]
B．命令 x=zeros(1,2);x>0，结果 ans=1
C．命令 y=diag(eye(3),1)&#,2)，结果 ans=[-5.0501 1.2311]
MATLAB数学实验练习题附答案PPT课件
%2、求100~999之间的水仙花数
clear all；
clc;
for n=100:999
n1=floor(n/100); %取出百位数字n1
ቤተ መጻሕፍቲ ባይዱ
n2=mod(floor(n/10),10); %取出十位数字n2
2
n3=mod(n,10) ; %取出个位数字n3
if n1^3+n2^3+n3^3==n
fprintf('%d是“水仙花数”\n', n)
让我们共同进步
2021/3/9
5
===============================================================
%5.作xoy面上的曲线(x-a).^2+y.^2=r.^2绕z轴旋转所得到的图形(0<r<a),取r=1,a=2.
clear all,clc;
x=1:0.1:3;
3
[x,y,z]=cylinder(sqrt(1-(x-2).^2));
surf(x,y,z)
2021/3/9
%6. 作x2/2-y2/3=2z图形. clear,clf,clc; [X,Y,Z]=meshgrid(-1:0.1:1); V=X.^2./2+Y.^2./2-2*Z; fv=isosurface(X,Y,Z,V,0);%显示V=0等值面 p=patch(fv); set(p,'FaceColor','blue','EdgeColor','none');%设置图形的相关属性 view(3) hold on camlight
end
202e1n/3d/9
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 limx 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 计算2211x 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项（要求将结果用向量的形式给出）。