matlab复习资料 16页

习题P771、计算y=x3%P77_3题a=3;A=4;b=a*a;B=b*b-1;c=a+A-2*B;C=a+2*B+c8、创建一个5*5随机阵并求其逆M文件程序：inv%P77_8题%创渣一a=rand(5)+(x-0.98)2/(x+1.35)3-5(x+l/x),当x=2和x=4时的M文件程序：%P77_1题x=[2,4];y=x.A3+(x-0.98).A2./(x+1.35).A3-5.*(x+l./x)2、计算cos6(T-V9-V^的值。M文件程序:%P77_2题cos(60)-(9-sqrt(2))A(l/3)3、已知a=3,A=4,b=a*a,B=b*b-1，c=a+A-2B，C=a+2B+c，求C.M文件程序：10、设A=[1，4,8,13;-3,6,-5,-9;2,-7,-12,-8];B=[5,4,3,—2;6,-2,3,一8;-1，3,—9,7];求C=A*B,D=A.*B。M文件程序:%P7710题\nCommandWindow???Errorusing==>mtimesInnermatrixdimensionsmustagree.1部矩阵维度必须一致Errorin==>P77」0at5C=A*B»(3Editor•E:VS^\MATLAB\新建文件夹\P77J0.mFileEditTextCellToolsDebugDesktopWindowHelpd运曇脚九趋粕邐鴒~奸77-10题2-clear;clc;3-A=[1J4J8,13;-3,6,-5,-9;2,-7,-12,-8];4-B=[5,4,3,-2:6,-2,3,-81,3,-9,7]:5-C=A*b|6-D=A.*B11、设y=cosx[0.5+Jsin:']，把x=0-2*pi区间分为125点，画出以x(l+r)为横坐标，y为纵坐标的曲线。M文件程序：%P77_11x=0:(2*pi+l)/125:2*pi;y=cos(x).氺(0.5+3.*sin(x)./(1+x.2));plot(x,y);xlabelCx’);ylabel(’y’)；gridon;12、设zsin3z,y=zcos3z，要求在z=-45-45区间内画出X，y，z三维曲线。M文件程序：%P77_12z=-45:1:45;x=z.本sin(3.*z);y=z.本cos(3.*z);plot3(x,y,z);xlabelCx’)；ylabel(’y’)；zlabel(’z’)；gridon;13、设z=x2,2”\求定义域又=[-2,2]，尸[-2,2]内的2值，并绘出三维曲M文件程序：\n%P77_13clear;closcall;clc;x=-2:0.1:2;y=x;[X,Y]=meshgrid(x,y);Kl.*exp(-X."2-Y.^2);surf(X,Y,Z);\nP123j、在某处测得海洋不同深度处的水温如下:深度(m)44671495014221634水温7.044.283.402.542.13M文件程序：%P122_1x:[500,1000,1500];shcn=[446,714,950^143271^34];利用分段线性插值函数，求在深度为500m、1000m、1500m处的水温wen:[7.04^-28X40,2.54,2.13];y=interpl(shen,wen,x,5氺linear’plot(shen,wen,’—k，，x,y,’氺r’一维插倍函数P1234、求代数方程3x5+4x4+7x3+2f+9%+12=0的根。M文件程序：%P123_4clear;clc;closeall;p=[3472912];r=roots(p)%r=roots(p)多项式求根。%验证：将根代入多项式，其值应该为0y=polyval(p,r)K=Polyder(p):多项式求导数X=Pby(A):按根组成多项式Polyint:多项式积、分Y=polyval(p，x):求多项式P中x处的值。P123_5、设方程的根为x=[-3,-5,-8,-9]，求他们对应的x多项式的M文件程序：%P123_5clc;clear;closeall;x=[-3,-5,-8,-9];y=poly(x)去均值：x-mean(x)归一化：x/max(abs(x))P123_8、设对称实矩阵a=[2,4,9;4,2,4;9,4,18]，求其特征根和特征向量M文件程序：%P123_8clc;clear;closeall;八二[2,4,9;4,2,3;9,4,18];多项式相乘：w=conv(a,b)多项式相除：[q,r]=deconv(a,b)\n[V，D]=eig(A)%V为A相应的特征向量，D为A的特征根。3.线性代数方程求解一般线性方程组：Qanxx+a,2x2+^^axnxxn^bx+a22X2+Ha2nxin-^2+am2X2+…+amA=bm□□□写成矩阵形式可表示为：AX=B或XA=B。其中系数矩阵A的阶数为mxn。在MATLAB中，引入矩阵除法求解。(1)求解方程AX#—►，很重要格式：X=A\Bt左除_条件：矩阵A与!行数b须相奪。(2)求解方程XA=pJ1棬式：X=B/AJE^^S<=B*inv(A))条件：矩阵A与矩阵须相等。P123J0、求函数f(x，y)=4(%-y)-x2->，的极值。M文件程序：%P123_10functionf=myfun(v)clear;clc;closeall;x=v(l);y=v(2);f=4.*(x-y)-x.八2-y.八2;f=-f;单变量函数求极小x=fminbnd(‘fun',xl,x2)多变量函数求极小rsx=fminsearch(‘fun’，v)单变量函数求零点：x=fzero(‘fun',xO)»v=[-2,2]；>>vmin=fminsearch(’myfun、v)vmin=FileEditTextCellToolsDebug□(多■|其1functionf=myfun(v)2%J?lear;clc;closeall;32.0000-2.0000»运行〉〉4-x=v(l);y=v(2):5-f=4.*(x-y)-x."2-y."26-f=-f;7P123_ll、设siny+ex-a^2=0，求一。dxM文件程序：\nSymsxy;F=sin(y)+exp(x)-x*yA2;\nFy=diff(F，y);Fx=diff(F，x);dydx=-Fx/FyP123_12、计算二重积分S=£n£2n(x-y)2sin2(x+oM文件程序：s=dblquad_(inline(’(x-y).八2.*(sin(x+y)).八2’)，0，pi，pi，2*pi)定义一个函数二重积分函数，13、设(x，Y)的概率密度为求E(X)，E(Y)，E(XY)。M文件程序：%pl23_13根据概率密度求期望clear;clc;closeall;f二inline(’12*y."2’)；delta=0.001;%当delta=0.001吋，误差很小x=0:delta:1-delta;y=x；len=length(x);z=zeros(len);fori=l:len%x下标forj=l:i%y卜•标，上下限：0:xz(j,i)=f(y(j));endendmesh(x，y,z)，xlabel(’x’)，ylabel(’y’)，zlabel(’f’)%验证sum(sum(f(x)dxdy))==1f_sum=sum(sum(z))^delta2%当delta=0.001时的结果为1.0000%求概率密度函数f(x)，对dy求积分\nf_x=sum(z)*de1ta;figure,plot(x,f_x)holdonplot(x,4*x."3,’*’)％与理论计算相比较%求期望E(X)E_X=sum(f_x.*x)*delta%理论值为0.8,当delta=0.001时的结果为0.7995%求概率密度函数f(y),对dx求积分f_y=(sum(z，2)*delta)’；figure,plot(y,f_y)holdonplot(y,12*y,2.*(l-y),’*’)%与理论计算相比较%求期望E(Y)E_Y=sum(f_y.*y)*delta%理论值为0.6，当delta=0.001时的结果为0.6000%求期望E(XY)[X,Y]=meshgrid(x,y);XY=X.*Y;E_XY=sum(sum(z.*XY))*delta"2%理论值为0.5，当delta=0.001时的结果为0.4997P182j、若4n)=cos是一个N=12的有限长序列，计算它的DFTN=12;tax=0.01;n=0:tax:N_tax;y=cos(n*pi/6);plot(n,y);\nP182一2、己知两序列x⑻=df11,h⑻=求两序列的线性卷积。M文件程序：%P182_3and4nx=12;m=5;n=O:nx~l;xl二0.8*ones(1,nx);subplot(3,1,1),stem(n,xl);nh=6;nl=0:nh—l;hl=ones(1,nh);subplot(3,1，2)，stcm(nl,hl);yl=conv(xl，hl);ny=(n⑴+nl(1)):(n(end)+nl(end));subplot(3,1,3),stem(ny,yl);xlabel(’n’)，ylabcl(’y(n)’),title(’x(n)*h(n)’);\n1若x(n)=cos3己知两序列为|«4.59、屯网去哦处押£一个例2的娜序列.计算它的D打并画出狐x(n)0.80^w^ll0其他h(n)10其他求两序列的线性卷积•4用FFT实现h题中两序列的线性卷积2+3z-15求传递函数"(z)6求传递函数H(z)1+0V+4?+15.6/+6:2+2.4:-6.43z4+2.4z;+6.3z2-11.4z+6的因式形式，并■出；极点图1826?n&b“对传递函数HC?)18”一6TiTiPTrjpTp进⑽分分式族开.8存一模拟城波器*⑽递峨如下，軸它的■和相報曲线.S2+0.4J+1设让个10阶的带通BuUerworth滤波器，它的通带范啪为⑽〜200%并®出它的单位冲激啗应.用双线性变换法没计一个Butterworth低通滤波器，要求狂通带截止頻率为100Hz，阻带极止频举为2(^Hz，通带哀减Rp<2dB，阻带哀减Rs>15dB.采样频率Fs=500Hz.0没计一个阶数为48.通带范围为0.35