课本中若干习题解答

课本中若干习题MATLAB解答祝大家元旦快乐、学业有成！代向家长问好，春节愉快、合家欢乐！曹青松6#-228，，87046160（O），，2012-12-27老曹较为了解的方向控制理论及应用传感器技术、无损检测技术振动、测试与诊断机电一体化日后相关问题，欢迎讨论；可以替老曹“全方位”宣传（赞与骂都可以），欢迎你认识的学弟学妹们来老曹这“受苦受累”！（远离223（老曹），珍爱生命）快乐脉冲与痛苦阶跃另解！>>sys=tf(1,[1,2,3]);>>impulse(sys)>>step(sys)围绕老曹发生的事某某某同学说：曹老师你也从学生过来的，你也知道，……一句话：老曹会身体力行！让研究生周一跑步，我会坚持；让同学8：30到实验室，我会带头！QQ日志：“不要挑战老曹，应该要挑战你自己”围绕考试发生的事情怎么考啊？难！课又上得差，听不懂，书又难死了，看不懂……这日子没办法过了……还考这么多题？给题不就得了！又不说答案？天哪，答案又不复印下来？天哪，还闭卷！直接复印在答题纸上更省事……围绕作业发生的事情1、电子邮件的那礼仪！2、作业值与不值？研究生的成长之路？？？A）科学研究相当于众人挖隧道！B）敢于怀疑权威，敢于胜利！C）不要轻言“难”！1）大学研究生课程阶段，精读精典教材！2）研究开始阶段（思考方向、博缆全书），不求质量但求数量、看目录、看摘要3）深入（成果）阶段，钻进去，出得来！4）螺旋上升式学习研究理论精读泛读精读6、登堂入室其他工作者权威团队自己外部关系%P33example2-2%P33example2-2a=1.0;b=1.5;

whileb-a>epsx=a+(b-a)/2iff_ex22(a)*f_ex22(x)>0a=x;elseb=x;endend%P36example2-4x=1.5fori=1:10x=(x+2)^(1/4)end

%P42example2-7x=1.5fori=1:10x=(3*x^4+2)/(4*x^3-1)endP42例2.8C=3x=0.5fori=1:10x=1/2*(x+3/x)endx=0.5000000000000003.2500000000000002.0865384615384621.7621632399858211.7323080932066351.7320508266751491.732050807568877例2.9%P54example2-9x01=1;x02=2;x0=[x01;x02];

fori=1:8x01=x0(1,1);x02=x0(2,1);

f01=x01+2*x02-2;f02=x01^2+4*x02^2-4;J01=[1,2];J02=[2*x01,8*x02];s0=inv([J01;J02])*[f01;f02]x1=x0-s0

x0=x1;endP68example3-4h=hilb(3);h_1=inv(h);r=max(sum(abs(h),2))r1=max(sum(abs(h_1),2))r*r1P70example3-5P73example3-6P75example3-7P79example3-8P80example3-9P88example3-12P94example3-14A=[5-1-1;-13-1;-1-15];A(:,1)=A(:,1)./sqrt(A(1,1));A(2,2)=sqrt(A(2,2)-A(2,1)^2);A(3,2)=(A(3,2)-A(3,1)^2)/A(2,2);A(3,3)=sqrt(A(3,3)-A(3,1)^2-A(3,2)^2);%P121example4-3A=[1031;2-103;1310];b=[14-514]';x0=[000]';

%JacobiAA=[0-1/10*3-1/10;1/503/10;-1/10-3/100];bb=[14/105/1014/10]';fori=1:9x(:,i)=AA*x0+bb;x0=x(:,i);endx%G-Sb=[14-514]';x0=[000]';D=[1000;0-100;0010];L=[000;-200;-1-30];U=[0-3-1;00-3;000];fori=1:9xx(:,i)=inv(D-L)*U*x0+inv(D-L)*b;x0=xx(:,i);endxx%SORb=[14-514]';x0=[000]';D=[1000;0-100;0010];L=[000;-200;-1-30];U=[0-3-1;00-3;000];w=1.1;fori=1:9x_sor(:,i)=inv(D./w-L)*((1/w-1)*D+U)*x0+inv(D./w-L)*b;x0=x_sor(:,i);endx_sorP131example4-7A=[32;26];b=[2-8]';x0=[-2-2]';

