资源描述:
《《MATLAB与仿真系统》上机测验09》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、《MATLAB与仿真系统》上机测验题号一二三四五六七八九合计分数一、简述题。(共15分)1.简述MATLAB语言特点。(5分)2.简述三维绘图中绘制曲线的最基本函数以及三维曲面的基本函数,并且阐述绘制三维曲面时采用不同函数时的区别。试通过实际例子加以说明。(5分)3.试举例说明矩阵的乘法运算和点乘运算的不同。(5分)答:A=[1,2,3;3,2,1;1,1,1],B=[101;011;110];C=A*BC=453435222D=A.*BD=103021110二、已知,求的特征值以及特征向量。(共10分)答:程序A=[-9,6,8;2,5,12;-8,8,5];[V,D]=eig(A
2、)结果:V=-0.1145+0.2123i-0.1145-0.2123i0.38930.73000.73000.8075-0.6364+0.0625i-0.6364-0.0625i0.4431D=-5.7747+1.6096i000-5.7747-1.6096i00012.5494三、已知(共15分)求(1).;(2).的逆矩阵(3).包括矩阵主对角元素的上三角矩阵(4).的行列式值答:所有程序:symsabcdefghi;P1=[010;100;001];P2=[100;010;101];A=[abc;def;ghi];%(1)B=P1*P2*A%(2)b=inv(B)%(3)c=
3、triu(B)%(4)d=det(B)运行结果:B=[d,e,f][a,b,c][a+g,b+h,c+i]b=[-(b*i-c*h)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b),(e*c+e*i-f*b-f*h)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b),-(e*c-f*b)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b)][(a*i-c*g)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b),(-d*c-d*i+f*a+f*g)/(-d*b*i+d*c*h
4、+a*e*i-a*f*h-g*e*c+g*f*b),-(-d*c+f*a)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b)][-(a*h-b*g)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b),-(-d*b-d*h+e*a+e*g)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b),(-d*b+e*a)/(-d*b*i+d*c*h+a*e*i-a*f*h-g*e*c+g*f*b)]c=[d,e,f][0,b,c][0,0,c+i]d=d*b*i-d*c*h-a*e*i+a*f*h+g*e*c-g
5、*f*b四、用符号法求下列积分(共10分)(1).答:程序:clearsymsx;f=1/(1+2*x^4+x^8);a=int(f,x)结果:a=1/4*x/(1+x^4)+3/32*2^(1/2)*log((x^2+x*2^(1/2)+1)/(x^2-x*2^(1/2)+1))+3/16*2^(1/2)*atan(x*2^(1/2)+1)+3/16*2^(1/2)*atan(x*2^(1/2)-1)(2).答:程序:clearsymsx;f=(x^2+1)/(x^4+1);y=int(f,x,-inf,0)结果:y=1/2*2^(1/2)*pi五、求级数之和(共10分)(1).答
6、:程序:clearsymsnx;f=n^2*x^(n-1);s=symsum(f,n,1,inf)结果:s=(-x-1)/(x-1)^3(2).答:程序:clearsymsn;f=n^2/5^n;s=symsum(f,n,1,inf)结果:s=15/32六、绘图:(共10分)绘制极坐标曲线并说明参数k,m,n对图形的影响。答:程序:cleark=input('k=');m=input('m=');n=input('n=');theta=0:0.01:2*pi;rho=k*sin(m*theta).*cos(n*theta);polar(theta,rho)分别输入k=1,m=1,n=
7、1,可得如下图形图1极坐标图k=1,m=1,n=1当输入k=1,m=3,n=1时图2k=1,m=3,n=1当输入k=1,m=1,n=3时图3k=1,m-1,n=3七、编程(共10分)对于向量A里面的6个数据进行加密后输出。规则如下:每个数字减去8,然后再除以7的余数代替该数字;把第一位与第二位交换,第三位和第五位交换,第四位和第六位交换。(10分)答:程序:A=[];fori=1:6a=input('A(i)');A=[A,a];endB=A-8;C=re