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程序设计与应用(第二版) 刘卫国课后实验报告

0 0
0 0
4.0000 0
0 9.0000
A,H,Trace,Rank,Norm 结果 H = -3.2509e+004 Trace =33 Rank =5 Norm =25.3208 4、 A=[-29,6,18;20,5,12;-8,8,5]; [V,D]=eig(A); 结果 A = -29 20 -8 5、 5 8 6 18 12 5
0 0 2、 H=hilb(5); P=pascal(5); Hh=det(H); Hp=det(P); Th=cond(H); Tp=cond(P); H,P,Hh,Hp,Th,Tp 结果 Hh = 3.7493e-012 Hp =1 Th = 4.7661e+005 Tp =8.5175e+003 3、 A=fix(10*rand(5)); H=det(A); Trace=trace(A); Rank=rank(A); Norm=norm(A);
b=?1 ans =1 x1 =0.6180 r10 =0.6180 r20 =-1.6180 4、 for i=1:100 if i==1 f(i)=1; elseif i==2 f(i)=0; elseif i==3 f(i)=1; else f(i)=f(i-1)-2*f(i-2)+f(i-3); end end a=max(f); b=min(f); c=sum(f); disp(['最大值',num2str(a)]); disp(['最小值',num2str(b)]); disp(['各数之和',num2str(c)]); zhengshu=length(find(f>0)); zeroshu=length(find(f==0)); fushu=length(find(f<0)); 结果最大值；4.377633e+011 最小值：-8.994121e+011 各数之和：-7.427456e+011 zhengshu =49 zeroshu = 2 fushu =49 5、 s=0;

Matlab程序设计与应用第二版刘卫国课后实验答案

