matlab作

文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.习题一12B1、试个MATLAB的工作空间中建立以下2个矩阵：A＝[12]34，求出矩阵A和B的乘积，并将结果赋给变量C。>>A=[12];B=[12;34];>>C=A*BC=7102、利用MATLAB提供的帮助信息，了解inv命令的调用。>>helpinv3、使用help命令查询函数plot的功能以及调用方法，然后利用plot命令绘制函数0xy=sin(x)的图形，其中.>>helpplot>>x=0:pi/1000:pi;>>y=sin(x);>>plot(x,y)>>xlabel('x');>>ylabel('y=sin(x)');>>title('正弦函数');4、试用不同的方法建立数组A＝[11.52.02.53.0]，了解怎样访问数组A的第二个元素，然后将其更换为4.0。建立数组:>>A＝[11.52.02.53.0]>>A=1:0.5:3.0访问元素：>>A(2)ans=1.5000更换元素：>>A(2)=4.0A=1.00004.00002.00002.50003.000012B5、已知矩阵34，试用MATLAB提供的关系运算命令将B中所有大于2的元素全改为0。第一种：>>B=[12;34];>>B(find(B>2))=0B=12001文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.第二种：>>B=[12;34];>>form=1:2forn=1:2ifB(m,n)>2B(m,n)=0;endendend>>BB=1200123A4566、已知矩阵789，试求矩阵A的左右翻转矩阵，上下翻转矩阵，然后在工作空间中利用size命令查看矩阵A的大小。(1)生成矩阵A>>A=[1:9];>>A=reshape(A,3,3)A=147258369(2)左右翻转：>>fliplr(A)ans=741852963(3)上下翻转：>>flipud(A)ans=369258147(4)查看矩阵A的大小：>>size(A)ans=332文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.12B7、已知矩阵34，试求其转置、逆、迹、特征值、特征向量和B对应的行列式的值。>>B=[1:4];B=reshape(B,2,2)%生成矩阵BB=1324>>B'%矩阵转置ans=1234>>inv(B)%矩阵的逆运算ans=-2.00001.50001.0000-0.5000a=trace(B)%矩阵的迹a=5>>[x,y]=eig(B)%矩阵的特征值y，特征向量xx=-0.9094-0.56580.4160-0.8246y=-0.3723005.3723>>V=det(B)%矩阵的行列式V=-28、分别建立一个33阶的单位阵、随机阵和魔方阵。>>A=eye(3)%生成3*3阶的单位矩阵AA=100010001>>B=rand(3)%生成3*3阶的随机矩阵BB=0.65960.64900.43240.51860.80030.82530.97300.45380.0835>>C=magic(3)%生成3*3阶的魔方矩阵CC=3文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.816357492ax()x22x2b()xx3x23x19、已知多项式，。试求两个多项式的和与乘积。>>A=[12-2];B=[11-31];>>a=poly2sym(A),b=poly2sym(B)a=x^2+2*x-2b=x^3+x^2-3*x+1>>c=a+bc=x^3+2*x^2-x-1>>d=conv(A,B);poly2sym(d)ans=x^5+3*x^4-3*x^3-7*x^2+8*x-2iz34i,z12i,z2e610.复数123表达，及计算>>zzz12z3。>>z1=3+4*i;z2=1+2*i;z3=2*exp(pi/6*i);>>z=z1*z2/z3z=0.3349+5.5801i11、产生1×5的均布随机数组，进行如下操作：1）寻访数组的第三个元素；2)寻访数组的第一、二、五个元素组成的子数组；3）寻访前三个元素组成的子数组；4）寻访除前2个元素外的全部其他元素。>>rand('state',0)%产生的随机数都与第一次运行产生的相同>>A=rand(1,5)%产生1×5的均布随机数组A=0.95010.23110.60680.48600.8913>>A(3)%寻访数组x的第三个元素ans=0.6068>>A([125])%寻访数组x的第一、二、五个元素组成的子数组ans=0.95010.23110.8913>>A(1:3)%寻访前三个元素组成的子数组ans=4文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.0.95010.23110.6068>>A(3:end)%寻访除前2个元素外的全部其他元素ans=0.60680.48600.891312、试用两种方法用MATLAB计算38第一种：>>nthroot(-8,3)ans=-2第二种：>>symsxsolve('x^3+8',x)ans=-21-3^(1/2)*i1+3^(1/2)*i第三种：>>x=-8;sign(x)*abs(x).^(1/3)ans=-2(s22)(s4)(s1)13、求s3s1的“商”及“余”多项式。>>a=[102];b=[14];c=[11];d=[1011];e=conv(a,b);f=conv(e,c);>>a=[102];b=[14];c=[11];d=[1011];>>e=conv(a,b);f=conv(e,c);>>[k,r]=deconv(f,d)%k是商，r是余式k=15r=0054314、求方程x^4+7x^3+9x-20=0的全部根。第一种方法：>>p=[1709-20];roots(p)ans=-7.2254-0.4286+1.5405i-0.4286-1.5405i5文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.1.0826第二种方法：>>compan(p)ans=-70-920100001000010>>eig(ans)ans=-7.2254-0.4286+1.5405i-0.4286-1.5405i1.082615、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。>>a=[100-21];b=[14-0.5];>>p1=poly2sym(a);p2=poly2sym(b);>>p=p1+p2p=x^4+x^2+2*x+1/2>>P=sym2poly(p);>>i=0:5;xi=0.2*i;>>polyval(P,xi)ans=0.50000.94161.48562.18963.14964.500016、已知一线性方程组如下所示：3xxx3.6123x2x4x2.1123x4x5x1.4123，试求其结果。>>A=[31-1;124;-145];B=[3.62.1-1.4];B/Aans=1.4818-0.46060.3848习题二1、编制一个函数，使得该函数能对输入的两个数值进行比较，并返回其中的最小值。functiony=min(m,n)ifm1)error('参数太多');endlen=length(x);fori=1:len/2tmp=x(i);x(i)=x(len-i+1);x(len-i+1)=tmp;endx3、编制一个m程序，计算阶乘n！=1×2×3×…×nfunctiony=jiecheng(N)ifN<0error('Nshouldbepositiveorzero');endifN==0|N==1y=1;elsey=N*jiecheng(N-1);end4.利用循环语句进行程序设计：假设定义m×n的矩阵A。判断矩阵A的第1列元素是否为0，若全为0，则从矩阵A中删除第1列functionB=liu(A)[m,n]=size(A);b=0;a=1;whilea<=mifA(a,1)==0b=b+1;enda=a+1;endifb==m7文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.A(:,1)=[];B=A;elseB=A;end5、利用循环语句进行程序设计：在区间[－2,－0.75]内，步长为0.25，对函数y=f(x)=1+1/x求值，并列表。将所得x值和y值分别存入向量r和s中。function[r,s]=funa=0;x=-2;whilex<=-0.75y=1+1/x;a=a+1;r(a)=x;s(a)=y;x=x+0.25;end运行结果：>>[r,s]=funr=-2.0000-1.7500-1.5000-1.2500-1.0000-0.7500s=0.50000.42860.33330.20000-0.333363k2i6、编程计算i0编制m程序：functionk=fun1(x)k=0;fori=0:xk=k+2^i;end计算：>>fun(63)ans=1.8447e+0197.写一个MATLAB小程序findN01.m，求出最小的n值，使得n!>realmax。请问n的值是多少？此时(n-1)!的值又是多少？functiony=findN01(realmax)maxN=1000;forn=1:maxN8文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.value=prod(1:n);ifvalue>realmaxbreak;endendfprintf('n=%d\n',n);fprintf('(n-1)!=%d\n',prod(1:n-1））；习题三1、用subplot命令在同一图形输出窗口中绘制以下4个函数的图形：yxx,[0,3]yxsinx,x[1,1]yx2,x[0,1.5]，，，ytanx,x[0,1.3]。>>clear>>subplot(2,2,1);>>x=1:0.001:3;>>plot(x,x,'-r*');>>subplot(2,2,2);>>x=-1:0.001:1;>>plot(x,x.*sin(x),'-r*');>>subplot(2,2,3);>>x=0:0.001:1.5;>>plot(x,x.^2,'-r*');>>subplot(2,2,4);>>x=0:0.001:1.3;>>plot(x,tan(x),'-r*');9文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.ye0.2xsinxx[0,5]2、绘制曲线在区间上的阶梯图。>>x=0:0.01*pi:5*pi;y=exp(-0.2*x).*sin(x);plot(x,y);axis([0,0.2*pi,0,0.2*pi]);holdonstairs(x,y,'r');51cos[0,8]3、试绘制以极坐标形式表示的图形：43，其中的范围为。>>x=0:0.001*pi:8*pi;>>p=cos(5*x/4)+1/3;>>polar(x,p,'-r')ttye3sin3tye3t4、画出衰减振荡曲线及其它的包络线0。的取值范围是[0,4]。>>t=0:0.001*pi:4*pi;10文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.y0=exp(-t/3);>>y=y0.*sin(3*t);>>plot(t,y0,'-r',t,y,'-b');>>legend('包络线','衰减振荡曲线');sin(x2y2)zx2y2x,y5、画出所表示的三维曲面。的取值范围是[8,8]。>>x=-8:0.001:8;>>y=-8:0.001:8;>>m=sqrt(x.^2+y.^2);>>z=sin(m)./m;>>plot3(x,y,z,'-r');6、在[02π]范围内绘制二维曲线图y=sin(x)*cos(5x)。>>x=0:0.01*pi:2*pi;>>y=sin(x).*cos(5*x);>>plot(x,y,'-rx')。7、在[02π]范围内绘制以Y轴为对数的二维曲线图。11文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.y=|1000sin(4x)|+1>>x=0:0.001*pi:2*pi;>>y=abs(1000*sin(4*x))+1;>>semilogy(x,y,'-*r')8、绘制z=sin(x)*cos(y)的三维网格和三维曲面图，x,y变化范围均为[02π]。>>x=[0:0.1:2*pi];>>y=[0:0.1:2*pi];>>z=sin(y')*cos(x);>>mesh(x,y,z);>>xlabel('x-axis'),ylabel('y-axis'),zlabel('z-label');>>x=0:0.1:2*pi;>>y=0:0.1:2*pi;>>[x,y]=meshgrid(x,y);>>z=sin(x).*cos(y);>>surf(x,y,z);>>xlabel('x-axis'),ylabel('y-axis'),zlabel('z-label');29、用简短的MATLAB命令在一个图上绘制在0≤x≤7范围内的sin(2x)、sinx和cosx2三条曲线，并将其一一标明。>>clear>>x=0:0.01*pi:7;>>y=sin(2*x);>>z=sin(x.^2);>>m=cos(x.^2);>>plot(x,y,'-r*')>>holdon>>plot(x,z,'-.bo')>>plot(x,m,'--')>>legend('sin(2x)','sin(x^2)','cos(x^2)')>>holdoff10、在极坐标系中绘制曲线>>t=0:0.01*pi:8*pi;>>r=exp(cos(t))-2*cos(4*t)+sin(t/12).^5;>>polar(t,r,'-r')习题四aaA1112aa1、求矩阵2122的行列式值、逆和特征根>>symsa11a12a21a22;>>A=[a11a12;a21a22];12文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.>>B=det(A)B=a11*a22-a12*a21>>inv(A)ans=[a22/(a11*a22-a12*a21),-a12/(a11*a22-a12*a21)][-a21/(a11*a22-a12*a21),a11/(a11*a22-a12*a21)]>>[x,y]=eig(A)x=[(a11/2+a22/2-(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)/2)/a21-a22/a21,(a11/2+a22/2+(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)/2)/a21-a22/a21][1,1]y=[a11/2+a22/2-(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)/2,0][0,a11/2+a22/2+(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)/2]sin/22AeitdtA/22、验证积分2。用matlab计算左边的：>>symsAwtm;>>Y=int('A*exp(-w*t*i)',t,-1*m/2,m/2)%不知道在matlab中如何输入在这里就用m表示Y=(2*A*sin((m*w)/2))/w右边化简即得证。t11(1)ktk3(2k1)2k3、求t0，k1>>symstk;>>s=symsum([t,k^3],t,0,t-1)s=[(t*(t-1))/2,k^3*t]>>symk;>>s=symsum([1/(2*k-1)^2(-1)^k/k],k,1,inf)s=13文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.[pi^2/8,-log(2)]2dat3dat324、求dxtcosxlnx、dttcosxlnx>>symsaxt;>>A=[a,t^3;t*cos(x),log(x)];>>diff(A,'x')%微分一次ans=[0,0][-t*sin(x),1/x]>>diff(A,'x',2)%对表达式A微分2次ans=[0,0][-t*cos(x),-1/x^2]axbx21dxsinx5.求x>>symsabx;>>A=[a*x,b*x.^2;1/x,sin(x)];>>int(A,x)ans=[(a*x^2)/2,(b*x^3)/3][log(x),-cos(x)]1xdt6、求0lnt>>symsx>>y=int(1/log(t)','t',1,x)>>diff(y,x)npndq,ndqpq10,dp,qpn8d17、求224线性方程组的解。>>A=[1,1/2,1/2,-1;11-11;1-1/4-11;-8-111];>>B=[01001]';>>A\Bans=1.000014文档收集于互联网，如有不妥请联系删除.文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.8.00008.00009.00008、求(x2)x2的解。>>symsx;y=(x+2).^x-2;S=solve(y,x)S=matrix([[0.1]])uy^2vzw0y,z9、求方程组，yzw0关于的解。>>symsuyvzw;>>m=u*y^2+v*z+w;>>n=y+z+w;>>q=solve(m,n,'y','z')q=y:[2x1sym]z:[2x1sym]>>q.yans=(v+2*u*w+(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)-w(v+2*u*w-(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)-w>>q.zans=-(v+2*u*w+(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)-(v+2*u*w-(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)15文档收集于互联网，如有不妥请联系删除.