MATLAB教程R2012a课后习题答案.pptxVIP

?;2;3;4;5;6;〖解答〗 A=sym([a11 a12 a13;a21 a22 a23;a31 a32 a33]) DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31);? 符号变量的限定性定义的作用。 〖解答〗 syms k; x=sym(x,positive); f_k=2/(2*k+1)*((x-1)/(x+1))^(2*k+1); s=simple(symsum(f_k,k,0,inf)) %结果与理论值 lnx 相符！ s = piecewise([abs(x - 1) x + 1, log(x)]) 〖注意〗 ? 解答中，条件 abs(x - 1) x + 1 意味着：;9;;11;;;14;;G1 G2 G3 G4 H3 H5 + G1 G3 G4 G5 H2 H4 + 1) （2）局部性验证 syms a b c d e f g y2u=subs(Y2U,[G1,G2,G3,G4,G5,G6,G7,H1,H2,H3,H4,H5],[a,e,f,1,b,c,0,g,0 ,0,0,d]); [nn,dd]=numden(y2u); dd=sort(dd); disp([blanks(5),局部性验证用的传递函数 y2u]) pretty(nn/dd) 局部性验证用的传递函数 y2u a (e f + b c f + d e f) --------------------------- d + e g + b c g + d e g + 1 此结果与胡寿松主编的“自动控制原理”例 2-21 一致。;Y ? ?c(I ? A)?1 b ? d ?U ? ?c(I ? A)?1 f ? g?W 据此可写出传递函数;;axis square grid on;;21;22;23;〖解答〗 %标量循环运算法 t=linspace(0,10,200); N=length(t); y1=zeros(size(t)); for k=1:N y1(k)=1-exp(-0.5*t(k))*cos(2*t(k)); end subplot(1,2,1),plot(t,y1),xlabel(t),ylabel(y1),grid on %数组运算法 y2=1-exp(-0.5*t).*cos(2*t); subplot(1,2,2),plot(t,y2),xlabel(t),ylabel(y2),grid on;25;;27;或数值梯度 gradient 指令计算 y?(t) ，然后把 y(t) 和 y?(t) 曲线绘制在 同 一 张 图 上 ， 观 察 数 值 求 导 的 后 果 。 （ 模 拟 数 据 从 prob_data401.mat 获得） 〖目的〗 强调：要非常慎用数值导数计算。 练习 mat 数据文件中数据的获取。 实验数据求导的后果 把两条曲线绘制在同一图上的一种方法。 〖解答〗 从数据文件获得数据的指令 假如 prob_data401.mat 文件在当前目录或有哪些信誉好的足球投注网站路径上 clear load prob_data401.mat 用 diff 求导的指?? dt=t(2)-t(1); yc=diff(y)/dt; %注意 yc 的长度将比 y 短 1 plot(t,y,b,t(2:end),yc,r) grid on;上。 ? 求导会使数据中原有的噪声放大。;〖提示〗 数值积分均可尝试。 符号积分的局限性。 〖目的〗 符号积分的局限性。 〖解答〗 dx=pi/2000; x=0:dx:pi; s=trapz(exp(sin(x).^3))*dx s =