当前位置:
文档之家› matlab程序设计教程 第二版
matlab程序设计教程 第二版
3 a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] b=max(a) c=min(a) b=0 c=0 for i=1:20 if a(i)>b b=a(i); end
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
Tp=cond(P) Tp = 8.5175e+003 条件数越趋近于 1,矩阵的性能越好,所以帕斯卡矩阵性能更好。 5. A=[-29,6,18;20,5,12;-8,8,5] A = -29 6 18 20 5 12 -8 8 5 [V,D]=eig(A) V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351 V 为 A 的特征向量,D 为 A 的特征值。它们之间满足 A*V=V*D
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
1 0 0 0 0 0
课后答案网:
ans = 0 0 1 1 0 0 0 isfinite(A) ans = 1 0 0 0 1 1 1 7.结构矩阵建立 P46 结构矩阵名.成员名=表达式 单元矩阵建立 P47 矩阵元素用大括号括起来。 8.用稀疏存储方式可以可以节省内存空间; 运算规则上,稀疏矩阵与普通矩阵一样,只是 矩阵的存储方式不同。在运算过程中,稀疏存储矩阵可以直接参与运算。当参与运算的对象 不全是稀疏存储矩阵时,所得结果一般是完全存储形式。 第三章 实验 1.提示: 先根据输入的 3 位整数分别求个、十、百位上的数字 a、b、c,反向输出整数: 100*a+10*b+c 2(1) m=input('请输入一个百分制的成绩: '); if m>=0&m<=100 switch fix(m/10) case {9,10} disp('这个成绩所属的成绩等级:A'); case {8} disp('这个成绩所属的成绩等级:B'); case {7} disp('这个成绩所属的成绩等级:C'); case {6} disp('这个成绩所属的成绩等级:D'); otherwise disp('这个成绩所属的成绩等级:E'); end else disp('错误:输入的成绩不合理'); end
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
Matlab 程序设计教程(第二版)课后参考答案
第一章 实验 1.实验一第 1 题 2.自己验证总结 Matlab 的主要优点 3.实验一第 2 题 4.实验一第 3 题 5.网站 思考练习 1. 启动见书 P5 退出见书 P6 2. Matlab 主要功能见书 P2 4 种功能 3. 分行输入行末尾加续行符,即三个点“…” 4. 见书 P11 5. 直接在命令窗口输入 fac 第二章 实验 1.(1) w=sqrt(2)*(1+0.34245*10^(-6)) w = 1.4142 (2) a=3.5; b=5; c=-9.8; x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) x =0.9829 (3) a=3.32; b=-7.9; y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) y = -128.4271 (4) t=[2,1-3i;5,-0.65]; z=0.5*exp(2*t)*log(t+sqrt(1+t.*t)) z = 1.0e+004 * 0.0048 + 0.0002i 0.0048 - 0.0034i 1.5899 2.0090 - 1.3580i 2.实验二 第 1 题 3.实验二 第 2 题 4. H=hilb(5); P=pascal(5); Hh=det(H) Hh = 3.7493e-012 Hp=det(P) Hp = 1 Th=cond(H) Th = 4.7661e+005
课后答案网:
课后答案网 您最真诚的朋友
网团队竭诚为学生服务,免费提供各门课后答案,不用积分,甚至不用注 册,旨在为广大学生提供自主学习的平台! 课后答案网: 视频教程网: PPT课件网:
4(1) A=1 for i=2:100 A=(-1)^(i+1)*1/(2*n-1)+A; end A (2) A=0 for i=1:1000 A=1/((4)^i)+A; end A
5 function [m,n]=chengji(A,B) m=A.*B; n=A*B; 然后在命令文件p3_12.m中调用chengji.m文件 a=[1 2 3;4 5 6;7 8 9]
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
function m=f1(n) m=n+10*log(n^2+5); f2.m 代码如下: function m=f2(n) m=0; for i=1:n m=m+i*(i+1); end 然后在命令文件p3_11.m中调用f1.m和f2.m p3_11.m代码如下: y1=f1(40)/(f1(30)+f1(20)); y2=f2(40)/(f2(30)+f2(20)); y1 y2 8
课后答案网:
if a(i)<c c=a(i); end end b c 4 a=[-3:0.1:3]; b=((exp(0.3.*a)-exp(-0.3.*a))./2) c=sin(a+0.3) y=b.*c+log((0.3+a)./2) for i=1:61 y=(exp(0.3*a(i))-exp(-0.3*a(i)))/2*sin(a(i)+0.3)+log((0.3+a(i))/2) end 5 (1) A=0; n=100 for i=1:n A=1/n/n+A; end (2) A=1; n=10000; for i=1:n A=(2*n*n*2)/((2*n-1)*(2*n+1))*A; end 6 A=rand(5,6) n=input('ÇëÊäÈëÒ»¸ön'); if n>5 n=5 A(n,n) else A(n) end 7 f1.m 代码如下:
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
polar(x,y) 4 t=0:pi/100:2*pi; x=exp(-t/20).*cos(t); y=exp(-t/20).*sin(t); z=t; plot3(x,y,z)
100
∑ k 函数文件factor1.m代码如下:
k =1
function f=factor1(k) if k==1 f=1; else f=factor1(k-1)+k; end
50
∑ k ^2 函数文件factor2.m代码如下:
k =1
function f=factor2(k) if k==1 f=1; else f=factor2(k-1)+^2; end
(2) [x,y]=meshgrid(-5:5); z=ones(11); surf(x,y,z);shading interp;
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
思考练习 1 务必自己看课本完成 2(1) x=1:0.1:10 a=exp((x.^2)/2) y=a.*1/2/pi plot(x,y) (2) t=-pi:0.1:2*pi X=t.*sin(t) Y=t.*cos(t) plot(X,Y) 3 t=0:0.00005:pi; x=sin(3*t).*cos(t); y1=2*x-0.5; y2=sin(3*t).*sin(t); plot(x,y2,x,y1); hold on k=find(abs(y2-y1)<1e-4); x1=x(k); y3=2*x1-0.5; plot(x1,y3,'rp')
10
∑ 1 / k 函数文件factor3.m代码如下:
k =1
function if k==1 f=1; else
f=factor3(k)
若侵犯了您的版权利益,敬请来信告知!
课 ww 后 w. 答 ha 案 ck 网 sh p. cn
课后答案网:
思考练习 1. 在 MATLAB 中 6+7i 是一个复数常量,6+7*i 则是一个表达式。i 是虚数单位,而 I 是单 位向量。 2. A.*B 表示 A 和 B 单个元素之间对应相乘,A*B 是按矩阵乘法的规则乘。 A./B 表示 A 矩阵除以 B 矩阵的对应元素,B.\A 等价于 A./B。 A/B 表示 A *inv(B)而 B\A 表示 inv(B)* A,对于矩阵运算一般两式不等。 3.(1) A(7)=[] (2) abs(x) (3) reshape(x,3,4) (4) k=find(t==0); t(k)=0 (5) ones(size(A)) (6) diag(A) B= diag(diag(A)) 4. 3+randn(25,20) 5.实验二 第 3 题 6. A=[34,NaN,Inf,-Inf,-pi,eps,0]; all(A) ans = 0 any(A) ans = 1 isnan(A) ans = 0 isinf(A)