fori=1:8r0=b-A*x0;a0=r0'*r0/(r0'*A*r0);x1=x0+a0*r0;x(:,i)=x1;x0=x1;endxP157example5-6a=[3,1;1,3];v0=[0,1]';fori=1:20v=a*v0;max(v)V=v/max(v);Eig=(V'*a*V)/(V'*V);v0=V;endP160，例5.8a=[3,1;1,3];v0=[0,1]';fori=1:10v=inv(a)*v0;max(v)V=v/max(v)v0=V;endP165example5-9P167example5-10P168example5-11P172example5-12%P202example6-6A=[11111;1.001.251.51.752.00]';f=[1.62921.75611.87642.00822.1353]';

x=inv(A'*A)*A'*f

%P205example6-8Sigma1=sqrt(sum(A(:,1).^2));v1=A(:,1)+Sigma1*[10000]';A1=A(:,1)-2*(v1'*A(:,1))./(v1'*v1)*v1A2=A(:,2)-2*(v1'*A(:,2))./(v1'*v1)*v1f2=f-2*(v1'*f)./(v1'*v1)*v1AA=[A1,A2]

A2(1)=0;Sigma2=-sqrt(sum(A2.^2))v2=A2+Sigma2*[01000]';f3=f2-2*(v2'*f2)./(v2'*v2)*v2A3=A2-2*(v2'*A2)./(v2'*v2)*v2A3(1)=AA(1,2)AAA=[A1,A3]

x=inv(AAA(1:2,1:2))*f3(1:2,1)P202example6.56.6A=[11111;1.001.251.51.752.00]';f=[1.62921.75611.87642.00822.1353]';

x=inv(A'*A)*A'*fP2136.9,P2176.10书本已较详细P218example6-11x=[0.40.550.650.80.91.05];f=[0.410750.578150.696750.888111.026521.25382];

diff_1=(f(2:6)-f(1:5))./(x(2:6)-x(1:5))

diff_2=(diff_1(2:5)-diff_1(1:4))./(x(3:6)-x(1:4))

diff_3=(diff_2(2:4)-diff_2(1:3))./(x(4:6)-x(1:3))

diff_4=(diff_3(2:3)-diff_3(1:2))./(x(5:6)-x(1:2))

diff_5=(diff_4(2)-diff_4(1))./(x(6)-x(1))P230exmaple6-13x=1:6;y=[161821171512];xi=0.75:0.05:6.25;p=polyfit(x,y,5);v1=polyval(p,xi);v2=interp1(x,y,xi,'linear');v3=interp1(x,y,xi,'pchip');v4=interp1(x,y,xi,'spline');subplot(2,2,1);plot(x,y,'o',xi,v1,'-');subplot(2,2,2);plot(x,y,'o',xi,v2,'-');subplot(2,2,3);plot(x,y,'o',xi,v3,'-');subplot(2,2,4);plot(x,y,'o',xi,v4,'-');P245example7-3I0=(1-0)*exp(-(1/2)^2)I1=(1-0)/2*(exp(-0^2)+exp(-1^2))I2=(1-0)/6*(exp(-0^2)+4*exp(-(1/2)^2)+exp(-1^2))P249example7-4x=0:1/8:1;f=sin(x)./x;f(1)=1;

n=8;h=(1-0)/n;T=1/2*h*(f(1)+f(9)+2*sum(f(2:8)))

n=4;h=(1-0)/n;S=h/6*(f(1)+f(9)+4*sum(f(2:2:8))+2*sum(f(3:2:7)))%P255example7-5%T01=1/2*(f(1)+f(9))T02=T01/2+1/2*(f(5))T03=T02/2+1/4*(f(3)+f(7))

%T11=1/3*(4*T02-T01)T12=1/3*(4*T03-T02)

%T21=1/15*(16*T12-T11)P255example7-6T0=0.5;h=1;fori=1:5a=1:2:2^i-1;

x=a/(2^i);sum1=sum(x.^(3/2));T=1/2*T0+h/2*sum1T0=T;

h=h/2;endTT=[0.5,TT]TT1=1/3*(4*TT(2:6)-TT(1:5))TT2=1/15*(16*TT1(2:5)-TT1(1:4))TT3=1/(4^3-1)*(4^3*TT2(2:4)-TT2(1:3))TT4=1/(4^4-1)*(4^4*TT3(2:3)-TT3(1:2))TT5=1/(4^5-1)*(4^5*TT4(2)-TT4(1))P263example7-9a=[0.23692690.23692690.47862870.47862870.5688889];x=[0.9061798-0.90617980.5384693-0.53846930];b=sin(0.5+0.5*x)./(0.5+0.5*x)c=a.*bsum(c)/2P267example7-10fori=1:12h=10^(-i);x=2;Dc=(sqrt(x+h)-sqrt(x-h))/(2*h)endP270example7.11h=0.1;x=0.5;x1=x+h;x2=x-h;Dc1=(x1^2*exp(-x1)-x2^2*exp(-x2))/(2*h);

h=0.05;x=0.5;x1=x+h;x2=x-h;Dc2=(x1^2*exp(-x1)-x2^2*exp(-x2))/(2*h);h=0.025;x=0.5;x1=x+h;x2=x-h;Dc3=(x1^2*exp(-x1)-x2^2*exp(-x2))/(2*h);D1c1=1/3*(4*Dc2-Dc1)D1c2=1/3*(4*Dc3-Dc2)

D2c1=1/15*(16*D1c2-D1c1)P280example8.4y0=1;h=0.1;t0=0;

fori=1:5

t=t0+h;y=y0+h*(t-y0+1)yy=t+exp(-t)

t0=t;y0=y;endP289Example8.7format