MATLAB程序设计与应用课后习题答案
MATLAB程序设计与应用课后实验答案
M A T L A B程序设计与应用(刘卫国编)课后实验答案(总66页)本页仅作为文档封面，使用时可以删除This document is for reference only-rar21year.March实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z 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 --+=++=-- (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0::2. 已知：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]解：运算结果：3. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ (1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
解：. 运算结果：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 +⎡⎤=⎢⎥⎣⎦。
MATLAB程序设计与应用(刘卫国编)课后实验答案
实验一 MATLAB 运算基础1. 先求以下表达式的值，然后显示MATLAB 工作空间的利用情形并保留全数变量。
(1) 0122sin851z 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 --+=++=-- (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:: 解：2. 已知：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]解：运算结果：3. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ (1) 求它们的乘积C 。
(2) 将矩阵C的右下角3×2子矩阵赋给D。
(3) 查看MATLAB工作空间的利用情形。
解：. 运算结果：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 +⎡⎤=⎢⎥⎣⎦。
matlab程序设计与应用第二版习题答案
matlab程序设计与应用第二版习题答案Matlab程序设计与应用第二版习题答案Matlab是一种强大的数学软件，广泛应用于科学计算、数据分析和工程设计等领域。
《Matlab程序设计与应用》是一本经典的教材，对于学习和掌握Matlab编程语言具有重要的意义。
本文将为大家提供《Matlab程序设计与应用第二版》中部分习题的答案，帮助读者更好地理解和应用Matlab。
第一章：Matlab基础1.1 基本操作1. a = 3; b = 4; c = sqrt(a^2 + b^2); disp(c);2. x = linspace(-pi, pi, 100); y = sin(x); plot(x, y);3. A = [1 2 3; 4 5 6; 7 8 9]; B = [9 8 7; 6 5 4; 3 2 1]; C = A + B; disp(C);1.2 控制结构1. for i = 1:10disp(i);end2. n = 0; sum = 0; while sum < 100n = n + 1;sum = sum + n;enddisp(n);3. x = 5; if x > 0disp('x is positive');elseif x < 0disp('x is negative');elsedisp('x is zero');end第二章：向量和矩阵运算2.1 向量运算1. A = [1 2 3]; B = [4 5 6]; C = A .* B; disp(C);2. A = [1 2 3]; B = [4 5 6]; C = A ./ B; disp(C);3. A = [1 2 3]; B = [4 5 6]; C = dot(A, B); disp(C);2.2 矩阵运算1. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A * B; disp(C);2. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = B * A; disp(C);3. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A .* B; disp(C); 第三章：函数和脚本文件3.1 函数1. function y = myfunc(x)y = x^2 + 3*x + 2;end2. function [y1, y2] = myfunc(x1, x2)y1 = x1^2 + 3*x1 + 2;y2 = x2^2 + 3*x2 + 2;end3. function [y1, y2] = myfunc(x)y1 = x^2 + 3*x + 2;y2 = sin(x);end3.2 脚本文件1. x = linspace(0, 2*pi, 100); y = sin(x); plot(x, y);2. x = linspace(-10, 10, 100); y = x.^2 + 3*x + 2; plot(x, y);3. x = linspace(0, 2*pi, 100); y1 = sin(x); y2 = cos(x); plot(x, y1, x, y2);通过以上习题的答案，读者可以对Matlab程序设计的基本语法和常用函数有一个初步的了解。
MATLAB程序设计与应用 刘卫国主编 高等教育出版社课后答案
y='A'; elseif x>=80
y='B'; elseif x>=70
y='C'; elseif x>=60
y='D'; else
y='E'; end y
3：
x=input('请输入工作小时数');
天天learn
if 84*x-700<=0 y=0;
1、 A=rand(100,300)
(1) mean(A) std(A,0,1);
(2) [B,C]=max(A) [E,F]=max(B)
[B1,C1]=min(A)
(3) Length(find(A>0.5))/30000
[E1,F2]=min(B1)
2、 P=45+(95-45)*rand(100,5) (1) [Y,U]=max(P) [Y1,U1]=min(P) (2) mean(P) std(P,0,1) (3) K=sum(P,2) [E,F]=max(K) [E1,F1]=min(K) (4) [zcj,xsxh]=-sort(-K,1)
break end end format long y=y-1/(2*k-1) k=k-1 3：a=input('a='); b=input('b='); n=1; x0=1; x1=a/(x0+b); while abs(x0-x1)>=10^(-5)|n<=500 k=a/(x1+b); x0=x1; x1=k; n=n+1; end x1
MATLAB程序设计及应用(第二版)课后实验答案
Matlab 课后实验题答案实验一 MATLAB 运算基础1.先求下列表达式的值，然后显示 MATLA 工作空间的使用情况并保存全部变量。
0.3a 0.3ae e0.3sin (a 0.3) ln(1)Zi2sin85 0 1 e 2Z 211n(x .1 x 2)，其中 x22i0.45Z 32t20 t 1(4) z t2 1 1 t 2，其中t =0::t22t 1 2 t 3解：2.已知：12 34 4 1 3 13.0, 2.9,川,2.9,3.0A 34 7 87 ,B 2 0 33 65 7 3 2 7求下列表达式的值:(1)A+6*B和A-B+I (其中I为单位矩阵)(2)A*B 和 A.*B(3)A A3 和 A.A3(4)A/B 及B\A(5)[A,B]和[A([1,3],:);BA2]解：AA3A.A3A/BB\A[A,B][A([1,3],:);BA2]3.设有矩阵A和B1 2 3 4 5 3 0 166 7 8 9 10 17 6 9A 11 12 13 14 15 ,B 0 23 416 17 18 19 20 9 7 021 22 23 24 25 4 13 11(1)求它们的乘积Co(2)将矩阵C的右下角3X 2子矩阵赋给D(3)查看MATLA工作空间的使用情况。
解：.运算结果：E=(reshape(1:1:25,5,5))';F二[3 0 16;17 -6 9;0 23 -4;9 7 0;4 1311];C= E*FH=C(3:5,2:3)C =93 150 77258 335 237423 520 397588 705 557753 890 717 H =520 397705 557890 7174.完成下列操作：(1)求［100,999］之间能被21整除的数的个数。
(2)建立一个字符串向量，删除其中的大写字母解: (1)结果：m=100:999;n=fin d(mod(m,21)==0);len gth( n)ans =43(2).建立一个字符串向量例如：ch二'ABC123d4e56Fg9：则要求结果是:ch二'ABC123d4e56Fg9：k=fi nd(ch>='A'&ch<='Z');ch(k)=[]ch =123d4e56g9实验二MATLAB 矩阵分析与处理1.设有分块矩阵AO33 R32'其中E R O S分别为单位矩阵、解：M 文件如下;随机矩阵、零矩阵和对角阵，试通过数值计算验证AE R2RSO S 2B Fditor - D TI T it L^d?tEBB2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P,且求其行列式的 值Hh 和Hp 以及它们的条件数Th 和Tp ,判断哪个矩阵性能更好。
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.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.1766Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.63510.6119Columns 57 through 600.5777 0.5327 0.4774 0.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.25002. 已知：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 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =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 7ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C。
MATLAB程序设计与应用课后习题答案Word版
1/ 4 1/ 5
xx12
0.95
0.67
1/ 4 1/ 5 1/ 6 x3
0.52
（1） 求方程的解；
（2）将方程右边向量元素 改为 0.53，在求解，并比较 的变化和解的相对 变化；
（3）计算系数矩阵 A 的条件数并分析结论。
（2） 变大，其解中，相对未变化前的 的解：x1 变大，x2 变小，x3 变大。 （3）
(2)当时 f (n)122334nn1 ，求 y 的值。
(1)
(2)
实验六
1.
设
y
0.5
3sin x 1 x2
cos
x
，在 x=0~2 区间取 101 点，绘制函数
的曲线。
4．绘制极坐标曲线 a sinbn ,
并分析参数 a、b、n 对曲线形状的影响。
以上五张截图分别是 a=1,b=1,n=1、2、
实验五
二、实验内容
4.设
f
(x)
(x2)120.1
(
x
1 3)40.01
，编写一个
MATLAB
函数
文件 fx.m，使得调用 f (x) 时,x 可用矩阵代入，得出的 f (x) 为同阶矩阵。
5.已知 y
f (40)
f (30) f (20)
(1)当 f (n)n10 ln n25 时，求 y 的值。
试用三次样条插值分别求出该日室内外 6:30~17:30 之间每隔 2h 各 点的近似温度（℃）。
5.有
3
P 个多项式， (x) 1
x4 2x3 4x2 5
， P2 (x)
x
2
，时进
行下列操作：
（1）求 P(x) P1 (x) P2 (x) P3 (x) 。
MATLAB程序设计及应用(第二版)课后实验答案
Matlab课后实验题答案《实验一 MATLAB运算基础1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)0 122sin851ze =+(2)21ln( 2z x=，其中2120.455i x+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0::^2. 已知：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])3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
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 OS +⎡⎤=⎢⎥⎣⎦。
MATLAB程序设计及应用(第二版)课后实验答案
Matlab 课后实验题答案实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z 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 --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:: 解：2. 已知：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]解：3. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ (1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
解：. 运算结果：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 +⎡⎤=⎢⎥⎣⎦。
MATLAB程序设计与应用课后习题答案
西安科技大学MATLAB程序设计专业：信息与计算科学班级：1001班学号：1008060129姓名：刘仲能2021年6月27日实验一2.已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=76538773443412A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=723302131B求以下表达式的值：（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]3.设有矩阵A 和B ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ，⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=11134079423096171603B （1） 求它们的乘积C 。
（2） 将矩阵C 的右下角3×2子矩阵赋给D 。
（3） 查看MATLAB 工作空间的利用情形(1)(2)（3）4.完成以下操作（1）求[100,999]之间能被21整除的数的个数。
（2）成立一个字符串向量，删除其中的大写字母。
(1)(2)实验二3.成立一个5×5矩阵，求它的行列式值、迹、秩和范数。
运行截图：A 矩阵的行列式值、迹、秩别离如下：范数如下：4.已知 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=5881252018629A 求A 的特点值及特点向量，并分析其数学意义。
运行截图：5.下面是一个线性方程组：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解；（2） 将方程右边向量元素b 3改成0.53，在求解，并比较b 3的转变和解的相对转变；（3） 计算系数矩阵A 的条件数并分析结论。
（2） b 3变大，其解中，相对未转变前的b 3的解：x1变大，x2变小，x3变大。
（3）由于A矩阵的条件数专门大，故当线性方程组中的b变大时，x也将发生专门大的转变，即数值稳固性较差。
MATLAB程序设计与应用 答案
第1章MATLAB概论1.1与其他计算机语言相比较，MA TLAB语言突出的特点是什么？MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。
1.2 MA TLAB系统由那些部分组成？MATLAB系统主要由开发环境、MA TLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装MA TLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。
第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。
1.4 MA TLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？在MA TLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。
1.5 如何启动M文件编辑/调试器？在操作桌面上选择‚建立新文件‛或‚打开文件‛操作时，M文件编辑/调试器将被启动。
在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。
1.6 存储在工作空间中的数组能编辑吗？如何操作？存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。
1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
MATLAB程序设计与应用课后实验答案
A-B+eye(3)
A*B
A.*B
A^3
A.^3
A/B
B\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.*B
A-B+eye(3)
A*B
A.*B
A^3
A.^3
A/B
B\A
0 0 0 0
ans =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
由ans,所以
2. 产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好。为什么？
解：M文件如下：
输出结果：
H =
P =
1 1 1 1 1
Columns 57 through 60
Column 61
z4 =
0 0
2. 已知：
求下列表达式的值：
(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];
实验四 循环结构程序设计
一、实验目的
1. 掌握利用for语句实现循环结构的方法。
2. 掌握利用while语句实现循环结构的方法。
MATLAB程序设计与应用课后习题答案
由于A矩阵的条件数很大，故当线性方程组中的b变大时，x也将发生很大的变化，即数值稳定性较差。
实验三
3.硅谷公司员工的工资计算方法如下：
（1）工作时数超过120小时者，超过部分加发15%；
（2）工作时数低于60小时者，扣发700元；
（3）其余按每小时84元计发。
试编程按输入的工号和该员工的工时数，计算应发工资。
3.设有矩阵A和B ，
（1）求它们的乘积C。
（2）将矩阵C的右下角3×2子矩阵赋给D。
（3）查看MATLAB工作空间的使用情况
(1) (2)
（3）
4.完成下列操作
（1）求[100,999]之间能被21整除的数的个数。
（2）建立一个字符串向量，删除其中的大写字母。
(1)(2)
实验二
3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
西安科技大学
MATLAB程序设计
专业：信息与计算科学
班级：1001班
学号：1008060129
姓名：刘仲能
20XX年6月27日
实验一
2.已知：
，
求下列表达式的值：
（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]
（2）最大元素和最小元素。
（3）大于0.5的随机数个数占总数的百分比。
(1)(2)（3）
2.某气象观测站测得某日6:00~18:00之间每隔2h的室内外温度（℃）如实验表1所示。
实验表1 室内外温度观测结果（℃）
时间h
6
8
10
MATLAB程序设计与应用第二版课后题答案
如对您有帮助，欢迎下载支持，谢谢！第二章3.设矩阵A为A=[24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21];（1）B=A(2:5,1:2:5)B =65 24 2134 98 218 42 12143 45 21（2）A(7)=[]A =24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21（3）A+30（4）size(A);ndims(A)（5）题目有误（6）reshape(x,3,4)（7）abs(x)（8）char(x)4. L1 =0 0 0 0 1 0 0 0 0L2 =1 1 1 1 1 0 0 0 0L3 =0 0 0 1 1 1 0 0 0L4 =4 5 65.(1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D =65.0000 5.0000如对您有帮助，欢迎下载支持，谢谢！0 32.000054.0000 3.1400E =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.2103(2)E<Dans =0 10 00 1E&Dans =1 10 11 1E|Dans =1 11 11 1~D|~Eans =0 01 00 0find(A>=10&A<25)ans =156.all(A)ans =any(A)ans =1isnan(A)ans =0 1 0 0 0 0 0isinf(A)ans =如对您有帮助，欢迎下载支持，谢谢！0 0 1 1 0 0 0isfinite(A)ans =1 0 0 0 1 1 17.A(1).x1=’学号’;A(1).x2=’姓名’;A(1).x3=’专业’;A(1).x4.x41=’成绩1’;……….A(2).x1=’学号’;A(2).x2=’姓名’;A(2).x3=’专业’;A(2).x4.x41=’成绩1’;……….A(3).x1=’学号’;A(3).x2=’姓名’;A(3).x3=’专业’;A(3).x4.x41=’成绩1’;……….A(4).x1=’学号’;A(4).x2=’姓名’;A(4).x3=’专业’;A(4).x4.x41=’成绩1’;……….A(5).x1=’学号’;A(5).x2=’姓名’;A(5).x3=’专业’;A(5).x4.x41=’成绩1’;……….8.(1)size(B)ans =2 2ndims(B)ans =2(2)B(2)ans =[3x3 doubleB(4)ans ={3x3 cell}(3)B(3)=[]B =[1] [3x3 double] {3x3 cell}B{3}=[]B =[1] [3x3 double] []第三章1.(1)A=eye(3)(2) A=100+100*rand(5,6)(3)A=1+sqrt(0.2)*randn(10,50)(4)B=ones(size(A))(5)A+30*eye(size(A))(6)B=diag(diag(A))2.B=rot90(A)C=rot90(A,-1)如对您有帮助，欢迎下载支持，谢谢！3.B=inv(A) ;A的逆矩阵C=det(A) ;A的行列式的值D=A*BE=B*AD=E 因此A与A-1是互逆的。
MATLAB程序设计与应用第二版课后题答案 (1)
第二章3.设矩阵A为A=[24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21];（1）B=A(2:5,1:2:5)B =65 24 2134 98 218 42 12143 45 21（2）A(7)=[]A =24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21（3）A+30（4）size(A);ndims(A)（5）题目有误（6）reshape(x,3,4)（7）abs(x)（8）char(x)4. L1 =0 0 0 0 1 0 0 0 0L2 =1 1 1 1 1 0 0 0 0L3 =0 0 0 1 1 1 0 0 0L4 =4 5 65.(1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D =65.0000 5.00000 32.000054.0000 3.1400E =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.2103(2)E<Dans =0 10 00 1E&Dans =1 10 11 1E|Dans =1 11 11 1~D|~Eans =0 01 00 0find(A>=10&A<25)ans =156.all(A)ans =any(A)ans =1isnan(A)ans =0 1 0 0 0 0 0isinf(A)ans =0 0 1 1 0 0 0isfinite(A)ans =1 0 0 0 1 1 17.A(1).x1=’学号’;A(1).x2=’姓名’;A(1).x3=’专业’;A(1).x4.x41=’成绩1’;………. A(2).x1=’学号’;A(2).x2=’姓名’;A(2).x3=’专业’;A(2).x4.x41=’成绩1’;………. A(3).x1=’学号’;A(3).x2=’姓名’;A(3).x3=’专业’;A(3).x4.x41=’成绩1’;………. A(4).x1=’学号’;A(4).x2=’姓名’;A(4).x3=’专业’;A(4).x4.x41=’成绩1’;………. A(5).x1=’学号’;A(5).x2=’姓名’;A(5).x3=’专业’;A(5).x4.x41=’成绩1’;……….8.(1)size(B)ans =2 2ndims(B)ans =2(2)B(2)ans =[3x3 doubleB(4)ans ={3x3 cell}(3)B(3)=[]B =[1] [3x3 double] {3x3 cell}B{3}=[]B =[1] [3x3 double] []第三章1.(1)A=eye(3)(2) A=100+100*rand(5,6)(3)A=1+sqrt(0.2)*randn(10,50)(4)B=ones(size(A))(5)A+30*eye(size(A))(6)B=diag(diag(A))2.B=rot90(A)C=rot90(A,-1)3.B=inv(A) ;A的逆矩阵C=det(A) ;A的行列式的值D=A*BE=B*AD=E 因此A与A-1是互逆的。
MATLAB程序设计和应用课后习题答案解析
西安科技大学MATLAB程序设计专业：信息与计算科学班级：1001班学号：1008060129姓名：刘仲能2012年6月27日实验一2.已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=76538773443412A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=723302131B 求下列表达式的值：（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]3.设有矩阵A 和B ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ，⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=11134079423096171603B （1） 求它们的乘积C 。
（2） 将矩阵C 的右下角3×2子矩阵赋给D 。
（3） 查看MATLAB 工作空间的使用情况(1) (2)（3）4.完成下列操作（1）求[100,999]之间能被21整除的数的个数。
（2）建立一个字符串向量，删除其中的大写字母。
(1)(2)实验二3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
运行截图：A 矩阵的行列式值、迹、秩分别如下：范数如下：4.已知 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=5881252018629A 求A 的特征值及特征向量，并分析其数学意义。
运行截图：5.下面是一个线性方程组：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解； （2） 将方程右边向量元素改为0.53，在求解，并比较的变化和解的相对变化；（3） 计算系数矩阵A 的条件数并分析结论。
（2） 变大，其解中，相对未变化前的的解：x1变大，x2变小，x3变大。
（3）由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变化，即数值稳定性较差。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
对，5 是亲密素数。求[2,50]区间内： （1）亲密数对的对数。 （2）与上述亲密数对对应的所有亲密素数之和。
实验五
二、实验内容
4.设
f
(x)
(x2)120.1
(
x
1 3)40.01
，编写一个
MATLAB
函数
文件，使得调用 f (x) 时,x 可用矩阵代入，得出的 f (x) 为同阶矩阵。
5.已知 y
（2） 如果迭代过程收敛于 r，那么 r 的准确值是 b b24a ，当(a， 2
b)的值取(1，1)、(8，3)、(10，时，分别对迭代结果和准确值进行比 较。
(1) (2)
5.若两个连续自然数的乘积减 1 是素数，则称这两个连续自然数是亲密数 对，该素数是亲密素数。例如，2×3—1＝5 是素数，所以 2 和 3 是亲密数
1
2
1
2
4. 已知
0 1 0
1 0 0
a b c
p 1
1
0
0
,
p 2
0
1
0 , A d
e
f
0 0 1
1 0 1
g h i
完成下列运算：
（1）
B
p1
•
p 2
•
A
。
（2）B 的逆矩阵并验证结果。
（2） 包括 B 矩阵主对角线元素的下三角阵。（4）B 的行列式值。
5. 用符号方法求下列极限或倒数。
(1)
(1)
(2)
实验二
3.建立一个 5×5 矩阵，求它的行列式值、迹、秩和范数。 运行截图：
A 矩阵的行列式值、迹、秩分别如下： 范数如下：
4.已知
29 6 18
A
20
5 12
8 8 5
求 A 的特征值及特征向量，并分析其数学意义。
运行截图：
1/ 2 5.下面是一个线性方程组： 1/ 3
1/ 3 1/ 4
y
0
y'(0) 0
8. 求微分方程组的数值解，并绘制解的曲线。
y y y
'
1
2
3
y' 2
y y 13
y
'
0.05 y
y
3
1
2
y (0) 1
0, y 2(0)
1,
y (0)1 3
实验十
1. 已知
，利用符号表达求
。
2. 分解因式。 (1)
3. 化简表达式。
(1) sin cos cos sin
f (40)
f (30) f (20)
(1)当 f (n)n10 ln n25 时，求 y 的值。
(2)当时 f (n)122334nn1 ，求 y 的值。
(1)
(2)
实验六
1.
设
y
0.5
3sin x 1 x2
cos
x
，在 x=0~2 区间取 101 点，绘制函数的
曲线。
4．绘制极坐标曲线 a sinbn ,
6. 用符号方法求下列积分。
(2)
实验十一
1. 计算 s
10 n1
1 2n 1
2. 将 ln x 在 x=1 处按 5 次多项式展开为泰勒级数。
3. 求下列方程的符号解。
（1） ln(1+x)
=2
4.求微分方程初值问题的符号解，并与数值解进行比较。
5.求微分方程组的通解。
$
（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]
1 2 3 4 5
3 0 16
6
7
8
9
10
17 6
9
3.设有矩阵 A 和 B A 11 12 13 14 15 ， B 0 23 4
表 1 所示。
实验表 1 室内外温度观测结果（℃）
时间 h
6
8
10
12
14
16
18
室内温度 t1
室外温度 t2
试用三次样条插值分别求出该日室内外 6:30~17:30 之间每隔 2h 各点的近似温度（℃）。
5.有 3
P 个多项式， (x) 1
x4 2x3 4x2 5
， P2 (x)
x
2
，时进
行下列操作：
（1）求 P(x) P1 (x) P2 (x) P3 (x) 。
（2）求 P(x) 的根。
（3）当 x 取矩阵 A 的每一元素时，求 P(x) 的值。其中：
1 1.2 1.4
A 0.75 2
3.5
0 5 2.5
（4） 当以矩阵 A 为自变量时，求 P(x) 的值。其中 A 的值与第（3）题相同。
9x y 4z u 13
3x
4y
2z
2u
1
3x 9 y 2u 11
直接解法：
LU 分解：
通解法：
4. 求非齐次线性方程组的通解。
32xx1 157xx2 223xx3 32xx4464
9x1
4x2
ห้องสมุดไป่ตู้x3
7x4
2
5. 求代数方程的数值解。
z x y （2） 在给定的初值
1，
0
1，
MATLAB 程序设计与应用课后习 题答案
西安科技大学
MATLAB 程序设计
专业： 班级： 学号： 姓名：
信息与计算科学 1001 班
29 刘仲能
2012 年 6 月 27 日
实验一
2.已知：
12 34 4
1 3 1
A 34
7
87
，
B
2
0
3
3 65 7
3 2 7
求下列表达式的值：
（3）
由于 A 矩阵的条件数很大，故当线性方程组中的 b 变大时，x 也将发生很大的变 化，即数值稳定性较差。
实验三
3.硅谷公司员工的工资计算方法如下： （1）工作时数超过 120 小时者，超过部分加发 15%； （2）工作时数低于 60 小时者，扣发 700 元； （3）其余按每小时 84 元计发。 试编程按输入的工号和该员工的工时数，计算应发工资。
16 17 18 19 20
9
7
0
21 22 23 24 25
4 13 11
（1） 求它们的乘积 C。
（2） 将矩阵 C 的右下角 3×2 子矩阵赋给 D。
（3） 查看 MATLAB 工作空间的使用情况
(1)
(2)
（3）
4.完成下列操作 （1）求[100,999]之间能被 21 整除的数的个数。 （2）建立一个字符串向量，删除其中的大写字母。
实验四
1.根据 2 1 1 1 1 ，求 的近似值。当 n 分别取 100、1000、
6 12 22 32
n2
10000 时，结果是多少要求：分别用循环结构和向量运算来实现。
向量运算：
3. 式：
其中 a、b 为正的常数。
考虑
以下迭代公 a。
xn1 b xn
x （1） 编写程序求迭代的结果，迭代的终止条件为，迭代初值 =，迭代次数 0 不超过 500 次。
并分析参数 a、b、n 对曲线形状的影响。
以上五张截图分别是 a=1,b=1,n=1、2、3、4、7 时的情况，不难发现，当 n 为奇数时画出的图有奇数个环，而当 n 为偶数时画出的图有该偶数的两倍个 环。参数 a 控制极坐标的半径，参数 b 可对图进行角度旋转。 6.绘制曲面图形，并进行插值着色处理
1. 利用 MATLAB 提供的 rand 函数生成 30000 个符合均匀分布的随机数，然后检 验随机数的性质：
（1） 均值和标准方差。 （2） 最大元素和最小元素。 （3） 大于的随机数个数占总数的百分比。
(1)
(2)
（3）
2. 某气象观测站测得某日 6:00~18:00 之间每隔 2h 的室内外温度（℃）如实验
1/ 4 1/ 5
xx12
0.95
0.67
1/ 4 1/ 5 1/ 6 x3
0.52
（1） 求方程的解；
（2） 将方程右边向量元素 改为，在求解，并比较 的变化和解的相对变化； （3） 计算系数矩阵 A 的条件数并分析结论。
（2） 变大，其解中，相对未变化前的 的解：x1 变大，x2 变小，x3 变大。
(1)
(2)
(2)
(3)
实验九
1. 求函数在指定点的数值导数。
x x2 x3 f (x) 1 2x x2 ， x 1,2,3
0 2 6x
2. 用数值方法求定积分。
I （1）
2
10
cost2 4sin(2t)2 1dt 的近似值。
3. 分别用 3 种不同的数值方法解线性方程组。
6x 5y 2z 5z 4
xcos S cos t
y
cos
S
sin
t
z sin S
0s ,0t 3
2
2
实验七
2. 利用曲面对象绘制曲面 v(x,t)10e0.01x sin(2000t 0.2x ) ，先
利用默认属性绘制曲线，然后通过图形句柄操作来改变曲线的颜色、线 型和线宽，并利用文字对象给曲线添加文字标注。
实验八
0
1 下，求方程组的数值解。
0
sin
x
y2
ln
x
7
0
3x 2y z3 1 0