MATLAB程序设计与应用课后习题标准答案
实验四
1．根据 2 1 1 1 1 ，求 的近似值。当 n 分别取 100、1
6 12 22 32
n2
０00、１0000 时，结果是多少？要求：分别用循环结构和向量运算来实现。
向量运算:
3．考
虑以下迭代公
式：
xn1
b
a
xn
。其
中 a、b 为正的常数。
x （1） 编写程序求迭代的结果，迭代的终止条件为，迭代初值 ＝1.０,迭 0
实验五
二、实验内容
4.设
f
(x)
(
x
1 2)2
0.1
1 (x3)40.01
，编写一个ＭATLＡB
函
数文件 fx.m，使得调用 f (x) 时,x 可用矩阵代入，得出的 f (x) 为同阶矩阵。
５.已知 y
f (40)
f (30) f (20)
（1)当 f (n)n10 ln n25 时,求 y 的值。
３、4、7 时的情况,不难发现,当 n 为奇数时画出的图有奇数个环,而当 n 为偶
数时画出的图有该偶数的两倍个环。参数ａ控制极坐标的半径,参数 b 可对图
进行角度旋转。
6.绘制曲面图形，并进行插值着色处理
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 ) ,
西安科技大学
MAＴLAB 程序设计
专业： 班级: 学号: 姓名：
信息与计算科学 1００１班
１008060129 刘仲能
201２年 6 月 27 日
实验一
2.已知:
12 34 4
1 3 1
A 34 7
87
，
B
2
0
3
3 65 7
3 2 7
ﻩ求下列表达式的值：
（1)Ａ+6*Ｂ和 A-B+I(其中 I 为单位矩阵）
代次数不超过 500 次。
（2） 如果迭代过程收敛于ｒ,那么 r 的准确值是 b b24a ，当 2
（a,b）的值取(1，1)、(8,3)、（10，0.1）时,分别对迭代结果和准确值
进行比较。
(１) (2)
5.若两个连续自然数的乘积减 1 是素数，则称这两个连续自然数是亲密数对, 该素数是亲密素数。例如,2×3—１＝5 是素数,所以 2 和 3 是亲密数对,5 是亲 密素数。求[２,５0]区间内： (１)亲密数对的对数。 （2)与上述亲密数对对应的所有亲密素数之和。
1/ 4 1/ 5
xx12
0.95 0.67
1/ 4 1/ 5 1/ 6 x3
0.52
（1） 求方程的解；
（2）将方程右边向量元素 改为０.５3，在求解，并比较 的变化和解的相 对变化；
（3）计算系数矩阵 A 的条件数并分析结论。
（2） 变大，其解中,相对未变化前的 的解:x１变大，ｘ2 变小,x３变大。
先利用默认属性绘制曲线,然后通过图形句柄操作来改变曲线的颜色、线 型和线宽，并利用文字对象给曲线添加文字标注。
实验八
1. 利用 MＡTLAB 提供的 rand 函数生成 30000 个符合均匀分布的随机 数，然后检验随机数的性质： （1） 均值和标准方差。 （2） 最大元素和最小元素。 （3） 大于 0.5 的随机数个数占总数的百分比。
（2）当时 f (n)122334nn1 ,求 y 的值。
(１）
（２)
实验六
1.
设
y 0.5
3sin x 1 x2
cos
x
,在 x=0~2 区间取１０1 点,绘制函数的
曲线。
4.绘制极坐标曲线 a sinbn ，
并分析参数 a、b、ｎ对曲线形状的影响。
以上五张截图分别是 a=1,ｂ=1,n=１、2、
（3)当 x 取矩阵 A 的每一元素时，求 P(x) 的值。其中:
1 1.2 1.4
A 0.75 2
3.5
0 5 2.5
（4） 当以矩阵Ａ为自变量时,求 P(x) 的值。其中 A 的值与第（3）
题相同。
(1)
（2)
(2)
(３)
实验九
9
7
0
21 22 23 24 25
4 13 11
（1） 求它们的乘积Ｃ。
（2） 将矩阵 C 的右下角 3×2 子矩阵赋给Ｄ。
（3） 查看 MATLAB 工作空间的使用情况
(１)
(2)
(3)
4.完成下列操作 (1）求[１0０,9９9］之间能被２1 整除的数的个数。 （2)建立一个字符串向量，删除其中的大写字母。
２
５.0
试用三次样条插值分别求出该日室内外 6:3０~
3
P 个多项式， (x) 1
x4 2x3 4x2 5
, P2 (x)
x
2
,时进行
下列操作:
(1）求 P(x) P1 (x) P2 (x) P3 (x) 。
（2)求 P(x) 的根。
(1)
(2)
(3）
2. 某气象观测站测得某日６:0０～18:00 之间每隔２h 的室内外温度
（℃）如实验表 1 所示。
实验表 1 室内外温度观测结果（℃）
时间 h
6
8
10
１2
14
1６
18
室内温度 t1 18．0 2０.0 22.0 ２5.0 ３0．0 ２8.0 2４.０
室外温度 t １ 19．0 24．0 28．0 ３4.０ ３2.0 30.0
(２）A*B 和 A.*Ｂ
(3）A^3 和 A.^３
(４)Ａ/B 及 B＼Ａ
（5)［A,B]和［A([1，3],:);B＾２]
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
16 17 18 19 20
(1)
(２)
实验二
3.建立一个 5×５矩阵，求它的行列式值、迹、秩和范数。 运行截图：
Ａ矩阵的行列式值、迹、秩分别如下：
范数如下：
4.已知
29 6 18
A
20
5 12
8 8 5
求 A 的特征值及特征向量，并分析其数学意义。
运行截图：
1/ 2 5.下面是一个线性方程组： 1/ 3
1/ 3 1/ 4
(3）
由于 A 矩阵的条件数很大,故当线性方程组中的ｂ变大时，x 也将发生很大的变化, 即数值稳定性较差。
实验三
3.硅谷公司员工的工资计算方法如下： (1）工作时数超过１20 小时者,超过部分加发 15%； (2)工作时数低于６0 小时者，扣发 70０元; (3）其余按每小时 8４元计发。 试编程按输入的工号和该员工的工时数，计算应发工资。
MATLAB程序设计与应用第二版课后题答案2(最新整理)
-0.6863 0.5621 0.4615
-0.0937 -0.6976 0.7103
D=
-0.0166
0
0
0 1.4801
0
0
0 2.5365
第四章
1.a=input('请输入一个 4 位数:');
while (a<1000|a>9999)
a=input('输入错误，请重新输入一个4位数:');
end
9
function f=factor(n,m); y=0; for k=1:n
y=y+k^m; end
10.(1)S=108 (2)x=4 12 20
y=2 4 6
1. (1) x=-10:0.1:10; y=100./(1+x.^2); plot(x,y)
第五章
(2) x=-10:0.1:10; y=1/(2*pi)*exp(-x.^2/2); plot(x,y)
第二章 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 21 34 98 21 8 42 121 43 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 = 000010000 L2 =
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程序设计与应用课后习题答案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)0 122sin851ze =+（2）21ln(2z x=，其中2120.455ix+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--（4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t=0：0.5:2。
52. 已知：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.（1） z1=2*sin(pi*85/180)/(1+exp(2))（2）x=[2,1+2i;-0.45,5];z2=log(x+sqrt(1+x.^2))/2（3）a=-3.0:0.1:3.0;z3=0.5*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+log((0.3+a)/2)（4）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)2. A=[12,34,-4;34,7,87;3,65,7];B=[1,3,-1;2,0,3;3,-2,7];（1）a=A+6*Bb=A-B+eye(size(A)) %I=eye(size(A))（2）c=A*Bd=A.*B（3）e=A^3f=A.^3（4）g=A/Bh=B\A（5）m=[A,B]n=[A([1,3],:);B^2]3. A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];（1）C=A*B（2）D=C(3:end,2:end)4.（1） a=100:999;b=rem(a,21)==0;c=find(b);d=length(c)（2）ch='Just as Bianhaiman said,Xiehong is ...';e=find(ch>='A'&ch<='Z');ch(e)=[]实验二 MATLAB矩阵分析与处理1.E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([1,2]);A=[E,R+R*S;O,S^2];a=A^22. H=hilb(5);P=pascal(5);Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)3. A=rand(5,5);a=det(A) %行列式的值b=trace(A) %行列式的迹c=rank(A) %行列式的秩d=norm(A) %行列式的范数4. A=[-29,6,18;20,5,12;-8,8,5];[V,D]=eig(A)5. A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];b=[0.95,0.67,0.52]';（1）x=inv(A)*b（2）b1=[0.95,0.67,0.53]';x=inv(A)*b1（3）c=cond(A)6. A=rand(3,3);B=sqrtm(A)C=sqrt(A)实验三选择结构程序设计1. x=input('请输入x的值:');if x<0&x~=-3y=x*x+x-6;elseif x>=0&x<5&x~=2&x~=3y=x*x-5*x+6;elsey=x*x-x-1;endy2.（if语句）score=input('请输入成绩:');if score>=0&score<60disp('E');elseif score>=60&score<70disp('D');elseif score>=70&score<80disp('C');elseif score>=80&score<90disp('B');elseif score>=90&score<=100disp('A');elsedisp('您输入的成绩有误，请重新输入！');end（switch语句）score=input('请输入成绩:');switch fix(score/10)case num2cell(0:5)disp('E');case {6}disp('D');case {7}disp('C');case {8}disp('B');case {9}disp('A');otherwisedisp('您输入的成绩有误，请重新输入！'); end3. a=input('请输入工号:','s');b=input('请输入该员工的工时数:');if b>120y=120*84+(b-120)*(1+15/100)*84;elseif b<60y=b*84-700;elsey=b*84;endy4．a=10+fix(89*rand());b=10+fix(89*rand());c=input('请输入一个运算符:','s'); if c=='+'disp(a+b);elseif c=='-'disp(a-b);elseif c=='*'disp(a*b);elsedisp(a/b);end5. A=rand(5,6);n=1:5;n=input('请输入n的值:');tryB=A(n,:);catchB=A(end,:);endBlasterr实验四循环结构程序设计1.（循环结构）n=input('请输入n的值:');y=0;for i=1:ny=y+1/i/i;pi=sqrt(6*y);endpi（向量运算）n=input('请输入n的值:');i=1:n;y=1./i.^2;f=sum(y);pi=sqrt(f*6)2. y=0;n=0;while 1n=n+1;y=y+1/(2*n-1);if y>=3breakendendny=y-1/(2*n-1)3. a=input('输入a的值:');b=input('输入b的值:');x0=1.0;n=1;x1=a/(b+x0);while abs(x1-x0)>=10^(-5)|n<=500 x1=a/(b+x0);x0=x1;n=n+1;endx14. f(1)=1;f(2)=0;f(3)=1;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3); enda=max(f) %最大值b=min(f) %最小值c=sum(f) %各数之和d=length(find(f>0)) %正数的个数e=length(find(f==0)) %零的个数f=length(find(f<0)) %负数的个数5. n=0;s=0;for k=2:49sn=k*(k+1)-1;m=fix(sqrt(sn));for i=2:mif rem(sn,i)==0breakelseendendif i==mn=n+1;s=s+sn;endendns实验五函数文件1.（函数文件）function [e,l,s,c]=f1(x)e=exp(x);l=log10(x);s=sin(x);c=cos(x);（命令文件）x=input('请输入一个复数:'); [e,l,s,c]=f1(x);elsc2.（函数文件）function [x]=f2(A,B)C=B';x=C*inv(A);（命令文件）m1=input('Enter m1:');m2=input('Enter m2:');theta=input('Enter theta:'); g=9.8;A=[m1*cos(theta),-m1,-sin(theta),0;m1*sin(theta),0,cos(theta),0;0,m2,-sin(theta),0;0,0,-cos(theta),1];B=[0,m1*g,0,m2*g]';[x]=f2(A,B);x3.（函数文件）function f=f3(x)f=1;for i=2:sqrt(x)if rem(x,i)==0f=0;breakendend（命令文件）for i=10:99j=10*rem(i,10)+fix(i/10);if f3(i)&f3(j)disp(i);endend4.（函数文件）function f=f4(x)f=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);（命令文件）x=[1,2,3;4,5,6];y=f4(x)5.（1）（函数文件）function f=f5(n)f=n+10*log(n^2+5);（命令文件）y=f5(40)/(f5(30)+f5(20))（2）（函数文件）function f=f6(i)g=0;for i=1:nf=g+i.*(i+1);end（命令文件）y=f6(40)/(f6(30)+f6(20))实验六高层绘图操作1. x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.*x));plot(x,y);title('实验六第一题');xlabel('X');ylabel('Y');text(3,0.6,'y=0.5+3sin(x)/(1+x^2)');2.（1）x=linspace(0,pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'b',x,y2,'g',x,y3,'r'); title('同一坐标系下得三条曲线');text(2,4,'y1=x^2');text(2,-0,8,'y2=cos(2x)');text(2,-2.2,'y=x^2*cos(2x)');xlabel('X');ylabel('Y');（2）x=linspace(0,pi,100); y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);plot(x,y1);title('y1=x^2');xlabel('X');ylabel('Y'); subplot(1,3,2);plot(x,y2);title('y2=cos(x)');xlabel('X');ylabel('Y') subplot(1,3,3);plot(x,y3);title('y3=x^2*cos(x)'); xlabel('X');ylabel('Y');（3）x=linspace(0,2*pi,20); y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(4,3,1);bar(x,y1);subplot(4,3,2);bar(x,y2);subplot(4,3,3);bar(x,y3);subplot(4,3,4);stairs(x,y1);subplot(4,3,5);stairs(x,y2);subplot(4,3,6);stairs(x,y3);subplot(4,3,7);stem(x,y1);subplot(4,3,8);stem(x,y2);subplot(4,3,9);stem(x,y3);subplot(4,3,10);fill(x,y1,'b');subplot(4,3,11);fill(x,y2,'r');subplot(4,3,12);fill(x,y3,'g');3. x=linspace(-5,5,100);y=[];for x0=xif x0>=-5&x<=0y=[y,(x0+sqrt(pi))/exp(2)];elsey=[y,log(x0+sqrt(1+x0.^2))/2];endendplot(x,y);title('实验六第三题');xlabel('X');ylabel('Y');4. a=input('请输入a:');b=input('请输入b:');n=input('请输入n:');theta=0:0.01:2*pi;rho=a.*sin(b+n.*theta);polar(theta,rho,'r');title('极坐标曲线');（如：当输入a=3,b=4,n=6时，图形为：）5. x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-sqrt(x.*x+y.*y)/4); subplot(2,1,1);surf(x,y,z);title('曲面图');xlabel('X');ylabel('Y');zlabel('Z'); subplot(2,1,2);surfc(x,y,z);title('等高线图');xlabel('X');ylabel('Y');zlabel('Z');6．（方法一）s=0:pi/100:pi/2;t=0:pi/100:3*pi/2;[s,t]=meshgrid(s,t);x=cos(s).*cos(t);y=cos(s).*sin(t);z=sin(s);surf(x,y,z);shading interp;title('实验六第六题');xlabel('X');ylabel('Y');zlabel('Z');（方法二）ezsurf('cos(s)*cos(t)','cos(s)*sin(t)','sin(s)',[0,0.5*pi,0,1.5*pi]); shading interp;实验七低层绘图操作1. hf=figure('MenuBar','figure','NumberTitle','off','Color','r',...'Name','图形窗口示例',...'WindowButtonDownFcn','disp(''Left Button Pressed.'')');2. x=-1:0.1:1;y=x.^2.*exp(2*x);h=line('XData',x,'YData',y);text(0.6,0.36*exp(1.2),'\leftarrow y=x^2*exp(2*x)');set(h,'Color','r','LineStyle','--','LineWidth',2);grid on;title('函数y=x^2*exp(2*x)');3. x=0:0.00001:0.001;[x,t]=meshgrid(x);v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi);axes('view',[-37.5,30]);hf=surface(x,t,v);grid on;xlabel('X');ylabel('Y');zlabel('Z');title('曲面v=10*exp(-0.01*x)*sin(2000*pi*t-0.2*x+pi)');set(hf,'EdgeColor','interp','FaceColor','r','LineStyle',':','LineWidth',2);4. x=linspace(0,2*pi,20);y1=sin(x);y2=cos(x);y3=sin(x)./(cos(x)+eps);y4=cos(x)./(sin(x)+eps);axes('Position',[0.15,0.1,0.2,0.4]);plot(x,y1);title('y1=sin(x)');axes('Position',[0.05,0.6,0.2,0.3]);plot(x,y2);title('y2=cos(x)');axes('Position',[0.45,0.6,0.3,0.3]);plot(x,y3);title('y3=tan(x)');axes('Position',[0.5,0.2,0.3,0.2]);plot(x,y4);title('y4=cot(x)');5. cylinder(3);light('Position',[1,2,4]); shading interp;lighting gouraud;material shiny实验八数据处理与多项式计算1．A=rand(1,30000);a=mean(A) %均值b=std(A,0,2) %标准方差c=max(A) %最大元素d=min(A) %最小元素n=0;for i=1:30000if A(i)>0.5n=n+1; %大于0.5的随机数的个数endendm=n/30000 %百分比2. P=45+fix(50*rand(100,5));[zgf,xh]=max(P) %每门课的最高分及相应的学号[zdf,xh]=min(P) %每门课的最低分及相应的学号pjf=mean(P) %每门课的平均分fc=std(P) %标准方差A=sum(P,2); %总分[zfzg,xh]=max(A) %总分最高分及学号[zfzd,xh]=min(A) %总分最低分及学号[zcj,xsxh]=sort(A,1,'descend')3. h=6:2:18;x=6.5:2:17.5;t1=[18.0 20.0 22.0 25.0 30.0 28.0 24.0];t2=[15.0 19.0 24.0 28.0 34.0 32.0 30.0];a=interp1(h,t1,x,'spline')b=interp1(h,t2,x,'spline')4. x=linspace(1,101,10);y=log10(x);p=polyfit(x,y,5);y1=polyval(p,x);plot(x,y,'r:o',x,y1,'-*')5. p1=[1,2,4,0,5];p2=[1,2];p3=[1,2,3];p4=conv(p2,p3);p41=[0,p4];p=p1+p41 %P（x）=P1(x)+P2(x)*P3(x) x=roots(p) %根A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];y1=polyval(p,A)y2=polyvalm(p,A)。
《MATLAB程序设计与应用》刘卫国高等教育出版社课后答案 已解锁-去水印 适合打印
SY404 clear all for n=1:4 if n==1 f1=1; elseif n==2 f2=0; elseif n==3 f3=1; else a=f3-2*f2+f1; b=a-2*f3+f2; c=b-2*a+f3; d=c-2*b+a; H=[1,0,1,a,b,c,d]; for m=8:4:99 a=d-2*c+b; b=a-2*d+c; c=b-2*a+d; d=c-2*b+a; H=[H,a,b,c,d]; end f100=d-2*c+b; end end max=max(H);
[f1,f2]=f(n); a=f1; b=f2; elseif n==30; [f1,f2]=f(n); c=f1; d=f2; else [f1,f2]=f(n); e=f1; f=f2; end end y1=e/(a+c); y2=f/(b+d); disp(['(1) y=',num2str(y1)]) disp(['(2) y=',num2str(y2)]) f function [f1,f2]=f(n) f1=n+10*log(n^2+5); x=0; for a=1:n b=a*(a+1); x=x+b; end f2=x; fushu function [e,l,s,c]=fushu(x) e=exp(x); l=log(x); s=sin(x); c=cos(x); disp(['复数e的指数是：',num2str(e)]) disp(['复数e的对数是：',num2str(l)]) disp(['复数e的正弦是：',num2str(s)]) disp(['复数e的余弦是：',num2str(c)])
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也将发生专门大的转变，即数值稳固性较差。