【免费下载】matlab及在电子信息课程中的应用课后习题答案第四章
《MATLAB_程序设计与应用》刘卫国高等教育出版社-第4章课后答案
第四章1.a=input('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位数:'); endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*b+100*c+10*d+e;disp(['加密后：',num2str(a)])2.逻辑表达式法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');x=0.5:1:5.5;x1=(x>=0.5&x<1.5);x2=(x>=1.5&x<3.5);x3=(x>=3.5&x<=5.5);y1=a.*(x.^2)+b.*x+c;y2=a*(sin(b)^c)+x;y3=log(abs(b+c./x));y=y1.*x1+y1.*x2+y3.*x3; disp(y)if语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5if x>=0.5 & x<1.5y=a.*(x.^2)+b.*x+celseif x>=1.5 & x<3.5y=a*(sin(b)^c)+xelseif x>=3.5 & x<5.5y=log(abs(b+c./x))endendswitch语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5switch floor(x/0.5)case {1,2}y=a.*(x.^2)+b.*x+c;case {3,4,5,6}y=a*(sin(b)^c)+x;case {7,8,9,10}y=log(abs(b+c./x)); enddisp(y)end3.x=fix(rand(1,20)*89)+10;x1=mean(x);n=find(rem(x,2)==0 & x<x1);disp(['小于平均数的偶数是：',num2str(x(n))]);4.(1)A=input('请输入20个数的一个行向量：');a=A(1);b=A(1);for m=Aif a>=ma=m;elseif b<=mb=m;endenddisp(['最小数是：',num2str(a)])disp(['最大数是：',num2str(b)])(2)A=input('请输入20个数的一个行向量：'); maxval=max(A)minval=min(A)5.s=0;for a=0:63c=2^a;s=s+c;enddisp(['2的0次方到63次方的和是：',num2str(s)])k=0:63n=2.^ks=sum(n)6.（1）sum1=0;for n=1:100x=(-1)^(n+1)*(1/n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=(-1)^(n+1)*(1/n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=(-1)^(n+1)*(1/n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（2）sum1=0;n1=0;for n=1:2:100x=(-1)^n1*(1/n);sum1=sum1+x;n1=n1+1;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;n2=0;for n=1:2:1000x=(-1)^n2*(1/n);sum2=sum2+x;n2=n2+1;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;n3=0;for n=1:2:10000x=(-1)^n3*(1/n);sum3=sum3+x;n3=n3+1;enddisp(['当n取10000时：sum=',num2str(sum3)])（3）sum1=0;for n=1:100x=1/(4^n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=1/(4^n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=1/(4^n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（4）sum1=1;for n=1:100x=4*n*n/(2*n-1)/(2*n+1);sum1=sum1*x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=1;for n=1:1000x=4*n*n/(2*n-1)/(2*n+1);sum2=sum2*x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=1;for n=1:10000x=4*n*n/(2*n-1)/(2*n+1);sum3=sum3*x;enddisp(['当n取10000时：sum=',num2str(sum3)])7.函数文件function f=fibnacci(n)if n==1 | n==2f=1;elsef=fibnacci(n-1)+fibnacci(n-2);end命令文件：shulie=[];for k=1:nshulie=[shulie fibnacci(k)];endshulie8.function [f1,f2]=juzhenji(x1,x2)f1=x1*x2;f2=x1.*x2;命令文件：clear alla=input('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)'); [f1,f2]=juzhenji(a,b);disp(f1)disp(f2)9.function sum=qiuhe(n,m)if n<=1sum=0;elsesum=n^m+qiuhe(n-1,m);end命令文件：clear ally=qiuhe(100,1)+qiuhe(50,2)+qiuhe(10,-1);disp(y)10.s=0；a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0 s=s+k(j);endendendss =108（2）global xx=1:2:5;y=2:2:6;sub(y);xyfunction fun=sub(z) global xz=3*x;x=x+z;x =4 12 20 y =2 4 6。
【免费下载】MATLAB程序设计与应用第二版第4章课后题答案
maxval=max(A) =min(A)
5. s=0; for a=0:63
c=2^a; s=s+c; end disp(['2的0次方到63次方的和是：',num2str(s)])
1
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内 纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
第4章MATLAB在信号处理中的应用_Matlab及在电子信息类课程中的应用
第 4 章 MATLAB在信号处理中的应用
4.1.3离散时间信号的表示
在 MATLAB 中，离散时间信号 x(n) 的表示：需用一个向量 x 表示序列幅 值，用另一个等长的定位时间变量n，才能完整地表示一个序列。 [例4-10] 绘制离散时间信号的棒状图。其中x(-1)=-1, x(0)=1, x(1)=2, x(2)=1, x(3)=0, x(4)=-1。MATLAB源程序为： n=-3:5; %定位时间变量 x=[0,0,-1,1,2,1,-1,0,0]; stem(n,x); grid; % 绘制棒状图 line([-3,5],[0,0]); %画x轴线 xlabel('n'); ylabel('x[n]') 运行结果如图4.11所示。
6. Chirp（产生线性调频信号函数chirp）
格式：y=chirp(t,f0,t1,f1) 功能：产生一个线性扫频（频率随时间线性变化）信号，其时间轴的设置由数组t定义。 时刻0的瞬时频率为f0，时刻t1的瞬时频率为f1。默认情况下，f0=0Hz，t1=1,f1=100Hz。 例4-7：绘制一个线性调频信号。
3．零－极点增益模型 连续系统： 离散系统：
H ( z) k
( z q1 )( z q2 )( z qM ) ( z p1 )( z p2 )( z pN )
第 4 章 MATLAB在信号处理中的应用
4．极点留数模型 连续系统： 离散系统：
H ( s) r r1 r 2 N s p1 s p 2 s pN
y
(N)
(t ) ai y (t ) bi x(i ) (t )
(i ) i 0 i 0
N 1
M
matlab课后习题答案第四章
第4章数值运算习题 4 及解答1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff或数值梯度gradient指令计算)(t y'，然后把)(t y和)(t y'曲线绘制在同一张图上，观察数值求导的后果。
（模拟数据从prob_data401.mat获得）〖目的〗●强调：要非常慎用数值导数计算。
●练习mat数据文件中数据的获取。
●实验数据求导的后果●把两条曲线绘制在同一图上的一种方法。
〖解答〗（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat（2）用diff求导的指令dt=t(2)-t(1);yc=diff(y)/dt; %注意yc的长度将比y短1plot(t,y,'b',t(2:end),yc,'r')（3）用gradent 求导的指令（图形与上相似）dt=t(2)-t(1);yc=gradient(y)/dt;plot(t,y,'b',t,yc,'r')grid on〖说明〗● 不到万不得已，不要进行数值求导。
● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级以上。
● 求导会使数据中原有的噪声放大。
2 采用数值计算方法，画出dt tt x y x ⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖提示〗● 指定区间内的积分函数可用cumtrapz 指令给出。
● )5.4(y 在计算要求不太高的地方可用find 指令算得。
〖目的〗● 指定区间内的积分函数的数值计算法和cumtrapz 指令。
● find 指令的应用。
〖解答〗dt=1e-4;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3)ii=find(t==4.5);s45=s(ii)s45 =1.65413 求函数x ex f 3sin )(=的数值积分⎰=π0 )(dx x f s ，并请采用符号计算尝试复算。
MATLAB及其应用课后习题解答
5
“MATLAB 及其应用”课程作业
图 4- 2 画出函数图像
第5章
1． 已知椭圆的长、短轴 a 4, b 2 ，用“小红点线”画如下图所示的椭圆
x a cos t 。 y b sin t
（提示：参量 t ；点的大小；axis equal） 答 ： 运 行 :clf; a=4;b=2;t=0:pi/80:2*pi;x=a*cos(t);y=b*sin(t);plot(x,y,’r’,’markersize’,15);axis equal;xlabel(‘x’);ylabel(‘y’);shg;结果如图 5-1 所示：
1
“MATLAB 及其应用”课程作业
a5=sin(sym(pi/4)+exp(sym(0.7+pi/3))); a6=sin(sym(pi/4)+sym(exp(0.7+pi/3))); a7=sin(sym(pi/4+exp(0.7+pi/3))); a8=sym(sin(pi/4+exp(0.7+pi/3))); da2 = vpa(a1-a2,30) da3 = vpa(a1-a3,30) da4 = vpa(a1-a4,30) da5 = vpa(a1-a5,30) da6 = vpa(a1-a6,30) da7 = vpa(a1-a7,30) da8 = vpa(a1-a8,30) 结果如图 2-2 所示： 只有 a2 是精准的，观察可知，这是由 于 sym('Num')和 sym(Num)的区别带来 图 2- 2 vpa 观察误差 的前者以字符串的形式传给符号运算内 核，可以保留完整的精度；而后者经过浮点数运算之后再转换为符号类型，存在精度损 失。 3．在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。 sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*th)') 答：分别为 w,a,z，具体图 2-3 所示。
MATLAB与在电子信息类课程中的应用复习题答案电子科大第二版
第二章1、x=[2,4];y=x八3+(x-0.98)八2./(x+1.35)八3-5*(x+1./x) 2、y=cos(pi/3)-(9-sqrt(2))A(1/3) 3、a=3;A=4; b=a.A2; B=b.A2-1; c=a+A-2*B;C=a+2*B+c4、x=[1 2 3;4 5 6;7 8 9]Desktop->Workspace, 双击变量x 5、clearx=magic(3) y=randn(3,3) xy=[x,y] yx=[x,y] z=xy(:,1:2) 6、clear x=eye(4,4)y=triu(x)7、clear x=rand(4,5) y=x([1,2],:)z=(y>=0.3).*y8、clear x=randn(5,5) y=inv(x)9、clear x=randn(5,5)z=x A510、clearA=[1,4,8;-3,6,-5;2,-7,-12];B=[5,4,3;6,-2,3;-1,3,-9];C=A*BD=A.*B11、clear x=linspace(0,2*pi,125);y=cos(x).*(0.5+3*sin(x)./(1+x.A2)); plot(x,y)12、z=-45:1:45;x=z.*sin(3*z);y=z.*cos(3*z);plot3(x,y,z);13、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2)surf(x,y,z);14、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2);surf(x,y,z);z1=0.05*x-0.05*y+0.1;hold on,mesh(x,y,z1);15、(1)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);(2)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);n=2;alfa=pi/3;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,2);plot(t,x,t,y);n=2;alfa=pi/2;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,3);plot(t,x,t,y);n=2;alfa=pi;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,4);plot(t,x,t,y);球体clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/100);y(i,j)=i/100*sin(j*2*pi/100);z(i,j)=sqrt(1.001-x(i,jF2-y(i,j)八2);endendsurf(x,y,z);hold on;surf(x,y,-z);axisequal棱锥clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/4);y(i,j)=i/100*sin(j*2*pi/4);z(i,j)=sqrt(x(i,j)A2+y(i,j)A2);endendsurf(x,y,-z);求最大值x=[67 87 56 99 43] max=0;for i=1:5 if max>=x(i) max=max; elsemax=x(i) endendmax求最小值x=[67 87 56 99 43] min=inf;for i=1:5if min<=x(i) min=min;else min=x(i) end end min求和x=[67 87 56 99 43]; sum=0; for i=1:5 sum=sum+x(i) end sum第三章1、h0=[446,714,950,1422,1634];t0=[7.04,4.28,3.40,2.54,2.13];t1=interp1(h0,t0,500,'linear') 2、x0=[1,0,-1];y0=[0,1,0]; p=polyfit(x0,y0,3); x=-1:0.1:1; y=polyval(p,x);plot(x,y,-x,-y),axis equal3、clear x0=[1.0,1.1,1.2,1.3,1.4];y0=[0.25,0.2268,0.2066,0.1890,0.1736];p=polyfit(x0,y0,3);p1=polyder(p); y=polyval(p1,[1.0,1.2])4、p=[3 4 7 2 9 12]; r=roots(p)5、r=[-3 -5 -8 -9]; p=poly(r)6、function ypie=fun1(x,y)ypie=x42./y-x.*cos(y);[x,y]=ode45( 'fun1' ,[0,5],1)plot(x,y)7、function y=fun2(x)y=x.A4-3*x A3+5*cos(x)+8;fplot( 'fun2' ,[1,5])hold on;x=0:5;y=0*x;plot(x,y)a=fzero(fun2,2)b=fzero(fun2,3)8、a=[2 4 9;4 2 4;9 4 18];[v,gama]=eig(a)9、function xpie=fun3(t,x)f=exp(-t)xpie=[0 1 0 0;1 0 0 -1;0 0 0 1;0 -1 -10]*x+[0;1;0;0]*f[t,x]=ode45( 'fun3' ,[0,5],[0 0 0 0]'); plot(t,x(:,1),'b' ,t,x(:,3), 'r' )10、function f=fun4(x,y)f=4*(x-y)-x八2-y八2;v=-2:0.2:2;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)求极小值 (有误) function f=fun5(x,y)x=v(1);y=v(2)f=-(4*(x-y)-x.A2-y.A2);v=-50:2:50;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)min=fminsearch('fun4',[0,0])11、function f=fun6(x,y)f=sin(y)+exp(x)-x.*y.A2ezplot( 'fun6' )12、function f=fun7(x,y)f=(x-y).A2.*(sin(x+y)).A2s=dblquad( 'fun7' ,pi,2*pi,0,pi)14、a=randn(4,4)[l,u]=lu(a)a=randn(4,4)[u,gama,v]=svd(a)15、b=[3 0.5 4];a=[3 4 5 6]; [r,p,k]=residue(b,a)b=[3 0.5 4];a=[3 4 5 6];[r,p,k]=residue(b,a)t=0:0.1:100;y=0*t;for i=1:3 y=y+r(i)*exp(p(i)*t) end plot(t,y) b=[3 0.5 4];a=[3 4 5 6]; impulse(b,a) b=[3 0.5 4];a=[3 4 5 6]; step(b,a)第四章1、n=1:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x)); subplot(2,1,2);stem(y);clear n=0:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x));subplot(2,1,2);stem(y);2、(有误)n=0:19;x=5*0.6.八n;subplot(3,1,1);stem(n,x);for i=-20:39xl(i+21)=x(mod(i+40,20)+1); end n1=-20:39; subplot(3,1,2);stem(n1,x1);x2=x1(10:29); subplot(3,1,3);stem(n,x2)新方法m=10; e=0:19;c=0.6;k=5;a=k*c.A e;a=a';b=circshift(a,m);L=length(a)-1;n=0:L;subplot(2,1,1);stem(n,a);axis([0,L,min(a),max(a)]);subplot(2,1,2);stem(n,b);axis([0,l,min(a),max(a)])3、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);stem(y)4、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);subplot(2,1,1);stem(y)xx=fft(x);hh=fft(h);yy=ifft(xx.*hh);subplot(2,1,2);stem(yy)5、b=[2 3];a=[1 0.4 1];[z p k]=tf2zp(b,a) 6、b=[4 15.6 6 2.4 -6.4];a=[3 2.4 6.3 -11.4 6];[z p k]=tf2zp(b,a)7、b=[18];a=[18 3 -4 -1];[r p k]=residue(b,a)8、b=[0.2 0.3 1];a=[1 0.4 1]; freqs(b,a)9、Fs=2000;wc=[100 200]/(Fs/2);N=10;[b,a]=butter(N,wc, 'z' ); %freqz(b,a)dimpulse(b,a,100)10、wp=100;ws=200;rp=2;rs=15;Fs=500;wp=wp/(Fs/2);ws=ws/(Fs/2);[N,wc]=buttord(wp,ws,rp,rs); [b,a]=butter(N,wc, 'z' ) freqz(b,a)b=fir1(48,[0.35 0.65]);a=1; freqz(b,a)12、b=fir1(37,0.3);a=1; freqz(b,a)13、低通：F=[0:1/56:1];A=[ones(1,50),zeros(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)11、带通：F=[0:1/56:1];A=[ones(1,25),zeros(1,25),ones(1,7)];b=fir2(56,F,A);a=1;freqz(b,a)高通F=[0:1/56:1];A=[zeros(1,50),ones(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)%滤波器设计(椭圆滤波器)wp1=300;wp2=700;ws1=301;ws2=699; rp=0.1;rs=50;Fs=2000; wp=[wp1]/(Fs/2);ws=[ws1]/(Fs/2);[N,wc]=ellipord(wp,ws,rp,rs,'z')[b,a]=ellip(N,rp,rs,wc,'high')%滤波器特性分析freqz(b,a)%信号采样及时频变换Fs=2000;t=0:1/Fs:10;yt=sin(2*pi*200*t)+sin(2*pi*500*t)+sin(2*pi*800*t); figure(2);subplot(3,1,1);plot(t,yt)figure(2);subplot(3,1,2);plot(t,yt);axis([0,0.1,-5,5]) yf=abs(fft(yt));f=(1:length(yf))/length(yf)*Fs;figure(2);subplot(3,1,3);plot(f,yf)%信号滤波及时频分析yt1=filter(b,a,yt);figure(3);subplot(3,1,1);plot(t,yt1)figure(3);subplot(3,1,2);plot(t,yt1);axis([0,0.1,-5,5]) yf1=abs(fft(yt1));f=(1:length(yf1))/length(yf1)*Fs; figure(3);subplot(3,1,3);plot(f,yf1)。
matlab课后习题答案1到6章
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：(2)(3)答：(4)答：即逆左乘3.(1)答：5])(2)答：(3)答：(4)答：(5)答：(6)将含有12个元素的向量x转换成34⨯矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0,0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0,0]L4的值为[4, 5, 6]6.当A=[34, NaN, Inf, -Inf, -pi, eps, 0]时，分析下列函数的执行结果：all(A)，any(A)，isnan(A)，isinf(A)，isfinite(A)。
答：all(A)的值为0any(A) 的值为1isnan(A) 的值为[ 0, 1, 0, 0, 0, 0, 0]isinf(A) 的值为[ 0, 0, 1, 1, 0, 0, 0]isfinite(A) 的值为[1, 0, 0, 0, 1, 1, 1]7.用结构体矩阵来存储5名学生的基本情况数据，每名学生的数据包括学号、姓名、专业和6门课程的成绩。
MATLAB原理及应用实验报告第四章答案
《MATLAB 原理及应用》实验报告三.课后练习题答案1.为 ⎪⎩⎪⎨⎧<-=>+=01001x x x xx x y 编写赋值程序。
程序如下：①建立如下的M 文件：x=input('x=');%让用户通过键盘输入数值、字符串或表达式if x>0y=x+1;elseif x==0y=x;else x<0y=x-1;e nd程序执行结果如下>> kh1 %在当前工作目录下，文件名为“kh.1.m ” x=1>> yy =22.使用for ... end循环的array向量编程求出1+3+5...+100 的值程序如下：sum=0;>> for k=1:2:100sum=sum+k;end>> sumsum =25003.计算1+3+5...+100 的值，当和大于1000时终止计算。
程序如下：sum=0;for m=1:2:100; %建立1 3 5….100的向量if sum<=1000 %如果sum小于1000则可以继续加sum=sum+m; %累加elsebreak; %若sum的结果不符合条件就跳出整个循环endend结果为：sum =1024k =653.1计算从1开始多少个自然数之和超过100。
程序如下：>> sum=0;n=0;>> while sum<=100n=n+1;sum=sum+n;end结果为：n =14sum =1054.求1！+2！+3！+……+8！的值程序如下：n=1;sum=1;for m=2:8; %循环7次使得得到各次阶乘n=n*m;sum=sum+n; %累加end结果为:sum =462335.写程序，判断一年是否为闰年，符合下面两条件之一：（1990~2014）A、能被4整除，不能被100整除B、能被400整除程序如下count=0;for y=1990:2014;if((rem(y,4)==0&rem(y,100)~=0)|(rem(y,4)==0&rem(y,400)~=0));count=count+1;endend结果为:count =5。
matlab课后习题答案1到6章
欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
(2) 删除矩阵A 的第7号元素。
答：A(7)=[](3) 将矩阵A 的每个元素值加30。
答：A=A+30;(4) 求矩阵A 的大小和维数。
答：size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。
答：t(find(t==0))=eps; (6) 将含有12个元素的向量 x 转换成34⨯矩阵。
答：reshape(x,3,4); (7) 求一个字符串的ASCII码。
答：abs(‘123’); 或E 。
答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。
答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，11⎡⎤⎢⎥答：student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。
B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少？答：size(B) 的值为2, 2。
ndims(B) 的值为2。
(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A，其元素为[100,200]范围内的随机整数。
MATLAB与在电子信息课程中的应用课后答案
1. 求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x程序A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];B=[4;-3;9;-8];X=A\B解：X =-1.4841-0.68160.5337-1.24292. 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B 求C1=AB’;C2=A’B;C3=A.*B,并求它们的逆阵。
程序 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];C1= A*B',C2 = A'*B, C3 = A.*Binv(C1),inv(C2),inv(C3)3. a. 列出2×2阶的单位矩阵I, 4×4阶魔方矩阵M 和4×2阶的全幺矩阵A,全零矩阵Bb. 将这些矩阵拼接为6×6阶的矩阵C:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=M B A I C ' c. 求出C 的第2,4,6行,组成3×6阶的矩阵C1,及第2,4,6,裂,组成6×3阶的矩阵C2,d. 求D=C1C2及D1=C2C1.程序 >>I=eye(2),A=ones(4,2),B=zeros(4,2),M=magic(4),C=[I,A';B,M]>> C1=C([2,4,6,],:),C2=C(:,[2,4,6,])>> D=C1*C2, D1=C2*C14．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线解：程序x=linspace(0,2*pi,101)y=cos(x).*(0.5+3*sin(x)./(1+x.^2));plot(x,y),grid5．求代数方程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。
MATALB第四章部分课后答案
习题4.1绘制图形)y t--=8t)(tsin(21e≤t并在横轴标注Time纵0≤轴标注“amplitude”，图形的标题为Decaying-oscillating Exponential t=0:0.1:8;y=1-2*exp(-t).*sin(t);plot(t,y,'k')grid ontitle('Decaying-oscillating Exponential')xlabel('Time');ylabel('amplitude')4.2绘制如下函数图形并加上适当的修饰3008.0)309.0cos(5)(22.0≤≤+︒-=--t e t e t y ttt=0:0.1:30;y=5*exp(-0.2*t).*cos(0.9*t-pi/6); plot(t,y,'r:')4.3在同一张图中绘制下列函数曲线100625.0)24083.2cos(23.1)(625.0)(≤≤+︒+==t t t z t yt=0:0.1:10;y=0.625;z=1.23*cos(2.83*t+4*pi/3)+0.625;plot(t,y,'r:')hold onplot(t,z,'b:')text(0,(0.625-1.23*0.5),'这是z(t=0)的位置')text(10,4*pi/3,'这是z(t=10)的位置')4.4对应于250≤≤t 区间内，在同一图中绘制下列函数曲线tt t te t e t y e t y et y ----+︒-===25.1)128554.0cos(02.2)(02.2)(25.1)(3.033.021t=0:0.1:25;y1=1.25*exp(-t);y2=2.02*exp(-0.3*t);y3=2.02*exp(-0.3*t).*cos(0.554*t-128)+1.25*exp(-t);plot(t,y1,'k')hold onplot(t,y2,'b')hold on plot(t,y3,'r')4.5已知椭圆的长、短轴分别为a=4，b=2，用“小红点线”画椭圆 {)cos()sin(t a x t b y ==t=0:pi/100:2*pi;a=4;b=2;x=a.*cos(t);y=b.*sin(t);plot(x,y,'r.')4.6在同一图形框中绘制两个子图，分别显示下列曲线x=0:0.01:10;y1=sin(2*x).*sin(3*x);y2=0.4*x;subplot(1,2,1),plot(x,y1,'r'),title('y1')subplot(1,2,2),plot(x,y2,'b'),title('y2')4.7试将图形分割成3个子图，并分别绘制lgx在100≤区间内对数坐0≤t标、x半对数坐标和y半对数坐标的函数曲线，并加上适当的图形修饰。
（完整版）matlab经典习题及解答
（完整版）matlab经典习题及解答第1章 MATLAB 概论1.1 与其他计算机语⾔相⽐较，MATLAB 语⾔突出的特点是什么？MATLAB 具有功能强⼤、使⽤⽅便、输⼊简捷、库函数丰富、开放性强等特点。
1.2 MATLAB 系统由那些部分组成？MATLAB 系统主要由开发环境、MATLAB 数学函数库、MATLAB 语⾔、图形功能和应⽤程序接⼝五个部分组成。
1.4 MATLAB 操作桌⾯有⼏个窗⼝？如何使某个窗⼝脱离桌⾯成为独⽴窗⼝？⼜如何将脱离出去的窗⼝重新放置到桌⾯上？在MATLAB 操作桌⾯上有五个窗⼝，在每个窗⼝的右上⾓有两个⼩按钮，⼀个是关闭窗⼝的Close 按钮，⼀个是可以使窗⼝成为独⽴窗⼝的Undock 按钮，点击Undock 按钮就可以使该窗⼝脱离桌⾯成为独⽴窗⼝，在独⽴窗⼝的view 菜单中选择Dock ……菜单项就可以将独⽴的窗⼝重新防⽌的桌⾯上。
1.5 如何启动M ⽂件编辑/调试器？在操作桌⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M ⽂件编辑/调试器将被启动。
在命令窗⼝中键⼊edit 命令时也可以启动M ⽂件编辑/调试器。
1.6 存储在⼯作空间中的数组能编辑吗？如何操作？存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。
1.7 命令历史窗⼝除了可以观察前⾯键⼊的命令外，还有什么⽤途？命令历史窗⼝除了⽤于查询以前键⼊的命令外，还可以直接执⾏命令历史窗⼝中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M ⽂件中。
1.8 如何设置当前⽬录和搜索路径，在当前⽬录上的⽂件和在搜索路径上的⽂件有什么区别？当前⽬录可以在当前⽬录浏览器窗⼝左上⽅的输⼊栏中设置，搜索路径可以通过选择操作桌⾯的file 菜单中的Set Path 菜单项来完成。
在没有特别说明的情况下，只有当前⽬录和搜索路径上的函数和⽂件能够被MATLAB 运⾏和调⽤，如果在当前⽬录上有与搜索路径上相同⽂件名的⽂件时则优先执⾏当前⽬录上的⽂件，如果没有特别说明，数据⽂件将存储在当前⽬录上。
matlab课后练习习题及答案详解
matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较，MATLAB语言突出的特色是什么？MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。
MATLAB系统由那些部分构成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
安装MATLAB时，在选择组件窗口中哪些部分一定勾选，没有勾选的部分此后怎样补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框能否被勾选来决定，能够依据自己的需要选择安装内容，但基本平台（即MATLAB选项）一定安装。
第一次安装没有选择的内容在补安装时只要依据安装的过程进行，不过在选择组件时只勾选要补装的组件或工具箱即可。
MATLAB操作桌面有几个窗口？怎样使某个窗口离开桌面成为独立窗口？又怎样将离开出去的窗口从头搁置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是封闭窗口的Close按钮，一个是能够使窗口成为独立窗口的Undock 按钮，点击Undock按钮就能够使该窗口离开桌面成为独立窗口，在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。
怎样启动M文件编写/调试器？在操作桌面上选择“成立新文件”或“翻开文件”操作时，M文件编写/调试器将被启动。
在命令窗口中键入edit命令时也能够启动M文件编写/调试器。
储存在工作空间中的数组能编写吗？怎样操作？储存在工作空间的数组能够经过数组编写器进行编写：在工作空间阅读器中双击要编写的数组名翻开数组编写器，再选中要改正的数据单元，输入改正内容即可。
命令历史窗口除了能够察看前方键入的命令外，还有什么用途？命令历史窗口除了用于查问从前键入的命令外，还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
怎样设置目前目录和搜寻路径，在目前目录上的文件和在搜寻路径上的文件有什么差别？目前目录能够在目前目录阅读器窗口左上方的输入栏中设置，搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。
Matlab语言及其在电子信息工程中的应用：第4讲V.ppt
[输出参数1, 输出参数2, …]=函数名（输入参数1，输入参数2，…）
例1： 给定数a,b,和一个正整数n，对k=1,2,…,n，计算
(a b)k , (a b)k
(1) 先建立一个函数文件power_ab.m function [out1,out2]=power_ab(a,b,n) %power.m是计算(a+b)^n 和(a-b)^n的函数文件
显示出错信息*****，并终止程序； ➢warning(‘****’)
显示告警信息****，程序继续运行。
许钢 2020/4/19
Matlab语言及其在电子信息工程中的应用
9
第四讲 MATLAB程序设计
三、字符与字符串
1、在MATLAB中字符串要用单引号来界定 注：单引号必须是在英文状态下输入 例： a='Hello' a=[a, ' World'] %字符串可以合并 >>a =
Matlab语言及其在 电子信息工程中的应用
主讲：许钢
（V1.2）
许钢 2020/4/19
Matlab语言及其在电子信息工程中的应用
1
第四讲 MATLAB程序设计
一、M文件 二、控制语句 三、字符与字符串 四、函数的调用 五、程序的优化
（ “逻辑判断和流程控制”其编程思想和C、 Basic等高级语言类似，具体语法请同学们自学！）
deblank 移去空字符
strcmp 比较字符串
num2str 将数字变为字符串
str2num 将字符串变为数字
str2mat 将字符号串变为文本矩阵
strrep 将一个字符串变为另一个字符串
许钢 2020/4/19
matlab课后习题答案第四章
第4章数值运算习题 4 及解答1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff或数值梯度gradient指令计算)(t y'，然后把)(t y和)(t y'曲线绘制在同一张图上，观察数值求导的后果。
（模拟数据从prob_data401.mat 获得）〖目的〗●强调：要非常慎用数值导数计算。
●练习mat数据文件中数据的获取。
●实验数据求导的后果●把两条曲线绘制在同一图上的一种方法。
〖解答〗（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat（2）用diff求导的指令dt=t(2)-t(1);yc=diff(y)/dt; %注意yc的长度将比y短1plot(t,y,'b',t(2:end),yc,'r')grid on（3）用gradent 求导的指令（图形与上相似）dt=t(2)-t(1);yc=gradient(y)/dt;plot(t,y,'b',t,yc,'r')grid on〖说明〗● 不到万不得已，不要进行数值求导。
● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级以上。
● 求导会使数据中原有的噪声放大。
2 采用数值计算方法，画出dt tt x y x ⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖提示〗● 指定区间内的积分函数可用cumtrapz 指令给出。
● )5.4(y 在计算要求不太高的地方可用find 指令算得。
〖目的〗● 指定区间内的积分函数的数值计算法和cumtrapz 指令。
● find 指令的应用。
〖解答〗dt=1e-4;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3)ii=find(t==4.5);s45=s(ii)s45 =1.65413 求函数x ex f 3sin )(=的数值积分⎰=π0 )(dx x f s ，并请采用符号计算尝试复算。
MATLAB及在电子信息课程中的应用.ppt
1.3.1 命令窗(Command Window) 单行命令执行方式
执行结果直接显示
1.3.2 图形窗(Figure Window) 用图形方式表示计算结果
1.3.3 文本编辑窗(File Editor) 多行命令组成语言组，可以文件方式存盘
•
下面就具体看一下MATLAB 的工作环境演示。
2020/8/16
• 1.1.1 MATLAB 概述
MATLAB是集数值计算、符号运算及图 形处理等强大功能于一体的科学计算语言，是 一种交互式的以矩阵为基础的系统计算平台,它 用于科学和工程的计算与可视化。它的优点在 于快速开发计算方法，而不在于计算速度。
MATLAB已成为一门高校必修的课程， 也是最为普遍的计算工具之一。
2020/8/16
8
1.2 MATLAB语言的特点
• 友好的工作平台和编程环境 • 简单易用的程序语言 • 强大的科学计算及数据处理能力 • 出色的图形处理功能 • 应用广泛的模块集和工具箱 • 实用的程序接口和发布平台 • 模块化的设计和系统级的仿真
2020/8/16
9
1.3 MATLAB的工作环境
2020/8/16
3
第1章 MATLAB语言概述
• 1.1 MATLAB语言的发展 • 1.2 MATLAB语言的特点 • 1.3 MATLAB的工作环境
1.3.1 命令窗 1.3.2 图形窗 1.3.3 文本编辑窗 • 1.4 演示程序 • 1.5 网络资源
2020/8/16
4
1.1 MATLAB语言的发展
t=t+dt; theta=theta0*cos(sqrt(g/l)*t); x=l*sin(theta);y=-l*cos(theta); set(head,'xdata',x,'ydata',y); set(body,'xdata',[0;x],'ydata',[0;y]); drawnow; end
matlab及在电子信息课程中的应用课后习题答案第四章
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),shgx3=4.*cos(t)+3.*t;y3=4.*sin(t)+3;plot(x3,y3),shghold 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),pauseaxis([-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:41for j=1:41if z(i,j)==z1(i,j)a=z(i,j);plot3(i,j,a,'r+');endendendhold off;disp('题解4.18');pause;type f1.mi=1;a(1)=fzero('f1',-2);for m=-2:1:8b=fzero('f1',m);if abs(a(i)-b)>0.00001i=i+1;a(i)=b;endendapause,echo on,fplot('f1',[-2 8]);shg;hold on,plot(a,0,'ro'),shg,pause,axis([-2 8 -20 20]),grid,hold off;echo off;disp('题解4.18第二问');pause; type f11.mi=1;a(1)=fzero('f11',-2);for m=-2:1:8b=fzero('f11',m);if abs(a(i)-b)>0.00001i=i+1;a(i)=b;endendapause,echo on,fplot('f11',[-2 8]);shg;hold on,plot(a,0,'ro'),shg,pause,axis([-2 8 -20 20]),grid,hold off;echo off;disp('题解4.19');echo on;x=[0:0.05:4];y=f31(x);plot(x,y)echo off;disp('题解420');pause;type f420.mi=1;a(1)=fzero('f420',0.1);for m=0.1:0.1:4b=fzero('f420',m);if abs(a(i)-b)>0.0001i=i+1;a(i)=b;endendapause,echo on,fplot('f420',[0 4]);shg;hold on,plot(a,f420(a),'ro'),shg,pause, axis([0 4 -20 20]),grid,hold off;pause,quad('f420',a(1),a(2)),echo off;disp('题解421');pause;type f421.mecho on[x,y]=ode23('f421',[0 5],1);plot(x,y),shg,pause,p=polyfit(x,y,9),xi=0:0.1:5;yi=polyval(p,xi);plot(x,y,'ro',xi,yi,'b')echo off,disp('题解422');pause;echo onpause,s='magic(3)',pause,eval(s);pause,y(3,:)='magic(3)';y(4,:)='magic(4)';y(5,:)='magic(5)';y,pausefor n=3:5eval(y(n,:));endecho offdisp('题解424');pause;echo one=exp(1);pause,st=sprintf('自然对数底数e=%21.20f',e);disp(st) echo offdisp('题解425');pause;echo on='Jone';student.birthplace='Shanghai';student(2).name='David';student(2).birthplace='Nanjing';student(3).name='Tom';student(3).birthplace='Hangzhou';student,pause,name={'Jone';'David';'Tom'};birthplace={'Shanghai';'Nanjing';'Hangzhou'}; studentt={name;birthplace};studentt,pause,student1=['Jone ';'Shanhai'],student2=['David ';'Nanjing'],student3=['Tom ';'Hongzhou'],echo off调用子程序f1.mfunction y=f1(x)y=x.^5-4*x.^4+3*x.^2-2*x+6;f11.mfunction y=f11(x)y=x.^5-4*x.^4+3*x.^2-2*x+6+x.*sin(x);f31.mfunction f1=f31(x)f1=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);f420.mfunction y=f420(x)y=x.^3-2*(x.^2).*sin(x)+5*x.*cos(x)+1./(x); f421.mfunction yp=f421(x,y)yp=(x.^2)./y-x.*cos(y);。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
disp('题解 4.18 第二问');pause; type f11.m i=1; a(1)=fzero('f11',-2); for m=-2:1:8
b=fzero('f11',m); if abs(a(i)-b)>0.00001
i=i+1; a(i)=b;
end end a pause,echo on, fplot('f11',[-2 8]);shg; hold on, plot(a,0,'ro'),shg,pause, axis([-2 8 -20 20]),grid, hold off;echo off; disp('题解 4.19'); echo on; x=[0:0.05:4]; y=f31(x); plot(x,y) echo off; disp('题解 420');pause; type f420.m i=1; a(1)=fzero('f420',0.1); for m=0.1:0.1:4
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
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内 纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
i=i+1; a(i)=b;
end end a pause,echo on, fplot('f1',[-2 8]);shg; hold on, plot(a,0,'ro'),shg,pause, axis([-2 8 -20 20]),grid, hold off;echo off;
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内 纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定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);
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
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
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内 纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。