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实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
“MATLAB”练习题要求:抄题、写出操作命令、运行结果，并根据要求，贴上运行图。
1、求230x e x -=的所有根.（先画图后求解）（要求贴图）>> solve('exp(x)—3*x^2',0)ans =—2*lambertw （—1/6＊3^（1/2））-2＊lambertw(—1,—1/6＊3^(1/2）)—2*lambertw （1/6*3^(1/2))3、求解下列各题：1）30sin lim x x x x ->->> sym x ；〉> limit((x-sin （x))/x^3)ans =1/62） (10)cos ,x y e x y =求>> sym x;>> diff （exp(x ）＊cos(x)，10）ans =（-32）＊exp(x)*sin （x)3）21/20(17x e dx ⎰精确到位有效数字）〉〉 sym x;〉〉 vpa((int(exp(x^2），x,0,1/2)），17)ans =0.544987104183622224）42254x dx x+⎰〉> sym x ；>〉 int （x^4/(25+x^2）,x)ans =125*atan （x/5) - 25＊x + x^3/35)求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。
〉> sym t;>> x=log(sqrt （1+t^2)）;y=atan(t)；〉> diff （y ，t ）/diff （x ，t)ans =1/t6)设函数y =f (x ）由方程xy +e y = e 所确定,求y ′(x ).>> syms x y ；f=x ＊y+exp(y ）—exp （1);〉> -diff(f,x ）/diff （f,y)ans =-y/（x + exp （y))7)0sin 2x e xdx +∞-⎰>〉 syms x ；>〉 y=exp(-x)＊sin(2*x ）；〉> int(y ，0,inf ）ans =2/58） 08x =展开（最高次幂为）〉> syms xf=sqrt(1+x);taylor(f,0，9)ans =— (429＊x^8）/32768 + (33*x^7）/2048 — (21＊x^6)/1024 +（7*x^5）/256 - (5＊x^4)/128 + x^3/16 - x^2/8 + x/2 + 19) 1sin (3)(2)x y e y =求〉> syms x y ；>〉 y=exp(sin （1/x)）;>〉 dy=subs （diff(y,3）,x ，2）dy =—0.582610)求变上限函数2x x ⎰对变量x 的导数.>> syms a t ；>〉 diff （int(sqrt(a+t),t,x ，x^2）)Warning: Explicit integral could not be found 。
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版韩明版)5.1,5.3,5.5,5.6部分答案
数学实验(M A T L A B版韩明版)5.1,5.3,5.5,5.6部分答案-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN练习5.1B的分布规律和分布函数的图形，通过观1、仿照本节的例子，分别画出二项分布()7.0,20察图形，进一步理解二项分布的性质。
解：分布规律编程作图：>> x=0:1:20;y=binopdf(x,20,0.7);>> plot(x,y,'*')图像：yx分布函数编程作图：>> x=0:0.01:20;>>y=binocdf(x,20,0.7)>> plot(x,y)图像：1x观察图像可知二项分布规律图像像一条抛物线，其分布函数图像呈阶梯状。
2、仿照本节的例子，分别画出正态分布()25,2N的概率密度函数和分布函数的图形，通过观察图形，进一步理解正态分布的性质。
解：概率密度函数编程作图：>> x=-10:0.01:10;>> y=normpdf(x,2,5);>> plot(x,y)图像：00.010.020.030.040.050.060.070.08x y分布函数编程作图：>> x=-10:0.01:10;>> y=normcdf(x,2,5);>> plot(x,y)图像：01x y观察图像可知正态分布概率密度函数图像像抛物线，起分布函数图像呈递增趋势。
3、设()1,0~N X ，通过分布函数的调用计算{}11<<-X P ,{}22<<-X P , {}33<<-X P .解：编程求解：>> x1=normcdf(1)-normcdf(-1),x2=normcdf(2)-normcdf(-2),x3=normcdf(3)-normcdf(-3) x1 = 0.6827x2 = 0.9545x3 = 0.9973即：{}6827.011=<<-X P ，{}9545.022=<<-X P ，{}9973.033=<<-X P .4、设()7.0,20~B X ,通过分布函数的调用计算{}10=X P 与{}10<X P .解：编程求解：>> x1=binopdf(10,20,0.7),x2=binocdf(10,20,0.7)-binopdf(10,20,0.7) x1 = 0.0308x2 = 0.0171即：{}0308.010==X P ，{}0171.010=<X P5、设()8~P X ,求：（1）{}4≤X P ;(2){}52≤<X P .解：（1）编程求解：>> p=poisscdf(4,8)p = 0.0996即：{}0996.04=≤X P（2）编程求解：>> p=poisscdf(5,8)-poisscdf(2,8)p = 0.1775即：{}1775.052=≤<X P6、（1）设()1,0~N X ,求01.0z ；（2）对2χ分布，求()8205.0χ；（3）对()1305.0t ；（4）对F 分布，求()10,1505.0F 。
MATLAB)课后实验答案
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e =+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=-- (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) t=0:0.5:2.5; z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)运算结果：z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044i2. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵）(2) A*B 和A.*B(3) A^3和A.^3(4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2]解： M 文件：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]运算结果：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1 ans =68 44 62309 -72 596154 -5 241 ans =12 102 468 0 2619 -130 49 ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000ans =109.4000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7 ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ (1) 求它们的乘积C 。
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
MATLAB数学实验100例题解
一元函数微分学实验1 一元函数的图形（基础实验）实验目的 通过图形加深对函数及其性质的认识与理解， 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法，建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧。
初等函数的图形2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势。
解：程序代码：>〉 x=linspace （0,2＊pi,600）; t=sin （x)。
/(cos （x ）+eps ）;plot(x ，t)；title （’tan （x ）')；axis (［0，2*pi ，-50,50]）; 图象:程序代码: 〉〉 x=linspace （0,2*pi,100); ct=cos （x)。
/（sin(x)+eps ）； plot(x,ct ）；title(’cot(x)')；axis (［0,2＊pi ，—50,50])； 图象：cot(x)4在区间]1,1[-画出函数xy 1sin =的图形。
解：程序代码:>> x=linspace （-1，1，10000）;y=sin(1。
/x ）; plot （x,y ）; axis （［-1，1，—2,2］） 图象：二维参数方程作图6画出参数方程⎩⎨⎧==t t t y tt t x 3cos sin )(5cos cos )(的图形：解：程序代码:>〉 t=linspace(0,2＊pi,100)； plot(cos(t ）.＊cos （5＊t ），sin(t ）。
*cos(3*t)); 图象：极坐标方程作图8 作出极坐标方程为10/t e r =的对数螺线的图形. 解：程序代码：〉〉 t=0:0.01:2*pi ； r=exp （t/10);polar(log(t+eps ）,log （r+eps)); 图象：90270分段函数作图10 作出符号函数x y sgn =的图形。
matlab数学实验考试题及答案
matlab数学实验考试题及答案一、选择题（每题2分，共10分）1. MATLAB中用于生成0到1之间均匀分布的随机数的函数是？A. randB. randiC. randnD. randperm答案：A2. 下列哪个命令可以计算矩阵的行列式？A. detB. rankC. eigD. inv答案：A3. MATLAB中用于求解线性方程组的命令是？A. solveB. linsolveC. fsolveD. ode45答案：A4. 在MATLAB中，如何创建一个3x3的单位矩阵？A. eye(3)B. ones(3)C. zeros(3)D. identity(3)答案：A5. MATLAB中用于绘制二维图形的函数是？A. plotB. surfC. meshD. contour答案：A二、填空题（每题3分，共15分）1. MATLAB中，使用________函数可以计算矩阵的迹。
答案：trace2. 若要在MATLAB中创建一个从1到10的向量，可以使用________函数。
答案：1:103. MATLAB中，使用________函数可以计算矩阵的特征值。
答案：eig4. 若要在MATLAB中绘制一个正弦波，可以使用________函数。
答案：sin5. MATLAB中，使用________函数可以计算矩阵的逆。
答案：inv三、简答题（每题10分，共20分）1. 描述MATLAB中如何使用循环结构来计算并打印1到100之间所有奇数的和。
答案：可以使用for循环结构，初始化一个变量sum为0，然后遍历1到100之间的每个数，使用模运算符判断是否为奇数，如果是，则将其加到sum上，最后打印sum的值。
2. 简述MATLAB中如何使用条件语句来检查一个数是否为素数，并打印出所有小于100的素数。
答案：可以使用for循环遍历2到99之间的每个数，对于每个数，使用一个while循环检查它是否有除1和它本身之外的因数，如果没有，则使用if语句判断该数是否为素数，如果是，则打印该数。
MATLAB全部实验及答案
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,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.^B5、已知+?=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。
数学实验(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实验题答案result5 =( 1 ) a = 1 : 2 : 5a =1 3 5( 2 ) b = [ a' , a' , a' ；a ]b =1 1 13 3 35 5 51 3 5( 3 ) c = a + b ( 2 , : )c =4 6 82、下列运算是否合法，为什么如合法，结果是多少>> result2=a*bError using *Inner matrix dimensions must agree. >> result3=a+b result3 =3 6 258 11>> result4=b*dresult4 =31 22 2240 49 1331 22 2240 49 13-5 -8 7>> result6=a.*b result6 =2 8 -3415 30>> result7=a./b result7 =>> result9=a.\b result9 =>> result10=a92result10 =1 4 916 25 36>> resultl 1=29aresult11 =2 4 816 32 64>>result5=[b;c']*d 3、⽤MATLAB求解下⾯的的⽅程组。
1、求以下变量的值，并在MATLAB^验证。
1 2 x13 2 x211 5 x32 13 x4>> 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>> 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)*B27 2 1 29 15 3 22 2 11 51 32 13(1)求矩阵A的秩(rank)(2)求矩阵 A 的⾏列式(determinant)(3)求矩阵 A 的逆(inverse)(4)求矩阵 A 的特征值及特征向量(eigenvalue and eigenvector)>> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3)n 10查看y 的值)m1=0;for m=-10:10 m仁m1+2^m;endm1m1 =6、求分段函数的值。
matlab数学实验答案第十章
1.>> f=[-3,4,-2,5];A=[1 1 3 -1;2 -3 1 -2];>> b=[14;-2];>> Aeq=[4 -1 2 -1];Beq=-2;lb=zeros(3,1);>> [x,feval]=linprog(f,A,b,Aeq,Beq,lb) Optimization terminated.x =0.00008.00000.0000-6.0000feval =2.00002.>> f=[5;4;8];A=[2 -1 0;5 3 0];b=[4;15];>> Aeq=[1 2 1];Beq=6;lb=zeros(3,1);>> [x,feval]=linprog(f,A,b,Aeq,Beq,lb) Optimization terminated.x =0.00003.00000.0000feval =12.00003.M文件>> A=[1 1 0 0 0 0 0 0 0 0 0;1 -0.06 1 1 0 0 0 0 0 0 0 ;-0.15 -0.06 1 1 1 1 1 1 0 0 0 ;-0.15 -0.06 -0.15 1 -0.06 1 1 -0.06 1 1 0 ;-0.15 -0.06 -0.15 1 -0.06 1 1 -0.06 1 -0.06 1;0 0 0 1 -0.06 -0.15 1 0 0 0 0;0 0 0 1 0 0 1 0 0 0 0];>> b=[10;10;10;10;10;3;4];>> f=[0;0;0;-1.40;0;0;-1.25;0;-1.14;0;-1.06];>> lb=[0;0;0;0;0;0;0;0;0;0;0];>> ub=[10;10;10;10;10;10;10;10;10;10;10];>> [x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb,ub) Optimization terminated.x =6.73163.26840.00003.46457.74140.00000.00000.00008.20580.00000.0000fval =-14.2049五年末拥有的资金金额总额最大为 14.20494.>> f=[-30;-40];A=[20 100;40 80;90 60];>> b=[200;200;360];>> [x,feval]=linprog(f,A,b,[],[],[0;0])Optimization terminated.x =3.50000.7500feval =-135.00005.>> f=[-7000;-12000];A=[9 5;4 5;3 10];b=[360;200;300]; >> [x,feval]=linprog(f,A,b)Optimization terminated.x =20.000024.0000feval =-4.2800e+0056.>> f=[10,5,6,4,8,12];>> A=[-1 -1 -1 0 0 0;0 0 0 -1 -1 -1; -1 0 0 -1 0 0 ; 0 -1 0 0 -1 0 ;0 0 -1 0 0 -1]; >> b=[-60;-100;-50;-70;-40];>> Aeq=[];>> beq=[]; vlb=zeros(6,1);vub=[];>> [x,fval]=linprog(f,A,b,Aeq,beq,vlb,vub)Optimization terminated.x =0.000020.000040.000050.000050.00000.0000fval =940.00007.>> f=[8400;7500;6700];>> A=[60 65 80;-3800 -3500 -2500];>> b=[70 ;-2900];>> Aeq=[1 1 1];>> beq=1;>> lb=[0 0 0];>> ub=[inf inf inf];>> [x,fval]=linprog(f,A,b,Aeq,beq,lb,ub)Optimization terminated.x =0.00000.66670.3333fval =7.2333e+0038.>> f=[-1;-2;1;-1;-4;2];>> A=[1 1 1 1 1 1;2 1 -2 1 0 0;0 0 1 1 2 1];>> b=[6;4;4];>> lb=zeros(6,1);>> [x,fval]=linprog(f,A,b,[],[],lb)Optimization terminated.x =0.00004.00000.00000.00002.00000.0000fval =-16.000010.建立M函数fun44.mfunction f=fun44(x)f=-(sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4))); 建立M函数myconl.mfunction [g,ceq]=mycon1(x)g(1)=x(1)-400;g(2)=1.1*x(1)+x(2)-440;g(3)=1.21*x(1)+1.1*x(2)+x(3)-484;g(4)=1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4;ceq=0command窗口>> x0=[1;1;1;1];>> vlb=[0;0;0;0];vub=[];>> A=[];b=[];Aeq=[];beq=[];>>[x,fval]=fmincon('fun44',x0,A,b,Aeq,beq,vlb,vub,'mycon1') 12.function f=fmb(x)f=-x(1)*x(2)*x(3);>> A=[1 -2 -2;1 2 2];b=[0;72];>> lb=[-inf;10;-inf];ub=[inf;20;inf];>> Aeq=[1 -1 0];Beq=10;>> [x,feval]=fmincon(@fmb,[20,10,0],A,b,Aeq,Beq,lb,ub,[])x =22.5850 12.5850 12.1225feval =-3.4456e+003。
MATLAB)课后实验答案[1](最新整理)
](https://img.taocdn.com/s3/m/2bd1298ccf84b9d529ea7a94.png)
5. 下面是一个线性方程组：
1112340.95x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥
实验三 选择结构程序设计
1. 求分段函数的值。
222603560523
1x 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。
其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。
要求：
(1) 分别用if语句和switch语句实现。
(2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
解：M文件如下
3. 硅谷公司员工的工资计算方法如下：
(1) 工作时数超过120小时者，超过部分加发15%。
(2) 工作时数低于60小时者，扣发700元。
(3) 其余按每小时84元计发。
试编程按输入的工号和该号员工的工时数，计算应发工资。
解：M文件下
运算结果如下；
5. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。
例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。
求[2,50]区间内：
(1) 亲密数对的对数。
(2) 与上述亲密数对对应的所有亲密素数之和。
解：
M文件：。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/d27ebaa958fb770bf78a55b8.png)
实验一 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 a z 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 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E 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 值。
MATLAB数学实验课后答案
数学实验MATLAB参考答案（重要部分）P20,ex1(5) 等于[exp(1),exp(2);exp(3),exp(4)](7) 3=1*3, 8=2*4(8) a为各列最小值，b为最小值所在的行号(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)P20, ex2(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码P20,ex3>> r=2;p=0.5;n=12;>> T=log(r)/n/log(1+0.01*p)T =11.5813P20,ex4>> x=-2:0.05:2;f=x.^4-2.^x;>> [fmin,min_index]=min(f)fmin =-1.3907 %最小值min_index =54 %最小值点编址>> x(min_index)ans =0.6500 %最小值点>> [f1,x1_index]=min(abs(f)) %求近似根--绝对值最小的点f1 =0.0328x1_index =24>> x(x1_index)ans =-0.8500>> x(x1_index)=[];f=x.^4-2.^x; %删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) %求另一近似根--函数绝对值次小的点f2 =0.0630x2_index =65>> x(x2_index)ans =1.2500P20,ex5>> z=magic(10)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6679 6 13 95 97 29 31 38 45 7210 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59>> sum(z)ans =505 505 505 505 505 505 505 505 505 505 >> sum(diag(z))ans =505>> z(:,2)/sqrt(3)ans =57.157746.188046.765450.229553.693613.85642.88683.46416.928210.3923>> z(8,:)=z(8,:)+z(3,:)z =92 99 1 8 15 67 74 51 58 40 98 80 7 14 16 73 55 57 64 41 4 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 34 17 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6683 87 101 115 119 83 87 101 115 11910 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59P 40 ex1先在编辑器窗口写下列M函数，保存为eg2_1.m function [xbar,s]=ex2_1(x)n=length(x);xbar=sum(x)/n;s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));例如>>x=[81 70 65 51 76 66 90 87 61 77];>>[xbar,s]=ex2_1(x)xbar =72.4000s =12.1124P 40 ex2s=log(1);n=0;while s<=100n=n+1;s=s+log(1+n);endm=n计算结果m=37P 40 ex3clear;F(1)=1;F(2)=1;k=2;x=0;e=1e-8; a=(1+sqrt(5))/2;while abs(x-a)>ek=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); enda,x,k计算至k=21可满足精度P 40 ex4clear;tic;s=0;for i=1:1000000s=s+sqrt(3)/2^i;ends,toctic;s=0;i=1;while i<=1000000s=s+sqrt(3)/2^i;i=i+1;ends,toctic;s=0;i=1:1000000;s=sqrt(3)*sum(1./2.^i);s,tocP 40 ex5t=0:24;c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...31 32 31 29 27 25 24 22 20 18 17 16];plot(t,c)P 40 ex6(1)clear;fplot('x^2*sin(x^2-x-2)',[-2,2])x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y)y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2]) (2)参数方法t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t); plot(x,y)(3)x=-3:0.1:3;y=x;[x,y]=meshgrid(x,y);z=x.^2+y.^2;surf(x,y,z)(4)x=-3:0.1:3;y=-3:0.1:13;[x,y]=meshgrid(x,y);z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;surf(x,y,z)(5)t=0:0.01:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z)(6)theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); [theta,fai]=meshgrid(theta,fai);x=2*sin(fai).*cos(theta);y=2*sin(fai).*sin(theta);z=2*cos(fai);surf(x,y,z)(7)x=linspace(0,pi,100);y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)page41, ex7x=-1.5:0.05:1.5;y=1.1*(x>1.1)+x.*(x<=1.1).*(x>=-1.1)-1.1*(x<-1.1);plot(x,y)page41,ex8分别使用which trapz, type trapz, dir C:\MATLAB7\toolbox\matlab\datafun\page41,ex9clear;close;x=-2:0.1:2;y=x;[x,y]=meshgrid(x,y);a=0.5457;b=0.7575;p=a*exp(-0.75*y.^2-3.75*x.^2-1.5*x).*(x+y>1);p=p+b*exp(-y.^2-6*x.^2).*(x+y>-1).*(x+y<=1);p=p+a*exp(-0.75*y.^2-3.75*x.^2+1.5*x).*(x+y<=-1);mesh(x,y,p)page41, ex10lookfor lyapunovhelp lyap>> A=[1 2 3;4 5 6;7 8 0];C=[2 -5 -22;-5 -24 -56;-22 -56 -16]; >> X=lyap(A,C)X =1.0000 -1.0000 -0.0000-1.0000 2.0000 1.0000-0.0000 1.0000 7.0000Chapter 3%Exercise 1>> a=[1,2,3];b=[2,4,3];a./b,a.\b,a/b,a\bans =0.5000 0.5000 1.0000ans =2 2 1ans =0.6552 %一元方程组x[2,4,3]=[1,2,3]的近似解ans =0 0 00 0 00.6667 1.3333 1.0000%矩阵方程[1,2,3][x11,x12,x13;x21,x22,x23;x31,x32,x33]=[2,4,3]的特解Exercise 2(1)>> A=[4 1 -1;3 2 -6;1 -5 3];b=[9;-2;1];>> rank(A), rank([A,b]) %[A,b]为增广矩阵ans =3ans =3 %可见方程组唯一解>> x=A\bx =2.38301.48942.0213Exercise 2(2)>> A=[4 -3 3;3 2 -6;1 -5 3];b=[-1;-2;1]; >> rank(A), rank([A,b])ans =3ans =3 %可见方程组唯一解>> x=A\bx =-0.4706-0.2941Exercise 2(3)>> A=[4 1;3 2;1 -5];b=[1;1;1];>> rank(A), rank([A,b])ans =2ans =3 %可见方程组无解>> x=A\bx =0.3311-0.1219 %最小二乘近似解Exercise 2(4)>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1 2 3]';%注意b的写法>> rank(a),rank([a,b])ans =3ans =3 %rank(a)==rank([a,b])<4说明有无穷多解>> a\bans =110 %一个特解Exercise 3>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1,2,3]'; >> x=null(a),x0=a\bx =-0.62550.6255-0.20850.4170x0 =11%通解kx+x0Exercise 4>> x0=[0.2 0.8]';a=[0.99 0.05;0.01 0.95]; >> x1=a*x, x2=a^2*x, x10=a^10*x >> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> x0=[0.8 0.2]';>> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> [v,e]=eig(a)v =0.9806 -0.70710.1961 0.7071e =1.0000 00 0.9400>> v(:,1)./xans =1.17671.1767 %成比例，说明x是最大特征值对应的特征向量Exercise 5%用到公式(3.11)(3.12)>> B=[6,2,1;2.25,1,0.2;3,0.2,1.8];x=[25 5 20]';>> C=B/diag(x)C =0.2400 0.4000 0.05000.0900 0.2000 0.01000.1200 0.0400 0.0900>> A=eye(3,3)-CA =0.7600 -0.4000 -0.0500-0.0900 0.8000 -0.0100-0.1200 -0.0400 0.9100>> D=[17 17 17]';x=A\Dx =37.569625.786224.7690%Exercise 6(1)>> a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a) ans =-94ans =0.2553 -0.0213 0.04260.1596 -0.1383 -0.22340.1809 -0.2234 -0.0532v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766%Exercise 6(2)>> a=[1 1 -1;0 2 -1;-1 2 0];det(a),inv(a),[v,d]=eig(a) ans =1ans =2.0000 -2.0000 1.00001.0000 -1.0000 1.00002.0000 -3.0000 2.0000v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i -0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i%Exercise 6(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> det(A),inv(A), [v,d]=eig(A)ans =1ans =68.0000 -41.0000 -17.0000 10.0000 -41.0000 25.0000 10.0000 -6.0000 -17.0000 10.0000 5.0000 -3.0000 10.0000 -6.0000 -3.0000 2.0000v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887%Exercise 6(4)、(以n=5为例)%关键是矩阵的定义%方法一（三个for）n=5;for i=1:n, a(i,i)=5;endfor i=1:(n-1),a(i,i+1)=6;endfor i=1:(n-1),a(i+1,i)=1;enda%方法二（一个for）n=5;a=zeros(n,n);a(1,1:2)=[5 6];for i=2:(n-1),a(i,[i-1,i,i+1])=[1 5 6];enda(n,[n-1 n])=[1 5];a%方法三（不用for）n=5;a=diag(5*ones(n,1));b=diag(6*ones(n-1,1));c=diag(ones(n-1,1));a=a+[zeros(n-1,1),b;zeros(1,n)]+[zeros(1,n);c,zeros(n-1,1)] %下列计算>> det(a)ans =665>> inv(a)ans =0.3173 -0.5865 1.0286 -1.6241 1.9489-0.0977 0.4887 -0.8571 1.3534 -1.62410.0286 -0.1429 0.5429 -0.8571 1.0286 -0.0075 0.0376 -0.1429 0.4887 -0.5865 0.0015 -0.0075 0.0286 -0.0977 0.3173 >> [v,d]=eig(a)v =-0.7843 -0.7843 -0.9237 0.9860 -0.9237 0.5546 -0.5546 -0.3771 -0.0000 0.3771 -0.2614 -0.2614 0.0000 -0.1643 0.0000 0.0924 -0.0924 0.0628 -0.0000 -0.0628 -0.0218 -0.0218 0.0257 0.0274 0.0257d =0.7574 0 0 0 00 9.2426 0 0 00 0 7.4495 0 00 0 0 5.0000 00 0 0 0 2.5505%Exercise 7(1)>> a=[4 1 -1;3 2 -6;1 -5 3];[v,d]=eig(a) v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766>> det(v)ans =-0.9255 %v行列式正常, 特征向量线性相关，可对角化>> inv(v)*a*v %验算ans =-3.0527 0.0000 -0.00000.0000 3.6760 -0.0000-0.0000 -0.0000 8.3766>> [v2,d2]=jordan(a) %也可用jordanv2 =0.0798 0.0076 0.91270.1886 -0.3141 0.1256-0.1605 -0.2607 0.4213 %特征向量不同d2 =8.3766 0 00 -3.0527 - 0.0000i 00 0 3.6760 + 0.0000i>> v2\a*v2ans =8.3766 0 0.00000.0000 -3.0527 0.00000.0000 0.0000 3.6760>> v(:,1)./v2(:,2) %对应相同特征值的特征向量成比例ans =2.44912.44912.4491%Exercise 7(2)>> a=[1 1 -1;0 2 -1;-1 2 0];[v,d]=eig(a)v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i>> det(v)ans =-5.0566e-028 -5.1918e-017i %v的行列式接近0, 特征向量线性相关，不可对角化>> [v,d]=jordan(a)v =1 0 11 0 01 -1 0d =1 1 00 1 10 0 1 %jordan标准形不是对角的，所以不可对角化%Exercise 7(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> [v,d]=eig(A)v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887>> inv(v)*A*vans =0.0102 0.0000 -0.0000 0.00000.0000 0.8431 -0.0000 -0.0000-0.0000 0.0000 3.8581 -0.0000-0.0000 -0.0000 0 30.2887%本题用jordan不行, 原因未知%Exercise 7(4)参考6(4)和7(1), 略%Exercise 8 只有(3)对称, 且特征值全部大于零, 所以是正定矩阵. %Exercise 9(1)>> a=[4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0]>> rank(a)ans =3>> rank(a(1:3,:))ans =2>> rank(a([1 2 4],:)) %1,2,4行为最大无关组ans =3>> b=a([1 2 4],:)';c=a([3 5],:)';>> b\c %线性表示的系数ans =0.5000 5.0000-0.5000 1.00000 -5.0000%Exercise 10>> a=[1 -2 2;-2 -2 4;2 4 -2]>> [v,d]=eig(a)v =0.3333 0.9339 -0.12930.6667 -0.3304 -0.6681-0.6667 0.1365 -0.7327d =-7.0000 0 00 2.0000 00 0 2.0000>> v'*vans =1.0000 0.0000 0.00000.0000 1.0000 00.0000 0 1.0000 %v确实是正交矩阵%Exercise 11%设经过6个电阻的电流分别为i1, ..., i6. 列方程组如下%20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0; %i1=i3+i4;i5=i2+i3;i6=i4+i5;%计算如下>> A=[1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0 0;1 0 -1 0 0 -3 0 0 0;1 -1 0 0 0 0 -4 0 0;0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3;0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;0 0 0 0 0 0 -1 -1 1];>>b=[20 5 0 0 0 0 0 0 0]'; A\b ans =13.34536.44018.54203.3274-1.18071.60111.72630.42042.1467>> A=[1 2 3;4 5 6;7 8 0];>> left=sum(eig(A)), right=sum(trace(A))left =6.0000right =6>> left=prod(eig(A)), right=det(A) %原题有错, (-1)^n应删去left =27.0000right =27>> fA=(A-p(1)*eye(3,3))*(A-p(2)*eye(3,3))*(A-p(3)*eye(3,3)) fA =1.0e-012 *0.0853 0.1421 0.02840.1421 0.1421 0-0.0568 -0.1137 0.1705>> norm(fA) %f(A)范数接近0ans =2.9536e-013roots([1 1 1])%Exercise 1(2)roots([3 0 -4 0 2 -1])%Exercise 1(3)p=zeros(1,24);p([1 17 18 22])=[5 -6 8 -5];roots(p)%Exercise 1(4)p1=[2 3];p2=conv(p1, p1);p3=conv(p1, p2);p3(end)=p3(end)-4; %原p3最后一个分量-4roots(p3)%Exercise 2fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x'); fzero(fun,2)】%Exercise 3fun=inline('x^4-2^x');fplot(fun,[-2 2]);grid on;fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5)%Exercise 4fun=inline('x*sin(1/x)','x');fplot(fun, [-0.1 0.1]);x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=[x,-x]%Exercise 5fun=inline('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);1 6*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2]','x');[a,b,c]=fsolve(fun,[0 0 0])%Exercise 6fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))]; [a,b,c]=fsolve(fun,[0.5 0.5])%Exercise 7clear; close; t=0:pi/100:2*pi;x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t);x2=3+sqrt(2)*cos(t); y2=6*sin(t);plot(x1,y1,x2,y2); grid on; %作图发现4个解的大致位置，然后分别求解y1=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[ 1.5,2])y2=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[ 1.8,-2])y3=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[ 3.5,-5])y4=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[ 4,-4])%Exercise 8(1)clear;fun=inline('x.^2.*sin(x.^2-x-2)');fplot(fun,[-2 2]);grid on; %作图观察x(1)=-2;x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2);fun2=inline('-x.^2.*sin(x.^2-x-2)');x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)%答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小，x(2)最大。
数学实验（MATLAB版韩明版）21-24部分答案
练习2.1画出下列常见曲线的图形（其中a=1,b=2,c=3)1.的图像：31x y =()55≤≤-x 编程：>> x=-5:0.1:5;>> y=x.^(1/3);>> plot(x,y)>> grid on>> xlabel('x');ylabel('y')>> legend('y=x.^(1/3)')>> title('y=x.^(1/3)')图像：-5-4-3-2-101234500.20.40.60.811.21.41.61.8xyy=x.(1/3)y=x.(1/3)2.的图像：e x y -=2()55≤≤-x 编程：>> x=-5:0.1:5;>> y=exp(-x.^2);>> plot(x,y)>> grid on>> xlabel('x');ylabel('y')>> legend('y=e^(-x^2)');>> title('y=e^(-x^2)')图像：-5-4-3-2-101234500.10.20.30.40.50.60.70.80.91xyy=e (-x 2)y=e (-x 2)3.的图像：，a=1⎪⎭⎫ ⎝⎛=++=+=axy a y atx yx t t t 313,1333222()55≤≤-x 编程：>> t=-5:0.1:5;>> x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2);>> plot(x,y)>> grid on>> xlabel('x');ylabel('y')>> title('x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2)')>> legend('x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2)')图像：-1.5-1-0.500.51 1.500.511.522.53xyx=3*t./(1+t.2);y=3*t.2./(1+t.2)x=3*t./(1+t.2);y=3*t.2./(1+t.2)4.的图像: ,a=1⎪⎪⎭⎫ ⎝⎛-=+=+=x a a y a x x y t t t t3223221,1()55≤≤-t 编程：>> t=-5:0.1:5;>> x=t.^2./(1+t.^2);y=t.^3./(1+t.^2);>> plot(x,y)>> xlabel('x');ylabel('y')>> title('x=t.^2./(1+t.^2);y=t.^3./(1+t.^2)')>> legend('x=t.^2./(1+t.^2);y=t.^3./(1+t.^2)')>> grid on 图像：0.10.20.30.40.50.60.70.80.91-5-4-3-2-1012345xyx=t.2./(1+t.2);y=t.3./(1+t.2)x=t.2./(1+t.2);y=t.3./(1+t.2)5.的图像：,a=1,b=2()()t b y t t a x cos 1,sin -=-=pi t pi *2*2≤≤-编程：>> t=-2*pi:0.1:2*pi;>> x=t-sin(t);y=2*(1-cos(t));>> plot(x,y)>> grid on>> xlabel('x');ylabel('y')>> legend('x=t-sin(t);y=2*(1-cos(t))')>> title('x=t-sin(t);y=2*(1-cos(t))')图像：-8-6-4-20246800.511.522.533.54xyx=t-sin(t);y=2*(1-cos(t))x=t-sin(t);y=2*(1-cos(t))6.的图像：,a=1⎪⎪⎭⎫ ⎝⎛=+==a yx t a y t a x 32323233sincos ,pi t pi *2*2≤≤-编程：>> t=-2*pi:0.1:2*pi;>> x=(cos(t)).^3;y=(sin(t)).^3;>> plot(x,y)>> grid on>> xlabel('x');ylabel('y')>> title('x=(cos(t)).^3;y=(sin(t)).^3')图像：-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81xyx=(cos(t)).3;y=(sin(t)).3x=(cos(t)).3;y=(sin(t)).37.的图像：ct z t b y t a x ===,sin ,cos ()pi t pi c b a *2*2,3,2,1≤≤-===编程：>> t=-2*pi:0.1:2*pi;>> x=cos(t);y=2*sin(t);z=3*t;>> plot3(x,y,z)>> xlabel('x');ylabel('y');zlabel('z')>> grid on>> legend('x=cos(t);y=2*sin(t);z=3*t')>> title('x=cos(t);y=2*sin(t);z=3*t')图像：-101-2-1012-20-101020xx=cos(t);y=2*sin(t);z=3*tyzx=cos(t);y=2*sin(t);z=3*t8.的图像：θa r =()pi a *20,1≤≤=θ编程：>> theta=0.0:0.1:2*pi;>> r=theta;>> polar(theta,r)>> grid on>> legend('r=theta')>> title('r=theta')图像：24 68302106024090270120300150330180r=theta r=theta9.的图像：e a r θ=()pi a *20,1≤≤=θ编程：>> theta=-2*pi:0.1:2*pi;>> r=exp(theta);>> polar(theta,r)>> grid on>> title('r=exp(theta)')>> legend('r=exp(theta)')图像：100200 300400 5003021060240902701203001503301800r=exp(theta)r=exp(theta)10.的图像：()⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-==+y x a y xar 22222222,2cos θ1=a 编程：>> theta=0:0.1:2*pi;>> r=sqrt(abs(cos(2*theta)));>> polar(theta,r)>> grid on>> title('r=sqrt(abs(cos(2*theta)))');>> legend('r=sqrt(abs(cos(2*theta)))')图像：0.20.4 0.60.8 13021060240902701203001503301800r=sqrt(abs(cos(2*theta)))r=sqrt(abs(cos(2*theta)))11.的图像：a=1()⎪⎭⎫⎝⎛==+xy a y xar 2222*222,2sin θ编程：>> theta=0:0.1:2*pi;>> r=sqrt(abs(sin(2*theta)));>> polar(theta,r)>> grid on>> title('r=sqrt(abs(sin(2*theta)))')>> legend('r=sqrt(abs(sin(2*theta)))')图像：0.20.4 0.60.8 13021060240902701203001503301800r=sqrt(abs(sin(2*theta)))r=sqrt(abs(sin(2*theta)))12.的图像：a=1)cos 1(θ+=a r 编程：>> theta=0:0.1:2*pi;>> r=1+cos(theta);>> polar(theta,r)>> grid on>> legend('r=1+cos(theta)')>> title('r=1+cos(theta)')图像：0.51 1.52302106024090270120300150330180r=1+cos(theta)r=1+cos(theta)练习2.21.求出下列极限值.（1）； （2）；nnn n3lim 3+∞→()n n n n ++-+∞→122lim（3）；（4）；x x x 2cot lim 0→⎪⎭⎫ ⎝⎛→x m xx coslim 0（5）； （6）.⎪⎪⎭⎫ ⎝⎛--→111lim 1e x x x ⎪⎭⎫ ⎝⎛-+∞→x x xx 2lim 解：（1）编程：>> syms n >> limit((n^3+3^n)^(1/n),n,inf) ans = 3（2）编程：>> syms n>> limit(sqrt(n+2)-2*sqrt(n+1)+sqrt(n),n,inf) ans =0（3）编程：>> syms x >> limit(x*cot(2*x),x,0) ans = 1/2（4）编程：>> syms x m >> limit((cos(m/x))^x,x,inf) ans =1（5）编程：>> syms x>> limit(1/x-1/(exp(x)-1),x,1) ans = (exp(1)-2)/(exp(1)-1)（6）编程：>> syms x >> limit(sqrt(x^2+x)-x,x,1) ans = 2^(1/2)-12.有个客户看中某套面积为180，每平方米7500元的房子。
数学实验（MATLAB版韩明版）2.1-2.4部分答案
数学实验（MATLAB版韩明版）2.1-2.4部分答案练习2.1画出下列常见曲线的图形（其中a=1,b=2,c=3)1.31xy =的图像：()55≤≤-x编程：>> x=-5:0.1:5; >> y=x.^(1/3); >> plot(x,y) >> grid on>> xlabel('x');ylabel('y') >> legend('y=x.^(1/3)') >> title('y=x.^(1/3)') 图像：-5-4-3-2-101234500.20.40.60.811.21.41.61.8xyy=x.(1/3)y=x.(1/3)2.e x y -=2的图像：()55≤≤-x 编程：>> x=-5:0.1:5; >> y=exp(-x.^2); >> plot(x,y) >> grid on>> xlabel('x');ylabel('y') >> legend('y=e^(-x^2)'); >> title('y=e^(-x^2)') 图像：-5-4-3-2-101234500.10.20.30.40.50.60.70.80.91xyy=e (-x 2)y=e (-x 2)=++=+=axy a y at x yx tt t313,1333222的图像：()55≤≤-x ，a=1编程：>> t=-5:0.1:5;>> x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2); >> plot(x,y) >> grid on>> xlabel('x');ylabel('y')>> title('x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2)') >> legend('x=3*t./(1+t.^2);y=3*t.^2./(1+t.^2)') 图像：-1.5-1-0.500.51 1.500.511.522.53xyx=3*t./(1+t.2);y=3*t.2./(1+t.2)4.?+=+=x a a y a x xyt tt t3223221,1的图像: ()55≤≤-t ,a=1 编程：>> t=-5:0.1:5;>> x=t.^2./(1+t.^2);y=t.^3./(1+t.^2); >> plot(x,y)>> xlabel('x');ylabel('y')>> title('x=t.^2./(1+t.^2);y=t.^3./(1+t.^2)') >> legend('x=t.^2./(1+t.^2);y=t.^3./(1+t.^2)') >> grid on 图像：0.10.20.30.40.50.60.70.80.91-5-4-3-2-1012345xyx=t.2./(1+t.2);y=t.3./(1+t.2)5.()()t b y t t a x cos 1,sin -=-=的图像：pi t pi *2*2≤≤-,a=1,b=2 编程：>> t=-2*pi:0.1:2*pi; >> x=t-sin(t);y=2*(1-cos(t)); >> plot(x,y) >> grid on>> xlabel('x');ylabel('y')>> legend('x=t-sin(t);y=2*(1-cos(t))') >> title('x=t-sin(t);y=2*(1-cos(t))') 图像：-8-6-4-20246800.511.522.533.54xyx=t-sin(t);y=2*(1-cos(t))x=t-sin(t);y=2*(1-cos(t))6.=+==ayx t a y t a x 32323233sincos ,的图像：pi t pi *2*2≤≤-,a=1 编程：>> t=-2*pi:0.1:2*pi; >> x=(cos(t)).^3;y=(sin(t)).^3; >> plot(x,y) >> grid on>> xlabel('x');ylabel('y')>> title('x=(cos(t)).^3;y=(sin(t)).^3') 图像：-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81xyx=(cos(t)).3;y=(sin(t)).3x=(cos(t)).3;y=(sin(t)).37.ct z t b y t a x ===,sin ,cos 的图像：()pi t pi c b a *2*2,3,2,1≤≤-=== 编程：>> t=-2*pi:0.1:2*pi; >> x=cos(t);y=2*sin(t);z=3*t; >> plot3(x,y ,z)>> xlabel('x');ylabel('y');zlabel('z') >> grid on>> legend('x=cos(t);y=2*sin(t);z=3*t') >> title('x=cos(t);y=2*sin(t);z=3*t') 图像：-101-2-1012-20-101020xx=cos(t);y=2*sin(t);z=3*tyzx=cos(t);y=2*sin(t);z=3*t8.θa r =的图像：()pi a *20,1≤≤=θ编程：>> theta=0.0:0.1:2*pi; >> r=theta; >> polar(theta,r) >> grid on>> legend('r=theta') >> title('r=theta') 图像：24 68302106024090270120300150330180r=theta r=theta9.e a r θ=的图像：()pi a *20,1≤≤=θ编程：>> theta=-2*pi:0.1:2*pi; >> r=exp(theta); >> polar(theta,r) >> grid on >> title('r=exp(theta)') >> legend('r=exp(theta)') 图像：100200 300400 5003021060240902701203001503301800r=exp(theta)r=exp(theta)10.()?-==+yx ayxar 22222222,2cos θ的图像：1=a 编程：>> theta=0:0.1:2*pi;>> r=sqrt(abs(cos(2*theta))); >> polar(theta,r) >> grid on>> title('r=sqrt(abs(cos(2*theta)))'); >> legend('r=sqrt(abs(cos(2*theta)))') 图像：0.20.4 0.60.8 13021060240902701203001503301800r=sqrt(abs(cos(2*theta)))11.()==+xy a yxar 2222*222,2sin θ的图像：a=1编程：>> theta=0:0.1:2*pi;>> r=sqrt(abs(sin(2*theta))); >> polar(theta,r) >> grid on>> title('r=sqrt(abs(sin(2*theta)))') >> legend('r=sqrt(abs(sin(2*theta)))') 图像：0.4 0.60.8 13021060240902701203001503301800r=sqrt(abs(sin(2*theta)))12.)cos 1(θ+=a r 的图像：a=1 编程：>> theta=0:0.1:2*pi; >> r=1+cos(theta); >> polar(theta,r) >> grid on >> legend('r=1+cos(theta)') >> title('r=1+cos(theta)') 图像：0.51 1.52302106024090270120300150330180r=1+cos(theta)r=1+cos(theta)练习2.21.求出下列极限值. （1）nnn n33+∞→；（2）()n n n n ++-+∞→122lim；（3）x x x 2cot lim→；（4）??? ?→x m xx cos lim 0；（5）--→111lim1e xx x ；（6）??-+∞→x x xx 2lim .解：（1）编程：>> syms n >> limit((n^3+3^n)^(1/n),n,inf) ans = 3（2）编程：>> syms n>> limit(sqrt(n+2)-2*sqrt(n+1)+sqrt(n),n,inf) ans =0（3）编程：>> syms x >> limit(x*cot(2*x),x,0) ans = 1/2（4）编程：>> syms x m >> limit((cos(m/x))^x,x,inf) ans =1（5）编程：>> syms x>> limit(1/x-1/(exp(x)-1),x,1) ans = (exp(1)-2)/(exp(1)-1) （6）编程：>> syms x >> limit(sqrt(x^2+x)-x,x,1) ans = 2^(1/2)-1 2.有个客户看中某套⾯积为180m 2，每平⽅⽶7500元的房⼦。