资源描述:
《东北大学 matlab实验.docx》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、《MATLAB语言与应用》上机实验作业第一部分2、用MATLAB语句输入矩阵和,前面给出的是矩阵,如果给出命令将得出什么结果?MATLAB结果:>>A=[1234;4321;2341;3241];>>B=[1+4j2+3j3+2j4+1j;4+1j3+2j2+3j1+4j;2+3j3+2j4+1j1+4j;3+2j2+3j4+1j1+4j];>>AA=1234432123413241>>BB=1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.000
2、0i4.0000+1.0000i3.0000+2.0000i2.0000+3.0000i1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.0000i1.0000+4.0000i3.0000+2.0000i2.0000+3.0000i4.0000+1.0000i1.0000+4.0000i>>A(5,6)=5A=1234004321002341003241000000053、假设已知矩阵,试给出相应的MATLAB命令,将其全部偶数行提取出来,赋给矩阵
3、,用命令生成矩阵,用上述命令检验一下结果是不是正确。MATLAB结果:>>A=magic(8);>>AA=64236160675795554121351501617474620214342244026273736303133323435292838392541232244451918484915145253111056858595462631>>B=A(2:2:end,:)B=955541213515016402627373630313341232244451918488585954626314、用
4、数值方法可以求出,试不采用循环的形式求出和式的数值解。由于数值方法是采用double形式进行计算的,难以保证有效位数字,所以结果不一定精确。试采用运算的方法求该和式的精确值。MATLAB结果:>>formatlong;>>sum(2.^[0:63])ans=1.0955e+0195、选择合适的步距绘制出下面的图形。(1),其中;(2),其中。MATLAB结果:(1),其中t=[-1:0.05:-0.2,-0.199:0.001:0.2,0.2:0.05:1];y=sin(1./t);Warning:
5、Dividebyzero.>>plot(t,y)(2),其中>>x=[-pi:0.05:-1.8,-1.7999:0.001:-1.2,-1.2:0.05:1.2,1.201:0.001:1.8,1.81:0.05:pi];>>y=sin(tan(x))-tan(sin(x));>>plot(x,y)6、试绘制出二元函数的三维图和三视图。MATLAB结果:三维图:xx=[-2:0.1:-1.2,-1.1:0.02:-0.9,-0.8:0.1:0.8,0.9:0.02:1.1,1.2:0.1:2];y
6、y=[-1:0.1:-0.2,-0.1:0.02:0.1,0.2:0.1:1];>>[x,y]=meshgrid(xx,yy);>>z=1./(sqrt((1-x).^2+y.^2))+1./(sqrt((1+x).^2+y.^2));Warning:Dividebyzero.Warning:Dividebyzero.>>surf(x,y,z)三视图:>>subplot(224),surf(x,y,z);>>subplot(221),surf(x,y,z),view(0,90);>>subplot(
7、222),surf(x,y,z),view(90,0);>>subplot(223),surf(x,y,z),view(0,0)7、试求出如下极限。(1);(2);(3)。MATLAB结果:(1)>>symsx;>>f=(3^x+9^x)^(1/x);>>L=limit(f,x,inf)L=9(2)>>symsxy;>>f=(x*y)/(((x*y)+1)^(1/2)-1);>>L=limit(limit(f,x,0),y,0)L=2(3)>>symsxy;f=(1-cos(x^2+y^2))/((
8、x^2+y^2)*exp(x^2+y^2));>>L=limit(limit(f,x,0),y,0)L=08、已知参数方程,试求出和MATLAB结果:>>symst;x=log(cos(t));y=cos(t)-t*sin(t);f1=diff(y,t)/diff(x,t)f1=-(-2*sin(t)-t*cos(t))/sin(t)*cos(t)>>f2=diff(y,t,2)/diff(x,t,2);subs(f2,t,sym(pi/3))ans=3/8-1/24