- 1、本文档共20页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB程序设计与应用(第二版)实验答案
MATLAB程序设计与应用(第二版)实验参考答案
%实验一 MATLAB运算基础
%第1题
%(1)
z1=2*sin(85*pi/180)/(1+exp(2))
%(2)
x=[2,1+2i;-0.45,5];
z2=0.5*log(x+sqrt(1+x.^2))
%(3)
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)
%(4)
t=0:0.5:2.5;
z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3)
%第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(size(A))
A*B
A.*B
A^3
A.^3
A/B
B\A
[A,B]
[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]
C=A*B
F=size(C)
D=C(F(1)-2:F(1),F(2)-1:F(2))
whos
%第4题
%(1):
A=100:999;
B=rem(A,21);
C=length(find(B==0))
%(2):
A=lsdhKSDLKklsdkl;
k=find(A=AA=Z);
A(k)=[]
%实验二 MATLAB矩阵分析与处理
%第1题
E=eye(3);
R=rand(3,2);
O=zeros(2,3);
S=diag([2,3]);
A=[E,R;O,S];
A^2
B=[E,(R+R*S);O,S^2]
%第2题
H=hilb(5)
P=pascal(5)
Hh=det(H)
Hp=det(P)
Th=cond(H)
Tp=cond(P)
%第3题
A=fix(10*rand(5))
H=det(A)
Trace=trace(A)
Rank=rank(A)
Norm=norm(A)
%第4题
A=[-29,6,18;20,5,12;-8,8,5]
[V,D]=eig(A)
%数学意义略
%第5题方法一
%(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)
%第5题方法二
A=hilb(4)
A(:,1)=[]
A(4,:)=[]
B=[0.95,0.67,0.52];
X=inv(A)*B
B1=[0.95,0.67,0.53];
X1=inv(A)*B1
N=cond(B)
N1=cond(B1)
Na=cond(A) %矩阵A为病态矩阵
%第6题
A=[1,4,9;16,25,36;49,64,81]
B=sqrtm(A)
C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算,sqrt函数是以矩阵中的元素进行计算
%实验三 选择程序结构设计
%第1题程序一
x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];
y=[]; %建立存放所有y值的矩阵
for x0=x
if x00x0~=-3
y=[y,x0*x0+x0-6];
elseif x0=0x05x0~=2x0~=3
y=[y,x0*x0-5*x0+6];
else
y=[y,x0*x0-x0-1];
end
end
x %输出所有x
y %输出所有y
%第1题程序二
x=[-5,-3,1,2,2.5,3,5];
y=[];
for a=1:7
if x(a)0x(a)~=-3
y=[y,(x(a))^2+x(a)-6];
elseif x(a)=0x(a)5x(a)~=2x(a)~=3
y=[y,(x(a))^2-5*x(a)+6];
else
y=[y,x(a)*x(a)-x(a)-1];
end
end
%第2题程序一
x=input(请输入一个百分制成绩:);
if x100|x0
disp(您输入的成绩不是百分制成绩,请重新输入。);
else
if x=100x=90
disp(A);
el
您可能关注的文档
- Liquid Art.ppt
- Liquid Computing A Real Effect.pdf
- Liquid Crystal Lasers_ Recent Advances.pdf
- liquid phase catalyst.pdf
- Liquid Polymer Nanocomposites PEGME-SnO2 and PEGME-TiO2.pdf
- Liquid Polymorphism and Density Anomaly in a Lattice Gas Model.pdf
- LIQUID SKIES.ppt
- Liquid sloshing in cylindrical tank with elastic space.pdf
- Liquid software a new paradigm for networked systems.pdf
- Liquid stability in a model for ortho-terphenyl.pdf
文档评论(0)