matlab与在电子信息课程中应用课后习题答案第四章.docVIP

clc,clear,clf,close all; disp(题解4.5,4.6,4.7);pause; echo on; a=[3 4 7 2 9 12], roots(a), pause, b=[1 0 0 0 0 -1], roots(b), pause, c=[-3 -5 -8 -9], poly(c), echo off; disp(题解4.8);pause; echo on; a=[1 2 5 4 3],b=[1], [r,p,k]=residue(b,a) pause, t=0:0.2:10; ys=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t); h=plot(t,ys);set(h,linewidth,2), echo off; disp(题解4.9,4.10);pause; echo on; R1=randn(8,6), mean(R1), std(R1), pause, mean(R1(:)), std(R1(:)), pause, R=ceil(16*rand(4,6)), a=R(:,[1,2,3,4]), inv(a), echo off; disp(题解4.11); echo on; t=[0:0.05:10]; x1=2.*cos(t)+3.*t;y1=2.*sin(t)+3; plot(x1,y1),shg,hold on; x2=3.*cos(t)+3.*t;y2=3.*sin(t)+3; plot(x2,y2),shg x3=4.*cos(t)+3.*t;y3=4.*sin(t)+3; plot(x3,y3),shg hold off; pause; figure; plot(x1,y1,x2,y2,x3,y3) echo off; disp(题解4.12); echo on; t=[-2*pi:0.05:2*pi]; a=0; x=sin(t); y1=sin(t+a); y2=sin(2.*t+a); y3=sin(3.*t+a); y4=sin(4.*t+a); subplot(2,2,1),plot(x,y1), subplot(2,2,2),plot(x,y2), subplot(2,2,3),plot(x,y3), subplot(2,2,4),plot(x,y4), pause; y11=sin(2.*t); y22=sin(2.*t+pi/3); y33=sin(2.*t+pi/2); y44=sin(2.*t+pi); figure; subplot(2,2,1),plot(x,y11), subplot(2,2,2),plot(x,y22), subplot(2,2,3),plot(x,y33), subplot(2,2,4),plot(x,y44), echo off; disp(题解4.13);pause; echo on; p=[1 -4 0 3 -2 6], x=linspace(-2,8,100); y=polyval(p,x); plot(x,y), pause axis([-2 8 -20 20]) pause, roots(p), echo off; disp(题解4.14) echo on; z=[0:0.05:10]; x=z.*sin(3.*z);y=z.*cos(3.*z); plot3(x,y,z); pause; figure,plot3(x,y,z,r:) pause; xlabel(x),ylabel(y),zlabel(z), echo off; disp(题解4.15 4.16); echo on; [x,y]=meshgrid(-2:0.1:2,-2:0.1:2); z=x.^2.*exp(-(x.^2+y.^2)); mesh(z);hold on; pause; z1=0.05.*x-0.05.*y+0.1; mesh(z1); echo off; for i=1:41 for j=1:41 if z(i,j)==z1(i,j) a=z(i,j); plot3(i,j,a,r+); end end end hold off; disp(题解4.18);pause; type f1.m i=1; a(1)=fzero(f1,-2); for m=-2:1:8 b=fzero(f1,m); if abs(a(i)-b)0.00001 i=i+1; a(i)=b;