实验一：T1：%%第一小题z1=2*sin(85*pi/180)/(1+exp(2))%%第二小题x=[2,1+2i;,5];z2=1/2*log(x+sqrt(1+x.^2));z2%%第三小题a=::;z3=1/2*(exp*a)-exp*a)).*sin(a++log(+a )/2)%%第四题t=0::z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t. ^2-1)+(t>=2&t<3).*(t.^2-2*t+1)T2：A=[12,34,-4;34,7,87;3,65,7]B=[1,3,-1;2,0,3;3,-2,7]disp ('A+6*B=');disp(A+6*B);disp('A-B+I=');disp(A-B+eye(3));disp('A*B=');disp(A*B);disp('A.*B=');disp(A.*B);disp('A^3=');disp(A^3);disp('A.^3=');disp(A.^3);disp('A/B=');disp(A/B);disp('B\A=');disp(B\A);disp('[A,B]=');disp([A,B]);disp('[A([1,3],:);B^2]=');disp([A([1,3],:);B^2]);T3:z=1:25;A=reshape(z,5,5)';B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3)T4-1:a=100:999;b=find(rem(a,21)==0);c=length(b)T4-2:a=input('请输入一个字符串：','s');b=find(a>='A'&a<='Z');a(b)=[];disp(a);实验二：T1：E=eye(3),R=rand(3,2),O=zeros(2,3),S=d iag([1,2]);A=[E,R;O,S]disp('A^2=');disp(A^2);disp('[E,R+RS;O,S^2]');B=[E,R+R*S;O,S^2]T2：H=hilb(5)P=pascal(5)Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)a=abs(Th-1);b=abs(Tp-1);if a>bdisp('帕萨卡矩阵P性能更好'); elseif a<bdisp('希尔伯特矩阵H性能更好'); elsedisp('两个矩阵性能相同');endT3:a=1:25;A=reshape(a,5,5)disp('行列式的值：');disp(det(A));disp('矩阵的秩：');disp(rank(A));disp('矩阵的迹：');disp(trace(A));disp('矩阵的范数：');disp(norm(A));T4:A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)T5:A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/ 6]B=[,,]'X1=A\BB(3)=X2=A\Bdisp('系数矩阵A的条件数：');disp(cond(A));T6:a=1:25;A=reshape(a,5,5)disp('B1=sqrtm(A)');B1=sqrtm(A)disp('B2=sqrt(A)');B2=sqrt(A)disp('B1*B1');B1*B1disp('B2.*B2');B2.*B2实验三：T1:x=:2:for i=1:length(x);if (x(i)<0)&(x(i)~=-3)y(i)=x(i)^2+x(i)-6;elseif(x(i)>=0)&(x(i)<5)&(x(i)~=2)&(x(i)~=3 );y(i)=x(i)^2-5*x(i)+6;elsey(i)=x(i)^2-x(i)-1;endendyT2-if:s=input('please enter the score:');while (s<0||s>100)disp('the score is not reasonable'); s=input('please enter the score:'); endif (s>=90&s<=100);disp('A');elseif(s>=80&s<90);disp('B');elseif(s>=70&s<80);disp('C');elseif(s>=60&s<70);disp('D');elsedisp('E');endT2-switch:s=input('please enter the score:'); while (s<0||s>100)disp('the score is not reasonable'); s=input('please enter the score:'); endswitch fix(s/10)case{9,10}disp('A');case{8}disp('B');case{7}disp('C');case{6}disp('D');case{0,1,2,3,4,5}disp('E');endT3:t=input('请输入工时time=');if t>120w=120*84+(t-120)**84;elseif t<60w=t*84-700;elsew=84*t;enddisp('应发工资为：');disp(w);T4:a=10+floor(rand(1)*89)b=10+floor(rand(1)*89)s=input('请输入+或-或*或/','s');while(s~='+'&&s~='-'&&s~='*'&&s~='/') disp('输入的符号错误，请重新输入'); s=input('请输入+或-或*或/','s'); endswitch scase{'+'}c=a+b;case{'-'}c=a-b;case{'*'}c=a*b;case{'/'}c=a/b;endcT5:A=rand(5,6)n=input('请输入n:');while (n<1)disp('输入的n有误，请重新输入'); n=input('请输入n:');endif n>5B=A(5,:)elseB=A(n,:)End实验四：T1:n=input('请输入n:');m=0;for i=1:n;m=m+1/i^2;endpii=sqrt(6*m)T1-2:n=input('请输入n:');m=1./(1:n).^2;s=sum(m);pii=sqrt(6*s)T2:n=1000;y=0;for i=1:n;y=y+1/(2*i-1);if y>=3disp('最大n值为：');i-1disp('对应的y值为：');y-1/(2*i-1)break;endendT3:format longa=input('请输入a:');b=input('请输入b:');x(1)=a/+b);i=1;x(2)=a/(x(i)+b);while (abs(x(i+1)-x(i))>&i<500)i=i+1;x(i+1)=a/(x(i)+b);endr1=(-b+sqrt(b^2+4*a))/2;r2=(-b-sqrt(b^2+4*a))/2;disp(['x',num2str(i+1),'=',num2str(x( i+1),8)]);disp(['x',num2str(i),'=',num2str(x(i) ,8)]);disp(['r1=',num2str(r1,8),',r2=',num2 str(r2,8)]);format shortT4:f(1)=1;f(2)=0;f(3)=1;i=4;while(i<=100)f(i)=f(i-1)-2*f(i-2)+f(i-3);i=i+1;endfmax=max(f)fmin=min(f)fsum=sum(f)fp=length(find(f>0))fn=length(find(f<0))f0=length(find(f==0))T5:j=0;s=0;for i=3:50m=i*(i-1)-1;if isprime(m)==1ms=s+m;j=j+1;endenddisp(['亲密数的对数：',num2str(j)]); disp(['所有亲密素数之和：',num2str(s)]);实验五：fn1:function f=fn1(n)f=n+10*log(n.^2+5);fn2:function f=fn2(n)f=0;for i=1:nf=f+i*(i+1);endfx:function f=fx(x)f=1./((x-2).^2++1./((x-3).^4+;T1:x=input('请输入一个复数：');fe=exp(x)flg=log(x)fsin=sin(x)fcos=cos(x)T2:m1=input('请输入m1:');m2=input('请输入m2:');t=input('请输入theta:');A=[m1*cos(t),-m1,-sin(t),0;m1*sin(t), 0,...cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1 ]B=[0,m1*,0,m2*]'disp('[a1,a2,N1,N2]=');disp([A\B]');T3:j=1;n=1;k=1;for i=10:99 % 挑选出10~99中所有的素数if isprime(i)==1m(j)=i;j=j+1;endendfor t=1:length(m);s(n)=10*rem(m(t),10)+fix(m(t)/10); %挑选出的素数进行位置交换if isprime(s(n))==1 %再挑选出交换绝对素数r(k)=m(t);k=k+1;endn=n+1;enddisp(['所有两位绝对素数：',num2str(r)])T5:y1=fn1(40)/(fn1(30)+fn1(20)) %函数调用y2=fn2(40)/(fn2(30)+fn2(20))实验六：T1:x=linspace(0,2*pi,101);y=+3*sin(x)./(1+x.^2)).*cos(x);plot(x,y);T2:x=linspace(0,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'y-',x,y3,'b-.'); subplot(2,2,2);plot(x,y1,'r:');subplot(2,2,3);plot(x,y2,'y-');subplot(2,2,4);plot(x,y3,'b-.');T2-3:x=linspace(0,2*pi,10);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(3,4,1);bar(x,y1,'r');axis([0,7,0,40]); subplot(3,4,2);stairs(x,y1,'y');axis([0,7,0,40]); subplot(3,4,3);stem(x,y1,'b');axis([0,7,0,40]); subplot(3,4,4);fill(x,y1,'g');axis([0,7,0,40]); subplot(3,4,5);bar(x,y2,'r');axis([-1,7,-1,1]); subplot(3,4,6);stairs(x,y2,'y');axis([-1,7,-1,1]); subplot(3,4,7);stem(x,y2,'b');axis([-1,7,-1,1]); subplot(3,4,8);fill(x,y2,'g');axis([-1,7,-1,1]); subplot(3,4,9);bar(x,y3,'r');axis([0,7,-30,50]); subplot(3,4,10);stairs(x,y3,'y');axis([0,7,-30,50]); subplot(3,4,11);stem(x,y3,'b');axis([0,7,-30,50]); subplot(3,4,12);fill(x,y3,'g');axis([0,7,-30,50]);T3:x=linspace(-5,5,100);y=[];for x0=x;if x0<=0;y=[y,(x0+sqrt(pi))/exp(2)]; elsey=[y,1/2*log(x0+sqrt(1+x0^2))];endendsubplot(1,2,1);plot(x,y);for i=1:length(x)if x(i)<=0y(i)=(x(i)+sqrt(pi))/exp(2); elsey(i)=1/2*log(x(i)+sqrt(1+x(i)^2)); endendsubplot(1,2,2);plot(x,y);T4:a=input('请输入a:');b=input('请输入b:');n=input('请输入n:');theta=0::2*pi;y=a.*sin(b+n.*theta);polar(theta,y);T5:x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2 )/4);subplot(1,2,1);surf(x,y,z);subplot(1,2,2);contour3(x,y,z);T6:s=0:pi/2;t=0: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;实验八：T1:A=rand(1,30000);disp('均值为：');mean(A)disp('标准方差为：');std(A)disp('最大元素为：');max(A)disp('最小元素为：');min(A)k=find(A>;b=length(k);disp('大于的随机数个数占总数的百分比：');sprintf('%f%%', b/300)T2:P=45+50*rand(100,5)disp('最高分及序号：');[Y,U]=max(P,[],1)disp('最低分及序号：');[Y1,U1]=min(P,[],1)disp('每门课平均分：');mean(P,1)disp('每门课的标准方差：');std(P,0,1)s=sum(P,2);disp('总分最高分：');[Y2,U2]=max(s,[],1)disp('总分最低分：');[Y3,U3]=min(s,[],1)[zcj,xsxh]=sort(s,1,'descend') %desc end降序T3:h=6:2:18;t1=[,,,,,,];t2=[,,,,,,];t=:2:;T1=interp1(h,t1,t,'spline')T2=interp1(h,t2,t,'spline')T4:x=1:10:101;y=log10(x);p=polyfit(x,y,5)x1=1::101;y1=log10(x1);y2=polyval(p,x1);plot(x1,y1,'r-',x1,y2,'b-.');T5:P1=[1,2,4,0,5];P2=[0,0,0,1,2];P3=[0,0,1,2,3];P4=conv(P2,P3);sp1=length(P1);sp2=length(P4);P=[zeros(1,sp2-sp1),P1]+P4x=roots(P)A=[-1,,;,2,;0,5,];Ax=polyval(P,A)AA=polyvalm(P,A)实验九：T1:f=inline('det([x x^2 x^3;1 2*x 3*x^2;0 2 6*x])');i=1;x=1;while x<=g(i)=f(x);i=i+1;x=x+;enddx=diff(g)/;dx(1)dx(101)dx(length(g)-1)T2:f1=inline('sqrt(cos(t.^2)+4*sin(2*t). ^2+1)')f2=inline('log(1+x)./(1+x.^2)')I1=quad(f1,0,2*pi)I2=quad(f2,0,2*pi)T3:A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0 ,2];B=[-4,13,1,11]';x1=A\B[L,U]=lu(A);x2=U\(L\B)[Q,R]=qr(A);x3=R\(Q\B)fun:function F=fun(X);x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y^2+log(z)-7;F(2)=3*x+2^y-z^3+1;F(3)=x+y+z-5;T5:f=inline('3*x+sin(x)-exp(x)');fzero(f,X=fsolve('fun',[1 1 1]',optimset('Display','off'))实验十：T1:x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y))T2:syms x y;t=sym('5135');factor(x^4-y^4)factor(t)T3:syms beta1 beta2 x;simplify(sin(beta1)*cos(beta2)-cos(be ta1)*sin(beta2))simplify((4*x^2+8*x+3)/(2*x+1))T4:syms a b c d e f g h k;p1=[0,1,0;1,0,0;0,0,1];p2=[1,0,0;0,1,0;1,0,1];A=[a,b,c;d,e,f;g,h,k];B=p1*p2*AB1=inv(B)B1*B;tril(B)det(B)T5:syms x t a yf1=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/sin(x)^3;limit(f1,x,0) %(1)f2=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1) ;limit(f2,x,-1,'right') %(2)f3=(1-cos(2*x))/x;diff(f3,x,1) %(3)diff(f3,x,2)A=[a^x,t^3;t*cos(x),log(x)];diff(A,x,1) &(4)diff(A,t,2)diff(diff(A,x,1),t,1)f5=(x^2-2*x)*exp(-x^2-y^2-x*y);yx=-diff(f5,x,1)/diff(f5,y,1) %(5) fxy=diff(diff(f5,x,1),y,1)x=sym('0');y=sym('1');eval(fxy)T6:x=sym('x');f1=1/(1+x^4+x^8);int(f1,x)f2=1/(asin(x))^2/sqrt(1-x^2);int(f2,x)f3=(x^2+1)/(x^4+1);int(f3,x,0,inf)f4=exp(x)*(1+exp(x))^2;int(f4,x,0,log(2))实验十一：T1：syms n x ;S1=symsum(1/(2*n-1),n,1,10)S2=symsum(n^2*x^(n-1),n,1,inf)S3=symsum(n^2/5^n,n,1,inf)T2：x=sym('x');f=log(x);taylor(f,x,6,1)T3：x1=solve('log(1+x)-5/(1+sin(x))=2')x2=solve('x^2+9*sqrt(x+1)-1=0')x3=solve('3*x*exp(x)+5*sin(x)=0')[x4,y4]=solve('sqrt(x^2+y^2)-100=0',' 3*x+5*y-8=0')T4：dsolve('D2y+4*Dy+29*y=0','y(0)=0','Dy (0)=15','x')T5：[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4 *x-5*y+3*z','Dz=4*x-4*y+2*z','t')。

MATLAB程序设计与应用第二、三章作业;

MATLAB程序设计与应用第二、三章作业
第一章
1.A/B和A\B的数学含义。

解：A/B等效于B矩阵的逆右乘A矩阵，也就是A*inv(B)，A、B矩阵的列数应该相同；A\B等效于A的逆左乘B矩阵，也就是inv(A)*B，A、B矩阵的行数应该相同。

2. ①将矩阵A的２～５行，１，３，５列元素赋值给Ｂ矩阵；
②将矩阵A中的0元素用机器零eps代替；
③将12元素的向量转换成3×4的矩阵。

第二章
1.利用冒号法分别按升序和降序生成字符串，要求字符串包括大写字母A到小写字母z的所有字符。

2.向量Z为1×5的全1向量，要求将Z扩充为4×5的全1矩阵，然后将其转换为2×10的矩阵。

3.编写程序，将当前目录中的a1.bmp，……，a9.bmp读入变量空间，调用命令I=imread （’filename’）。

4.解线性方程
4x1+2x2−x3=2 3x1−x2+2x3=10 12x1+3x2=8。

《MATLAB_程序设计与应用》刘卫国第2与4章课后答案-推荐下载

matlab程序设计与应用刘卫国第2与4章课后答案
西安建筑科技大学陈
3.写出完成下列操作的命令假设： A=rand(5,5) B=A(2:5,[1,3,5]) A(7)=[] A+30 size(A), ndims(A) t=0:0.1:1,t(find(t==0))=eps x=1:12, reshape(x,3,4) abs('b') char(89)
西安建筑科技大学陈
isinf(A)
ans = 0011000
isfinite(A)
ans = 1000111
7. student(1).number=123;student(1).name='lili'; student(1).specialty='computer'; student(1).grade=[87 68 70 70 97 95]; student(2).number=123;student(2).name='lili'; student(2).specialty='computer'; student(2).grade=[87 68 70 70 97 95]; student(3).number=123;student(3).name='lili'; student(3).specialty='computer'; student(3).grade=[87 68 70 70 97 95]; student(4).number=123;student(4).name='lili'; student(4).specialty='computer'; student(4).grade=[87 68 70 70 97 95]; student(5).number=123;student(5).name='lili'; student(5).specialty='computer'; student(5).grade=[87 68 70 70 97 95]; student

Matlab矩阵分析与处理-刘卫国第二版-实验二

程序：
运行结果：
本题中条件数比较大所以跟别元素的微小变动会因起解向量的很大变化，验证了（2）的结果。
6．比较F1=sqrtm（A）和F2=sqrt（A）
程序：
运行结果：
通过分析可得：sqrtm()为对矩阵的平方根。Sqrt（）为对矩阵中每个元素的平方根。
（3）计算系数矩阵A的条件数并分析结论。
6，建立A矩阵，是比较sqrtm（A）和sqrt（A），分析他们的区别。
三、实验结果
1．程序：
运行结果：
D等于E即矩阵H和5阶帕斯卡矩阵P
程序：
运行结果：
（2）其行列式的值Hh和Hp
程序：
运行结果：
（3）条件数Th和Tp
程序：
运行结果：
因为帕斯卡矩阵的条件数较小，且其值为1所以综合考虑帕斯卡矩阵的性能较好。
3。在本题中我建立5阶的帕斯卡矩阵。
求值C1，迹C2，秩C3和范数C4
程序：
运行结果：
4．A的特征值D1及特征向量D2；
程序：
运行结果：
求得三个特征值为：-3595/142，-3755/357和4697/279。
各特征值对应的特征向量为D2的各列构成的向量。
分析数学意义通过验证结果，D3=A*D2和D4=D2*D1的值相等：
程序：
运行结果：
5．（1）求方程的解。
程序：
运行结果：
（2）将方程右边向量元素改为0.53，再求解，并比较的变化和解的相对变化。
程序：
运行结果：
（3）计算系数矩阵A的条件数E并分析结论。
数学实验报告
实验名称：Matlab矩阵分析与处理
实验项目：Matlab矩阵分析与处理
专业班级：

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

实验一 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=0:0.5:2.5解：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]解：运算结果：3. 设有矩阵A和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 +⎡⎤=⎢⎥⎣⎦。

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

