会计学4.1特征值与特征向量4.1.1实验目的1.掌握MATLAB软件线性变换的几何含义；2.掌握MATLAB软件求解(qiújiě)特征值和特征向量的方法。4.1.2实验(shíyàn)指导4.1.3实验(shíyàn)内容附录程序(chéngxù)1.绘制向量程序(chéngxù)：%绘制二维向量函数functiondrawvec(u)plot([0;u(1)],[0;u(2)]);%画向量线段holdontheta=acos(u(1)/norm(u));%计算夹角if(u(2)<0)theta=2*pi-theta;%当向量在第四象限endfill([u(1)-0.5*cos(theta+pi/12),u(1),u(1)-0.5*cos(theta-pi/12)],[u(2)-0.5*sin(theta+pi/12),u(2),u(2)-0.5*sin(theta-pi/12)],'black');%用线段填充holdoff2．线性变换程序：%线性变换的几何(jǐhé)含义x=[2;1];A1=[-1,0;0,1];A2=[1,0;0,-1];A3=[0.5,0;0,2];A4=[cos(pi/2),sin(pi/2);-sin(pi/2),cos(pi/2)];y1=A1*x;y2=A2*x;y3=A3*x;y4=A4*x;subplot(2,2,1);drawvec(x);holdon;drawvec(y1);axisequal;axis([-3,3,-2,2]);text(x(1),x(2)+0.5,'x');text(y1(1),y1(2)+0.5,'y_1');title('y_1=A_1x');gridon;subplot(2,2,2);drawvec(x);holdon;drawvec(y2);axisequal;axis([-3,3,-2,2]);text(x(1),x(2)+0.5,'x');text(y2(1),y2(2)+0.5,'y_2');title('y_2=A_2x');gridon;subplot(2,2,3);drawvec(x);holdon;drawvec(y3);axisequal;axis([-3,3,-2,2]);text(x(1),x(2)+0.5,'x');text(y3(1)-1,y3(2)-0.2,'y_3');title('y_3=A_3x');gridon;subplot(2,2,4);drawvec(x);holdon;drawvec(y4);axisequal;axis([-3,3,-2,2]);text(x(1),x(2)+0.5,'x');text(y4(1),y4(2)+0.5,'y_4');title('y_4=A_4x');gridon;方法二：用MATLAB命令：>>A=[0.50.25;0.250.5];%输入矩阵(jǔzhèn)A>>[Q,d]=eig(A)%求矩阵(jǔzhèn)A的特征值与特征向量Q=-0.70710.70710.70710.7071d=%注意在对角线上的值才是特征值0.2500000.75004.3.3实验(shíyàn)内容4.4MATLAB在微积分中的应用4.3.1实验目的1.了解MATLAB软件在微积分中应用；2.掌握MATLAB软件求泰勒展开式、不定积分和定积分的方法；3.进一步理解MATLAB软件在数学学习(xuéxí)中的重要作用。4.3.2实验(shíyàn)指导4.4.1实验内容泰勒(tàilè)多项式逼近泰勒(tàilè)多项式逼近//例4.4.1将计算(jìsuàn)关于x-1展式的前3项,(展开至(x-1)2)例4.4.2关于(guānyú)>>ezplot('sin(x)');>>holdon>>ezplot('x');>>holdon>>ezplot('x-x^3/6');>>holdon>>ezplot('x^9/362880-x^7/5040+x^5/120-…x^3/6+x');>>holdon>>ezplot('-x^11/39916800+x^9/362880-…x^7/5040+x^5/120-x^3/6+x');>>gridon/4.4.2实验内容符号(fúhào)积分int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间(qūjiān)[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。例4.4.2已知>>int(f,x,0,2)ans=(8*a)/3+2*b+2*c>>int(f,a)ans=a