MATLAB程序设计与应用课后习题标准答案

合集下载

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=' 输入的分数有错' elseif x>=90
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也将发生专门大的转变，即数值稳固性较差。

MATLAB程序设计与应用课后实验答案

A+6.*B
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程序设计与应用课后习题答案

（3）
由于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程序设计与应用(刘卫国编)课后实验答案

真验一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 tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t解：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 =z2 =z3 =Columns 1 through 4Columns 5 through 8Columns 9 through 12Columns 13 through 16Columns 17 through 20Columns 21 through 24Columns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 Columns 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.6351 0.6119Columns 57 through 600.5777 0.5327 0.4774 0.4126Column 610.3388z4 =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]; 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 118820ans =1728 39304 -6439304 343 65850327 274625 343ans =ans =ans =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战B(1) 供它们的乘积C.(2) 将矩阵C的左下角3×2子矩阵赋给D.(3) 查看MATLAB处事空间的使用情况.解：. 运算截止：E=(reshape(1:1:25,5,5))';F=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; C= E*FH=C(3:5,2:3)C =93 150 77258 335 237423 520 397588 705 557753 890 717H =520 397705 557890 7174. 完毕下列支配：(1) 供[100,999]之间能被21整除的数的个数.(2) 修坐一个字符串背量，简略其中的大写字母.解：(1) 截止：m=100:999;n=find(mod(m,21)==0);length(n)ans =43(2). 修坐一个字符串背量比圆：ch='ABC123d4e56Fg9';则央供截止是：ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch =真验二MATLAB矩阵分解与处理1. 设有分块矩阵33322322E RAO S⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E、R、O、S分别为单位矩阵、随机矩阵、整矩阵战对于角阵，试通过数值估计考证22E R RS AO S+⎡⎤=⎢⎥⎣⎦.解: M文献如下；输出截止：S =1 00 2A =0 0 0 1.0000 0a =0 0 0 1.0000 0 ans =0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0由ans,所以22E R RS AO S+⎡⎤=⎢⎥⎣⎦2. 爆收5阶希我伯特矩阵H战5阶帕斯卡矩阵P，且供其止列式的值Hh战Hp以及它们的条件数Th战Tp，推断哪个矩阵本能更佳.为什么？解：M文献如下：输出截止：H =P =1 1 1 1 1123d4e56g91 2 3 4 51 3 6 10 151 4 10 20 351 5 15 35 70Hh =Hp =1Th =4.7661e+005Tp =8.5175e+003果为它们的条件数Th>>Tp,所以pascal矩阵本能更佳.3. 修坐一个5×5矩阵，供它的止列式值、迹、秩战范数.解： M文献如下：输出截止为：A =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9d =5070000t =65c1 =c2 =cinf =4. 已知供A的特性值及特性背量，并分解其数教意思.解：M文献如图：输出截止为：V =D =-25.3169 0 00 -10.5182 0数教意思：V的3个列背量是A的特性背量，D的主对于角线上3个是A的特性值，特别的，V的3个列背量分别是D的3个特性值的特性背量.5. 底下是一个线性圆程组：(1) 供圆程的解.(2) 将圆程左边背量元素b3改为0.53再供解，并比较b3的变更妥协的相对于变更.(3) 估计系数矩阵A的条件数并分解论断.解： M文献如下：输出截止：X =X2 =C =1.3533e+003由截止，X战X2的值一般，那表示b的微弱变更对于圆程解也做用较小，而A的条件数算得较小，所以数值宁静性较佳，A是较佳的矩阵.6. 修坐A矩阵，试比较sqrtm(A)战sqrt(A)，分解它们的辨别.解：M文献如下：运止截止有：A =16 6 1820 5 129 8 5b1 =b2 =b =分解截止知：sqrtm(A)是类似A的数值仄圆根（那可由b1*b1=A的截止瞅出），而sqrt(A)则是对于A中的每个元素启根号，二则辨别便正在于此.真验三采用结构步调安排一、真验手段1. 掌握修坐战真止M文献的要领.2. 掌握利用if语句真止采用结构的要领.3. 掌握利用switch语句真止多分支采用结构的要领.4. 掌握try语句的使用.二、真验真质1. 供分段函数的值.用if语句真止，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值.解：M文献如下：运算截止有：f(-5)y =14>> f(-3)y =11>> f(1)y =2>> f(2)y =1>> f(2.5)y =>> f(3)y =5>> f(5)y =192. 输进一个百分造结果，央供输出结果等第A、B、C、D、E.其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E.央供：(1) 分别用if语句战switch语句真止.(2) 输进百分造结果后要推断该结果的合理性，对于分歧理的结果应输出堕落疑息.解：M文献如下试算截止：score=88grade =Bscore=123过失：输进的结果没有是百分造结果3. 硅谷公司职工的人为估计要领如下：(1) 处事时数超出120小时者，超出部分加收15%.(2) 处事时数矮于60小时者，扣收700元.(3) 其余按每小时84元计收.试编程按输进的工号战该号职工的工时数，估计应收人为.解：M文献下4. 安排步调，完毕二位数的加、减、乘、除四则运算，即爆收二个二位随机整数，再输进一个运算标记，干相映的运算，并隐现相映的截止.解：M文献如下；运算截止例：a =38b =33输进一个运算符:^c =falsea =92b =40输进一个运算符:+c =1325. 修坐5×6矩阵，央供输出矩阵第n止元素.当n值超出矩阵的止数时，自动转为输出矩阵末尾一止元素，并给出堕落疑息.解：M文献如下：运算截止如下：输进一个5止6列矩阵A=[1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7]输进一正整数n=411 2 3 9 7 3输进一个5止6列矩阵A=[1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7]输进一正整数n=62 3 4 5 6 7ans =Error using ==> dispToo many input arguments.真验四循环结构步调安排一、真验手段1. 掌握利用for语句真止循环结构的要领.2. 掌握利用while语句真止循环结构的要领.3. 认识利用背量运算去代替循环支配的要领.二、真验真质1. 根据2222211116123nπ=++++，供π的近似值.当n分别与100、1000、10000时，截止是几？央供：分别用循环结媾战背量运算（使用sum函数）去真止.解：M文献如下：运止截止如下：K>> %循环结构估计pi值y=0;n=input('n=');for i=1:ny=y+1/i/i;endpi=sqrt(6*y)n=100pi =n=1000pi =n=10000pi =%背量要领估计Pi值n=input('n=');i=1./(1:n).^2;s=sum(i);pi=sqrt(6*s)n=100pi =n=1000pi =n=100002. 根据11113521yn=++++-，供：(1) y<3时的最大n值.(2) 与(1)的n值对于应的y值.解：M—文献如下：运止截止如下：K>> y=0;n=0;while y<3n=n+1;y=y+1/(2*n-1);endynif y>3n=n-1;endny =n =57n =563. 思量以下迭代公式：其中a、b为正的教数.(1) 编写步调供迭代的截止，迭代的末止条件为|x n+1-x n|≤10-5，迭代初值x0=1.0，迭代次数没有超出500次.(2) 如果迭代历程支敛于r，那么r(a,b)的值与(1,1)、(8,3)、(10,0.1)时，分别对于迭代截止战准确值举止比较.解：M文献如下：运算截止如下；请输进正数a=1请输进正数b=1x =r =r =s =请输进正数a=8请输进正数b=3x =r =r =s =0.0请输进正数a=10x =r =r =s =4. 已知供f1~f100中：(1) 最大值、最小值、各数之战.(2) 正数、整、背数的个数.解：M—文献以下是运算截止：max(f)=437763282635min(f)=-899412113528sum(f)=-742745601951c1=49c2=2c3=495. 若二个连绝自然数的乘积减1是素数，则称那二个边陲自然数是亲稀数对于，该素数是亲稀素数.比圆，2×3-1=5，由于5是素数，所以2战3是亲稀数，5是亲稀素数.供[2,50]区间内：(1) 亲稀数对于的对于数.(2) 与上述亲稀数对于对于应的所有亲稀素数之战.解：M文献：运算截止为：j =29s =23615真验五函数文献一、真验手段1. 明白函数文献的观念.2. 掌握定义战调用MATLAB函数的要领.二、真验真质1. 定义一个函数文献，供给定复数的指数、对于数、正弦战余弦，并正在下令文献中调用该函数文献.解：M文献如下：函数fushu.M文献：function [e,l,s,c] = fushu(z)%fushu 复数的指数，对于数，正弦，余弦的估计%e 复数的指数函数值%l 复数的对于数函数值%s 复数的正弦函数值%c 复数的余弦函数值e=exp(z);l=log(z);s=sin(z);c=cos(z);下令文献M：z=input('请输进一个复数z=');[a,b,c,d]=fushu(z)运算截止如下：z=input('请输进一个复数z=');[a,b,c,d]=fushu(z)请输进一个复数z=1+ia =b =c =d =2. 一物理系统可用下列圆程组去表示：从键盘输进m1、m2战θ的值，供a1、a2、N1战N2的值.其中g与9.8，输进θ时以角度为单位.央供：定义一个供解线性圆程组AX=B的函数文献，而后正在下令文献中调用该函数文献.解： M文献函数fc.M文献:function X= fc(A,B)%fc fc是供解线性圆程的函数%A A是已知矩阵的系数矩阵X=A\B；下令M文献：clc;m1=input('输进m1=');m2=input('输进m2=');theta=input('输进theta=');x=theta*pi/180;g=9.8;A=[m1*cos(x) -m1 -sin(x) 0m1*sin(x) 0 cos(x) 00 m2 -sin(x) 00 0 -cos(x) 1];B=[0;m1*g;0;m2*g];X=fc(A,B)运算截止：输进m1=1输进m2=1输进theta=30X =3. 一个自然数是素数，且它的数字位子通过任性对于换后仍为素数.比圆13是千万于素数.试供所有二位千万于素数.央供：定义一个推断素数的函数文献.解：M文献：function [p] = prime(p)% 输进p的范畴，找出其中的素数m=p(length(p));for i=2:sqrt(m)n=find(rem(p,i)==0&p~=i);p(n)=[];%将p中能被i整除，而却没有等于i的元素，即下标为n的元素剔除，其余的即为素数endp;下令文献：clc;p=10:99;p=prime(p); %找出10到99内的所有素数p=10*rem(p,10)+(p-rem(p,10))/10;%将p素数矩阵每个元素个位十位变更程序p=prime(p)%再对于对于换后的素数矩阵找出所有的素数运算截止：。

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:0.5:2.52. 已知：求下列表达式的值：(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(1) 求它们的乘积C。

(2) 将矩阵C的右下角3×2子矩阵赋给D。

(3) 查看MATLAB工作空间的使用情况。

解：. 运算结果：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

解: M 文件如下；输出结果：由ans,所以22E R RS A OS +⎡⎤=⎢⎥⎣⎦ 2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。

为什么？解：M 文件如下：3. 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。

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程序设计与应用（第二版）实验参考答案%实验一MATLAB运算基础%第一题ftp://192.168.0.143/%（1）z1=2*sin(85*pi/180)/(1+exp(2))ftp://192.168.0.143/%（2）x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))ftp://192.168.0.143/%(3)a=-3.0:0.1:3.0;ftp://192.168.0.143/z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)%(4)t=0:0.5:2.5;z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)%第二题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(size(A))A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]%第三题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]C=A*BF=size(C)D=C(F(1)-2:F(1),F(2)-1:F(2))whos%第四题%(1):A=100:999;B=rem(A,21);C=length(find(B==0))%(2):A='lsdhKSDLKklsdkl';k=find(A>='A'&A<='Z');A(k)=[]%实验二MATLAB矩阵分析与处理%第一题E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([2,3]);A=[E,R;O,S];A^2B=[E,(R+R*S);O,S^2]%第二题H=hilb(5)P=pascal(5)Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)%第三题:A=fix(10*rand(5))H=det(A)Trace=trace(A)Rank=rank(A)Norm=norm(A)%第四题:A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)%数学意义略%第五题方法一:%(1):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]';x=inv(A)*b%(2):B=[0.95,0.67,0.53]';x=inv(A)*B%(3):cond(A)%第五题方法二:A=hilb(4)A(:,1)=[]A(4,:)=[]B=[0.95,0.67,0.52]';X=inv(A)*BB1=[0.95,0.67,0.53]';X1=inv(A)*B1N=cond(B)N1=cond(B1)Na=cond(A) %矩阵A为病态矩阵%第六题A=[1,4,9;16,25,36;49,64,81]B=sqrtm(A)C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算%实验三选择程序结构设计%第一题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x0<0&x0~=-3y=[y,x0*x0+x0-6];elseif x0>=0&x0<5&x0~=2&x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx %输出所有xy %输出所有y%第一题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)<0&x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)>=0&x(a)<5&x(a)~=2&x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第二题程序一:x=input('请输入一个百分制成绩：');if x>100|x<0disp('您输入的成绩不是百分制成绩，请重新输入。

MATLAB程序设计与应用课后习题标准答案

MATLAB程序设计与应用第二版课后题答案2(最新整理)

matlab程序设计与应用第二版习题答案

MATLAB程序设计与应用 刘卫国主编 高等教育出版社课后答案

MATLAB程序设计与应用课后习题答案Word版

MATLAB程序设计及应用(第二版)课后实验答案

MATLAB程序设计及应用(第二版)课后实验答案

MATLAB程序设计与应用(刘卫国)》(第二版)-课后实验答案

《MATLAB程序设计与应用》刘卫国高等教育出版社课后答案 已解锁-去水印 适合打印

MATLAB程序设计及应用(第二版)课后实验答案

MATLAB程序设计与应用课后习题答案

MATLAB程序设计与应用课后实验答案

MATLAB程序设计与应用课后习题答案

MATLAB程序设计与应用(刘卫国编)课后实验答案

MATLAB程序设计与应用课后实验答案

MATLAB程序设计与应用第二版课后题答案 (1)

MATLAB程序设计与应用(第二版)实验参考答案

MATLAB程序设计与应用刘卫国主编高等教育出版社课后答案

《MATLAB程序设计与应用》刘卫国高等教育出版社课后答案已解锁-去水印适合打印