为什么？解：M文件如下：输出结果：因为它们的条件数Th>>Tp,所以pascal矩阵性能更好。

最新MATLAB程序设计与应用第二版刘卫国实验部分

M A T L A B程序设计与应用第二版刘卫国实验部分实验一1.（1）z1=2*sin(85*pi/180)/(1+exp(2))（2）>> x=[2,1+2i;-0.45,5];>> z2=(1/2)*log(x+sqrt(1+x*x))(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.*t).*(t>=0&t<1)+(t.*t-1).*(t>=1&t<2)+(t.*t-2.*t+1).*(t>=2&t<3)2.(1)>> A=[12,34,-4;34,7,87;3,65,7];>> B=[1,3,-1;2,0,3;3,-2,7];>> A+6*B(2)>> A*B(3)>> A^3(4)>> A/B(5)>> [A,B]3.(1)>>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(2)>> D=C(3:end,2:end)(3)Name value classA <5×5double> doubleB <5×3double> doubleC <5×3double> doubleD [520 397;705 557;…] double4.(1)>> X=100:999;>> length(find(~rem(X,21)))(2)>> y='ahHGJjskjJjsdsaSa';>> M=y(find(y<='z'&y>='a'))实验二实验三1 x=input('请输入x的值')if x<0&x~=-3y=x^2+x-6;elseif x>=0&x<5&x~=2&x~=3y=x^2-5*x+6;elsey=x^2-x-1;endy2（1）x=input('请输入成绩')if x>=90&x<=100disp('等级为A')elseif x>=80&x<=89disp('等级为B')elseif x>=70&x<=79disp('等级为C')elseif x>=60&x<=69disp('等级为D')elseif x>=0&x<=59disp('等级为E')elsedisp('成绩有误')end(2)x=input('请输入成绩')switch(x)case{90,100}disp('等级为A')case{80,89}disp('等级为B')case{70,79}disp('等级为C')case{60,69}disp('等级为D')case{0,59}disp('等级为E')otherwisedisp('成绩错误')end>>3Num=input('输入员工的工号:');Time=input('该工号员工的工作时数:'); salary=input('该员工底薪:');if Time>=120salary=salary*(1+0.15);elseif Time<=60salary=salary-700;elseif Time>60&Time<120salary=Time*84;endsalary=salary实验四1：程序：n=input('输入n:');for a=1:n;b(a)=1/(a^2);endpi=sqrt(6*sum(b))2：.程序： sum=0;for n=1:200;b(n)=1/(2*n-1);sum=sum+b(n);if(sum>=3)break;endendn-14.程序：f(1)=1;max=0;f(2)=0;min=0;f(3)=1;c=0;zhengshu=0;fushu=0;zero=0;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3);endfor m=1:100if(f(m)>max)max=f(m);endif(f(m)<min)min=f(m);endc=sum(f);if(f>0)zhengshu=zhengshu+1;endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;endif(f(m)>0)zhengshu=zhengshu+1;endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;enddisp('最大值')maxdisp('最小值')mindisp('和')cdisp('正数的个数')zhengshudisp('负数的个数')fushudisp('零的个数')zero实验五1，函数文件：文件名：Yunsuan程序：function[z,d,s,c]=yunsuan(a,b)z=exp(a+b*i);d=log(a+b*i);s=sin(a+b*i);c=cos(a+b*i);disp(['复数',num2str(a),'+',num2str(b),'i的指数为:',num2str(z),'，对数为：',num2str(d),'，正弦为：',num2str(s),'，余弦为：',num2str(c)])函数调用：程序：clearx=input('请输入复数实部的值');y=input('请输入复数虚部的值');[x,y]=yunsuan(x,y);2函数文件：文件名：zhu程序：function[p,q,w]=zhu(p,q,w)A=[p*cos(w),-p,-sin(w),0;p*sin(w),0,cos(w),0;0,q,-sin(w),0;0,0,-cos(w),1];B=[0;9.8*p;0;9.8*q];X=A\B函数调用：程序：clearm1=input('m1=');m2=input('m2=');th=input('theta=');theta=pi*th/180;[m1,m2,theta]=zhu(m1,m2,theta);4函数文件：文件名：hanshu程序：function[y]=hanshu(x)y=1/((x-2).^2-0.1)+1./((x-3).^4+0.01);实验六1、程序：、x=0:0.02*pi:2*pi;y=(0.5+3*sin(x)./(1+x.*x)).*cos(x)plot(x,y)2、（1）程序：x=(0:0.02*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'y-'),hold onplot(x,y2,'c-'),hold onplot(x,y3,'k-'),hold off（2）程序：x=(0:0.02*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1)stairs(x,y1);plot(x,y1,'r-'),hold onsubplot(2,2,2)stairs(x,y2);plot(x,y2,'b-'),hold onsubplot(2,2,3)stairs(x,y3);plot(x,y3,'k-'),hold off （3）（1）绘制条形图程序：x=(0:0.2*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1)stairs(x,y1);bar(x,y1)subplot(2,2,2)stairs(x,y2);bar(x,y1)subplot(2,2,3)stairs(x,y3);bar(x,y3)（2）绘制梯形图程序：x=(0:0.05*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1)stairs(x,y1);subplot(2,2,2)stairs(x,y2);subplot(2,2,3)stairs(x,y3);（4）杆图：程序：x=(0:0.05*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1)stem(x,y1);subplot(2,2,2)stem(x,y2);subplot(2,2,3)stem(x,y3);（5）填充图程序：x=(0:0.05*pi:2*pi)';y1=x.*x;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1)fill(x,y1,'r');subplot(2,2,2)fill(x,y2,'b');subplot(2,2,3)fill(x,y3,'g');3、：程序：x=linspace(-5,5,1000);y=[];for x0=xif x0<=0y=[y,(x0+sqrt(pi))./exp(2)];elsey=[y,0.5*log(x0+sqrt(1+x0.*x0))];endendplot(x,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)5、程序：x=-linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y)z=cos(x).*cos(y).*exp(-0.25*sqrt(x.^2+y.^2))subplot(1,2,1);title('surf(x,y,z)');surf(x,y,z);subplot(1,2,2);title('surfc(x,y,z)');surfc(x,y,z);实验八1、程序：x=rand(1,30000);EX=mean(x) %返回向量x的算术平均值S=std(x) %返回向量x的标准方差MAX=max(x) %返回向量x的最大元素MIN=min(x) %返回向量x的最小元素p=sum(x>0.5)/30000 %满足条件的百分比2、程序：p=100*rand(100,5);[MAXOBJECT,U]=max(p) %求每门课的最高分,U记录最大元素的行号[MAXOBJECT,V]=min(p) %求每门课的最低分,V记录最小元素的行号EX=mean(p) %每门课的平均分S=std(p) %每门课的标准方差SUM_1=sum(p,2) ; %求每名学生的总成绩[MAXSUM,W]=max(SUM) %求五门课总分的最高分,W记录最高分学生的序号[MINSUM,Y]=min(SUM) %求五门课总分的最低分，Y记录最低分学生的序号[zcj,xsxh]=sort(SUM_1);zcj(end:-1:1);xsxh(end:-1:1);实验九2、程序：I1=quad('sqrt(cos(t.^2)+4*sin(2*t).^2+1)',0,2*pi)I2=quad('(log(1+x))./(1+x.*x)',0,1)3、程序：A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2]b=[-4;13;1;11]x=inv(A)*b实验十1、程序：x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y))2、分解因式（1）程序：syms x y;A=x^4-y^4;factor(A)（2）程序：factor(sym('5135'))3、化简表达式（1）程序：syms beta1 beta2y=sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2) simple(y)（2）程序：syms xy=(4*x^2+8*x+3)/(2*x+1)simple(y)5、用符号方法求下列极限或导数（1）程序：syms xf=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(sin(x)) limit(f)（2）程序：syms xy=(sqrt(pi)-sqrt(acos(x)))/(sqrt(x+1));limit(f,x,-1,'right')（3）程序：syms xy=(1-cos(2*x))/x;y1=diff(y)y2=diff(y,x,2)6、用符号方法求下列积分（1）程序：syms xf=1/(1+x^4+x^8)int(f)（2）程序：syms xf=1/(((asin(x))^2)*sqrt(1-x^2)) int(f)（3）程序：syms xf=(x^2+1)/(x^4+1)int(f,x,0,inf)（4）程序：syms xf=exp(x)*(1+exp(x))^2y=int(f,x,0,log(2))double(y)。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

（6）.将含有12个元素的向量x转换为矩阵
MATLAB程序;
(7求一个字符串的ASC值
MATLAB程序;
(8).求一个ASC值对应的字符
MATLAB程序;
4.下列命令执行后L1，L2，L3,L4的值分别是多少？
MATLAB程序：
5.已知矩阵，完成下列操作。
（1）。取出A的前3列构成矩阵B，前两列构成矩阵C，右下角3乘2子矩阵构成矩阵D，B与C的乘积构成矩阵E
MATLAB程序：
(2).分别求E<D,E&D,E|D,~E|,~D和find(A>=10&A<25)
MATLAB程序;ቤተ መጻሕፍቲ ባይዱ
第二章作业
3.写出完成下列操作的命令。
（1）将矩阵A第2~5行中第1，3，5列元素赋给矩阵B。
MATLAB程序：
（2）删除矩阵A的第七号元素。
MATLAB程序：
(3).将矩阵A每个元素值加30
MATLAB程序;
(4).求矩阵A的大小和维数
MATLAB程序：
（5）.将向量t的0元素用机器零来代替。
MATLAB程序：