当前位置:文档之家› 中南大学系统仿真实验报告

中南大学系统仿真实验报告

实验一MATLAB 中矩阵与多项式的基本运算实验任务1. 了解MATLAB命令窗口和程序文件的调用。

2 •熟悉如下MATLAB的基本运算:①矩阵的产生、数据的输入、相关元素的显示;②矩阵的加法、乘法、左除、右除;③特殊矩阵:单位矩阵、“ 1 ”矩阵、“0”矩阵、对角阵、随机矩阵的产生和运算;④多项式的运算:多项式求根、多项式之间的乘除。

基本命令训练1、>> eye(2)ans =1 00 1>> eye(4)ans =1 0 0 00 1 0 00 0 1 00 0 0 12、>> ones(2)1 1ans =1 1>> ones(4)ans =1 1 1 11 1 1 11 1 1 11 1 1 1 >> ones(2,2)ans =1 11 1>> ones(2,3)ans =1 1 11 1 1>> ones(4,3)ans =1 1 11 1 11 1 11 1 1 3、>> zeros(2)ans =0 00 0>> zeros(4)ans =0 0 0 00 0 0 00 0 0 00 0 0 0>> zeros(2,2)ans =0 00 0>> zeros(2,3)ans =0 0 00 0 0>> zeros(3,2)ans =0 00 0004、随机阵>> rand(2,3) ans =0.2785 0.9575 0.15760.5469 0.9649 0.9706 >> rand(2,3)ans =0.9572 0.8003 0.42180.4854 0.1419 0.9157 5、>> diag(5)ans =5>> diag(5,5) ans =0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> diag(2,3)ans =0 0 0 20 0 0 00 0 0 00 0 0 0 6、(inv (A)为求A的逆矩阵) >> B=[5 3 1;2 3 8;1 1 1],inv(B)5 3 12 3 81 1 1ans =0.6250 0.2500 -2.6250-0.7500 -0.5000 4.75000.1250 0.2500 -1.1250>> A=[2 3;4 4],B=[5 3;3 8],inv(A),inv(B);A\B,A/B,inv(A)*B,B*inv(A) A =2 34 4B =5 33 8ans =-1.0000 0.75001.0000 -0.5000ans =-2.7500 3.00003.5000 -1.0000ans =0.2258 0.29030.6452 0.2581ans =-2.7500 3.00003.5000 -1.0000 ans =-2.0000 2.25005.0000 -1.75007、>> p =[1,-6,-72,-27], roots(p)p =1 -6 -72 -27ans =12.1229-5.7345-0.3884>> p=[2,3,6],roots(p)p =2 3 6ans =-0.7500 + 1.5612i-0.7500 - 1.5612i8、( A 为n*n 的方阵)>> A=[0 1 0;-4 4 0;-2 1 2],poly(A),B=sym(A),poly(B) A =0 1 0-4 4 0-2 1 2 ans =1 -6 12 -8B =[ 0, 1,0][ -4, 4, 0][ -2, 1, 2]ans =x A3-6*x A2+12*x-89, 、( conv 是多项式相乘,deconv 是多项式相除) >> u=[1 2 4 6 ],v=[5 0 0 -6 7],conv(u,v)u =1 2 4 6v =5 0 0 -67ans =5 10 20 24 -5 -10 -8 42 >> v=[1 2 46 ],u=[5 0 0 -6 7],deconv(u,v) v =1 2 4 6u =5 0 0 -67 ans =5 -1010、(点乘是数组的运算,没有点的乘是矩阵运算)>> a = [2 5;3 4], b =[3 1;4 7],a.*b,a*ba =2 53 4b =3 14 7ans =6 512 28ans =26 3725 31>> a = [2 3]; b = [4 7];a.*b = [8 21];a*b %错误a*b' = 29;11、(who 可以看到你用过的一些变量,来了)>> whowhos 是把该变量及所存储的大小等信息都显示出Your variables are:A B a ans b p u>> whosName Size BytesClass Attributes2x2 32 doubleB 2x2 32 doublea 1x2 16 d oubleans 1x2 16 d oubleb 1x2 16 d oublep1x3 24 d oubleu 1x5 40 d oublev 1x4 32 d ouble12、>> A=[2 5 3;6 5 4],disp(A),size(A),length(A)A =2 5 36 5 42 5 36 5 4ans =2 3ans =3实验二MATLAB 绘图命令实验任务熟悉MATLAB基本绘图命令,掌握如下绘图方法:1.坐标系的选择、图形的绘制;2 •图形注解(题目、标号、说明、分格线)的加入;3 •图形线型、符号、颜色的选取基本命令训练1、>>t=[0:pi/360:2*pi];x=cos(t)+ cos(t*4);y=si n( t)+ sin (t*4);xlabel('x 轴');ylabel('y 轴');plot(y,x),grid;2、>>t=0:0.1:100; x=3*t;y=4*t;z=si n( 2*t);■15 i 0 5 0 05 1 1 5 2 plot3(x,y,z, 'g:')3、>>x = linspace(-2*pi,2*pi,40);y=si n( x);stairs(x,y)4、>>t=[0:pi/360:2*pi]; x=cos(t)+cos(t*4) + sin (t*4);y=si n( t)+ si n( t*4);plot(y,x, 'r:');xlabel('x 轴');ylabel('y 轴');6、>>th=[0:pi/20:2*pi];x=exp(j*th);plot(real(x),imag(x),'r-.');grid; text(0,0,'中心');5、>>th=[0:pi/1000:2*pi]';r=cos(2*th);polar(th,r);title( '四叶草图')四叶草图107、>>x=-2:0.01:2; 8、y=-2:0.01:2; 9、[X,Y] = meshgrid(x,y); Z = Y.*exp(-X.A 2-Y.A 2); [C,h] = con tour(X, YZ); set(h,'ShowText','o n','TextStep',get(h,'LevelStep')*2)_1O S•I,52■n.s o o.s8、>>x = 0:0.2:10; y = 2*x+3; subplot(411);plot(x,y); grid;title('y 的原函数');subplot(412) ;semilogy(x,y); grid;title('对 y 取对数'); 丫的原画数40 ----------- 1 ------------------------------------------- 1 --------------------- 1 ----------------------1 ------------------------------------------- 1 --------------------- 1---------------------1 --------l|Ip|il■|i|I九 ____ 1-___ — ___ I ____ L ___ :……: ___ J _ _______ L ___ ui|I|lI, il _ -■」 " jI■I__ h- ____________ II■Q 」【I1FIIIIIIIIIT 101 □123456r69 10, 对y 取对数对弋观对数subplot(413) ;semilogx(x,y);10 10 id101lZ10w10104020 0 对好对数grid;title('对x 取对数');subplot(414) ;loglog(x,y);grid;title('对xy均取对数');9、>>x = -3:0.3:3;bar(x,exp(-x.*x),'g')实验三MATLAB程序设计实验任务1 •熟悉MATLAB 程序设计的方法和思路;2 •掌握循环、分支语句的编写,学会使用look for、help命令程序举例1、>>f=[1,1];i=1;while f(i)+f(i+1)<1000f(i+2)=f(i)+f(i+1);i=i+1;endf,if =Columns 1 through 141 123 5 8 13 21 34 55 89 144 233 377 Columns 15 through 16610 987i =152、>>m=3;n=4;for i=1:mfor j=1:na(i,j)=1/(i+j-1);endendformat rat1 1/2 1/3 1/2 1/3 1/4 1/31/41/5 (分数格式形式。

相关主题